Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/594/A116 into FORTRAN code for loading all data files into arrays.
Note that special values are assigned to unknown or unspecified
numbers (also called NULL numbers);
when necessary, the coordinate components making up the right ascension
and declination are converted into floating-point numbers
representing these angles in degrees.
program load_ReadMe C============================================================================= C F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2018-Apr-25 C============================================================================= * This code was generated from the ReadMe file documenting a catalogue * according to the "Standard for Documentation of Astronomical Catalogues" * currently in use by the Astronomical Data Centers (CDS, ADC, A&A) * (see full documentation at URL http://vizier.u-strasbg.fr/doc/catstd.htx) * Please report problems or questions to C============================================================================= implicit none * Unspecified or NULL values, generally corresponding to blank columns, * are assigned one of the following special values: * rNULL__ for unknown or NULL floating-point values * iNULL__ for unknown or NULL integer values real*4 rNULL__ integer*4 iNULL__ parameter (rNULL__=--2147483648.) ! NULL real number parameter (iNULL__=(-2147483647-1)) ! NULL int number integer idig ! testing NULL number C============================================================================= Cat. J/A+A/594/A116 HI4PI spectra and column density maps (HI4PI team+, 2016) *================================================================================ *HI4PI: a full-sky HI survey based on EBHIS and GASS. * HI4PI collaboration, Ben Bekhti N., Floeer L., Keller R., Kerp J., Lenz D., * Winkel B., Bailin J., Calabretta M.R., Dedes L., Ford H.A., Gibson B.K., * Haud U., Janowiecki S., Kalberla P. M. W., Lockman F.J., * McClure-Griffiths N.M., Murphy T., Nakanishi H., Pisano D.J., * Staveley-Smith L. * <Astron. Astrophys. 594, A116 (2016)> * =2016A&A...594A.116H (SIMBAD/NED BibCode) C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'nhi_eq.dat' ! List of all-sky NHI FITS images, Equatorial integer*4 nr__ parameter (nr__=44) ! Number of records character*31 ar__ ! Full-size record character*3 WCSproj (nr__) ! WCS projection code integer*4 b_RVel (nr__) ! (km/s) Lower radial velocity limit (LSR, Radio) integer*4 B_RVel_1 (nr__) ! (km/s) Upper radial velocity limit (LSR, Radio) character*17 FileName (nr__) ! Name of FITS file in subdirectory NHI/EQ2000 c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'nhi_gal.dat' ! List of all-sky NHI FITS images, Galactic integer*4 nr__1 parameter (nr__1=44) ! Number of records character*31 ar__1 ! Full-size record character*3 WCSproj_1 (nr__1) ! WCS projection code integer*4 b_RVel_2 (nr__1) ! (km/s) Lower radial velocity limit (LSR, Radio) integer*4 B_RVel_3 (nr__1) ! (km/s) Upper radial velocity limit (LSR, Radio) character*17 FileName_1 (nr__1) ! Name of FITS file in subdirectory NHI/GAL c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'allsky_eq.dat' ! List of all-sky FITS data cubes, Equatorial integer*4 nr__2 parameter (nr__2=4) ! Number of records character*12 ar__2 ! Full-size record character*3 WCSproj_2 (nr__2) ! WCS projection code character*8 FileName_2 (nr__2) ! Name of FITS file in subdirectory ALLSKY/EQ2000 c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'allsky_gal.dat' ! List of all-sky FITS data cubes, Galactic integer*4 nr__3 parameter (nr__3=4) ! Number of records character*12 ar__3 ! Full-size record character*3 WCSproj_3 (nr__3) ! WCS projection code character*8 FileName_3 (nr__3) ! Name of FITS file in subdirectory ALLSKY/GAL c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'cubes_eq.dat' ! List of 20dx20d FITS data cubes, Equatorial integer*4 nr__4 parameter (nr__4=534) ! Number of records character*28 ar__4 ! Full-size record C J2000 position composed of: RAdeg DEdeg real*4 RAdeg (nr__4) ! (deg) Right ascension (J2000) of data cube center real*4 DEdeg (nr__4) ! (deg) Declination (J2000) of data cube center character*3 WCSproj_4 (nr__4) ! WCS projection code character*12 FileName_4 (nr__4) ! Name of FITS file in subdirectory CUBES/EQ2000 c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'cubes_gal.dat' ! List of 20dx20d FITS data cubes, Galactic integer*4 nr__5 parameter (nr__5=534) ! Number of records character*28 ar__5 ! Full-size record C Position composed of: GLON GLAT real*4 GLON (nr__5) ! (deg) Galactic longitude of data cube center real*4 GLAT (nr__5) ! (deg) Galactic latitude of data cube center character*3 WCSproj_5 (nr__5) ! WCS projection code character*12 FileName_5 (nr__5) ! Name of FITS file in subdirectory CUBES/GAL c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'nhi_hpx.dat' ! HI column densities integer*4 nr__6 parameter (nr__6=12582912) ! Number of records character*71 ar__6 ! Full-size record C J2000 position composed of: RAdeg DEdeg integer*4 HPX (nr__6) ! HPX index (HPXINDEX) real*8 RAdeg_1 (nr__6) ! (deg) Right ascension (J2000) real*8 DEdeg_1 (nr__6) ! (deg) Declination (J2000) real*8 GLON_1 (nr__6) ! (deg) Galactic longitude (J2000) real*8 GLAT_1 (nr__6) ! (deg) Galactic latitude (J2000) real*4 NHI (nr__6) ! (cm-2) HI column density c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'hpxspec.dat' ! List of HPX-grid FITS binary tables integer*4 nr__7 parameter (nr__7=192) ! Number of records character*52 ar__7 ! Full-size record C J2000 position composed of: RAdeg DEdeg real*8 GLON_2 (nr__7) ! (deg) Galactic longitude real*8 GLAT_2 (nr__7) ! (deg) Galactic latitude real*8 RAdeg_2 (nr__7) ! (deg) Right ascension (J2000) real*8 DEdeg_2 (nr__7) ! (deg) Declination (J2000) integer*4 HPXidx (nr__7) ! HealPIX index (nside=4, RING scheme,0-based) character*12 FileName_6 (nr__7) ! Name of the FITS file in subdirectory HPX C============================================================================= C Loading file 'nhi_eq.dat' ! List of all-sky NHI FITS images, Equatorial C Format for file interpretation 1 format(A3,1X,I4,1X,I4,1X,A17) C Effective file loading open(unit=1,status='old',file= +'nhi_eq.dat') write(6,*) '....Loading file: nhi_eq.dat' do i__=1,44 read(1,'(A31)')ar__ read(ar__,1) + WCSproj(i__),b_RVel(i__),B_RVel_1(i__),FileName(i__) c ..............Just test output........... write(6,1) + WCSproj(i__),b_RVel(i__),B_RVel_1(i__),FileName(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'nhi_gal.dat' ! List of all-sky NHI FITS images, Galactic C Format for file interpretation 2 format(A3,1X,I4,1X,I4,1X,A17) C Effective file loading open(unit=1,status='old',file= +'nhi_gal.dat') write(6,*) '....Loading file: nhi_gal.dat' do i__=1,44 read(1,'(A31)')ar__1 read(ar__1,2) + WCSproj_1(i__),b_RVel_2(i__),B_RVel_3(i__),FileName_1(i__) c ..............Just test output........... write(6,2) + WCSproj_1(i__),b_RVel_2(i__),B_RVel_3(i__),FileName_1(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'allsky_eq.dat' ! List of all-sky FITS data cubes, Equatorial C Format for file interpretation 3 format(A3,1X,A8) C Effective file loading open(unit=1,status='old',file= +'allsky_eq.dat') write(6,*) '....Loading file: allsky_eq.dat' do i__=1,4 read(1,'(A12)')ar__2 read(ar__2,3)WCSproj_2(i__),FileName_2(i__) c ..............Just test output........... write(6,3)WCSproj_2(i__),FileName_2(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'allsky_gal.dat' ! List of all-sky FITS data cubes, Galactic C Format for file interpretation 4 format(A3,1X,A8) C Effective file loading open(unit=1,status='old',file= +'allsky_gal.dat') write(6,*) '....Loading file: allsky_gal.dat' do i__=1,4 read(1,'(A12)')ar__3 read(ar__3,4)WCSproj_3(i__),FileName_3(i__) c ..............Just test output........... write(6,4)WCSproj_3(i__),FileName_3(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'cubes_eq.dat' ! List of 20dx20d FITS data cubes, Equatorial C Format for file interpretation 5 format(F5.1,1X,F5.1,1X,A3,1X,A12) C Effective file loading open(unit=1,status='old',file= +'cubes_eq.dat') write(6,*) '....Loading file: cubes_eq.dat' do i__=1,534 read(1,'(A28)')ar__4 read(ar__4,5) + RAdeg(i__),DEdeg(i__),WCSproj_4(i__),FileName_4(i__) c ..............Just test output........... write(6,5) + RAdeg(i__),DEdeg(i__),WCSproj_4(i__),FileName_4(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'cubes_gal.dat' ! List of 20dx20d FITS data cubes, Galactic C Format for file interpretation 6 format(F5.1,1X,F5.1,1X,A3,1X,A12) C Effective file loading open(unit=1,status='old',file= +'cubes_gal.dat') write(6,*) '....Loading file: cubes_gal.dat' do i__=1,534 read(1,'(A28)')ar__5 read(ar__5,6) + GLON(i__),GLAT(i__),WCSproj_5(i__),FileName_5(i__) c ..............Just test output........... write(6,6) + GLON(i__),GLAT(i__),WCSproj_5(i__),FileName_5(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'nhi_hpx.dat' ! HI column densities C Format for file interpretation 7 format(I8,1X,F9.5,1X,F9.5,1X,F9.5,1X,F9.5,1X,E22.15) C Effective file loading open(unit=1,status='old',file= +'nhi_hpx.dat') write(6,*) '....Loading file: nhi_hpx.dat' do i__=1,12582912 read(1,'(A71)')ar__6 read(ar__6,7) + HPX(i__),RAdeg_1(i__),DEdeg_1(i__),GLON_1(i__),GLAT_1(i__), + NHI(i__) c ..............Just test output........... write(6,7) + HPX(i__),RAdeg_1(i__),DEdeg_1(i__),GLON_1(i__),GLAT_1(i__), + NHI(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'hpxspec.dat' ! List of HPX-grid FITS binary tables C Format for file interpretation 8 format(F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,I3,1X,A12) C Effective file loading open(unit=1,status='old',file= +'hpxspec.dat') write(6,*) '....Loading file: hpxspec.dat' do i__=1,192 read(1,'(A52)')ar__7 read(ar__7,8) + GLON_2(i__),GLAT_2(i__),RAdeg_2(i__),DEdeg_2(i__),HPXidx(i__), + FileName_6(i__) c ..............Just test output........... write(6,8) + GLON_2(i__),GLAT_2(i__),RAdeg_2(i__),DEdeg_2(i__),HPXidx(i__), + FileName_6(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end