Conversion of standardized ReadMe file for
file /./ftp/cats/J/PASJ/59/335 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-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/PASJ/59/335 Properties of 160 F-K disk dwarfs/subgiants (Takeda+, 2007) *================================================================================ *Fundamental parameters and elemental abundances of 160 F-G-K stars based on *OAO spectrum database. * Takeda Y. * <Publ. Astron. Soc. Jap., 59, 335-356 (2007)> * =2007PASJ...59..335T C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table1.dat' ! Fundamental stellar parameters of program stars integer*4 nr__ parameter (nr__=160) ! Number of records character*87 ar__ ! Full-size record integer*4 HD ! HD number of the star integer*4 HR ! ? HR number of the star integer*4 HIP ! ? HIP number of the star character*11 Name ! Star name character*8 SpType ! MK spectral type integer*4 Teff ! (K) Effective temperature (1) real*4 logg ! ([cm/s2]) Surface gravity (1) real*4 vt ! (km/s) Microturbulent velocity (1) real*4 v_Fe_H_ ! ([Sun]) Metallicity (1) real*4 logL ! ([solLum]) Stellar luminosity derived as -0.4(M_bol_-4.75) real*4 Mass ! (solMass) Mass (2) real*4 logAge ! ([yr]) Age (2) character*1 q_logAge ! [ABC] Quality flag showing the reliability of * age (see Sect. 2.1 for more details) (3) real*4 V_LSR ! (km/s) Space speed relative to the LSR (G1) character*1 n_HD ! [P] 'P' when star harbors planet(s) *Note (1): Determined from FeI and FeII lines taken from table 1 of Paper II. *Note (2): Mass and age derived by comparing the (logTeff, logL) position on * the HR diagram with Girardi et al.'s (2000, Cat. J/A+AS/141/371) * theoretical evolutionary tracks *Note (3): Quality flag as follows: * A = high * B = moderate * C = poor c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'kinepara.dat' ! Kinematic parameters of program stars (tablee1) integer*4 nr__1 parameter (nr__1=160) ! Number of records character*109 ar__1 ! Full-size record C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs (Epoch=1991.25,) real*8 RAdeg ! (deg) Right Ascension J2000, Epoch=1991.25, real*8 DEdeg ! (deg) Declination J2000, Epoch=1991.25, C ---------------------------------- ! (position vector(s) in degrees) character*1 n_HD_1 ! [P] 'P' when star harbors planet(s) integer*4 HD_1 ! HD number integer*4 RAh ! (h) Right ascension (J2000, Ep=1991.25, HIP) integer*4 RAm ! (min) Right ascension (J2000, Ep=1991.25, HIP) real*4 RAs ! (s) Right ascension (J2000, Ep=1991.25, HIP) character*1 DE_ ! Declination sign (J2000, Ep=1991.25, HIP integer*4 DEd ! (deg) Declination (J2000, Ep=1991.25, HIP) integer*4 DEm ! (arcmin) Declination (J2000, Ep=1991.25, HIP) real*4 DEs ! (arcsec) Declination (J2000, Ep=1991.25, HIP) real*4 plx ! (mas) Parallax (1) real*8 pmRA ! (mas/yr) Proper motion in alpha direction (1) real*8 pmDE ! (mas/yr) Proper motion in delta direction (1) real*4 HRV ! (km/s) Heliocentric radial velocity (from table 2 of * Paper I, Cat. J/PASJ/57/13) real*4 Rm ! (kpc) Mean galactocentric distance real*4 zm ! (kpc) Maximum separation from the galactic plane real*4 e ! Galactocentric orbital eccentricity real*4 Ulsr ! (km/s) Space velocity in radial direction relative * to LSR real*4 Vlsr ! (km/s) Space velocity in tangential direction * relative to LSR real*4 Wlsr ! (km/s) Space velocity in vertical direction * relative to LSR real*4 V_LSR_1 ! (km/s) speed relative to LSR (G1) *Note (1): From HIP catalog c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'starpara.dat' ! Fundamental stellar parameters of 160 stars (tablee2) integer*4 nr__2 parameter (nr__2=160) ! Number of records character*105 ar__2 ! Full-size record character*1 n_HD_2 ! [P] 'P' when star harbors planet(s) integer*4 HD_2 ! HD number real*4 Teff2 ! (K) Effective temperature (1) real*4 logg2 ! ([cm/s+2]) logarithmic surface gravity (1) real*4 vt2 ! (km/s) microturbulent velocity (1) real*4 v_Fe_H_2 ! ([Sun]) Fe abundance relative to the Sun (1) real*4 VMAG2 ! (mag) Absolute visual magnitude (2) real*4 BC2 ! (mag) Bolometric correction (2) real*4 Mbol2 ! (mag) Bolometric magnitude (2) real*4 logL2 ! ([solLum]) logarithmic luminosity real*4 M_trk ! (solMass) Stellar mass derived from evolutionary * track (3) real*4 E_M_trk ! (solMass) Possible error in M(trk) on the upper side real*4 e_M_trk_1 ! (solMass) Possible error in M(trk) on the lower side real*4 M_iso ! (solMass) Stellar mass derived from isochrones real*4 logg_TLM ! ([cm/s+2]) Theorerical logg derived from * Teff, L, and M(trk) real*4 logAge_trk ! ([yr]) log(stellar age) derived from * evolutionary track (3) real*4 E_logAge_trk ! ([yr]) Possible error in logAge(trk) on the * upper side real*4 e_logAge_trk_1 ! ([yr]) Possible error in logAge(trk) on the * lower side character*1 q_logAge_trk ! Reliability of age(trk) (4) real*4 logAge_iso ! ([yr]) log(stellar age) derived from isochrones *Note (1): Teff, log g, vt, and [Fe/H]) spectroscopically established in * Paper II, Cat. J/PASJ/57/27 *Note (2): From table1 of paper II, Cat. J/PASJ/57/27 *Note (3): Age and mass revisited in this paper by comparing the positions * on the HR diagram with theoretically evolutionary tracks/isochrones. *Note (4): Quality flag as follows: * A = high * B = moderate * C = poor c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'xhtable.dat' ! [X/H] (relative-to-Sun abundance for 19 species) (from tablee3) integer*4 nr__3 parameter (nr__3=160) ! Number of records character*141 ar__3 ! Full-size record integer*4 HD_3 ! HD number integer*4 Teff_1 ! (K) Effective temperature real*4 logg_1 ! ([cm/s2]) logarithmic surface gravity real*4 vt_1 ! (km/s) Microturbulent velocity real*4 v_Fe_H__1 ! ([Sun]) Fe abundance relative to the Sun real*4 v_Na_H_ ! ([Sun]) Na abundance relative to the Sun real*4 v_Mg_H_ ! ([Sun]) Mg abundance relative to the Sun * (MgI 8806 fit) real*4 v_Al_H_ ! ([Sun]) ?=99.99 Al abundance relative to the Sun real*4 v_Si_H_ ! ([Sun]) Si abundance relative to the Sun real*4 v_S_H_ ! ([Sun]) ?=99.99 S abundance relative to the Sun * (SI 8684 fit) real*4 v_Ca_H_ ! ([Sun]) Ca abundance relative to the Sun real*4 v_ScI_H_ ! ([Sun]) ?=99.99 Sc abundance relative to the Sun * (from ScI) real*4 v_ScII_H_ ! ([Sun]) Sc abundance relative to the Sun (from ScII) real*4 v_TiI_H_ ! ([Sun]) Ti abundance relative to the Sun (from TiI) real*4 v_TiII_H_ ! ([Sun]) Ti abundance relative to the Sun (from TiII) real*4 v_VI_H_ ! ([Sun]) ?=99.99 V abundance relative to the Sun * (from VI) real*4 v_VII_H_ ! ([Sun]) ?=99.99 V abundance relative to the Sun * (from VII) real*4 v_CrI_H_ ! ([Sun]) Cr abundance relative to the Sun (from CrI) real*4 v_CrII_H_ ! ([Sun]) Cr abundance relative to the Sun (from CrII) real*4 v_Mn_H_ ! ([Sun]) ?=99.99 Mn abundance relative to the Sun real*4 v_Co_H_ ! ([Sun]) Co abundance relative to the Sun real*4 v_Ni_H_ ! ([Sun]) Ni abundance relative to the Sun real*4 v_Cu_H_ ! ([Sun]) ?=99.99 Cu abundance relative to the Sun real*4 v_Zn_H_ ! ([Sun]) ?=99.99 Zn abundance relative to the Sun c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'mgsabund.dat' ! Profile-fit abundance for Mg and S and computed equivalence widhts (from tablee3) integer*4 nr__4 parameter (nr__4=161) ! Number of records character*61 ar__4 ! Full-size record integer*4 HD_4 ! HD number (999999 for the Sun) integer*4 Teff_2 ! (K) Effective temperature real*4 logg_2 ! ([cm/s2]) logarithmic surface gravity real*4 vt_2 ! (km/s) Microturbulent velocity real*4 v_Fe_H__2 ! ([Sun]) Fe abundance relative to the Sun real*4 A_Mg ! Mg abundance from 8806 fit * (normalisation of H=12.00) real*4 EW_Mg ! (0.1pm) EW(MgI 8806.8) inversely computed from A(Mg) real*4 A_S ! S abundance from 8693-5 fit * (normalisation of H=12.00) real*4 EW_S1 ! (0.1pm) ?=-99.9 EW(SI 8693.1) inversely computed * from A(S) real*4 EW_S2 ! (0.1pm) ?=-99.9 EW(SI 8693.9) inversely computed * from A(S) real*4 EW_S3 ! (0.1pm) ?=-99.9 EW(SI 8694.6) inversely computed * from A(S) c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'lgf.dat' ! *Wavelength and the loggf data of each hyper-fine splitting (hfs) component (from tablee4) integer*4 nr__5 parameter (nr__5=365) ! Number of records character*35 ar__5 ! Full-size record integer*4 lambda ! (0.1nm) Line wavelength integer*4 o_Ion ! Number of hyperfine components real*4 Ion ! Ion code (G2) (1) real*8 lam_hf ! (0.1nm) Wavelength of the hyperfine components real*4 chi ! (eV) ? lower excitation potential real*4 loggf ! log of (statistical weight times oscillator * strength) *Note (1): 21.01, 23.00, 25.00, 27.00 for ScII, VI, MnI, CoI respectively c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'dif.dat' ! Results of abundance corrections caused by inclusion of the hyper-fine splitting (hfs) effect for the relevant line (from tablee4) integer*4 nr__6 parameter (nr__6=2824) ! Number of records character*74 ar__6 ! Full-size record real*4 Ion_1 ! Ion code (G2) (1) integer*4 lambda_1 ! (0.1nm) Line wavelength integer*4 HD_5 ! HD number (0 for Sun) real*4 Teff_3 ! (K) Effective temperature real*4 logg_3 ! ([m/s2]) Logarithmic surface gravity real*4 vt_3 ! (km/s) Microturbulent velocity real*4 v_Fe_H__3 ! Fe abundance relative to the Sun real*4 EW ! (0.1pm) Equivalent width of the line real*4 A_hfs ! Abundance obtained with hyperfine splitting real*4 A_nohf ! Abundance obtained without hyperfine splitting real*4 delta ! hfs correction defined as A(hfs)-A(nohf) *Note (1): 21.01, 23.00, 25.00, 27.00 for ScII, VI, MnI, CoI respectively c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'xhcor.dat' ! Corrections on the final (line-averaged) results of [ScII/H], [VI/H], [MnI/H] and [CoI/H] caused by inclusion of the hyper-fine splitting effect for each of the relevant lines (from tablee4) integer*4 nr__7 parameter (nr__7=569) ! Number of records character*66 ar__7 ! Full-size record integer*4 HD_6 ! HD number real*4 Teff_4 ! (K) Effective temperature real*4 logg_4 ! ([cm/s2]) logarithmic surface gravity real*4 vt_4 ! (km/s) Microturbulent velocity real*4 v_Fe_H__4 ! ([Sun]) Fe abundance relative to the Sun real*4 Ion_2 ! Ion code (G2) (1) real*4 v_X_H__hfs ! final [X/H] value obtained with hfs real*4 v_X_H__nohf ! final [X/H] value obtained without hfs real*4 delta_1 ! hfs correction ([X/H](hfs)-[X/H](nohf)) *Note (1): 21.01, 23.00, 25.00, 27.00 for ScII, VI, MnI, CoI respectively c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'mabund.dat' ! Average abundances (from tablee3) integer*4 nr__8 parameter (nr__8=2552) ! Number of records character*45 ar__8 ! Full-size record integer*4 HD_7 ! HD number integer*4 Nline ! Number of lines for the species integer*4 o__delA_ ! Number of lines actually adopted for computing * the averaged abundances real*4 Ion_3 ! Ion code (G2) real*4 v_delA_ ! Average differential abundance (star-Sun) real*4 e__delA_ ! rms uncertainty on <delA> character*11 FileName ! Name of the file with individual measurements * in subdirectory cmb c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'cmb/*' ! Individual abundances (from tablee3) integer*4 nr__9 parameter (nr__9=160) ! Number of records character*79 ar__9 ! Full-size record real*4 Ion_4 ! Ion code (G2) real*8 lambda_2 ! (0.1nm) Wavelength real*4 chi_1 ! (eV) Lower excitation potential real*4 loggf_1 ! log(gf) value real*4 EW_star ! (0.1pm) Equivalent width of the star real*4 vt_star ! (km/s) Adopted microturbulence of the star real*4 A_star ! Abundance of the ion in the star real*4 EW_Sun ! (0.1pm) Solar flux equivalent width real*4 vt_Sun ! (km/s) Adopted microturbulence of the Sun real*4 A_Sun ! Abundance of the ion in the Sun real*4 delA ! Differential abundance (star-Sun) character*1 Note ! [x] x for rejection of inadequate data * according to Chauvenet's criterion C============================================================================= C Loading file 'table1.dat' ! Fundamental stellar parameters of program stars C Format for file interpretation 1 format( + I6,1X,I4,1X,I6,1X,A11,1X,A8,1X,I4,1X,F4.2,1X,F4.2,1X,F5.2,1X, + F5.2,1X,F4.2,1X,F5.2,1X,A1,1X,F5.1,1X,A1) C Effective file loading open(unit=1,status='old',file= +'table1.dat') write(6,*) '....Loading file: table1.dat' do i__=1,160 read(1,'(A87)')ar__ read(ar__,1) + HD,HR,HIP,Name,SpType,Teff,logg,vt,v_Fe_H_,logL,Mass,logAge, + q_logAge,V_LSR,n_HD if(ar__(8:11) .EQ. '') HR = iNULL__ if(ar__(13:18) .EQ. '') HIP = iNULL__ c ..............Just test output........... write(6,1) + HD,HR,HIP,Name,SpType,Teff,logg,vt,v_Fe_H_,logL,Mass,logAge, + q_logAge,V_LSR,n_HD c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'kinepara.dat' ! Kinematic parameters of program stars (tablee1) C Format for file interpretation 2 format( + A1,1X,I6,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,F4.1,1X,F6.2, + 1X,F8.2,1X,F8.2,1X,F6.1,1X,F6.3,1X,F5.3,1X,F5.3,1X,F6.1,1X, + F6.1,1X,F5.1,1X,F5.1) C Effective file loading open(unit=1,status='old',file= +'kinepara.dat') write(6,*) '....Loading file: kinepara.dat' do i__=1,160 read(1,'(A109)')ar__1 read(ar__1,2) + n_HD_1,HD_1,RAh,RAm,RAs,DE_,DEd,DEm,DEs,plx,pmRA,pmDE,HRV,Rm, + zm,e,Ulsr,Vlsr,Wlsr,V_LSR_1 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,2) + n_HD_1,HD_1,RAh,RAm,RAs,DE_,DEd,DEm,DEs,plx,pmRA,pmDE,HRV,Rm, + zm,e,Ulsr,Vlsr,Wlsr,V_LSR_1 write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'starpara.dat' ! Fundamental stellar parameters of 160 stars * (tablee2) C Format for file interpretation 3 format( + A1,1X,I6,1X,F6.1,1X,F5.3,1X,F4.2,1X,F6.3,1X,F4.2,1X,F5.2,1X, + F4.2,1X,F6.3,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F5.2, + 1X,F4.2,1X,F4.2,1X,A1,1X,F5.2) C Effective file loading open(unit=1,status='old',file= +'starpara.dat') write(6,*) '....Loading file: starpara.dat' do i__=1,160 read(1,'(A105)')ar__2 read(ar__2,3) + n_HD_2,HD_2,Teff2,logg2,vt2,v_Fe_H_2,VMAG2,BC2,Mbol2,logL2, + M_trk,E_M_trk,e_M_trk_1,M_iso,logg_TLM,logAge_trk, + E_logAge_trk,e_logAge_trk_1,q_logAge_trk,logAge_iso c ..............Just test output........... write(6,3) + n_HD_2,HD_2,Teff2,logg2,vt2,v_Fe_H_2,VMAG2,BC2,Mbol2,logL2, + M_trk,E_M_trk,e_M_trk_1,M_iso,logg_TLM,logAge_trk, + E_logAge_trk,e_logAge_trk_1,q_logAge_trk,logAge_iso c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'xhtable.dat' ! [X/H] (relative-to-Sun abundance for 19 species) * (from tablee3) C Format for file interpretation 4 format( + I6,1X,I4,1X,F4.2,1X,F4.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X, + F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2, + 1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X, + F5.2) C Effective file loading open(unit=1,status='old',file= +'xhtable.dat') write(6,*) '....Loading file: xhtable.dat' do i__=1,160 read(1,'(A141)')ar__3 read(ar__3,4) + HD_3,Teff_1,logg_1,vt_1,v_Fe_H__1,v_Na_H_,v_Mg_H_,v_Al_H_, + v_Si_H_,v_S_H_,v_Ca_H_,v_ScI_H_,v_ScII_H_,v_TiI_H_,v_TiII_H_, + v_VI_H_,v_VII_H_,v_CrI_H_,v_CrII_H_,v_Mn_H_,v_Co_H_,v_Ni_H_, + v_Cu_H_,v_Zn_H_ c ..............Just test output........... write(6,4) + HD_3,Teff_1,logg_1,vt_1,v_Fe_H__1,v_Na_H_,v_Mg_H_,v_Al_H_, + v_Si_H_,v_S_H_,v_Ca_H_,v_ScI_H_,v_ScII_H_,v_TiI_H_,v_TiII_H_, + v_VI_H_,v_VII_H_,v_CrI_H_,v_CrII_H_,v_Mn_H_,v_Co_H_,v_Ni_H_, + v_Cu_H_,v_Zn_H_ c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'mgsabund.dat' ! Profile-fit abundance for Mg and S and computed * equivalence widhts (from tablee3) C Format for file interpretation 5 format( + I6,1X,I4,1X,F4.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.1,1X,F4.2,1X, + F5.1,1X,F5.1,1X,F5.1) C Effective file loading open(unit=1,status='old',file= +'mgsabund.dat') write(6,*) '....Loading file: mgsabund.dat' do i__=1,161 read(1,'(A61)')ar__4 read(ar__4,5) + HD_4,Teff_2,logg_2,vt_2,v_Fe_H__2,A_Mg,EW_Mg,A_S,EW_S1,EW_S2, + EW_S3 c ..............Just test output........... write(6,5) + HD_4,Teff_2,logg_2,vt_2,v_Fe_H__2,A_Mg,EW_Mg,A_S,EW_S1,EW_S2, + EW_S3 c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'lgf.dat' ! *Wavelength and the loggf data of each hyper-fine * splitting (hfs) component (from tablee4) C Format for file interpretation 6 format(I4,1X,I2,1X,F5.2,1X,F8.3,1X,F5.3,1X,F6.3) C Effective file loading open(unit=1,status='old',file= +'lgf.dat') write(6,*) '....Loading file: lgf.dat' do i__=1,365 read(1,'(A35)')ar__5 read(ar__5,6)lambda,o_Ion,Ion,lam_hf,chi,loggf if(ar__5(24:28) .EQ. '') chi = rNULL__ c ..............Just test output........... write(6,6)lambda,o_Ion,Ion,lam_hf,chi,loggf c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'dif.dat' ! Results of abundance corrections caused by * inclusion of the hyper-fine splitting (hfs) * effect for the relevant line (from tablee4) C Format for file interpretation 7 format( + F5.2,1X,I4,1X,I6,2X,F6.1,2X,F5.3,2X,F4.2,1X,F6.3,2X,F5.1,1X, + F6.3,1X,F6.3,2X,F6.3) C Effective file loading open(unit=1,status='old',file= +'dif.dat') write(6,*) '....Loading file: dif.dat' do i__=1,2824 read(1,'(A74)')ar__6 read(ar__6,7) + Ion_1,lambda_1,HD_5,Teff_3,logg_3,vt_3,v_Fe_H__3,EW,A_hfs, + A_nohf,delta c ..............Just test output........... write(6,7) + Ion_1,lambda_1,HD_5,Teff_3,logg_3,vt_3,v_Fe_H__3,EW,A_hfs, + A_nohf,delta c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'xhcor.dat' ! Corrections on the final (line-averaged) results * of [ScII/H], [VI/H], [MnI/H] and [CoI/H] caused * by inclusion of the hyper-fine splitting effect * for each of the relevant lines (from tablee4) C Format for file interpretation 8 format( + I6,2X,F6.1,2X,F5.3,2X,F4.2,1X,F6.3,2X,F5.2,3X,F6.3,1X,F6.3,3X, + F6.3) C Effective file loading open(unit=1,status='old',file= +'xhcor.dat') write(6,*) '....Loading file: xhcor.dat' do i__=1,569 read(1,'(A66)')ar__7 read(ar__7,8) + HD_6,Teff_4,logg_4,vt_4,v_Fe_H__4,Ion_2,v_X_H__hfs, + v_X_H__nohf,delta_1 c ..............Just test output........... write(6,8) + HD_6,Teff_4,logg_4,vt_4,v_Fe_H__4,Ion_2,v_X_H__hfs, + v_X_H__nohf,delta_1 c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'mabund.dat' ! Average abundances (from tablee3) C Format for file interpretation 9 format(3X,I6,2X,I2,1X,I2,1X,F5.2,1X,F5.2,1X,F4.2,1X,A11) C Effective file loading open(unit=1,status='old',file= +'mabund.dat') write(6,*) '....Loading file: mabund.dat' do i__=1,2552 read(1,'(A45)')ar__8 read(ar__8,9) + HD_7,Nline,o__delA_,Ion_3,v_delA_,e__delA_,FileName c ..............Just test output........... write(6,9) + HD_7,Nline,o__delA_,Ion_3,v_delA_,e__delA_,FileName c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'cmb/*' ! Individual abundances (from tablee3) C Format for file interpretation 10 format( + 1X,F5.2,F9.3,2X,F5.3,1X,F6.3,2X,F5.1,2X,F4.2,2X,F5.3,3X,F4.1, + 2X,F4.2,2X,F5.3,3X,F5.2,1X,A1) C Effective file loading open(unit=1,status='old',file= +'cmb/*') write(6,*) '....Loading file: cmb/*' do i__=1,160 read(1,'(A79)')ar__9 read(ar__9,10) + Ion_4,lambda_2,chi_1,loggf_1,EW_star,vt_star,A_star,EW_Sun, + vt_Sun,A_Sun,delA,Note c ..............Just test output........... write(6,10) + Ion_4,lambda_2,chi_1,loggf_1,EW_star,vt_star,A_star,EW_Sun, + vt_Sun,A_Sun,delA,Note c .......End.of.Just test output........... end do close(1) C============================================================================= stop end