Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/358/257 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/358/257 Molecular clouds SO and CS obs. II. (Nilsson+, 2000) *================================================================================ *SO and CS observations of molecular clouds. *II. Analysis and modelling of the abundance ratios - probing O_2_/CO with SO/CS? * Nilsson A., Hjalmarson A., Bergman P., Millar T.J. * <Astron. Astrophys. 358, 257 (2000)> * =2000A&A...358..257N C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table2.dat' ! Maxima and minima of SO/CS integrated intensity ratios integer*4 nr__ parameter (nr__=22) ! Number of records character*77 ar__ ! Full-size record C J2000.0 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg ! (deg) Right Ascension J2000.0 real*8 DEdeg ! (deg) Declination J2000.0 C ---------------------------------- ! (position vector(s) in degrees) character*15 Source ! Source name character*11 n_Source ! Note (isotopomers or Residual) integer*4 RAh ! (h) Right ascension (J2000.0) integer*4 RAm ! (min) Right ascension (J2000.0) real*4 RAs ! (s) Right ascension (J2000.0) character*1 DE_ ! Declination sign (J2000.0) integer*4 DEd ! (deg) Declination (J2000.0) integer*4 DEm ! (arcmin) Declination (J2000.0) integer*4 DEs ! (arcsec) Declination (J2000.0) integer*4 ORAMax ! (arcsec) Maximum offset in right ascension integer*4 ODEMax ! (arcsec) Maximum offset in declination real*4 RatioMax ! Maximum ratio of SO/CS integer*4 ORAMin ! (arcsec) Minimum offset in right ascension integer*4 ODEMin ! (arcsec) Minimum offset in declination real*4 RatioMin ! Minimum ratio of SO/CS c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table3.dat' ! Integrated intensity ratios (A-E) and estimated abundance ratios (F) at certain positions integer*4 nr__1 parameter (nr__1=92) ! Number of records character*80 ar__1 ! Full-size record character*15 Source_1 ! Source name integer*4 ORA ! (arcsec) Offset in right ascension from source integer*4 ODE ! (arcsec) Offset in declination from source real*4 RatioA ! ? SO(3_2_-2_1_)/^34^SO(3_2_-2_1_) ratio real*4 e_RatioA ! ? rms uncertainty on RatioA (1) real*4 RatioB ! ? CS(2-1)/C^34^S(2-1) ratio real*4 e_RatioB ! ? rms uncertainty on RatioB (1) real*4 RatioC ! ? SO(3_2_-2_1_)/CS(2-1) ratio real*4 e_RatioC ! ? rms uncertainty on RatioC (1) real*4 RatioD ! ? ^34^SO(3_2_-2_1_)/C^34^S(2-1) ratio real*4 e_RatioD ! ? rms uncertainty on RatioD (1) real*4 RatioE ! ? Adopted SO/CS ratio (2) real*4 e_RatioE ! ? rms uncertainty on RatioE (1) character*1 Note ! [abcd] Note (3) real*4 RadioF ! Ratio [SO]/[CS] *Note (1): An error of (0.00) means 1{sigma} error <0.005 *Note (2): Adopted SO/CS ratio, calculated for optically thin lines, * T_ex_=20K (see text). *Note (3): Explanations of the adopted SO/CS ratio: * a: ^34^SO(3_2_-2_1_)/C^34^S(2-1) * b: SO(3_2_-2_1_)/(22 x C^34^S(2-1)) * c: 22 x ^34^SO(3_2_-2_1_)/CS(2-1) * d: SO(3_2_-2_1_)/CS(2-1) C============================================================================= C Loading file 'table2.dat' ! Maxima and minima of SO/CS integrated intensity * ratios C Format for file interpretation 1 format( + A15,1X,A11,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,1X,I4,1X, + I4,1X,F4.2,1X,I3,1X,I4,1X,F4.2) C Effective file loading open(unit=1,status='old',file= +'table2.dat') write(6,*) '....Loading file: table2.dat' do i__=1,22 read(1,'(A77)')ar__ read(ar__,1) + Source,n_Source,RAh,RAm,RAs,DE_,DEd,DEm,DEs,ORAMax,ODEMax, + RatioMax,ORAMin,ODEMin,RatioMin 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) + Source,n_Source,RAh,RAm,RAs,DE_,DEd,DEm,DEs,ORAMax,ODEMax, + RatioMax,ORAMin,ODEMin,RatioMin write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table3.dat' ! Integrated intensity ratios (A-E) and estimated * abundance ratios (F) at certain positions C Format for file interpretation 2 format( + A15,2X,I4,1X,I4,1X,F4.1,1X,F3.1,1X,F4.1,1X,F3.1,1X,F4.2,1X, + F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,A1,1X,F3.1) C Effective file loading open(unit=1,status='old',file= +'table3.dat') write(6,*) '....Loading file: table3.dat' do i__=1,92 read(1,'(A80)')ar__1 read(ar__1,2) + Source_1,ORA,ODE,RatioA,e_RatioA,RatioB,e_RatioB,RatioC, + e_RatioC,RatioD,e_RatioD,RatioE,e_RatioE,Note,RadioF if(ar__1(28:31) .EQ. '') RatioA = rNULL__ if(ar__1(33:35) .EQ. '') e_RatioA = rNULL__ if(ar__1(37:40) .EQ. '') RatioB = rNULL__ if(ar__1(42:44) .EQ. '') e_RatioB = rNULL__ if(ar__1(46:49) .EQ. '') RatioC = rNULL__ if(ar__1(51:54) .EQ. '') e_RatioC = rNULL__ if(ar__1(56:59) .EQ. '') RatioD = rNULL__ if(ar__1(61:64) .EQ. '') e_RatioD = rNULL__ if(ar__1(66:69) .EQ. '') RatioE = rNULL__ if(ar__1(71:74) .EQ. '') e_RatioE = rNULL__ c ..............Just test output........... write(6,2) + Source_1,ORA,ODE,RatioA,e_RatioA,RatioB,e_RatioB,RatioC, + e_RatioC,RatioD,e_RatioD,RatioE,e_RatioE,Note,RadioF c .......End.of.Just test output........... end do close(1) C============================================================================= stop end