Conversion of standardized ReadMe file for
file /./ftp/cats/J/PASJ/59/335 into FORTRAN code for loading all data files into arrays.
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 (nr__) ! HD number of the star integer*4 HR (nr__) ! ? HR number of the star integer*4 HIP (nr__) ! ? HIP number of the star character*11 Name (nr__) ! Star name character*8 SpType (nr__) ! MK spectral type integer*4 Teff (nr__) ! (K) Effective temperature (1) real*4 logg (nr__) ! ([cm/s2]) Surface gravity (1) real*4 vt (nr__) ! (km/s) Microturbulent velocity (1) real*4 v_Fe_H_ (nr__) ! ([Sun]) Metallicity (1) real*4 logL (nr__) ! ([solLum]) Stellar luminosity derived as -0.4(M_bol_-4.75) real*4 Mass (nr__) ! (solMass) Mass (2) real*4 logAge (nr__) ! ([yr]) Age (2) character*1 q_logAge (nr__) ! [ABC] Quality flag showing the reliability of * age (see Sect. 2.1 for more details) (3) real*4 V_LSR (nr__) ! (km/s) Space speed relative to the LSR (G1) character*1 n_HD (nr__) ! [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 (nr__1) ! (deg) Right Ascension J2000, Epoch=1991.25, real*8 DEdeg (nr__1) ! (deg) Declination J2000, Epoch=1991.25, C ---------------------------------- ! (position vector(s) in degrees) character*1 n_HD_1 (nr__1) ! [P] 'P' when star harbors planet(s) integer*4 HD_1 (nr__1) ! HD number integer*4 RAh (nr__1) ! (h) Right ascension (J2000, Ep=1991.25, HIP) integer*4 RAm (nr__1) ! (min) Right ascension (J2000, Ep=1991.25, HIP) real*4 RAs (nr__1) ! (s) Right ascension (J2000, Ep=1991.25, HIP) character*1 DE_ (nr__1) ! Declination sign (J2000, Ep=1991.25, HIP integer*4 DEd (nr__1) ! (deg) Declination (J2000, Ep=1991.25, HIP) integer*4 DEm (nr__1) ! (arcmin) Declination (J2000, Ep=1991.25, HIP) real*4 DEs (nr__1) ! (arcsec) Declination (J2000, Ep=1991.25, HIP) real*4 plx (nr__1) ! (mas) Parallax (1) real*8 pmRA (nr__1) ! (mas/yr) Proper motion in alpha direction (1) real*8 pmDE (nr__1) ! (mas/yr) Proper motion in delta direction (1) real*4 HRV (nr__1) ! (km/s) Heliocentric radial velocity (from table 2 of * Paper I, Cat. J/PASJ/57/13) real*4 Rm (nr__1) ! (kpc) Mean galactocentric distance real*4 zm (nr__1) ! (kpc) Maximum separation from the galactic plane real*4 e (nr__1) ! Galactocentric orbital eccentricity real*4 Ulsr (nr__1) ! (km/s) Space velocity in radial direction relative * to LSR real*4 Vlsr (nr__1) ! (km/s) Space velocity in tangential direction * relative to LSR real*4 Wlsr (nr__1) ! (km/s) Space velocity in vertical direction * relative to LSR real*4 V_LSR_1 (nr__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 (nr__2) ! [P] 'P' when star harbors planet(s) integer*4 HD_2 (nr__2) ! HD number real*4 Teff2 (nr__2) ! (K) Effective temperature (1) real*4 logg2 (nr__2) ! ([cm/s+2]) logarithmic surface gravity (1) real*4 vt2 (nr__2) ! (km/s) microturbulent velocity (1) real*4 v_Fe_H_2 (nr__2) ! ([Sun]) Fe abundance relative to the Sun (1) real*4 VMAG2 (nr__2) ! (mag) Absolute visual magnitude (2) real*4 BC2 (nr__2) ! (mag) Bolometric correction (2) real*4 Mbol2 (nr__2) ! (mag) Bolometric magnitude (2) real*4 logL2 (nr__2) ! ([solLum]) logarithmic luminosity real*4 M_trk (nr__2) ! (solMass) Stellar mass derived from evolutionary * track (3) real*4 E_M_trk (nr__2) ! (solMass) Possible error in M(trk) on the upper side real*4 e_M_trk_1 (nr__2) ! (solMass) Possible error in M(trk) on the lower side real*4 M_iso (nr__2) ! (solMass) Stellar mass derived from isochrones real*4 logg_TLM (nr__2) ! ([cm/s+2]) Theorerical logg derived from * Teff, L, and M(trk) real*4 logAge_trk (nr__2) ! ([yr]) log(stellar age) derived from * evolutionary track (3) real*4 E_logAge_trk(nr__2) ! ([yr]) Possible error in logAge(trk) on the * upper side real*4 e_logAge_trk_1(nr__2) ! ([yr]) Possible error in logAge(trk) on the * lower side character*1 q_logAge_trk(nr__2) ! Reliability of age(trk) (4) real*4 logAge_iso (nr__2) ! ([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 (nr__3) ! HD number integer*4 Teff_1 (nr__3) ! (K) Effective temperature real*4 logg_1 (nr__3) ! ([cm/s2]) logarithmic surface gravity real*4 vt_1 (nr__3) ! (km/s) Microturbulent velocity real*4 v_Fe_H__1 (nr__3) ! ([Sun]) Fe abundance relative to the Sun real*4 v_Na_H_ (nr__3) ! ([Sun]) Na abundance relative to the Sun real*4 v_Mg_H_ (nr__3) ! ([Sun]) Mg abundance relative to the Sun * (MgI 8806 fit) real*4 v_Al_H_ (nr__3) ! ([Sun]) ?=99.99 Al abundance relative to the Sun real*4 v_Si_H_ (nr__3) ! ([Sun]) Si abundance relative to the Sun real*4 v_S_H_ (nr__3) ! ([Sun]) ?=99.99 S abundance relative to the Sun * (SI 8684 fit) real*4 v_Ca_H_ (nr__3) ! ([Sun]) Ca abundance relative to the Sun real*4 v_ScI_H_ (nr__3) ! ([Sun]) ?=99.99 Sc abundance relative to the Sun * (from ScI) real*4 v_ScII_H_ (nr__3) ! ([Sun]) Sc abundance relative to the Sun (from ScII) real*4 v_TiI_H_ (nr__3) ! ([Sun]) Ti abundance relative to the Sun (from TiI) real*4 v_TiII_H_ (nr__3) ! ([Sun]) Ti abundance relative to the Sun (from TiII) real*4 v_VI_H_ (nr__3) ! ([Sun]) ?=99.99 V abundance relative to the Sun * (from VI) real*4 v_VII_H_ (nr__3) ! ([Sun]) ?=99.99 V abundance relative to the Sun * (from VII) real*4 v_CrI_H_ (nr__3) ! ([Sun]) Cr abundance relative to the Sun (from CrI) real*4 v_CrII_H_ (nr__3) ! ([Sun]) Cr abundance relative to the Sun (from CrII) real*4 v_Mn_H_ (nr__3) ! ([Sun]) ?=99.99 Mn abundance relative to the Sun real*4 v_Co_H_ (nr__3) ! ([Sun]) Co abundance relative to the Sun real*4 v_Ni_H_ (nr__3) ! ([Sun]) Ni abundance relative to the Sun real*4 v_Cu_H_ (nr__3) ! ([Sun]) ?=99.99 Cu abundance relative to the Sun real*4 v_Zn_H_ (nr__3) ! ([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 (nr__4) ! HD number (999999 for the Sun) integer*4 Teff_2 (nr__4) ! (K) Effective temperature real*4 logg_2 (nr__4) ! ([cm/s2]) logarithmic surface gravity real*4 vt_2 (nr__4) ! (km/s) Microturbulent velocity real*4 v_Fe_H__2 (nr__4) ! ([Sun]) Fe abundance relative to the Sun real*4 A_Mg (nr__4) ! Mg abundance from 8806 fit * (normalisation of H=12.00) real*4 EW_Mg (nr__4) ! (0.1pm) EW(MgI 8806.8) inversely computed from A(Mg) real*4 A_S (nr__4) ! S abundance from 8693-5 fit * (normalisation of H=12.00) real*4 EW_S1 (nr__4) ! (0.1pm) ?=-99.9 EW(SI 8693.1) inversely computed * from A(S) real*4 EW_S2 (nr__4) ! (0.1pm) ?=-99.9 EW(SI 8693.9) inversely computed * from A(S) real*4 EW_S3 (nr__4) ! (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 (nr__5) ! (0.1nm) Line wavelength integer*4 o_Ion (nr__5) ! Number of hyperfine components real*4 Ion (nr__5) ! Ion code (G2) (1) real*8 lam_hf (nr__5) ! (0.1nm) Wavelength of the hyperfine components real*4 chi (nr__5) ! (eV) ? lower excitation potential real*4 loggf (nr__5) ! 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 (nr__6) ! Ion code (G2) (1) integer*4 lambda_1 (nr__6) ! (0.1nm) Line wavelength integer*4 HD_5 (nr__6) ! HD number (0 for Sun) real*4 Teff_3 (nr__6) ! (K) Effective temperature real*4 logg_3 (nr__6) ! ([m/s2]) Logarithmic surface gravity real*4 vt_3 (nr__6) ! (km/s) Microturbulent velocity real*4 v_Fe_H__3 (nr__6) ! Fe abundance relative to the Sun real*4 EW (nr__6) ! (0.1pm) Equivalent width of the line real*4 A_hfs (nr__6) ! Abundance obtained with hyperfine splitting real*4 A_nohf (nr__6) ! Abundance obtained without hyperfine splitting real*4 delta (nr__6) ! 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 (nr__7) ! HD number real*4 Teff_4 (nr__7) ! (K) Effective temperature real*4 logg_4 (nr__7) ! ([cm/s2]) logarithmic surface gravity real*4 vt_4 (nr__7) ! (km/s) Microturbulent velocity real*4 v_Fe_H__4 (nr__7) ! ([Sun]) Fe abundance relative to the Sun real*4 Ion_2 (nr__7) ! Ion code (G2) (1) real*4 v_X_H__hfs (nr__7) ! final [X/H] value obtained with hfs real*4 v_X_H__nohf(nr__7) ! final [X/H] value obtained without hfs real*4 delta_1 (nr__7) ! 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 (nr__8) ! HD number integer*4 Nline (nr__8) ! Number of lines for the species integer*4 o__delA_ (nr__8) ! Number of lines actually adopted for computing * the averaged abundances real*4 Ion_3 (nr__8) ! Ion code (G2) real*4 v_delA_ (nr__8) ! Average differential abundance (star-Sun) real*4 e__delA_ (nr__8) ! rms uncertainty on <delA> character*11 FileName (nr__8) ! 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 (nr__9) ! Ion code (G2) real*8 lambda_2 (nr__9) ! (0.1nm) Wavelength real*4 chi_1 (nr__9) ! (eV) Lower excitation potential real*4 loggf_1 (nr__9) ! log(gf) value real*4 EW_star (nr__9) ! (0.1pm) Equivalent width of the star real*4 vt_star (nr__9) ! (km/s) Adopted microturbulence of the star real*4 A_star (nr__9) ! Abundance of the ion in the star real*4 EW_Sun (nr__9) ! (0.1pm) Solar flux equivalent width real*4 vt_Sun (nr__9) ! (km/s) Adopted microturbulence of the Sun real*4 A_Sun (nr__9) ! Abundance of the ion in the Sun real*4 delA (nr__9) ! Differential abundance (star-Sun) character*1 Note (nr__9) ! [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(i__),HR(i__),HIP(i__),Name(i__),SpType(i__),Teff(i__), + logg(i__),vt(i__),v_Fe_H_(i__),logL(i__),Mass(i__), + logAge(i__),q_logAge(i__),V_LSR(i__),n_HD(i__) if(ar__(8:11) .EQ. '') HR(i__) = iNULL__ if(ar__(13:18) .EQ. '') HIP(i__) = iNULL__ c ..............Just test output........... write(6,1) + HD(i__),HR(i__),HIP(i__),Name(i__),SpType(i__),Teff(i__), + logg(i__),vt(i__),v_Fe_H_(i__),logL(i__),Mass(i__), + logAge(i__),q_logAge(i__),V_LSR(i__),n_HD(i__) 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(i__),HD_1(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__), + DEd(i__),DEm(i__),DEs(i__),plx(i__),pmRA(i__),pmDE(i__), + HRV(i__),Rm(i__),zm(i__),e(i__),Ulsr(i__),Vlsr(i__),Wlsr(i__), + V_LSR_1(i__) RAdeg(i__) = rNULL__ DEdeg(i__) = rNULL__ c Derive coordinates RAdeg and DEdeg from input data c (RAdeg and DEdeg are set to rNULL__ when unknown) if(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15. if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4. if(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240. if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__) if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60. if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600. if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__) c ..............Just test output........... write(6,2) + n_HD_1(i__),HD_1(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__), + DEd(i__),DEm(i__),DEs(i__),plx(i__),pmRA(i__),pmDE(i__), + HRV(i__),Rm(i__),zm(i__),e(i__),Ulsr(i__),Vlsr(i__),Wlsr(i__), + V_LSR_1(i__) write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__) 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(i__),HD_2(i__),Teff2(i__),logg2(i__),vt2(i__), + v_Fe_H_2(i__),VMAG2(i__),BC2(i__),Mbol2(i__),logL2(i__), + M_trk(i__),E_M_trk(i__),e_M_trk_1(i__),M_iso(i__), + logg_TLM(i__),logAge_trk(i__),E_logAge_trk(i__), + e_logAge_trk_1(i__),q_logAge_trk(i__),logAge_iso(i__) c ..............Just test output........... write(6,3) + n_HD_2(i__),HD_2(i__),Teff2(i__),logg2(i__),vt2(i__), + v_Fe_H_2(i__),VMAG2(i__),BC2(i__),Mbol2(i__),logL2(i__), + M_trk(i__),E_M_trk(i__),e_M_trk_1(i__),M_iso(i__), + logg_TLM(i__),logAge_trk(i__),E_logAge_trk(i__), + e_logAge_trk_1(i__),q_logAge_trk(i__),logAge_iso(i__) 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(i__),Teff_1(i__),logg_1(i__),vt_1(i__),v_Fe_H__1(i__), + v_Na_H_(i__),v_Mg_H_(i__),v_Al_H_(i__),v_Si_H_(i__), + v_S_H_(i__),v_Ca_H_(i__),v_ScI_H_(i__),v_ScII_H_(i__), + v_TiI_H_(i__),v_TiII_H_(i__),v_VI_H_(i__),v_VII_H_(i__), + v_CrI_H_(i__),v_CrII_H_(i__),v_Mn_H_(i__),v_Co_H_(i__), + v_Ni_H_(i__),v_Cu_H_(i__),v_Zn_H_(i__) c ..............Just test output........... write(6,4) + HD_3(i__),Teff_1(i__),logg_1(i__),vt_1(i__),v_Fe_H__1(i__), + v_Na_H_(i__),v_Mg_H_(i__),v_Al_H_(i__),v_Si_H_(i__), + v_S_H_(i__),v_Ca_H_(i__),v_ScI_H_(i__),v_ScII_H_(i__), + v_TiI_H_(i__),v_TiII_H_(i__),v_VI_H_(i__),v_VII_H_(i__), + v_CrI_H_(i__),v_CrII_H_(i__),v_Mn_H_(i__),v_Co_H_(i__), + v_Ni_H_(i__),v_Cu_H_(i__),v_Zn_H_(i__) 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(i__),Teff_2(i__),logg_2(i__),vt_2(i__),v_Fe_H__2(i__), + A_Mg(i__),EW_Mg(i__),A_S(i__),EW_S1(i__),EW_S2(i__),EW_S3(i__) c ..............Just test output........... write(6,5) + HD_4(i__),Teff_2(i__),logg_2(i__),vt_2(i__),v_Fe_H__2(i__), + A_Mg(i__),EW_Mg(i__),A_S(i__),EW_S1(i__),EW_S2(i__),EW_S3(i__) 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(i__),o_Ion(i__),Ion(i__),lam_hf(i__),chi(i__), + loggf(i__) if(ar__5(24:28) .EQ. '') chi(i__) = rNULL__ c ..............Just test output........... write(6,6) + lambda(i__),o_Ion(i__),Ion(i__),lam_hf(i__),chi(i__), + loggf(i__) 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(i__),lambda_1(i__),HD_5(i__),Teff_3(i__),logg_3(i__), + vt_3(i__),v_Fe_H__3(i__),EW(i__),A_hfs(i__),A_nohf(i__), + delta(i__) c ..............Just test output........... write(6,7) + Ion_1(i__),lambda_1(i__),HD_5(i__),Teff_3(i__),logg_3(i__), + vt_3(i__),v_Fe_H__3(i__),EW(i__),A_hfs(i__),A_nohf(i__), + delta(i__) 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(i__),Teff_4(i__),logg_4(i__),vt_4(i__),v_Fe_H__4(i__), + Ion_2(i__),v_X_H__hfs(i__),v_X_H__nohf(i__),delta_1(i__) c ..............Just test output........... write(6,8) + HD_6(i__),Teff_4(i__),logg_4(i__),vt_4(i__),v_Fe_H__4(i__), + Ion_2(i__),v_X_H__hfs(i__),v_X_H__nohf(i__),delta_1(i__) 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(i__),Nline(i__),o__delA_(i__),Ion_3(i__),v_delA_(i__), + e__delA_(i__),FileName(i__) c ..............Just test output........... write(6,9) + HD_7(i__),Nline(i__),o__delA_(i__),Ion_3(i__),v_delA_(i__), + e__delA_(i__),FileName(i__) 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(i__),lambda_2(i__),chi_1(i__),loggf_1(i__),EW_star(i__), + vt_star(i__),A_star(i__),EW_Sun(i__),vt_Sun(i__),A_Sun(i__), + delA(i__),Note(i__) c ..............Just test output........... write(6,10) + Ion_4(i__),lambda_2(i__),chi_1(i__),loggf_1(i__),EW_star(i__), + vt_star(i__),A_star(i__),EW_Sun(i__),vt_Sun(i__),A_Sun(i__), + delA(i__),Note(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end