FORTRAN Generation
(/./ftp/cats/J/ApJ/814/133)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/814/133 into FORTRAN code for reading data files line by line.

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 2024-Mar-28
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/ApJ/814/133     Herschel Galactic plane survey of [NII]     (Goldsmith+, 2015)
*================================================================================
*Herschel Galactic plane survey of [NII] fine structure emission.
*    Goldsmith P.F., Yildiz U.A., Langer W.D., Pineda J.L.
*   <Astrophys. J., 814, 133 (2015)>
*   =2015ApJ...814..133G    (SIMBAD/NED BibCode)
C=============================================================================

C  Internal variables

      integer*4 i__

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table2.dat'	! Observed PACS Integrated Intensities for
                                [N II] 122 and 205 {mu}m Transitions are
                                Tabulated Together with the Inferred Electron
                                Volume Densities and N^+^ Column Densities

      integer*4 nr__
      parameter (nr__=70)	! Number of records
      character*91 ar__   	! Full-size record

C  Position composed of: GLON GLAT
      character*10  ID          ! Line of sight (LOS) label (GLLL.l+B.b)
      real*8        GLON        ! (deg) Galactic longitude
      real*4        GLAT        ! (deg) Galactic latitude
      integer*8     OBSID       ! Herschel OBSID
      integer*4     R           ! [0/3]? Detection status (1)
      character*1   l_NII122    ! [<] Limit flag on NII122
      real*4        NII122      ! (10-8W/m2/sr) ? [N II] 122 {mu}m intensity averaged over
*                                      the 25 PACS spaxels
      real*4        e_NII122    ! (10-8W/m2/sr) ? Uncertainty in NII122
      real*4        RMS         ! (10-8) ? RMS variation among the 25 spaxels at
*                                      122 {mu}m wavelength
      integer*4     R1          ! [0/3]? Detection status (1)
      character*1   l_NII205    ! [<] Limit flag on NII205
      real*4        NII205      ! (10-8W/m2/sr) ? [N II] 205 {mu}m intensity averaged over
*                                      the 25 PACS spaxels
      real*4        e_NII205    ! (10-8W/m2/sr) ? Uncertainty in NII205
      real*4        n_e         ! (cm-3) ? Electron density
      real*4        e_n_e       ! (cm-3) ? Uncertainty in n(e)
      real*4        N_N_        ! (10+16cm-2) ? N^+^ column density
      real*4        e_N_N_      ! (10+16cm-2) ? Uncertainty in N(N+)
*Note (1): Detection status as follows:
*  0 = Not observed;
*  1 = Non-detection;
*  2 = Detection in individual spaxel or spaxels;
*  3 = Only the average of all spaxels provides a detection.
*  For the non-detections, we give the 3{sigma} upper limits on the intensities.

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table3.dat'	! PACS Results

      integer*4 nr__1
      parameter (nr__1=79)	! Number of records
      character*91 ar__1  	! Full-size record

C  Position composed of: GLON GLAT
      character*10  ID_1        ! Line of sight (LOS) label (GLLL.l+B.b)
      real*8        GLON_1      ! (deg) Galactic longitude
      real*4        GLAT_1      ! (deg) Galactic latitude
      integer*8     OBSID_1     ! Herschel OBSID
      integer*4     R_1         ! [0/3]? Detection status (1)
      character*1   l_NII122_1  ! [<] Limit flag on NII122
      real*4        NII122_1    ! (10-8W/m2/sr) ? [N II] 122 {mu}m intensity averaged over
*                                      the 25 PACS spaxels
      real*4        e_NII122_1  ! (10-8W/m2/sr) ? Uncertainty in NII122
      real*4        RMS_1       ! (10-8) ? RMS variation among the 25 spaxels at
*                                      122 {mu}m wavelength
      integer*4     R1_1        ! [0/3]? Detection status (1)
      character*1   l_NII205_1  ! [<] Limit flag on NII205
      real*4        NII205_1    ! (10-8W/m2/sr) ? [N II] 205 {mu}m intensity averaged over
*                                      the 25 PACS spaxels
      real*4        e_NII205_1  ! (10-8W/m2/sr) ? Uncertainty in NII205
      real*4        n_e_1       ! (cm-3) ? Electron density
      real*4        e_n_e_1     ! (cm-3) ? Uncertainty in n(e)
      real*4        N_N__1      ! (10+16cm-2) ? N^+^ column density
      real*4        e_N_N__1    ! (10+16cm-2) ? Uncertainty in N(N+)
