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