FORTRAN Generation
(/./ftp/cats/J/A_A/504/853)

Conversion of standardized ReadMe file for file /./ftp/cats/J/A_A/504/853 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-29
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/504/853      Molecular lines in 5 massive dense cores     (Herpin, 2009)
*================================================================================
*S-bearing molecules in massive dense cores
*     Herpin F., Marseille M., Wakelam V., Bontemps S., Lis D.C.
*    <Astron. Astrophys. 504, 853 (2009)>
*    =2009A&A...504..853H
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table6.dat'	! Observed line emissions parameters for W43MM1

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

      character*9   Species     ! Species
      character*17  Trans       ! Observed line transition
      real*4        dV          ! (km/s) ? Linewidth
      real*4        e_dV        ! (km/s) ? rms uncertainty on dV
      character*1   u_dV        ! [?] Uncertainty flag on dV
      real*4        Tmb         ! (K) ? Main beam temperature
      character*1   u_Tmb       ! [?] Uncertainty flag on Tmb
      integer*4     e_Tmb       ! (mK) ? rms uncertainty on Tmb
      character*1   l_Flux      ! Limit flag on Flux
      real*4        Flux        ! (K.km/s) ? Integrated flux
      real*4        e_Flux      ! (K.km/s) ? rms uncertainty on Flux
      character*1   u_e_Flux    ! [?] Uncertainty flag on Flux
      real*4        dV3s        ! (km/s) ? linewidth at 3 sigma
      character*1   u_dV3s      ! [?] Uncertainty flag on dV3s
      real*4        dV1_2       ! (km/s) ? half-power linewidth
      character*1   u_dV1_2     ! [?] Uncertainty flag on dV1/2
      character*1   Flag        ! [abcdfg] Flag on profile (1)
      character*14  Com         ! Comments
      character*13  FileName    ! Name of the fits file in subdirectory fits
*Note (1): Flags as follows:
*      a = self-absorption at dV_source
*      b = line blend with H_2_CO line from image band
*      c = unidentified line at 34.87km/s (216.716GHz)
*      d = + other line at 29.1km/s (135.78223GHz),
*            Tmb=0.14K, {delta}=1.1km/s, Tmb{delta}=0.17
*      f = detected at 2{sigma} with two peaks at -15.57 and -16.21km/s
*      g = unidentified lines at 216.701026 and 216.71.6GHz (25 and 4km/s)

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

C  Declarations for 'table7.dat'	! Observed line emissions parameters for IRAS18264-1152

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

      character*9   Species_1   ! Species
      character*17  Trans_1     ! Observed line transition
      real*4        dV_1        ! (km/s) ? Linewidth
      real*4        e_dV_1      ! (km/s) ? rms uncertainty on dV
      character*1   u_dV_1      ! [?] Uncertainty flag on dV
      real*4        Tmb_1       ! (K) ? Main beam temperature
      character*1   u_Tmb_1     ! [?] Uncertainty flag on Tmb
      integer*4     e_Tmb_1     ! (mK) ? rms uncertainty on Tmb
      character*1   l_Flux_1    ! Limit flag on Flux
      real*4        Flux_1      ! (K.km/s) ? Integrated flux
      real*4        e_Flux_1    ! (K.km/s) ? rms uncertainty on Flux
      character*1   u_e_Flux_1  ! [?] Uncertainty flag on Flux
      real*4        dV3s_1      ! (km/s) ? linewidth at 3 sigma
      character*1   u_dV3s_1    ! [?] Uncertainty flag on dV3s
      real*4        dV1_2_1     ! (km/s) ? half-power linewidth
      character*1   u_dV1_2_1   ! [?] Uncertainty flag on dV1/2
      character*1   Flag_1      ! [abcdfg] Flag on profile (1)
      character*14  Com_1       ! Comments
      character*13  FileName_1  ! Name of the fits file in subdirectory fits
*Note (1): Flags as follows:
*      a = self-absorption at dV_source
*      b = line blend with H_2_CO line from image band
*      c = unidentified line at 34.87km/s (216.716GHz)
*      d = + other line at 29.1km/s (135.78223GHz),
*            Tmb=0.14K, {delta}=1.1km/s, Tmb{delta}=0.17
*      f = detected at 2{sigma} with two peaks at -15.57 and -16.21km/s
*      g = unidentified lines at 216.701026 and 216.71.6GHz (25 and 4km/s)

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

