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-Jul-18
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