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