C  Declarations for 'table8.dat'	! Observed line emissions parameters for IRAS05358+3543

      integer*4 nr__2
      parameter (nr__2=25)	! Number of records
      character*115 ar__2  	! Full-size record

      character*9   Species_2   ! Species
      character*17  Trans_2     ! Observed line transition
      real*4        dV_2        ! (km/s) ? Linewidth
      real*4        e_dV_2      ! (km/s) ? rms uncertainty on dV
      character*1   u_dV_2      ! [?] Uncertainty flag on dV
      real*4        Tmb_2       ! (K) ? Main beam temperature
      character*1   u_Tmb_2     ! [?] Uncertainty flag on Tmb
      integer*4     e_Tmb_2     ! (mK) ? rms uncertainty on Tmb
      character*1   l_Flux_2    ! Limit flag on Flux
      real*4        Flux_2      ! (K.km/s) ? Integrated flux
      real*4        e_Flux_2    ! (K.km/s) ? rms uncertainty on Flux
      character*1   u_e_Flux_2  ! [?] Uncertainty flag on Flux
      real*4        dV3s_2      ! (km/s) ? linewidth at 3 sigma
      character*1   u_dV3s_2    ! [?] Uncertainty flag on dV3s
      real*4        dV1_2_2     ! (km/s) ? half-power linewidth
      character*1   u_dV1_2_2   ! [?] Uncertainty flag on dV1/2
      character*1   Flag_2      ! [abcdfg] Flag on profile (1)
      character*14  Com_2       ! Comments
      character*13  FileName_2  ! Name of the fits file in subdirectory fits
*Note (1): Flags as follows:
*      a = self-absorption at dV_source
*      b = line blend with H_2_CO line from image band
*      c = unidentified line at 34.87km/s (216.716GHz)
*      d = + other line at 29.1km/s (135.78223GHz),
*            Tmb=0.14K, {delta}=1.1km/s, Tmb{delta}=0.17
*      f = detected at 2{sigma} with two peaks at -15.57 and -16.21km/s
*      g = unidentified lines at 216.701026 and 216.71.6GHz (25 and 4km/s)

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

C  Declarations for 'table9.dat'	! Observed line emissions parameters for IRAS18162-2048

      integer*4 nr__3
      parameter (nr__3=25)	! Number of records
      character*115 ar__3  	! Full-size record

      character*9   Species_3   ! Species
      character*17  Trans_3     ! Observed line transition
      real*4        dV_3        ! (km/s) ? Linewidth
      real*4        e_dV_3      ! (km/s) ? rms uncertainty on dV
      character*1   u_dV_3      ! [?] Uncertainty flag on dV
      real*4        Tmb_3       ! (K) ? Main beam temperature
      character*1   u_Tmb_3     ! [?] Uncertainty flag on Tmb
      integer*4     e_Tmb_3     ! (mK) ? rms uncertainty on Tmb
      character*1   l_Flux_3    ! Limit flag on Flux
      real*4        Flux_3      ! (K.km/s) ? Integrated flux
      real*4        e_Flux_3    ! (K.km/s) ? rms uncertainty on Flux
      character*1   u_e_Flux_3  ! [?] Uncertainty flag on Flux
      real*4        dV3s_3      ! (km/s) ? linewidth at 3 sigma
      character*1   u_dV3s_3    ! [?] Uncertainty flag on dV3s
      real*4        dV1_2_3     ! (km/s) ? half-power linewidth
      character*1   u_dV1_2_3   ! [?] Uncertainty flag on dV1/2
      character*1   Flag_3      ! [abcdfg] Flag on profile (1)
      character*14  Com_3       ! Comments
      character*13  FileName_3  ! Name of the fits file in subdirectory fits
*Note (1): Flags as follows:
*      a = self-absorption at dV_source
*      b = line blend with H_2_CO line from image band
*      c = unidentified line at 34.87km/s (216.716GHz)
*      d = + other line at 29.1km/s (135.78223GHz),
*            Tmb=0.14K, {delta}=1.1km/s, Tmb{delta}=0.17
*      f = detected at 2{sigma} with two peaks at -15.57 and -16.21km/s
*      g = unidentified lines at 216.701026 and 216.71.6GHz (25 and 4km/s)

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

C  Declarations for 'table10.dat'	! Opacities derived from modeling the molecular
                            emission of the sources

      integer*4 nr__4
      parameter (nr__4=24)	! Number of records
      character*55 ar__4  	! Full-size record

      character*9   Species_4   ! Species
      character*17  Trans_4     ! Observed line transition
      real*4        Opac1       ! ? Line opacity for source W43MM1
      character*2   n_Opac1     ! [ND] ND for no data
      real*4        Opac2       ! Line opacity for source IRAS18264-1152
      real*4        Opac3       ! ? Line opacity for source IRAS05358+3543
      character*2   n_Opac3     ! [ND] ND for no data
      real*4        Opac4       ! ? Line opacity for source IRAS18162-2048
      character*2   n_Opac4     ! [ND] ND for no data

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

