Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/631/A77 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 2026-Oct-05
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/631/A77 Low-mass stellar models and isochrones (Amard+, 2019)
*================================================================================
*First grids of low-mass stellar models and isochrones with self-consistent
*treatment of rotation.
*From 0.2 to 1.5 M_{sun}_ at seven metallicities from PMS to TAMS.
* Amard L., Palacios A., Charbonnel C., Gallet F., Georgy C., Lagarde N.,
* Siess L.
* <Astron. Astrophys. 631, A77 (2019)>
* =2019A&A...631A..77A (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'list.dat' ! List of all the tracks and isochrone files
integer*4 nr__
parameter (nr__=1568) ! Number of records
character*67 ar__ ! Full-size record
real*4 Mini ! (Msun) [0.20/1.50]? Initial mass (for tracks only) (1)
real*4 Zini ! [0.0002/0.0260] Initial metallicity (2)
real*4 Vini ! [0.00/0.60] Initial rotation rate (3)
real*4 logT ! ([yr]) [6.0/10.1]? Initial age (for isochrones only) (4)
character*39 FileName ! Name of the file
*Note (1): Initial masses from 0.2 to 1.5, by 0.1 steps.
*Note (2): Initial metallicities are
* 0.0020, 0.0060, 0.0080, 0.0100, 0.0130, 0.0190 and 0.0260.
*Note (3): Initial rotation rates:
* 0.00 = standard
* 0.20 = slow rotators
* 0.40 = median rotators
* 0.60 = fast rotators
*Note (4): Initial log(age) from 6.0 to 10.1, by 0.1 steps.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tracks/*' ! *Evolution track of a model of mass (M), a
metallicity (Z), and an initial rotation rate (V)
integer*4 nr__1
parameter (nr__1=392) ! Number of records
character*1631 ar__1 ! Full-size record
integer*4 Line ! Line in the file (line)
real*4 Time ! (yr) Age (time)
real*8 logTeff ! ([K]) log10 of Effective temperature (logTeff)
real*8 logL ! ([Lsun]) log10 of Surface luminosity (logL)
real*8 logLgrav ! ([Lsun]) log10 of Gravothermal luminosity
* (logLgrav)
real*8 M ! (Msun) Mass (M)
real*8 R ! (Rsun) Photospheric radius (R)
real*4 logg ! ([cm/s2]) log10 effective gravity (logg)
real*4 rhophot ! (g/cm3) Photospheric density (rho_phot)
real*8 logMdot ! ([Ms/yr]) Mass-loss rate (logMdot)
real*8 logTc ! ([K]) log10 of central temperature (logTc)
real*8 logPc ! ([10-7J/cm3]) log10 of central pressure, in erg/cm3
* (logPc)
real*8 logrhoc ! ([g/cm3]) log10 of central density (logrhoc)
real*8 logTmax ! ([K]) log10 of maximal temperature (logTmax)
real*8 MrTmax ! (Msun) Mass coordinate of Tmax (Mr_Tmax)
real*8 logrhomax ! ([g/cm3]) log10 of maximal density (logrhomax)
real*4 epsnucl ! (10-7J/g/s) Nuclear energy production rate,
* in erg/g/s (eps_nucl)
real*4 epsgrav ! (10-7J/g/s) Gravothermal energy production rate,
* in erg/g/s (eps_grav)
real*4 epsnu ! (10-7J/g/s) Plasma neutrino energy loss rate,
* in erg/g/s (eps_nu)
real*8 MrbCE ! (Msun) Mass coordinate of the base of CE
* (Mr_b_CE)
real*4 normRbCE ! Normalised radius of the base of CE
* (normR_b_CE)
real*8 logTb ! ([K]) log10 of Temp. at the base of CE
* (logT_b)
real*8 logrhob ! ([g/cm3]) log10 of density at the base of CE
* (logrho_b)
real*8 MrtCC ! (Msun) Mass coordinate of the top of CC
* (Mr_t_CC)
real*4 normRtCC ! Normalised radius of the top of CC
* (normR_t_CC)
real*8 logTt ! ([K]) log10 of Temp. at the top of CC (logT_t)
real*8 logrhot ! ([g/cm3]) log10 of density at the top of CC
* (logrho_t)
real*4 taumax ! (yr) Maximal turnover timescale in the CE
* (tau_max)
real*4 Romax ! Corresponding Rossby number (Ro_max)
real*4 taug ! (yr) Integrated turnover timescale over CE
* (tau_g)
real*4 Rog ! Corresponding Rossby number (Ro_g)
real*4 tauHp2 ! (yr) Turnover timescale at Hp/2 above CE bot.
* (tau_Hp2)
real*4 RoHp2 ! Corresponding Rossby number (Ro_Hp2)
real*4 tauHp ! (yr) Turnover timescale at Hp above CE bot.
* (tau_Hp)
real*4 RoHp ! Corresponding Rossby number (Ro_Hp)
real*4 tauR2 ! (yr) Turnover timescale at mid CE in radius
* (tau_R2)
real*4 RoR2 ! Corresponding Rossby number (Ro_R2)
real*4 tauM2 ! (yr) Turnover timescale at mid CE in mass
* (tau_M2)
real*4 RoM2 ! Corresponding Rossby number (Ro_M2)
real*4 taumaxc ! (yr) Maximal turnover timescale in the CC
* (tau_max_c)
real*4 Romaxc ! Corresponding Rossby number (Ro_max_c)
real*4 taugc ! (yr) Integrated turnover timescale over CC
* (tau_g_c)
real*4 Rogc ! Corresponding Rossby number (Ro_g_c)
real*4 tauHp2c ! (yr) Turnover timescale at Hp/2 below CC top
* (tau_Hp2_c)
real*4 RoHp2c ! Corresponding Rossby number (Ro_Hp2_c)
real*4 tauHpc ! (yr) Turnover timescale at Hp below CC top
* (tau_Hp_c)
real*4 RoHpc ! Corresponding Rossby number (Ro_Hp_c)
real*4 tauR2c ! (yr) Turnover timescale at mid CC in radius
* (tau_R2_c)
real*4 RoR2c ! Corresponding Rossby number (Ro_R2_c)
real*4 tauM2c ! (yr) Turnover timescale at mid CC in mass
* (tau_M2_c)
real*4 RoM2c ! Corresponding Rossby number (Ro_M2_c)
real*4 k2conv ! (yr) Tot reduced moment of inertia of the CE
* (k2_conv)
real*4 k2rad ! Tot reduced moment of inertia below CE
* (k2_rad)
real*4 Omegas ! (rad/s) Surface angular velocity (Omega_s)
real*4 Omegac ! (rad/s) Mean angular velocity of radiative zone
* (Omega_c)
real*4 Vsurf ! (km/s) Surface rotation velocity (Vsurf)
real*4 Prot ! (d) Rotation period (Prot)
real*4 Jact ! (g.cm2/s) Total specific Ang. momentum of the star
* (J_act)
real*4 Jcore ! (g.cm2/s) Specific AM of the core (J_core) (G2)
real*4 OOc ! Ratio of the angular and critical veloc.
* (OOc)
real*4 Vcrit ! (km/s) Break-up surface velocity (Vcrit)
real*4 torque ! (g.cm2/s2) Total AM variation due to stellar winds
* (torque)
real*4 Bequi ! (G) Magnetic field assuming equipartition
* (B_equi) (G3)
real*4 Dnu ! (muHz) Large separation from asympt. relation
* (D_nu)
real*4 Dnuech ! (muHz) Large separation from scaling relation
* (D_nu_ech)
real*4 Dnuerr ! relative error on large separation
* (D_nu_err)
real*4 numax ! (muHz) Frequency with maximum amplitude (nu_max)
real*4 DPiasym ! (s) Asymptotic period spacing of g-modes
* (DPi_asym.)
real*4 Racctot ! (s) Total acoustic radius (R_acc_tot)
real*4 RaccBCE ! (s) Acoustic radius at the base of CE
* (R_acc_BCE)
real*4 RaccHe ! (s) Acoustic radius at the He 2nd ionisation
* (R_acc_He)
real*8 Mbol ! (mag) Bolometric magnitude (Mbol)
real*8 BC ! (mag) Bolometric corrections (BC)
real*8 U_B ! (mag) U-B colour index (U-B)
real*8 B_V ! (mag) B-V colour index (B-V)
real*8 V_R ! (mag) V-R colour index (V-R)
real*8 V_I ! (mag) V-I colour index (V-I)
real*8 J_K ! (mag) J-K colour index (J-K)
real*8 H_K ! (mag) H-K colour index (H-K)
real*8 V_K ! (mag) V-K colour index (V-K)
real*8 G_V ! (mag) G-V colour index (G-V)
real*8 Gbp_V ! (mag) Gbp-V colour index (Gbp-V)
real*8 Grp_V ! (mag) Grp-V colour index (Grp-V)
real*8 UMAG ! (mag) U absolute magnitude (M_U)
real*8 BMAG ! (mag) B absolute magnitude (M_B)
real*8 VMAG ! (mag) V absolute magnitude (M_V)
real*8 RMAG ! (mag) R absolute magnitude (M_R)
real*8 IMAG ! (mag) I absolute magnitude (M_I)
real*8 HMAG ! (mag) H absolute magnitude (M_H)
real*8 JMAG ! (mag) J absolute magnitude (M_J)
real*8 KMAG ! (mag) K absolute magnitude (M_K)
real*8 GMAG ! (mag) G absolute magnitude (M_G)
real*8 GbpMAG ! (mag) Gbp absolute magnitude (M_Gbp)
real*8 GrpMAG ! (mag) Grp absolute magnitude (M_Grp)
real*4 H1s ! Surface abundance of H (H1s) (G4)
real*4 H2s ! Surface abundance of Deuterium (H2s) (G4)
real*4 He3s ! Surface abundance of 3He (He3s) (G4)
real*4 He4s ! Surface abundance of 4He (He4s) (G4)
real*4 Li6s ! Surface abundance of 6Li (Li6s) (G4)
real*4 Li7s ! Surface abundance of 7Li (Li7s) (G4)
real*4 Be7s ! Surface abundance of 7Be (Be7s) (G4)
real*4 Be9s ! Surface abundance of 9Be (Be9s) (G4)
real*4 B10s ! Surface abundance of 10B (B10s) (G4)
real*4 B11s ! Surface abundance of 11B (B11s) (G4)
real*4 C12s ! Surface abundance of 12C (C12s) (G4)
real*4 C13s ! Surface abundance of 13C (C13s) (G4)
real*4 C14s ! Surface abundance of 14C (C14s) (G4)
real*4 N14s ! Surface abundance of 14N (N14s) (G4)
real*4 N15s ! Surface abundance of 15N (N15s) (G4)
real*4 O16s ! Surface abundance of 16O (O16s) (G4)
real*4 O17s ! Surface abundance of 17O (O17s) (G4)
real*4 O18s ! Surface abundance of 18O (O18s) (G4)
real*4 F19s ! Surface abundance of 19F (F19s) (G4)
real*4 Ne20s ! Surface abundance of 20Ne (Ne20s) (G4)
real*4 Ne21s ! Surface abundance of 21Ne (Ne21s) (G4)
real*4 Ne22s ! Surface abundance of 22Ne (Ne22s) (G4)
real*4 Na23s ! Surface abundance of 23Na (Na23s) (G4)
real*4 Mg24s ! Surface abundance of 24Mg (Mg24s) (G4)
real*4 Mg25s ! Surface abundance of 25Mg (Mg25s) (G4)
real*4 Mg26s ! Surface abundance of 26Mg (Mg26s) (G4)
real*4 Al26s ! Surface abundance of 26Al (Al26s) (G4)
real*4 Al27s ! Surface abundance of 27Al (Al27s) (G4)
real*4 Si28s ! Surface abundance of 28Si (Si28s) (G4)
real*4 H1c ! Central abundance of H (H1c) (G4)
real*4 H2c ! Central abundance of Deuterium (H2c) (G4)
real*4 He3c ! Central abundance of 3He (He3c) (G4)
real*4 He4c ! Central abundance of 4He (He4c) (G4)
real*4 C12c ! Central abundance of 12C (C12c) (G4)
real*4 C13c ! Central abundance of 13C (C13c) (G4)
real*4 C14c ! Central abundance of 14C (C14c) (G4)
real*4 N14c ! Central abundance of 14N (N14c) (G4)
real*4 N15c ! Central abundance of 15N (N15c) (G4)
real*4 O16c ! Central abundance of 16O (O16c) (G4)
real*4 O17c ! Central abundance of 17O (O17c) (G4)
real*4 O18c ! Central abundance of 18O (O18c) (G4)
real*4 F19c ! Central abundance of 19F (F19c) (G4)
real*4 Ne20c ! Central abundance of 20Ne (Ne20c) (G4)
real*4 Ne21c ! Central abundance of 21Ne (Ne21c) (G4)
real*4 Ne22c ! Central abundance of 22Ne (Ne22c) (G4)
real*4 Na23c ! Central abundance of 23Na (Na23c) (G4)
real*4 Mg24c ! Central abundance of 24Mg (Mg24c) (G4)
real*4 Mg25c ! Central abundance of 25Mg (Mg25c) (G4)
real*4 Mg26c ! Central abundance of 26Mg (Mg26c) (G4)
real*4 Al26c ! Central abundance of 26Al (Al26c) (G4)
real*4 Al27c ! Central abundance of 27Al (Al27c) (G4)
real*4 Si28c ! Central abundance of 28Si (Si28c) (G4)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'iso/*' ! *Isochrones for each (Z) and (V) at ages (lgT)
integer*4 nr__2
parameter (nr__2=1176) ! Number of records
character*1635 ar__2 ! Full-size record
real*8 Mini_1 ! (Msun) Mass (M_ini)
real*8 Zini_1 ! Metallicity (Z_ini)
real*4 OmOcini ! Initial velocity (OmOc_ini) (G1)
real*8 logTeff_1 ! ([K]) log10 of Effective temperature (logTeff)
real*8 logL_1 ! ([Lsun]) log10 of Surface luminosity (logL)
real*8 logLgrav_1 ! ([Lsun]) log10 of Gravothermal luminosity
* (logLgrav)
real*8 M_1 ! (Msun) Mass (M)
real*8 R_1 ! (Rsun) Photospheric radius (R)
real*4 logg_1 ! ([cm/s2]) log10 effective gravity (logg)
real*4 rhophot_1 ! (g/cm3) Photospheric density (rho_phot)
real*8 logdM_dt ! ([Ms/yr]) Mass-loss rate (logMdot)
real*8 logTc_1 ! ([K]) log10 of central temperature (logTc)
real*8 logPc_1 ! ([10-7J/cm3]) log10 of central pressure, in erg/cm3
* (logPc)
real*8 logrhoc_1 ! ([g/cm3]) log10 of central density (logrhoc)
real*8 logTmax_1 ! ([K]) log10 of maximal temperature (logTmax)
real*8 MrTmax_1 ! (Msun) Mass coordinate of Tmax (Mr_Tmax)
real*8 logrhomax_1 ! ([g/cm3]) log10 of maximal density (logrhomax)
real*4 epsnucl_1 ! (10-7J/g/s) Nuclear energy production rate,
* in erg/g/s (eps_nucl)
real*4 epsgrav_1 ! (10-7J/g/s) Gravothermal energy production rate,
* in erg/g/s (eps_grav)
real*4 epsnu_1 ! (10-7J/g/s) Plasma neutrino energy loss rate,
* in erg/g/s (eps_nu)
real*8 MrbCE_1 ! (Msun) Mass coordinate of the base of CE
* (Mr_b_CE)
real*4 normRbCE_1 ! Normalised radius of the base of CE
* (normR_b_CE)
real*8 logTb_1 ! ([K]) log10 of Temp. at the base of CE (logT_b)
real*8 logrhob_1 ! ([g/cm3]) log10 of density at the base of CE
* (logrho_b)
real*8 MrtCC_1 ! (Msun) Mass coordinate of the top of CC
* (Mr_t_CC)
real*4 normRtCC_1 ! Normalised radius of the top of CC
* (normR_t_CC)
real*8 logTt_1 ! ([K]) log10 of Temp. at the top of CC (logT_t)
real*8 logrhot_1 ! ([g/cm3]) log10 of density at the top of CC
* (logrho_t)
real*4 taumax_1 ! (yr) Maximal turnover timescale in the CE
* (tau_max)
real*4 Romax_1 ! Corresponding Rossby number (Ro_max)
real*4 taug_1 ! (yr) Integrated turnover timescale over CE
* (tau_g)
real*4 Rog_1 ! Corresponding Rossby number (Ro_g)
real*4 tauHp2_1 ! (yr) Turnover timescale at Hp/2 above CE bot.
* (tau_Hp2)
real*4 RoHp2_1 ! Corresponding Rossby number (Ro_Hp2)
real*4 tauHp_1 ! (yr) Turnover timescale at Hp above CE bot.
* (tau_Hp)
real*4 RoHp_1 ! Corresponding Rossby number (Ro_Hp)
real*4 tauR2_1 ! (yr) Turnover timescale at mid CE in radius
* (tau_R2)
real*4 RoR2_1 ! Corresponding Rossby number (Ro_R2)
real*4 tauM2_1 ! (yr) Turnover timescale at mid CE in mass
* (tau_M2)
real*4 RoM2_1 ! Corresponding Rossby number (Ro_M2)
real*4 taumaxc_1 ! (yr) Maximal turnover timescale in the CC
* (tau_max_c)
real*4 Romaxc_1 ! Corresponding Rossby number (Ro_max_c)
real*4 taugc_1 ! (yr) Integrated turnover timescale over CC
* (tau_g_c)
real*4 Rogc_1 ! Corresponding Rossby number (Ro_g_c)
real*4 tauHp2c_1 ! (yr) Turnover timescale at Hp/2 below CC top
* (tau_Hp2_c)
real*4 RoHp2c_1 ! Corresponding Rossby number (Ro_Hp2_c)
real*4 tauHpc_1 ! (yr) Turnover timescale at Hp below CC top
* (tau_Hp_c)
real*4 RoHpc_1 ! Corresponding Rossby number (Ro_Hp_c)
real*4 tauR2c_1 ! (yr) Turnover timescale at mid CC in radius
* (tau_R2_c)
real*4 RoR2c_1 ! Corresponding Rossby number (Ro_R2_c)
real*4 tauM2c_1 ! (yr) Turnover timescale at mid CC in mass
* (tau_M2_c)
real*4 RoM2c_1 ! Corresponding Rossby number (Ro_M2_c)
real*4 k2conv_1 ! (yr) Tot reduced moment of inertia of the CE
* (k2_conv)
real*4 k2rad_1 ! Tot reduced moment of inertia below CE
* (k2_rad)
real*4 Omegas_1 ! (rad/s) Surface angular velocity (Omega_s)
real*4 Omegac_1 ! (rad/s) Mean angular velocity of radiative zone
* (Omega_c)
real*4 Vsurf_1 ! (km/s) Surface rotation velocity (Vsurf)
real*4 Prot_1 ! (d) Rotation period (Prot)
real*4 Jact_1 ! (g.cm2/s) Total specific Ang. momentum of the star
* (J_act)
real*4 Jcore_1 ! (g.cm2/s) Specific AM of the core (J_core) (G2)
real*4 OOc_1 ! Ratio of the angular and break-up veloc.
* (OOc)
real*4 Vcrit_1 ! (km/s) Break-up surface velocity (Vcrit)
real*4 torque_1 ! (g.cm2/s2) Total AM variation due to stellar winds
* (torque)
real*4 Bequi_1 ! (G) Magnetic field assuming equipartition
* (B_equi) (G3)
real*4 Dnu_1 ! (muHz) Large separation from asympt. relation
* (D_nu)
real*4 Dnuech_1 ! (muHz) Large separation from scaling relation
* (D_nu_ech)
real*4 Dnuerr_1 ! relative error on large separation
* (D_nu_err)
real*4 numax_1 ! (muHz) Frequency with maximum amplitude (nu_max)
real*4 DPiasym_1 ! (s) Asymptotic period spacing of g-modes
* (DPi_asym.)
real*4 Racctot_1 ! (s) Total acoustic radius (R_acc_tot)
real*4 RaccBCE_1 ! (s) Acoustic radius at the base of CE
* (R_acc_BCE)
real*4 RaccHe_1 ! (s) Acoustic radius at the He 2nd ionisation
* (R_acc_He)
real*8 Mbol_1 ! (mag) Bolometric magnitude (Mbol)
real*8 BC_1 ! (mag) Bolometric corrections (BC)
real*8 U_B_1 ! (mag) U-B colour index (U-B)
real*8 B_V_1 ! (mag) B-V colour index (B-V)
real*8 V_R_1 ! (mag) V-R colour index (V-R)
real*8 V_I_1 ! (mag) V-I colour index (V-I)
real*8 J_K_1 ! (mag) J-K colour index (J-K)
real*8 H_K_1 ! (mag) H-K colour index (H-K)
real*8 V_K_1 ! (mag) V-K colour index (V-K)
real*8 G_V_1 ! (mag) G-V colour index (G-V)
real*8 Gbp_V_1 ! (mag) Gbp-V colour index (Gbp-V)
real*8 Grp_V_1 ! (mag) Grp-V colour index (Grp-V)
real*8 UMAG_1 ! (mag) Absolute U magnitude (M_U)
real*8 BMAG_1 ! (mag) Absolute B magnitude (M_B)
real*8 VMAG_1 ! (mag) Absolute V magnitude (M_V)
real*8 RMAG_1 ! (mag) Absolute R magnitude (M_R)
real*8 IMAG_1 ! (mag) Absolute I magnitude (M_I)
real*8 HMAG_1 ! (mag) Absolute H magnitude (M_H)
real*8 JMAG_1 ! (mag) Absolute J magnitude (M_J)
real*8 KMAG_1 ! (mag) Absolute K magnitude (M_K)
real*8 GMAG_1 ! (mag) Absolute G magnitude (M_G)
real*8 GbpMAG_1 ! (mag) Absolute Gbp magnitude (M_Gbp)
real*8 GrpMAG_1 ! (mag) Absolute Grp magnitude (M_Grp)
real*4 H1s_1 ! Surface abundance of H (H1s) (G4)
real*4 H2s_1 ! Surface abundance of Deuterium (H2s) (G4)
real*4 He3s_1 ! Surface abundance of ^3^He (He3s) (G4)
real*4 He4s_1 ! Surface abundance of ^4^He (He4s) (G4)
real*4 Li6s_1 ! Surface abundance of ^6^Li (Li6s) (G4)
real*4 Li7s_1 ! Surface abundance of ^7^Li (Li7s) (G4)
real*4 Be7s_1 ! Surface abundance of ^7^Be (Be7s) (G4)
real*4 Be9s_1 ! Surface abundance of ^9^Be (Be9s) (G4)
real*4 B10s_1 ! Surface abundance of ^10^B (B10s) (G4)
real*4 B11s_1 ! Surface abundance of ^11^B (B11s) (G4)
real*4 C12s_1 ! Surface abundance of ^12^C (C12s) (G4)
real*4 C13s_1 ! Surface abundance of ^13^C (C13s) (G4)
real*4 C14s_1 ! Surface abundance of ^14^C (C14s) (G4)
real*4 N14s_1 ! Surface abundance of ^14^N (N14s) (G4)
real*4 N15s_1 ! Surface abundance of ^15^N (N15s) (G4)
real*4 O16s_1 ! Surface abundance of ^16^O (O16s) (G4)
real*4 O17s_1 ! Surface abundance of ^17^O (O17s) (G4)
real*4 O18s_1 ! Surface abundance of ^18^O (O18s) (G4)
real*4 F19s_1 ! Surface abundance of ^19^F (F19s) (G4)
real*4 Ne20s_1 ! Surface abundance of ^20^Ne (Ne20s) (G4)
real*4 Ne21s_1 ! Surface abundance of ^21^Ne (Ne21s) (G4)
real*4 Ne22s_1 ! Surface abundance of ^22^Ne (Ne22s) (G4)
real*4 Na23s_1 ! Surface abundance of ^23^Na (Na23s) (G4)
real*4 Mg24s_1 ! Surface abundance of ^24^Mg (Mg24s) (G4)
real*4 Mg25s_1 ! Surface abundance of ^25^Mg (Mg25s) (G4)
real*4 Mg26s_1 ! Surface abundance of ^26^Mg (Mg26s) (G4)
real*4 Al26s_1 ! Surface abundance of ^26^Al (Al26s) (G4)
real*4 Al27s_1 ! Surface abundance of ^27^Al (Al27s) (G4)
real*4 Si28s_1 ! Surface abundance of ^28^Si (Si28s) (G4)
real*4 H1c_1 ! Central abundance of H (H1c) (G4)
real*4 H2c_1 ! Central abundance of Deuterium (H2c) (G4)
real*4 He3c_1 ! Central abundance of ^3H^e (He3c) (G4)
real*4 He4c_1 ! Central abundance of ^4H^e (He4c) (G4)
real*4 C12c_1 ! Central abundance of ^12^C (C12c) (G4)
real*4 C13c_1 ! Central abundance of ^13^C (C13c) (G4)
real*4 C14c_1 ! Central abundance of ^14^C (C14c) (G4)
real*4 N14c_1 ! Central abundance of ^14^N (N14c) (G4)
real*4 N15c_1 ! Central abundance of ^15^N (N15c) (G4)
real*4 O16c_1 ! Central abundance of ^16^O (O16c) (G4)
real*4 O17c_1 ! Central abundance of ^17^O (O17c) (G4)
real*4 O18c_1 ! Central abundance of ^18^O (O18c) (G4)
real*4 F19c_1 ! Central abundance of ^19^F (F19c) (G4)
real*4 Ne20c_1 ! Central abundance of ^20^Ne (Ne20c) (G4)
real*4 Ne21c_1 ! Central abundance of ^21^Ne (Ne21c) (G4)
real*4 Ne22c_1 ! Central abundance of ^22^Ne (Ne22c) (G4)
real*4 Na23c_1 ! Central abundance of ^23^Na (Na23c) (G4)
real*4 Mg24c_1 ! Central abundance of ^24^Mg (Mg24c) (G4)
real*4 Mg25c_1 ! Central abundance of ^25^Mg (Mg25c) (G4)
real*4 Mg26c_1 ! Central abundance of ^26^Mg (Mg26c) (G4)
real*4 Al26c_1 ! Central abundance of ^26^Al (Al26c) (G4)
real*4 Al27c_1 ! Central abundance of ^27^Al (Al27c) (G4)
real*4 Si28c_1 ! Central abundance of ^28^Si (Si28c) (G4)
C=============================================================================
C Loading file 'list.dat' ! List of all the tracks and isochrone files
C Format for file interpretation
1 format(1X,F4.2,2X,F6.4,2X,F4.2,2X,F6.3,1X,A39)
C Effective file loading
open(unit=1,status='old',file=
+'list.dat')
write(6,*) '....Loading file: list.dat'
do i__=1,1568
read(1,'(A67)')ar__
read(ar__,1)Mini,Zini,Vini,logT,FileName
if(ar__(2:5) .EQ. '') Mini = rNULL__
if(ar__(22:27) .EQ. '') logT = rNULL__
c ..............Just test output...........
write(6,1)Mini,Zini,Vini,logT,FileName
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tracks/*' ! *Evolution track of a model of mass (M), a
* metallicity (Z), and an initial rotation rate (V)
C Format for file interpretation
2 format(
+ I4,2X,E16.10,1X,F9.6,1X,F9.6,1X,F9.6,4X,F10.8,3X,F9.5,1X,F5.2,
+ 2X,E9.3,1X,F10.6,2X,F8.6,1X,F9.6,1X,F9.6,1X,F9.6,1X,F9.6,1X,
+ F9.6,1X,E11.4,1X,E11.4,1X,E11.4,2X,F9.7,2X,E9.3,2X,F7.5,1X,
+ F8.5,2X,F9.7,2X,E9.3,1X,F8.5,1X,F8.5,2X,E10.4,1X,E11.4,1X,
+ E11.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,
+ 2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,
+ E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,
+ 2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E8.2,
+ 2X,E10.4,2X,E8.2,2X,E8.2,2X,E10.4,2X,E8.2,1X,E11.4,2X,E8.2,1X,
+ E11.4,1X,E11.4,1X,E11.4,1X,E11.4,1X,E11.4,1X,E11.4,1X,E11.4,
+ 1X,E11.4,2X,F7.4,2X,F7.4,2X,F7.4,2X,F7.4,2X,F7.4,2X,F7.4,2X,
+ F7.4,2X,F7.4,2X,F7.4,2X,F7.4,2X,F7.4,2X,F7.4,2X,F7.4,2X,F7.4,
+ 2X,F7.4,2X,F7.4,2X,F7.4,2X,F7.4,2X,F7.4,2X,F7.4,2X,F7.4,2X,
+ F7.4,2X,F7.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,
+ E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,
+ 2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,
+ E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,
+ 2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,
+ E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,
+ 2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,
+ E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,
+ 2X,E10.4)
C Effective file loading
open(unit=1,status='old',file=
+'tracks/*')
write(6,*) '....Loading file: tracks/*'
do i__=1,392
read(1,'(A1631)')ar__1
read(ar__1,2)
+ Line,Time,logTeff,logL,logLgrav,M,R,logg,rhophot,logMdot,
+ logTc,logPc,logrhoc,logTmax,MrTmax,logrhomax,epsnucl,epsgrav,
+ epsnu,MrbCE,normRbCE,logTb,logrhob,MrtCC,normRtCC,logTt,
+ logrhot,taumax,Romax,taug,Rog,tauHp2,RoHp2,tauHp,RoHp,tauR2,
+ RoR2,tauM2,RoM2,taumaxc,Romaxc,taugc,Rogc,tauHp2c,RoHp2c,
+ tauHpc,RoHpc,tauR2c,RoR2c,tauM2c,RoM2c,k2conv,k2rad,Omegas,
+ Omegac,Vsurf,Prot,Jact,Jcore,OOc,Vcrit,torque,Bequi,Dnu,
+ Dnuech,Dnuerr,numax,DPiasym,Racctot,RaccBCE,RaccHe,Mbol,BC,
+ U_B,B_V,V_R,V_I,J_K,H_K,V_K,G_V,Gbp_V,Grp_V,UMAG,BMAG,VMAG,
+ RMAG,IMAG,HMAG,JMAG,KMAG,GMAG,GbpMAG,GrpMAG,H1s,H2s,He3s,He4s,
+ Li6s,Li7s,Be7s,Be9s,B10s,B11s,C12s,C13s,C14s,N14s,N15s,O16s,
+ O17s,O18s,F19s,Ne20s,Ne21s,Ne22s,Na23s,Mg24s,Mg25s,Mg26s,
+ Al26s,Al27s,Si28s,H1c,H2c,He3c,He4c,C12c,C13c,C14c,N14c,N15c,
+ O16c,O17c,O18c,F19c,Ne20c,Ne21c,Ne22c,Na23c,Mg24c,Mg25c,Mg26c,
+ Al26c,Al27c,Si28c
c ..............Just test output...........
write(6,2)
+ Line,Time,logTeff,logL,logLgrav,M,R,logg,rhophot,logMdot,
+ logTc,logPc,logrhoc,logTmax,MrTmax,logrhomax,epsnucl,epsgrav,
+ epsnu,MrbCE,normRbCE,logTb,logrhob,MrtCC,normRtCC,logTt,
+ logrhot,taumax,Romax,taug,Rog,tauHp2,RoHp2,tauHp,RoHp,tauR2,
+ RoR2,tauM2,RoM2,taumaxc,Romaxc,taugc,Rogc,tauHp2c,RoHp2c,
+ tauHpc,RoHpc,tauR2c,RoR2c,tauM2c,RoM2c,k2conv,k2rad,Omegas,
+ Omegac,Vsurf,Prot,Jact,Jcore,OOc,Vcrit,torque,Bequi,Dnu,
+ Dnuech,Dnuerr,numax,DPiasym,Racctot,RaccBCE,RaccHe,Mbol,BC,
+ U_B,B_V,V_R,V_I,J_K,H_K,V_K,G_V,Gbp_V,Grp_V,UMAG,BMAG,VMAG,
+ RMAG,IMAG,HMAG,JMAG,KMAG,GMAG,GbpMAG,GrpMAG,H1s,H2s,He3s,He4s,
+ Li6s,Li7s,Be7s,Be9s,B10s,B11s,C12s,C13s,C14s,N14s,N15s,O16s,
+ O17s,O18s,F19s,Ne20s,Ne21s,Ne22s,Na23s,Mg24s,Mg25s,Mg26s,
+ Al26s,Al27s,Si28s,H1c,H2c,He3c,He4c,C12c,C13c,C14c,N14c,N15c,
+ O16c,O17c,O18c,F19c,Ne20c,Ne21c,Ne22c,Na23c,Mg24c,Mg25c,Mg26c,
+ Al26c,Al27c,Si28c
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'iso/*' ! *Isochrones for each (Z) and (V) at ages (lgT)
C Format for file interpretation
3 format(
+ F7.3,2X,F8.6,2X,F5.3,3X,F9.6,1X,F9.6,1X,F9.6,4X,F10.8,3X,F9.5,
+ 1X,F5.2,2X,E9.3,1X,F10.6,2X,F8.6,1X,F9.6,1X,F9.6,1X,F9.6,1X,
+ F9.6,1X,F9.6,1X,E11.4,1X,E11.4,1X,E11.4,2X,F9.7,2X,E9.3,2X,
+ F7.5,1X,F8.5,2X,F9.7,2X,E9.3,1X,F8.5,1X,F8.5,2X,E10.4,2X,
+ E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,
+ 2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,
+ E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,
+ 2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,
+ E10.4,2X,E8.2,2X,E10.4,2X,E8.2,2X,E8.2,2X,E10.4,2X,E8.2,1X,
+ E11.4,2X,E8.2,1X,E11.4,1X,E11.4,1X,E11.4,1X,E11.4,1X,E11.4,1X,
+ E11.4,1X,E11.4,1X,E11.4,2X,F7.4,2X,F7.4,2X,F7.4,2X,F7.4,2X,
+ F7.4,2X,F7.4,2X,F7.4,2X,F7.4,2X,F7.4,2X,F7.4,2X,F7.4,2X,F7.4,
+ 2X,F7.4,2X,F7.4,2X,F7.4,2X,F7.4,2X,F7.4,2X,F7.4,2X,F7.4,2X,
+ F7.4,2X,F7.4,2X,F7.4,2X,F7.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,
+ E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,
+ 2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,
+ E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,
+ 2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,
+ E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,
+ 2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,
+ E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,2X,E10.4,
+ 2X,E10.4,2X,E10.4,2X,E10.4)
C Effective file loading
open(unit=1,status='old',file=
+'iso/*')
write(6,*) '....Loading file: iso/*'
do i__=1,1176
read(1,'(A1635)')ar__2
read(ar__2,3)
+ Mini_1,Zini_1,OmOcini,logTeff_1,logL_1,logLgrav_1,M_1,R_1,
+ logg_1,rhophot_1,logdM_dt,logTc_1,logPc_1,logrhoc_1,logTmax_1,
+ MrTmax_1,logrhomax_1,epsnucl_1,epsgrav_1,epsnu_1,MrbCE_1,
+ normRbCE_1,logTb_1,logrhob_1,MrtCC_1,normRtCC_1,logTt_1,
+ logrhot_1,taumax_1,Romax_1,taug_1,Rog_1,tauHp2_1,RoHp2_1,
+ tauHp_1,RoHp_1,tauR2_1,RoR2_1,tauM2_1,RoM2_1,taumaxc_1,
+ Romaxc_1,taugc_1,Rogc_1,tauHp2c_1,RoHp2c_1,tauHpc_1,RoHpc_1,
+ tauR2c_1,RoR2c_1,tauM2c_1,RoM2c_1,k2conv_1,k2rad_1,Omegas_1,
+ Omegac_1,Vsurf_1,Prot_1,Jact_1,Jcore_1,OOc_1,Vcrit_1,torque_1,
+ Bequi_1,Dnu_1,Dnuech_1,Dnuerr_1,numax_1,DPiasym_1,Racctot_1,
+ RaccBCE_1,RaccHe_1,Mbol_1,BC_1,U_B_1,B_V_1,V_R_1,V_I_1,J_K_1,
+ H_K_1,V_K_1,G_V_1,Gbp_V_1,Grp_V_1,UMAG_1,BMAG_1,VMAG_1,RMAG_1,
+ IMAG_1,HMAG_1,JMAG_1,KMAG_1,GMAG_1,GbpMAG_1,GrpMAG_1,H1s_1,
+ H2s_1,He3s_1,He4s_1,Li6s_1,Li7s_1,Be7s_1,Be9s_1,B10s_1,B11s_1,
+ C12s_1,C13s_1,C14s_1,N14s_1,N15s_1,O16s_1,O17s_1,O18s_1,
+ F19s_1,Ne20s_1,Ne21s_1,Ne22s_1,Na23s_1,Mg24s_1,Mg25s_1,
+ Mg26s_1,Al26s_1,Al27s_1,Si28s_1,H1c_1,H2c_1,He3c_1,He4c_1,
+ C12c_1,C13c_1,C14c_1,N14c_1,N15c_1,O16c_1,O17c_1,O18c_1,
+ F19c_1,Ne20c_1,Ne21c_1,Ne22c_1,Na23c_1,Mg24c_1,Mg25c_1,
+ Mg26c_1,Al26c_1,Al27c_1,Si28c_1
c ..............Just test output...........
write(6,3)
+ Mini_1,Zini_1,OmOcini,logTeff_1,logL_1,logLgrav_1,M_1,R_1,
+ logg_1,rhophot_1,logdM_dt,logTc_1,logPc_1,logrhoc_1,logTmax_1,
+ MrTmax_1,logrhomax_1,epsnucl_1,epsgrav_1,epsnu_1,MrbCE_1,
+ normRbCE_1,logTb_1,logrhob_1,MrtCC_1,normRtCC_1,logTt_1,
+ logrhot_1,taumax_1,Romax_1,taug_1,Rog_1,tauHp2_1,RoHp2_1,
+ tauHp_1,RoHp_1,tauR2_1,RoR2_1,tauM2_1,RoM2_1,taumaxc_1,
+ Romaxc_1,taugc_1,Rogc_1,tauHp2c_1,RoHp2c_1,tauHpc_1,RoHpc_1,
+ tauR2c_1,RoR2c_1,tauM2c_1,RoM2c_1,k2conv_1,k2rad_1,Omegas_1,
+ Omegac_1,Vsurf_1,Prot_1,Jact_1,Jcore_1,OOc_1,Vcrit_1,torque_1,
+ Bequi_1,Dnu_1,Dnuech_1,Dnuerr_1,numax_1,DPiasym_1,Racctot_1,
+ RaccBCE_1,RaccHe_1,Mbol_1,BC_1,U_B_1,B_V_1,V_R_1,V_I_1,J_K_1,
+ H_K_1,V_K_1,G_V_1,Gbp_V_1,Grp_V_1,UMAG_1,BMAG_1,VMAG_1,RMAG_1,
+ IMAG_1,HMAG_1,JMAG_1,KMAG_1,GMAG_1,GbpMAG_1,GrpMAG_1,H1s_1,
+ H2s_1,He3s_1,He4s_1,Li6s_1,Li7s_1,Be7s_1,Be9s_1,B10s_1,B11s_1,
+ C12s_1,C13s_1,C14s_1,N14s_1,N15s_1,O16s_1,O17s_1,O18s_1,
+ F19s_1,Ne20s_1,Ne21s_1,Ne22s_1,Na23s_1,Mg24s_1,Mg25s_1,
+ Mg26s_1,Al26s_1,Al27s_1,Si28s_1,H1c_1,H2c_1,He3c_1,He4c_1,
+ C12c_1,C13c_1,C14c_1,N14c_1,N15c_1,O16c_1,O17c_1,O18c_1,
+ F19c_1,Ne20c_1,Ne21c_1,Ne22c_1,Na23c_1,Mg24c_1,Mg25c_1,
+ Mg26c_1,Al26c_1,Al27c_1,Si28c_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end