FORTRAN Generation
(/./ftp/cats/J/AZh/79/610)

Conversion of standardized ReadMe file for file /./ftp/cats/J/AZh/79/610 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/AZh/79/610  Methyl acetylene observations of warm clouds       (Alakoz+, 2002)
*================================================================================
*Parameters of warm molecular clouds from methyl acetylene observations
*    Alakoz A.V., Kalenskii S.V., Promislov V.G., Johansson L.E.B., Winnberg A.
*   <Astron. Zh. 79, 610 (2002)>
*   =2002AZh....79..610A
*   =2002ARep...46..551A
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table2.dat'	! Gaussian parameters for detected CH_3_CCH lines
                            with 1{sigma} errors

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

C  1950 position composed of: RAh RAm RAs DE- DEd DEm DEs
      real*8        RAdeg       ! (deg) Right Ascension 1950
      real*8        DEdeg       ! (deg)     Declination 1950
C  ---------------------------------- ! (position vector(s) in degrees)

      character*11  Name        ! Name of the source
      integer*4     RAh         ! (h) Right ascension (1950)
      integer*4     RAm         ! (min) Right ascension (1950)
      real*4        RAs         ! (s) Right ascension (1950)
      character*1   DE_         ! Declination sign (1950)
      integer*4     DEd         ! (deg) Declination (1950)
      integer*4     DEm         ! (arcmin) Declination (1950)
      real*4        DEs         ! (arcsec) Declination (1950)
      real*4        I102K0      ! (K.km/s) CH3CCH intensity integrated over the whole
*                                    line at 102GHz, 6_0-5_0 transition (K=0) (1)
      real*4        e_I102K0    ! (K.km/s) Uncertainty in I102K0
      real*4        I102K1      ! (K.km/s) CH3CCH intensity integrated over the whole
*                                    line at 102GHz, 6_1-5_1 transition (K=1) (1)
      real*4        e_I102K1    ! (K.km/s) Uncertainty in I102K1
      character*1   l_I102K2    ! [<] Limit symbol on I102K2
      real*4        I102K2      ! (K.km/s) ? CH3CCH intensity integrated over the whole
*                                    line at 102GHz, 6_2-5_2 transition (K=2) (1)
      real*4        e_I102K2    ! (K.km/s) ? Uncertainty in I102K2
      character*1   l_I102K3    ! [<] Limit symbol on I102K3
      real*4        I102K3      ! (K.km/s) ? CH3CCH intensity integrated over the whole
*                                    line at 102GHz, 6_3-5_3 transition (K=3) (1)
      real*4        e_I102K3    ! (K.km/s) ? Uncertainty in I102K3
      character*1   l_I102K4    ! [<] Limit symbol on I102K4
      real*4        I102K4      ! (K.km/s) ? CH3CCH intensity integrated over the whole
*                                    line at 102GHz, 6_4-5_4 transition (K=4) (1)
      real*4        e_I102K4    ! (K.km/s) ? Uncertainty in I102K4
      real*8        VLSR        ! (km/s) LSR velocity of 102GHz line
      real*4        e_VLSR      ! (km/s) Uncertainty in VLSR
      real*4        FWHM        ! (km/s) FWHM of the 102GHz line
      real*4        e_FWHM      ! (km/s) Uncertainty in FWHM
      character*1   n_I85K0     ! [*] *: 85GHz line isn't observed
      real*4        I85K0       ! (K.km/s) ? CH3CCH intensity integrated over the whole
*                                    line at 85GHz, 5_0-4_0 transition (K=0) (1)
      real*4        e_I85K0     ! (K.km/s) ? Uncertainty in I85K0
      real*4        I85K1       ! (K.km/s) ? CH3CCH intensity integrated over the whole
*                                    line at 85GHz, 5_1-4_1 transition (K=1) (1)
      real*4        e_I85K1     ! (K.km/s) ? Uncertainty in I85K1
      character*1   l_I85K2     ! [<] Limit symbol on I85K2
      real*4        I85K2       ! (K.km/s) ? CH3CCH intensity integrated over the whole
*                                    line at 85GHz, 5_2-4_2 transition (K=2) (1)
      real*4        e_I85K2     ! (K.km/s) ? Uncertainty in I85K2
      character*1   l_I85K3     ! [<] Limit symbol on I85K3
      real*4        I85K3       ! (K.km/s) ? CH3CCH intensity integrated over the whole