C  Loading file 'table6.dat'	! Observed line emissions parameters for W43MM1

C  Format for file interpretation

    1 format(
     +  A9,1X,A17,2X,F4.2,1X,F4.2,1X,A1,1X,F4.2,A1,1X,I3,2X,A1,F6.3,
     +  1X,F4.2,1X,A1,2X,F5.2,A1,2X,F4.2,A1,A1,5X,A14,1X,A13)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table6.dat')
      write(6,*) '....Loading file: table6.dat'
      do i__=1,25
        read(1,'(A115)')ar__
        read(ar__,1)
     +  Species,Trans,dV,e_dV,u_dV,Tmb,u_Tmb,e_Tmb,l_Flux,Flux,e_Flux,
     +  u_e_Flux,dV3s,u_dV3s,dV1_2,u_dV1_2,Flag,Com,FileName
        if(ar__(30:33) .EQ. '') dV = rNULL__
        if(ar__(35:38) .EQ. '') e_dV = rNULL__
        if(ar__(42:45) .EQ. '') Tmb = rNULL__
        if(ar__(48:50) .EQ. '') e_Tmb = iNULL__
        if(ar__(54:59) .EQ. '') Flux = rNULL__
        if(ar__(61:64) .EQ. '') e_Flux = rNULL__
        if(ar__(69:73) .EQ. '') dV3s = rNULL__
        if(ar__(77:80) .EQ. '') dV1_2 = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  Species,Trans,dV,e_dV,u_dV,Tmb,u_Tmb,e_Tmb,l_Flux,Flux,e_Flux,
     +  u_e_Flux,dV3s,u_dV3s,dV1_2,u_dV1_2,Flag,Com,FileName
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table7.dat'	! Observed line emissions parameters for IRAS18264-1152

C  Format for file interpretation

    2 format(
     +  A9,1X,A17,2X,F4.2,1X,F4.2,1X,A1,1X,F4.2,A1,1X,I3,2X,A1,F6.3,
     +  1X,F4.2,1X,A1,2X,F5.2,A1,2X,F4.2,A1,A1,5X,A14,1X,A13)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table7.dat')
      write(6,*) '....Loading file: table7.dat'
      do i__=1,25
        read(1,'(A115)')ar__1
        read(ar__1,2)
     +  Species_1,Trans_1,dV_1,e_dV_1,u_dV_1,Tmb_1,u_Tmb_1,e_Tmb_1,
     +  l_Flux_1,Flux_1,e_Flux_1,u_e_Flux_1,dV3s_1,u_dV3s_1,dV1_2_1,
     +  u_dV1_2_1,Flag_1,Com_1,FileName_1
        if(ar__1(30:33) .EQ. '') dV_1 = rNULL__
        if(ar__1(35:38) .EQ. '') e_dV_1 = rNULL__
        if(ar__1(42:45) .EQ. '') Tmb_1 = rNULL__
        if(ar__1(48:50) .EQ. '') e_Tmb_1 = iNULL__
        if(ar__1(54:59) .EQ. '') Flux_1 = rNULL__
        if(ar__1(61:64) .EQ. '') e_Flux_1 = rNULL__
        if(ar__1(69:73) .EQ. '') dV3s_1 = rNULL__
        if(ar__1(77:80) .EQ. '') dV1_2_1 = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  Species_1,Trans_1,dV_1,e_dV_1,u_dV_1,Tmb_1,u_Tmb_1,e_Tmb_1,
     +  l_Flux_1,Flux_1,e_Flux_1,u_e_Flux_1,dV3s_1,u_dV3s_1,dV1_2_1,
     +  u_dV1_2_1,Flag_1,Com_1,FileName_1
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table8.dat'	! Observed line emissions parameters for IRAS05358+3543

