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 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 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/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    (nr__) ! Species
      character*17  Trans      (nr__) ! Observed line transition
      real*4        dV         (nr__) ! (km/s) ? Linewidth
      real*4        e_dV       (nr__) ! (km/s) ? rms uncertainty on dV
      character*1   u_dV       (nr__) ! [?] Uncertainty flag on dV
      real*4        Tmb        (nr__) ! (K) ? Main beam temperature
      character*1   u_Tmb      (nr__) ! [?] Uncertainty flag on Tmb
      integer*4     e_Tmb      (nr__) ! (mK) ? rms uncertainty on Tmb
      character*1   l_Flux     (nr__) ! Limit flag on Flux
      real*4        Flux       (nr__) ! (K.km/s) ? Integrated flux
      real*4        e_Flux     (nr__) ! (K.km/s) ? rms uncertainty on Flux
      character*1   u_e_Flux   (nr__) ! [?] Uncertainty flag on Flux
      real*4        dV3s       (nr__) ! (km/s) ? linewidth at 3 sigma
      character*1   u_dV3s     (nr__) ! [?] Uncertainty flag on dV3s
      real*4        dV1_2      (nr__) ! (km/s) ? half-power linewidth
      character*1   u_dV1_2    (nr__) ! [?] Uncertainty flag on dV1/2
      character*1   Flag       (nr__) ! [abcdfg] Flag on profile (1)
      character*14  Com        (nr__) ! Comments
      character*13  FileName   (nr__) ! 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  (nr__1) ! Species
      character*17  Trans_1    (nr__1) ! Observed line transition
      real*4        dV_1       (nr__1) ! (km/s) ? Linewidth
      real*4        e_dV_1     (nr__1) ! (km/s) ? rms uncertainty on dV
      character*1   u_dV_1     (nr__1) ! [?] Uncertainty flag on dV
      real*4        Tmb_1      (nr__1) ! (K) ? Main beam temperature
      character*1   u_Tmb_1    (nr__1) ! [?] Uncertainty flag on Tmb
      integer*4     e_Tmb_1    (nr__1) ! (mK) ? rms uncertainty on Tmb
      character*1   l_Flux_1   (nr__1) ! Limit flag on Flux
      real*4        Flux_1     (nr__1) ! (K.km/s) ? Integrated flux
      real*4        e_Flux_1   (nr__1) ! (K.km/s) ? rms uncertainty on Flux
      character*1   u_e_Flux_1 (nr__1) ! [?] Uncertainty flag on Flux
      real*4        dV3s_1     (nr__1) ! (km/s) ? linewidth at 3 sigma
      character*1   u_dV3s_1   (nr__1) ! [?] Uncertainty flag on dV3s
      real*4        dV1_2_1    (nr__1) ! (km/s) ? half-power linewidth
      character*1   u_dV1_2_1  (nr__1) ! [?] Uncertainty flag on dV1/2
      character*1   Flag_1     (nr__1) ! [abcdfg] Flag on profile (1)
      character*14  Com_1      (nr__1) ! Comments
      character*13  FileName_1 (nr__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  (nr__2) ! Species
      character*17  Trans_2    (nr__2) ! Observed line transition
      real*4        dV_2       (nr__2) ! (km/s) ? Linewidth
      real*4        e_dV_2     (nr__2) ! (km/s) ? rms uncertainty on dV
      character*1   u_dV_2     (nr__2) ! [?] Uncertainty flag on dV
      real*4        Tmb_2      (nr__2) ! (K) ? Main beam temperature
      character*1   u_Tmb_2    (nr__2) ! [?] Uncertainty flag on Tmb
      integer*4     e_Tmb_2    (nr__2) ! (mK) ? rms uncertainty on Tmb
      character*1   l_Flux_2   (nr__2) ! Limit flag on Flux
      real*4        Flux_2     (nr__2) ! (K.km/s) ? Integrated flux
      real*4        e_Flux_2   (nr__2) ! (K.km/s) ? rms uncertainty on Flux
      character*1   u_e_Flux_2 (nr__2) ! [?] Uncertainty flag on Flux
      real*4        dV3s_2     (nr__2) ! (km/s) ? linewidth at 3 sigma
      character*1   u_dV3s_2   (nr__2) ! [?] Uncertainty flag on dV3s
      real*4        dV1_2_2    (nr__2) ! (km/s) ? half-power linewidth
      character*1   u_dV1_2_2  (nr__2) ! [?] Uncertainty flag on dV1/2
      character*1   Flag_2     (nr__2) ! [abcdfg] Flag on profile (1)
      character*14  Com_2      (nr__2) ! Comments
      character*13  FileName_2 (nr__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  (nr__3) ! Species
      character*17  Trans_3    (nr__3) ! Observed line transition
      real*4        dV_3       (nr__3) ! (km/s) ? Linewidth
      real*4        e_dV_3     (nr__3) ! (km/s) ? rms uncertainty on dV
      character*1   u_dV_3     (nr__3) ! [?] Uncertainty flag on dV
      real*4        Tmb_3      (nr__3) ! (K) ? Main beam temperature
      character*1   u_Tmb_3    (nr__3) ! [?] Uncertainty flag on Tmb
      integer*4     e_Tmb_3    (nr__3) ! (mK) ? rms uncertainty on Tmb
      character*1   l_Flux_3   (nr__3) ! Limit flag on Flux
      real*4        Flux_3     (nr__3) ! (K.km/s) ? Integrated flux
      real*4        e_Flux_3   (nr__3) ! (K.km/s) ? rms uncertainty on Flux
      character*1   u_e_Flux_3 (nr__3) ! [?] Uncertainty flag on Flux
      real*4        dV3s_3     (nr__3) ! (km/s) ? linewidth at 3 sigma
      character*1   u_dV3s_3   (nr__3) ! [?] Uncertainty flag on dV3s
      real*4        dV1_2_3    (nr__3) ! (km/s) ? half-power linewidth
      character*1   u_dV1_2_3  (nr__3) ! [?] Uncertainty flag on dV1/2
      character*1   Flag_3     (nr__3) ! [abcdfg] Flag on profile (1)
      character*14  Com_3      (nr__3) ! Comments
      character*13  FileName_3 (nr__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  (nr__4) ! Species
      character*17  Trans_4    (nr__4) ! Observed line transition
      real*4        Opac1      (nr__4) ! ? Line opacity for source W43MM1
      character*2   n_Opac1    (nr__4) ! [ND] ND for no data
      real*4        Opac2      (nr__4) ! Line opacity for source IRAS18264-1152
      real*4        Opac3      (nr__4) ! ? Line opacity for source IRAS05358+3543
      character*2   n_Opac3    (nr__4) ! [ND] ND for no data
      real*4        Opac4      (nr__4) ! ? Line opacity for source IRAS18162-2048
      character*2   n_Opac4    (nr__4) ! [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(i__),Trans(i__),dV(i__),e_dV(i__),u_dV(i__),Tmb(i__),
     +  u_Tmb(i__),e_Tmb(i__),l_Flux(i__),Flux(i__),e_Flux(i__),
     +  u_e_Flux(i__),dV3s(i__),u_dV3s(i__),dV1_2(i__),u_dV1_2(i__),
     +  Flag(i__),Com(i__),FileName(i__)
        if(ar__(30:33) .EQ. '') dV(i__) = rNULL__
        if(ar__(35:38) .EQ. '') e_dV(i__) = rNULL__
        if(ar__(42:45) .EQ. '') Tmb(i__) = rNULL__
        if(ar__(48:50) .EQ. '') e_Tmb(i__) = iNULL__
        if(ar__(54:59) .EQ. '') Flux(i__) = rNULL__
        if(ar__(61:64) .EQ. '') e_Flux(i__) = rNULL__
        if(ar__(69:73) .EQ. '') dV3s(i__) = rNULL__
        if(ar__(77:80) .EQ. '') dV1_2(i__) = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  Species(i__),Trans(i__),dV(i__),e_dV(i__),u_dV(i__),Tmb(i__),
     +  u_Tmb(i__),e_Tmb(i__),l_Flux(i__),Flux(i__),e_Flux(i__),
     +  u_e_Flux(i__),dV3s(i__),u_dV3s(i__),dV1_2(i__),u_dV1_2(i__),
     +  Flag(i__),Com(i__),FileName(i__)
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(i__),Trans_1(i__),dV_1(i__),e_dV_1(i__),u_dV_1(i__),
     +  Tmb_1(i__),u_Tmb_1(i__),e_Tmb_1(i__),l_Flux_1(i__),
     +  Flux_1(i__),e_Flux_1(i__),u_e_Flux_1(i__),dV3s_1(i__),
     +  u_dV3s_1(i__),dV1_2_1(i__),u_dV1_2_1(i__),Flag_1(i__),
     +  Com_1(i__),FileName_1(i__)
        if(ar__1(30:33) .EQ. '') dV_1(i__) = rNULL__
        if(ar__1(35:38) .EQ. '') e_dV_1(i__) = rNULL__
        if(ar__1(42:45) .EQ. '') Tmb_1(i__) = rNULL__
        if(ar__1(48:50) .EQ. '') e_Tmb_1(i__) = iNULL__
        if(ar__1(54:59) .EQ. '') Flux_1(i__) = rNULL__
        if(ar__1(61:64) .EQ. '') e_Flux_1(i__) = rNULL__
        if(ar__1(69:73) .EQ. '') dV3s_1(i__) = rNULL__
        if(ar__1(77:80) .EQ. '') dV1_2_1(i__) = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  Species_1(i__),Trans_1(i__),dV_1(i__),e_dV_1(i__),u_dV_1(i__),
     +  Tmb_1(i__),u_Tmb_1(i__),e_Tmb_1(i__),l_Flux_1(i__),
     +  Flux_1(i__),e_Flux_1(i__),u_e_Flux_1(i__),dV3s_1(i__),
     +  u_dV3s_1(i__),dV1_2_1(i__),u_dV1_2_1(i__),Flag_1(i__),
     +  Com_1(i__),FileName_1(i__)
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(i__),Trans_2(i__),dV_2(i__),e_dV_2(i__),u_dV_2(i__),
     +  Tmb_2(i__),u_Tmb_2(i__),e_Tmb_2(i__),l_Flux_2(i__),
     +  Flux_2(i__),e_Flux_2(i__),u_e_Flux_2(i__),dV3s_2(i__),
     +  u_dV3s_2(i__),dV1_2_2(i__),u_dV1_2_2(i__),Flag_2(i__),
     +  Com_2(i__),FileName_2(i__)
        if(ar__2(30:33) .EQ. '') dV_2(i__) = rNULL__
        if(ar__2(35:38) .EQ. '') e_dV_2(i__) = rNULL__
        if(ar__2(42:45) .EQ. '') Tmb_2(i__) = rNULL__
        if(ar__2(48:50) .EQ. '') e_Tmb_2(i__) = iNULL__
        if(ar__2(54:59) .EQ. '') Flux_2(i__) = rNULL__
        if(ar__2(61:64) .EQ. '') e_Flux_2(i__) = rNULL__
        if(ar__2(69:73) .EQ. '') dV3s_2(i__) = rNULL__
        if(ar__2(77:80) .EQ. '') dV1_2_2(i__) = rNULL__
c    ..............Just test output...........
        write(6,3)
     +  Species_2(i__),Trans_2(i__),dV_2(i__),e_dV_2(i__),u_dV_2(i__),
     +  Tmb_2(i__),u_Tmb_2(i__),e_Tmb_2(i__),l_Flux_2(i__),
     +  Flux_2(i__),e_Flux_2(i__),u_e_Flux_2(i__),dV3s_2(i__),
     +  u_dV3s_2(i__),dV1_2_2(i__),u_dV1_2_2(i__),Flag_2(i__),
     +  Com_2(i__),FileName_2(i__)
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(i__),Trans_3(i__),dV_3(i__),e_dV_3(i__),u_dV_3(i__),
     +  Tmb_3(i__),u_Tmb_3(i__),e_Tmb_3(i__),l_Flux_3(i__),
     +  Flux_3(i__),e_Flux_3(i__),u_e_Flux_3(i__),dV3s_3(i__),
     +  u_dV3s_3(i__),dV1_2_3(i__),u_dV1_2_3(i__),Flag_3(i__),
     +  Com_3(i__),FileName_3(i__)
        if(ar__3(30:33) .EQ. '') dV_3(i__) = rNULL__
        if(ar__3(35:38) .EQ. '') e_dV_3(i__) = rNULL__
        if(ar__3(42:45) .EQ. '') Tmb_3(i__) = rNULL__
        if(ar__3(48:50) .EQ. '') e_Tmb_3(i__) = iNULL__
        if(ar__3(54:59) .EQ. '') Flux_3(i__) = rNULL__
        if(ar__3(61:64) .EQ. '') e_Flux_3(i__) = rNULL__
        if(ar__3(69:73) .EQ. '') dV3s_3(i__) = rNULL__
        if(ar__3(77:80) .EQ. '') dV1_2_3(i__) = rNULL__
c    ..............Just test output...........
        write(6,4)
     +  Species_3(i__),Trans_3(i__),dV_3(i__),e_dV_3(i__),u_dV_3(i__),
     +  Tmb_3(i__),u_Tmb_3(i__),e_Tmb_3(i__),l_Flux_3(i__),
     +  Flux_3(i__),e_Flux_3(i__),u_e_Flux_3(i__),dV3s_3(i__),
     +  u_dV3s_3(i__),dV1_2_3(i__),u_dV1_2_3(i__),Flag_3(i__),
     +  Com_3(i__),FileName_3(i__)
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(i__),Trans_4(i__),Opac1(i__),n_Opac1(i__),
     +  Opac2(i__),Opac3(i__),n_Opac3(i__),Opac4(i__),n_Opac4(i__)
        if(ar__4(29:34) .EQ. '') Opac1(i__) = rNULL__
        if(ar__4(43:46) .EQ. '') Opac3(i__) = rNULL__
        if(ar__4(50:53) .EQ. '') Opac4(i__) = rNULL__
c    ..............Just test output...........
        write(6,5)
     +  Species_4(i__),Trans_4(i__),Opac1(i__),n_Opac1(i__),
     +  Opac2(i__),Opac3(i__),n_Opac3(i__),Opac4(i__),n_Opac4(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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