*                                    line at 85GHz, 5_3-4_3 transition (K=3) (1)
      real*4        e_I85K3     ! (K.km/s) ? Uncertainty in I85K3
      character*1   l_I85K4     ! [<] Limit symbol on I85K4
      real*4        I85K4       ! (K.km/s) ? CH3CCH intensity integrated over the whole
*                                    line at 85GHz, 5_4-4_4 transition (K=4) (1)
      real*8        Vlsr85      ! (km/s) ? LSR velocity of 85GHz line
      real*4        e_Vlsr85    ! (km/s) ? Uncertainty in Vlsr85
      real*4        FWHM85      ! (km/s) ? FWHM of the 85GHz line
      real*4        e_FWHM85    ! (km/s) ? Uncertainty in FWHM85
*Note (1): Intensitiies in K.km/s, T*A scale

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

C  Declarations for 'table3.dat'	! Rotational temperatures and column densities
                            derived from the methyl acetylene observations

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

      character*11  Name_1      ! Name of the source
      character*1   l_Trot102   ! [<] Limit symbol on Trot102
      integer*4     Trot102     ! (K) Rotational temperature at 102GHz
      real*4        e_Trot102   ! (K) ? Uncertainty in Trot102
      character*1   l_Nch3ccn102 ! [<] Limit symbol on Nch3ccn102
      real*4        Nch3ccn102  ! (10+14cm-2) ? Column density of CH3CCH at 102GHz
      real*4        e_Nch3ccn102 ! (10+14cm-2) ? Uncertainty in Nch3ccn102
      character*1   l_Trot85    ! [<] Limit symbol on Trot85
      real*4        Trot85      ! (K) ? Rotational temperature at 85GHz
      real*4        e_Trot85    ! (K) ? Uncertainty in Trot85
      character*1   l_Nch3ccn85 ! [<] Limit symbol on Nch3ccn85
      real*4        Nch3ccn85   ! (10+14cm-2) ? Column density of CH3CCH at 85GHz
      real*4        e_Nch3ccn85 ! (10+14cm-2) ? Uncertainty in Nch3ccn85
      real*4        N_CH3CCH__N_H2 ! (10-9) ? CH3CCH abundance

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

C  Loading file 'table2.dat'	! Gaussian parameters for detected CH_3_CCH lines
*                            with 1{sigma} errors