C  Format for file interpretation

    3 format(
     +  A9,1X,A17,2X,F4.2,1X,F4.2,1X,A1,1X,F4.2,A1,1X,I3,2X,A1,F6.3,
     +  1X,F4.2,1X,A1,2X,F5.2,A1,2X,F4.2,A1,A1,5X,A14,1X,A13)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table8.dat')
      write(6,*) '....Loading file: table8.dat'
      do i__=1,25
        read(1,'(A115)')ar__2
        read(ar__2,3)
     +  Species_2,Trans_2,dV_2,e_dV_2,u_dV_2,Tmb_2,u_Tmb_2,e_Tmb_2,
     +  l_Flux_2,Flux_2,e_Flux_2,u_e_Flux_2,dV3s_2,u_dV3s_2,dV1_2_2,
     +  u_dV1_2_2,Flag_2,Com_2,FileName_2
        if(ar__2(30:33) .EQ. '') dV_2 = rNULL__
        if(ar__2(35:38) .EQ. '') e_dV_2 = rNULL__
        if(ar__2(42:45) .EQ. '') Tmb_2 = rNULL__
        if(ar__2(48:50) .EQ. '') e_Tmb_2 = iNULL__
        if(ar__2(54:59) .EQ. '') Flux_2 = rNULL__
        if(ar__2(61:64) .EQ. '') e_Flux_2 = rNULL__
        if(ar__2(69:73) .EQ. '') dV3s_2 = rNULL__
        if(ar__2(77:80) .EQ. '') dV1_2_2 = rNULL__
c    ..............Just test output...........
        write(6,3)
     +  Species_2,Trans_2,dV_2,e_dV_2,u_dV_2,Tmb_2,u_Tmb_2,e_Tmb_2,
     +  l_Flux_2,Flux_2,e_Flux_2,u_e_Flux_2,dV3s_2,u_dV3s_2,dV1_2_2,
     +  u_dV1_2_2,Flag_2,Com_2,FileName_2
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table9.dat'	! Observed line emissions parameters for IRAS18162-2048

C  Format for file interpretation

    4 format(
     +  A9,1X,A17,2X,F4.2,1X,F4.2,1X,A1,1X,F4.2,A1,1X,I3,2X,A1,F6.3,
     +  1X,F4.2,1X,A1,2X,F5.2,A1,2X,F4.2,A1,A1,5X,A14,1X,A13)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table9.dat')
      write(6,*) '....Loading file: table9.dat'
      do i__=1,25
        read(1,'(A115)')ar__3
        read(ar__3,4)
     +  Species_3,Trans_3,dV_3,e_dV_3,u_dV_3,Tmb_3,u_Tmb_3,e_Tmb_3,
     +  l_Flux_3,Flux_3,e_Flux_3,u_e_Flux_3,dV3s_3,u_dV3s_3,dV1_2_3,
     +  u_dV1_2_3,Flag_3,Com_3,FileName_3
        if(ar__3(30:33) .EQ. '') dV_3 = rNULL__
        if(ar__3(35:38) .EQ. '') e_dV_3 = rNULL__
        if(ar__3(42:45) .EQ. '') Tmb_3 = rNULL__
        if(ar__3(48:50) .EQ. '') e_Tmb_3 = iNULL__
        if(ar__3(54:59) .EQ. '') Flux_3 = rNULL__
        if(ar__3(61:64) .EQ. '') e_Flux_3 = rNULL__
        if(ar__3(69:73) .EQ. '') dV3s_3 = rNULL__
        if(ar__3(77:80) .EQ. '') dV1_2_3 = rNULL__
c    ..............Just test output...........
        write(6,4)
     +  Species_3,Trans_3,dV_3,e_dV_3,u_dV_3,Tmb_3,u_Tmb_3,e_Tmb_3,
     +  l_Flux_3,Flux_3,e_Flux_3,u_e_Flux_3,dV3s_3,u_dV3s_3,dV1_2_3,
     +  u_dV1_2_3,Flag_3,Com_3,FileName_3
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table10.dat'	! Opacities derived from modeling the molecular
*                            emission of the sources

C  Format for file interpretation

    5 format(A9,1X,A17,1X,F6.2,A2,1X,F4.2,1X,F4.2,A2,1X,F4.2,A2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table10.dat')
      write(6,*) '....Loading file: table10.dat'
      do i__=1,24
        read(1,'(A55)')ar__4
        read(ar__4,5)
     +  Species_4,Trans_4,Opac1,n_Opac1,Opac2,Opac3,n_Opac3,Opac4,
     +  n_Opac4
        if(ar__4(29:34) .EQ. '') Opac1 = rNULL__
        if(ar__4(43:46) .EQ. '') Opac3 = rNULL__
        if(ar__4(50:53) .EQ. '') Opac4 = rNULL__
c    ..............Just test output...........
        write(6,5)
     +  Species_4,Trans_4,Opac1,n_Opac1,Opac2,Opac3,n_Opac3,Opac4,
     +  n_Opac4
c    .......End.of.Just test output...........
      end do
      close(1)

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