*Note (1): Detection status as follows:
*  0 = Not observed;
*  1 = Non-detection;
*  2 = Detection in individual spaxel or spaxels;
*  3 = Only the average of all spaxels provides a detection.
*  For the non-detections, we give the 3{sigma} upper limits on the intensities.

C=============================================================================

C  Loading file 'table2.dat'	! Observed PACS Integrated Intensities for
*                                [N II] 122 and 205 {mu}m Transitions are
*                                Tabulated Together with the Inferred Electron
*                                Volume Densities and N^+^ Column Densities

C  Format for file interpretation

    1 format(
     +  A10,1X,F8.4,1X,F3.1,1X,I10,1X,I1,1X,A1,F5.2,1X,F5.3,1X,F5.2,
     +  1X,I1,1X,A1,F5.2,1X,F5.3,1X,F5.1,1X,F4.1,1X,F4.1,1X,F4.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,70
        read(1,'(A91)')ar__
        read(ar__,1)
     +  ID,GLON,GLAT,OBSID,R,l_NII122,NII122,e_NII122,RMS,R1,l_NII205,
     +  NII205,e_NII205,n_e,e_n_e,N_N_,e_N_N_
        if(ar__(36:36) .EQ. '') R = iNULL__
        if(ar__(39:43) .EQ. '') NII122 = rNULL__
        if(ar__(45:49) .EQ. '') e_NII122 = rNULL__
        if(ar__(51:55) .EQ. '') RMS = rNULL__
        if(ar__(57:57) .EQ. '') R1 = iNULL__
        if(ar__(60:64) .EQ. '') NII205 = rNULL__
        if(ar__(66:70) .EQ. '') e_NII205 = rNULL__
        if(ar__(72:76) .EQ. '') n_e = rNULL__
        if(ar__(78:81) .EQ. '') e_n_e = rNULL__
        if(ar__(83:86) .EQ. '') N_N_ = rNULL__
        if(ar__(88:91) .EQ. '') e_N_N_ = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  ID,GLON,GLAT,OBSID,R,l_NII122,NII122,e_NII122,RMS,R1,l_NII205,
     +  NII205,e_NII205,n_e,e_n_e,N_N_,e_N_N_
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table3.dat'	! PACS Results

C  Format for file interpretation

    2 format(
     +  A10,1X,F8.4,1X,F3.1,1X,I10,1X,I1,1X,A1,F5.2,1X,F5.3,1X,F5.2,
     +  1X,I1,1X,A1,F5.2,1X,F5.3,1X,F5.1,1X,F4.1,1X,F4.1,1X,F4.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,79
        read(1,'(A91)')ar__1
        read(ar__1,2)
     +  ID_1,GLON_1,GLAT_1,OBSID_1,R_1,l_NII122_1,NII122_1,e_NII122_1,
     +  RMS_1,R1_1,l_NII205_1,NII205_1,e_NII205_1,n_e_1,e_n_e_1,
     +  N_N__1,e_N_N__1
        if(ar__1(36:36) .EQ. '') R_1 = iNULL__
        if(ar__1(39:43) .EQ. '') NII122_1 = rNULL__
        if(ar__1(45:49) .EQ. '') e_NII122_1 = rNULL__
        if(ar__1(51:55) .EQ. '') RMS_1 = rNULL__
        if(ar__1(57:57) .EQ. '') R1_1 = iNULL__
        if(ar__1(60:64) .EQ. '') NII205_1 = rNULL__
        if(ar__1(66:70) .EQ. '') e_NII205_1 = rNULL__
        if(ar__1(72:76) .EQ. '') n_e_1 = rNULL__
        if(ar__1(78:81) .EQ. '') e_n_e_1 = rNULL__
        if(ar__1(83:86) .EQ. '') N_N__1 = rNULL__
        if(ar__1(88:91) .EQ. '') e_N_N__1 = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  ID_1,GLON_1,GLAT_1,OBSID_1,R_1,l_NII122_1,NII122_1,e_NII122_1,
     +  RMS_1,R1_1,l_NII205_1,NII205_1,e_NII205_1,n_e_1,e_n_e_1,
     +  N_N__1,e_N_N__1
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================
      stop
      end