C  Format for file interpretation

    1 format(
     +  A11,1X,I2,1X,I2,1X,F5.2,2X,A1,I2,1X,I2,1X,F4.1,2X,F4.2,1X,
     +  F4.2,2X,F4.2,1X,F4.2,2X,A1,F4.2,1X,F4.2,2X,A1,F4.2,1X,F4.2,2X,
     +  A1,F4.2,1X,F4.2,2X,F7.3,1X,F4.2,2X,F5.2,1X,F4.2,1X,A1,3X,F4.2,
     +  1X,F4.2,2X,F4.2,1X,F4.2,1X,A1,F4.2,1X,F4.2,1X,A1,F4.2,1X,F4.2,
     +  2X,A1,F4.2,4X,F7.3,1X,F4.2,1X,F5.2,1X,F4.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,44
        read(1,'(A201)')ar__
        read(ar__,1)
     +  Name,RAh,RAm,RAs,DE_,DEd,DEm,DEs,I102K0,e_I102K0,I102K1,
     +  e_I102K1,l_I102K2,I102K2,e_I102K2,l_I102K3,I102K3,e_I102K3,
     +  l_I102K4,I102K4,e_I102K4,VLSR,e_VLSR,FWHM,e_FWHM,n_I85K0,
     +  I85K0,e_I85K0,I85K1,e_I85K1,l_I85K2,I85K2,e_I85K2,l_I85K3,
     +  I85K3,e_I85K3,l_I85K4,I85K4,Vlsr85,e_Vlsr85,FWHM85,e_FWHM85
        if(ar__(62:65) .EQ. '') I102K2 = rNULL__
        if(ar__(67:70) .EQ. '') e_I102K2 = rNULL__
        if(ar__(74:77) .EQ. '') I102K3 = rNULL__
        if(ar__(79:82) .EQ. '') e_I102K3 = rNULL__
        if(ar__(86:89) .EQ. '') I102K4 = rNULL__
        if(ar__(91:94) .EQ. '') e_I102K4 = rNULL__
        if(ar__(126:129) .EQ. '') I85K0 = rNULL__
        if(ar__(131:134) .EQ. '') e_I85K0 = rNULL__
        if(ar__(137:140) .EQ. '') I85K1 = rNULL__
        if(ar__(142:145) .EQ. '') e_I85K1 = rNULL__
        if(ar__(148:151) .EQ. '') I85K2 = rNULL__
        if(ar__(153:156) .EQ. '') e_I85K2 = rNULL__
        if(ar__(159:162) .EQ. '') I85K3 = rNULL__
        if(ar__(164:167) .EQ. '') e_I85K3 = rNULL__
        if(ar__(171:174) .EQ. '') I85K4 = rNULL__
        if(ar__(179:185) .EQ. '') Vlsr85 = rNULL__
        if(ar__(187:190) .EQ. '') e_Vlsr85 = rNULL__
        if(ar__(192:196) .EQ. '') FWHM85 = rNULL__
        if(ar__(198:201) .EQ. '') e_FWHM85 = rNULL__
        RAdeg = rNULL__
        DEdeg = rNULL__
c  Derive coordinates RAdeg and DEdeg from input data
c  (RAdeg and DEdeg are set to rNULL__ when unknown)
        if(RAh .GT. -180) RAdeg=RAh*15.
        if(RAm .GT. -180) RAdeg=RAdeg+RAm/4.
        if(RAs .GT. -180) RAdeg=RAdeg+RAs/240.
        if(DEd .GE. 0) DEdeg=DEd
        if(DEm .GE. 0) DEdeg=DEdeg+DEm/60.
        if(DEs .GE. 0) DEdeg=DEdeg+DEs/3600.
        if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg
c    ..............Just test output...........
        write(6,1)
     +  Name,RAh,RAm,RAs,DE_,DEd,DEm,DEs,I102K0,e_I102K0,I102K1,
     +  e_I102K1,l_I102K2,I102K2,e_I102K2,l_I102K3,I102K3,e_I102K3,
     +  l_I102K4,I102K4,e_I102K4,VLSR,e_VLSR,FWHM,e_FWHM,n_I85K0,
     +  I85K0,e_I85K0,I85K1,e_I85K1,l_I85K2,I85K2,e_I85K2,l_I85K3,
     +  I85K3,e_I85K3,l_I85K4,I85K4,Vlsr85,e_Vlsr85,FWHM85,e_FWHM85
        write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table3.dat'	! Rotational temperatures and column densities
*                            derived from the methyl acetylene observations

C  Format for file interpretation

    2 format(
     +  A11,5X,A1,I2,2X,F4.1,1X,A1,F4.1,1X,F3.1,3X,A1,F4.1,1X,F4.1,2X,
     +  A1,F4.1,1X,F3.1,2X,F4.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,40
        read(1,'(A65)')ar__1
        read(ar__1,2)
     +  Name_1,l_Trot102,Trot102,e_Trot102,l_Nch3ccn102,Nch3ccn102,
     +  e_Nch3ccn102,l_Trot85,Trot85,e_Trot85,l_Nch3ccn85,Nch3ccn85,
     +  e_Nch3ccn85,N_CH3CCH__N_H2
        if(ar__1(22:25) .EQ. '') e_Trot102 = rNULL__
        if(ar__1(28:31) .EQ. '') Nch3ccn102 = rNULL__
        if(ar__1(33:35) .EQ. '') e_Nch3ccn102 = rNULL__
        if(ar__1(40:43) .EQ. '') Trot85 = rNULL__
        if(ar__1(45:48) .EQ. '') e_Trot85 = rNULL__
        if(ar__1(52:55) .EQ. '') Nch3ccn85 = rNULL__
        if(ar__1(57:59) .EQ. '') e_Nch3ccn85 = rNULL__
        if(ar__1(62:65) .EQ. '') N_CH3CCH__N_H2 = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  Name_1,l_Trot102,Trot102,e_Trot102,l_Nch3ccn102,Nch3ccn102,
     +  e_Nch3ccn102,l_Trot85,Trot85,e_Trot85,l_Nch3ccn85,Nch3ccn85,
     +  e_Nch3ccn85,N_CH3CCH__N_H2
c    .......End.of.Just test output...........
      end do
      close(1)

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