Conversion of standardized ReadMe file for
file /./ftp/cats/J/AZh/84/997 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-Apr-19 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/84/997 Abundances of Sr, Y, Zr, Ce and Ba (Mashonkina+, 2007) *================================================================================ *Neutron capture elements in halo, thick disk and thin disk stars: strontium, *yttrium, zirconium, cerium. * Mashonkina L.I., Vinogradova A.B., Ptitsyn D.A., Khokhlova V.S., * Chernetsova T.A. * <Astron. Zh. 84, 997 (2007)> * =2007AZh....84..997M * =2007ARep...51..903M C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'stars.dat' ! Objects with coordinates integer*4 nr__ parameter (nr__=74) ! Number of records character*53 ar__ ! Full-size record C 2000 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg ! (deg) Right Ascension 2000 real*8 DEdeg ! (deg) Declination 2000 C ---------------------------------- ! (position vector(s) in degrees) character*10 HD_BD ! HD/BD number integer*4 RAh ! (h) Right Ascension (2000) (hours) integer*4 RAm ! (min) Right Ascension (2000) (minutes) real*8 RAs ! (s) Right Ascension (2000) (seconds) character*1 DE_ ! Declination (2000) (sign) integer*4 DEd ! (deg) Declination (2000) (degrees) integer*4 DEm ! (arcmin) Declination (2000) (minutes) real*4 DEs ! (arcsec) Declination (2000) (seconds) character*15 OName ! Other name c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table1.dat' ! Atomic data and solar abundances for investigated lines obtained from solar line profile fitting integer*4 nr__1 parameter (nr__1=15) ! Number of records character*66 ar__1 ! Full-size record character*5 Ion ! Species real*8 lambda ! (0.1nm) Wavelength real*4 Eexc ! (eV) Excitation energy (1) integer*4 Mult ! Multiplicity number real*4 loggf ! ([-]) log of oscillation strength integer*4 Ref ! Reference (2) real*4 logC6 ! C_6_ abundance character*1 n_logC6 ! [1-3] Origin of logC6 (3) real*4 log_eps ! ([-]) Solar abundance (H=12) character*13 Rem ! Remarks *Note (1): According to Moore, 1972, A multiplet table of astrophysical * interest, NSRD-NBS 40 (Washington, D.C.). *Note (2): References as follows: * 13 = Kurucz, 1994, CD-ROM Nos 18, 19 * 23 = Hannaford et al., 1982ApJ...261..736H * 26 = Brage et al., 1998ApJ...496.1051B * 27 = Ljung et al., 2006A&A...456.1181L * 28 = Palmeri et al., 2000, Phys. Scr. 61, 323 * 29 = Reader et al., Wavelengths and Transaction Probabilities for * atoms and atomic ions, PartII, NSRD-NBS 68 (Washington, D.C., 1980) *Note (3): Note as follows: * 1 = based on profile analysis for the solar spectrum * 2 = data of Kurucz, CD-ROM Nos 18, 19 * 3 = data of Barklem et O'Mara, 1998MNRAS.300..863B c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table2.dat' ! Stellar parameters and element abundance ratios of the stars with [Fe/H]>-0.9 integer*4 nr__2 parameter (nr__2=54) ! Number of records character*70 ar__2 ! Full-size record character*10 HD_BD_1 ! HD number integer*4 Teff ! (K) Effective temperature real*4 logg ! Surface gravity real*4 v_Fe_H_ ! Iron abundance real*4 VT ! (km/s) Microturbulence velocity integer*4 Code ! [0/3] Stellar population type (1) real*4 v_Sr_Fe_ ! ? Strontium to Fe abundance real*4 v_Y_Fe_ ! ? Yttrium to Fe abundance real*4 v_Zr_Fe_ ! ? Zirconium to Fe abundance real*4 v_Ba_Fe_ ! Barium to Fe abundance real*4 v_Ce_Fe_ ! ? Cerium to Fe abundance *Note (1): population is coded as follows: * 0 = thin disk * 1 = thick disk * 2 = halo * 3 = transition phase between thin and thick disks c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table3.dat' ! Stellar parameters and element abundance ratios of the stars with [Fe/H]<-0.9 integer*4 nr__3 parameter (nr__3=20) ! Number of records character*70 ar__3 ! Full-size record character*10 HD_BD_2 ! HD number integer*4 Teff_1 ! (K) Effective temperature real*4 logg_1 ! Surface gravity real*4 v_Fe_H__1 ! Iron abundance real*4 VT_1 ! (km/s) Microturbulence velocity integer*4 Code_1 ! [0/3] Stellar population type (1) real*4 v_Sr_Fe__1 ! ? Strontium to Fe abundance real*4 v_Y_Fe__1 ! ? Yttrium to Fe abundance real*4 v_Zr_Fe__1 ! ? Zirconium to Fe abundance real*4 v_Ba_Fe__1 ! Barium to Fe abundance real*4 v_Ce_Fe__1 ! ? Cerium to Fe abundance *Note (1): population is coded as follows: * 0 = thin disk * 1 = thick disk * 2 = halo * 3 = transition phase between thin and thick disks C============================================================================= C Loading file 'stars.dat' ! Objects with coordinates C Format for file interpretation 1 format(A10,I2,1X,I2,1X,F7.4,1X,A1,I2,1X,I2,1X,F6.3,1X,A15) C Effective file loading open(unit=1,status='old',file= +'stars.dat') write(6,*) '....Loading file: stars.dat' do i__=1,74 read(1,'(A53)')ar__ read(ar__,1)HD_BD,RAh,RAm,RAs,DE_,DEd,DEm,DEs,OName 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)HD_BD,RAh,RAm,RAs,DE_,DEd,DEm,DEs,OName write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table1.dat' ! Atomic data and solar abundances for investigated * lines obtained from solar line profile fitting C Format for file interpretation 2 format( + A5,2X,F7.2,2X,F5.3,2X,I2,1X,F6.3,2X,I2,2X,F6.2,1X,A1,1X,F4.2, + 2X,A13) C Effective file loading open(unit=1,status='old',file= +'table1.dat') write(6,*) '....Loading file: table1.dat' do i__=1,15 read(1,'(A66)')ar__1 read(ar__1,2) + Ion,lambda,Eexc,Mult,loggf,Ref,logC6,n_logC6,log_eps,Rem c ..............Just test output........... write(6,2) + Ion,lambda,Eexc,Mult,loggf,Ref,logC6,n_logC6,log_eps,Rem c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table2.dat' ! Stellar parameters and element abundance ratios * of the stars with [Fe/H]>-0.9 C Format for file interpretation 3 format( + A10,1X,I4,2X,F4.2,2X,F5.2,2X,F3.1,2X,I1,2X,F5.2,2X,F5.2,1X, + F5.2,2X,F5.2,2X,F5.2) C Effective file loading open(unit=1,status='old',file= +'table2.dat') write(6,*) '....Loading file: table2.dat' do i__=1,54 read(1,'(A70)')ar__2 read(ar__2,3) + HD_BD_1,Teff,logg,v_Fe_H_,VT,Code,v_Sr_Fe_,v_Y_Fe_,v_Zr_Fe_, + v_Ba_Fe_,v_Ce_Fe_ if(ar__2(39:43) .EQ. '') v_Sr_Fe_ = rNULL__ if(ar__2(46:50) .EQ. '') v_Y_Fe_ = rNULL__ if(ar__2(52:56) .EQ. '') v_Zr_Fe_ = rNULL__ if(ar__2(66:70) .EQ. '') v_Ce_Fe_ = rNULL__ c ..............Just test output........... write(6,3) + HD_BD_1,Teff,logg,v_Fe_H_,VT,Code,v_Sr_Fe_,v_Y_Fe_,v_Zr_Fe_, + v_Ba_Fe_,v_Ce_Fe_ c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table3.dat' ! Stellar parameters and element abundance ratios * of the stars with [Fe/H]<-0.9 C Format for file interpretation 4 format( + A10,1X,I4,2X,F4.2,2X,F5.2,2X,F3.1,2X,I1,2X,F5.2,2X,F5.2,1X, + F5.2,2X,F5.2,2X,F5.2) C Effective file loading open(unit=1,status='old',file= +'table3.dat') write(6,*) '....Loading file: table3.dat' do i__=1,20 read(1,'(A70)')ar__3 read(ar__3,4) + HD_BD_2,Teff_1,logg_1,v_Fe_H__1,VT_1,Code_1,v_Sr_Fe__1, + v_Y_Fe__1,v_Zr_Fe__1,v_Ba_Fe__1,v_Ce_Fe__1 if(ar__3(39:43) .EQ. '') v_Sr_Fe__1 = rNULL__ if(ar__3(46:50) .EQ. '') v_Y_Fe__1 = rNULL__ if(ar__3(52:56) .EQ. '') v_Zr_Fe__1 = rNULL__ if(ar__3(66:70) .EQ. '') v_Ce_Fe__1 = rNULL__ c ..............Just test output........... write(6,4) + HD_BD_2,Teff_1,logg_1,v_Fe_H__1,VT_1,Code_1,v_Sr_Fe__1, + v_Y_Fe__1,v_Zr_Fe__1,v_Ba_Fe__1,v_Ce_Fe__1 c .......End.of.Just test output........... end do close(1) C============================================================================= stop end