FORTRAN Generation
(/./ftp/cats/J/A_A/594/A116)

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-Jan-20
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