Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/605/A102 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-Apr-20
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/605/A102 Stellar models. 0.85<M<6, Z=0.0001-0.014 (Charbonnel+, 2017)
*================================================================================
*The magnetic strip(s) in the advanced phases of stellar evolution.
*Theoretical convective turnover timescale and Rossby number for low- and
*intermediate-mass stars up to the AGB at various metallicities.
* Charbonnel C., Decressin T., Lagarde N., Gallet F., Palacios A.,
* Auriere M., Konstantinova-Antova R., Mathis S., Anderson R.I., Dintrans B.
* <Astron. Astrophys., 605, A102 (2017)>
* =2017A&A...605A.102C (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'z0001.dat' ! Grid of stellar models for Z=0.0001
integer*4 nr__
parameter (nr__=7252) ! Number of records
character*900 ar__ ! Full-size record
real*4 Z (nr__) ! Initial Z abundance
character*1 T (nr__) ! [rs] s: standard prescriptions or
* r: th+rot prescriptions
real*4 Mass (nr__) ! (Msun) [0.85/6] Initial mass
* (0.85, 1.0, 1.5, 2.0, 2.5, 3.0, 4.0, 6.0)
integer*4 Mod (nr__) ! Model number
real*8 L (nr__) ! (Lsun) Surface luminosity
real*4 Reff (nr__) ! (Rsun) Photospheric radius
integer*4 Teff (nr__) ! (K) Effective temperature
real*4 rhoeff (nr__) ! (g/cm3) Photospheric density
real*4 log_g (nr__) ! ([cm/s2]) Photospheric gravity
real*4 Mloss (nr__) ! (Msun/yr) Mass-loss rate
real*8 M (nr__) ! (Msun) Stellar mass
real*4 t_1 (nr__) ! (yr) Age
real*4 tchp_2 (nr__) ! (d) Convective turnover timescale calculated at
* half of the pressure scale height
real*4 tchp (nr__) ! (d) Convective turnover timescale calculated at
* the pressure scale height
real*4 tcrconv_2 (nr__) ! (d) Convective turnover timescale calculated at
* half of the radius of the convective envelope
real*4 tcmconv_2 (nr__) ! (d) Convective turnover timescale calculated at
* half of the mass of the convective envelope
real*4 tcmax (nr__) ! (d) Maximum convective turnover timescale
real*4 tg (nr__) ! (d) Global convective turnover timescale
real*4 rhp_2 (nr__) ! (Rsun) Radius at which tc_hp/2 was estimated
real*4 rhp (nr__) ! (Rsun) Radius at which tc_hp was estimated
real*4 rrconv_2 (nr__) ! (Rsun) Radius at which tc_rconv/2 was estimated
real*4 rmconv_2 (nr__) ! (Rsun) Radius at which tc_mconv/2 was estimated
real*4 rmax (nr__) ! (Rsun) Radius at which tc_max was estimated
real*4 Rohp_2 (nr__) ! Rossby number associated to tc_hp/2
real*4 Rohp (nr__) ! Rossby number associated to tc_hp
real*4 Rorconv_2 (nr__) ! Rossby number associated to tc_rconv/2
real*4 Romconv_2 (nr__) ! Rossby number associated to tc_mconv/2
real*4 Romax (nr__) ! Rossby number associated to tc_max
real*4 Rog (nr__) ! Rossby number associated to tg
real*4 Tc (nr__) ! (K) Central temperature
real*4 Tmax (nr__) ! (K) Maximum of temperature
real*4 MrTmax (nr__) ! (Msun) Mass coordinate of Tmax
real*8 rhoc (nr__) ! (g/cm3) Central density
real*4 rhomax (nr__) ! (g/cm3) Density at the location of Tmax
real*4 Pc (nr__) ! (dPa) Central pressure (cgs)
real*4 Mbenv (nr__) ! (Msun) Mass at the base of convective envelope
real*8 XHc (nr__) ! Central abundance of H (mass fraction)
real*4 XHe3c (nr__) ! Central abundance of 3He (mass fraction)
real*4 XHe4c (nr__) ! Central abundance of 4He (mass fraction)
real*4 XC12c (nr__) ! Central abundance of 12C (mass fraction)
real*4 XC13c (nr__) ! Central abundance of 13C (mass fraction)
real*4 XC14c (nr__) ! Central abundance of 14C (mass fraction)
real*4 XN14c (nr__) ! Central abundance of 14N (mass fraction)
real*4 XO16c (nr__) ! Central abundance of 16O (mass fraction)
real*4 XO17c (nr__) ! Central abundance of 17O (mass fraction)
real*4 XO18c (nr__) ! Central abundance of 18O (mass fraction)
real*4 XF19c (nr__) ! Central abundance of 19F (mass fraction)
real*4 XNe20c (nr__) ! Central abundance of 20Ne (mass fraction)
real*4 XNe21c (nr__) ! Central abundance of 21Ne (mass fraction)
real*4 XNe22c (nr__) ! Central abundance of 22Ne (mass fraction)
real*4 XNa23c (nr__) ! Central abundance of 23Na (mass fraction)
real*4 XMg24c (nr__) ! Central abundance of 24Mg (mass fraction)
real*4 XMg25c (nr__) ! Central abundance of 25Mg (mass fraction)
real*4 XMg26c (nr__) ! Central abundance of 26Mg (mass fraction)
real*4 XAl26c (nr__) ! Central abundance of 26Al (mass fraction)
real*4 XAl27c (nr__) ! Central abundance of 27Al (mass fraction)
real*4 XSi28c (nr__) ! Central abundance of 28Si (mass fraction)
real*4 XHs (nr__) ! Surface abundance of H (mass fraction)
real*4 XH2s (nr__) ! Surface abundance of 2H (mass fraction)
real*4 XHe3s (nr__) ! Surface abundance of 3He (mass fraction)
real*4 XHe4s (nr__) ! Surface abundance of 4He (mass fraction)
real*4 XLi6s (nr__) ! Surface abundance of 6Li (mass fraction)
real*4 XLi7s (nr__) ! Surface abundance of 7Li (mass fraction)
real*4 XBe7s (nr__) ! Surface abundance of 7Be (mass fraction)
real*4 XBe9s (nr__) ! Surface abundance of 9Be (mass fraction)
real*4 XB10s (nr__) ! Surface abundance of 10B (mass fraction)
real*4 XB11s (nr__) ! Surface abundance of 11B (mass fraction)
real*4 XC12s (nr__) ! Surface abundance of 12C (mass fraction)
real*4 XC13s (nr__) ! Surface abundance of 13C (mass fraction)
real*4 XC14s (nr__) ! Surface abundance of 14C (mass fraction)
real*4 XN14s (nr__) ! Surface abundance of 14N (mass fraction)
real*4 XO16s (nr__) ! Surface abundance of 16O (mass fraction)
real*4 XO17s (nr__) ! Surface abundance of 17O (mass fraction)
real*4 XO18s (nr__) ! Surface abundance of 18O (mass fraction)
real*4 XF19s (nr__) ! Surface abundance of 19F (mass fraction)
real*4 XNe20s (nr__) ! Surface abundance of 20Ne (mass fraction)
real*4 XNe21s (nr__) ! Surface abundance of 21Ne (mass fraction)
real*4 XNe22s (nr__) ! Surface abundance of 22Ne (mass fraction)
real*4 XNa23s (nr__) ! Surface abundance of 23Na (mass fraction)
real*4 XMg24s (nr__) ! Surface abundance of 24Mg (mass fraction)
real*4 XMg25s (nr__) ! Surface abundance of 25Mg (mass fraction)
real*4 XMg26s (nr__) ! Surface abundance of 26Mg (mass fraction)
real*4 XAl26s (nr__) ! Surface abundance of 26Al (mass fraction)
real*4 XAl27s (nr__) ! Surface abundance of 27Al (mass fraction)
real*4 XSi28s (nr__) ! Surface abundance of 28Si (mass fraction)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'z002.dat' ! Grid of stellar models for Z=0.002
integer*4 nr__1
parameter (nr__1=7227) ! Number of records
character*900 ar__1 ! Full-size record
real*4 Z_1 (nr__1) ! Initial Z abundance
character*1 T_2 (nr__1) ! [rs] s: standard prescriptions or
* r: th+rot prescriptions
real*4 Mass_1 (nr__1) ! (Msun) [0.85/6] Initial mass
* (0.85, 1.0, 1.5, 2.0, 2.5, 3.0, 4.0, 6.0)
integer*4 Mod_1 (nr__1) ! Model number
real*8 L_1 (nr__1) ! (Lsun) Surface luminosity
real*4 Reff_1 (nr__1) ! (Rsun) Photospheric radius
integer*4 Teff_1 (nr__1) ! (K) Effective temperature
real*4 rhoeff_1 (nr__1) ! (g/cm3) Photospheric density
real*4 log_g_1 (nr__1) ! ([cm/s2]) Photospheric gravity
real*4 Mloss_1 (nr__1) ! (Msun/yr) Mass-loss rate
real*8 M_1 (nr__1) ! (Msun) Stellar mass
real*4 t_3 (nr__1) ! (yr) Age
real*4 tchp_2_1 (nr__1) ! (d) Convective turnover timescale calculated at
* half of the pressure scale height
real*4 tchp_1 (nr__1) ! (d) Convective turnover timescale calculated at
* the pressure scale height
real*4 tcrconv_2_1(nr__1) ! (d) Convective turnover timescale calculated at
* half of the radius of the convective envelope
real*4 tcmconv_2_1(nr__1) ! (d) Convective turnover timescale calculated at
* half of the mass of the convective envelope
real*4 tcmax_1 (nr__1) ! (d) Maximum convective turnover timescale
real*4 tg_1 (nr__1) ! (d) Global convective turnover timescale
real*4 rhp_2_1 (nr__1) ! (Rsun) Radius at which tc_hp/2 was estimated
real*4 rhp_1 (nr__1) ! (Rsun) Radius at which tc_hp was estimated
real*4 rrconv_2_1 (nr__1) ! (Rsun) Radius at which tc_rconv/2 was estimated
real*4 rmconv_2_1 (nr__1) ! (Rsun) Radius at which tc_mconv/2 was estimated
real*4 rmax_1 (nr__1) ! (Rsun) Radius at which tc_max was estimated
real*4 Rohp_2_1 (nr__1) ! Rossby number associated to tc_hp/2
real*4 Rohp_1 (nr__1) ! Rossby number associated to tc_hp
real*4 Rorconv_2_1(nr__1) ! Rossby number associated to tc_rconv/2
real*4 Romconv_2_1(nr__1) ! Rossby number associated to tc_mconv/2
real*4 Romax_1 (nr__1) ! Rossby number associated to tc_max
real*4 Rog_1 (nr__1) ! Rossby number associated to tg
real*4 Tc_1 (nr__1) ! (K) Central temperature
real*4 Tmax_1 (nr__1) ! (K) Maximum of temperature
real*4 MrTmax_1 (nr__1) ! (Msun) Mass coordinate of Tmax
real*8 rhoc_1 (nr__1) ! (g/cm3) Central density
real*4 rhomax_1 (nr__1) ! (g/cm3) Density at the location of Tmax
real*4 Pc_1 (nr__1) ! (dPa) Central pressure (cgs)
real*4 Mbenv_1 (nr__1) ! (Msun) Mass at the base of convective envelope
real*8 XHc_1 (nr__1) ! Central abundance of H (mass fraction)
real*4 XHe3c_1 (nr__1) ! Central abundance of 3He (mass fraction)
real*4 XHe4c_1 (nr__1) ! Central abundance of 4He (mass fraction)
real*4 XC12c_1 (nr__1) ! Central abundance of 12C (mass fraction)
real*4 XC13c_1 (nr__1) ! Central abundance of 13C (mass fraction)
real*4 XC14c_1 (nr__1) ! Central abundance of 14C (mass fraction)
real*4 XN14c_1 (nr__1) ! Central abundance of 14N (mass fraction)
real*4 XO16c_1 (nr__1) ! Central abundance of 16O (mass fraction)
real*4 XO17c_1 (nr__1) ! Central abundance of 17O (mass fraction)
real*4 XO18c_1 (nr__1) ! Central abundance of 18O (mass fraction)
real*4 XF19c_1 (nr__1) ! Central abundance of 19F (mass fraction)
real*4 XNe20c_1 (nr__1) ! Central abundance of 20Ne (mass fraction)
real*4 XNe21c_1 (nr__1) ! Central abundance of 21Ne (mass fraction)
real*4 XNe22c_1 (nr__1) ! Central abundance of 22Ne (mass fraction)
real*4 XNa23c_1 (nr__1) ! Central abundance of 23Na (mass fraction)
real*4 XMg24c_1 (nr__1) ! Central abundance of 24Mg (mass fraction)
real*4 XMg25c_1 (nr__1) ! Central abundance of 25Mg (mass fraction)
real*4 XMg26c_1 (nr__1) ! Central abundance of 26Mg (mass fraction)
real*4 XAl26c_1 (nr__1) ! Central abundance of 26Al (mass fraction)
real*4 XAl27c_1 (nr__1) ! Central abundance of 27Al (mass fraction)
real*4 XSi28c_1 (nr__1) ! Central abundance of 28Si (mass fraction)
real*4 XHs_1 (nr__1) ! Surface abundance of H (mass fraction)
real*4 XH2s_1 (nr__1) ! Surface abundance of 2H (mass fraction)
real*4 XHe3s_1 (nr__1) ! Surface abundance of 3He (mass fraction)
real*4 XHe4s_1 (nr__1) ! Surface abundance of 4He (mass fraction)
real*4 XLi6s_1 (nr__1) ! Surface abundance of 6Li (mass fraction)
real*4 XLi7s_1 (nr__1) ! Surface abundance of 7Li (mass fraction)
real*4 XBe7s_1 (nr__1) ! Surface abundance of 7Be (mass fraction)
real*4 XBe9s_1 (nr__1) ! Surface abundance of 9Be (mass fraction)
real*4 XB10s_1 (nr__1) ! Surface abundance of 10B (mass fraction)
real*4 XB11s_1 (nr__1) ! Surface abundance of 11B (mass fraction)
real*4 XC12s_1 (nr__1) ! Surface abundance of 12C (mass fraction)
real*4 XC13s_1 (nr__1) ! Surface abundance of 13C (mass fraction)
real*4 XC14s_1 (nr__1) ! Surface abundance of 14C (mass fraction)
real*4 XN14s_1 (nr__1) ! Surface abundance of 14N (mass fraction)
real*4 XO16s_1 (nr__1) ! Surface abundance of 16O (mass fraction)
real*4 XO17s_1 (nr__1) ! Surface abundance of 17O (mass fraction)
real*4 XO18s_1 (nr__1) ! Surface abundance of 18O (mass fraction)
real*4 XF19s_1 (nr__1) ! Surface abundance of 19F (mass fraction)
real*4 XNe20s_1 (nr__1) ! Surface abundance of 20Ne (mass fraction)
real*4 XNe21s_1 (nr__1) ! Surface abundance of 21Ne (mass fraction)
real*4 XNe22s_1 (nr__1) ! Surface abundance of 22Ne (mass fraction)
real*4 XNa23s_1 (nr__1) ! Surface abundance of 23Na (mass fraction)
real*4 XMg24s_1 (nr__1) ! Surface abundance of 24Mg (mass fraction)
real*4 XMg25s_1 (nr__1) ! Surface abundance of 25Mg (mass fraction)
real*4 XMg26s_1 (nr__1) ! Surface abundance of 26Mg (mass fraction)
real*4 XAl26s_1 (nr__1) ! Surface abundance of 26Al (mass fraction)
real*4 XAl27s_1 (nr__1) ! Surface abundance of 27Al (mass fraction)
real*4 XSi28s_1 (nr__1) ! Surface abundance of 28Si (mass fraction)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'z004.dat' ! Grid of stellar models for Z=0.004
integer*4 nr__2
parameter (nr__2=5811) ! Number of records
character*900 ar__2 ! Full-size record
real*4 Z_2 (nr__2) ! Initial Z abundance
character*1 T_4 (nr__2) ! [rs] s: standard prescriptions or
* r: th+rot prescriptions
real*4 Mass_2 (nr__2) ! (Msun) [0.85/6] Initial mass
* (0.85, 1.0, 1.5, 2.0, 2.5, 3.0, 4.0, 6.0)
integer*4 Mod_2 (nr__2) ! Model number
real*8 L_2 (nr__2) ! (Lsun) Surface luminosity
real*4 Reff_2 (nr__2) ! (Rsun) Photospheric radius
integer*4 Teff_2 (nr__2) ! (K) Effective temperature
real*4 rhoeff_2 (nr__2) ! (g/cm3) Photospheric density
real*4 log_g_1 (nr__2) ! ([cm/s2]) Photospheric gravity
real*4 Mloss_2 (nr__2) ! (Msun/yr) Mass-loss rate
real*8 M_2 (nr__2) ! (Msun) Stellar mass
real*4 t_5 (nr__2) ! (yr) Age
real*4 tchp_2_2 (nr__2) ! (d) Convective turnover timescale calculated at
* half of the pressure scale height
real*4 tchp_2 (nr__2) ! (d) Convective turnover timescale calculated at
* the pressure scale height
real*4 tcrconv_2_2(nr__2) ! (d) Convective turnover timescale calculated at
* half of the radius of the convective envelope
real*4 tcmconv_2_2(nr__2) ! (d) Convective turnover timescale calculated at
* half of the mass of the convective envelope
real*4 tcmax_2 (nr__2) ! (d) Maximum convective turnover timescale
real*4 tg_2 (nr__2) ! (d) Global convective turnover timescale
real*4 rhp_2_2 (nr__2) ! (Rsun) Radius at which tc_hp/2 was estimated
real*4 rhp_2 (nr__2) ! (Rsun) Radius at which tc_hp was estimated
real*4 rrconv_2_2 (nr__2) ! (Rsun) Radius at which tc_rconv/2 was estimated
real*4 rmconv_2_2 (nr__2) ! (Rsun) Radius at which tc_mconv/2 was estimated
real*4 rmax_2 (nr__2) ! (Rsun) Radius at which tc_max was estimated
real*4 Rohp_2_2 (nr__2) ! Rossby number associated to tc_hp/2
real*4 Rohp_2 (nr__2) ! Rossby number associated to tc_hp
real*4 Rorconv_2_2(nr__2) ! Rossby number associated to tc_rconv/2
real*4 Romconv_2_2(nr__2) ! Rossby number associated to tc_mconv/2
real*4 Romax_2 (nr__2) ! Rossby number associated to tc_max
real*4 Rog_2 (nr__2) ! Rossby number associated to tg
real*4 Tc_2 (nr__2) ! (K) Central temperature
real*4 Tmax_2 (nr__2) ! (K) Maximum of temperature
real*4 MrTmax_2 (nr__2) ! (Msun) Mass coordinate of Tmax
real*8 rhoc_2 (nr__2) ! (g/cm3) Central density
real*4 rhomax_2 (nr__2) ! (g/cm3) Density at the location of Tmax
real*4 Pc_2 (nr__2) ! (dPa) Central pressure (cgs)
real*4 Mbenv_2 (nr__2) ! (Msun) Mass at the base of convective envelope
real*8 XHc_2 (nr__2) ! Central abundance of H (mass fraction)
real*4 XHe3c_2 (nr__2) ! Central abundance of 3He (mass fraction)
real*4 XHe4c_2 (nr__2) ! Central abundance of 4He (mass fraction)
real*4 XC12c_2 (nr__2) ! Central abundance of 12C (mass fraction)
real*4 XC13c_2 (nr__2) ! Central abundance of 13C (mass fraction)
real*4 XC14c_2 (nr__2) ! Central abundance of 14C (mass fraction)
real*4 XN14c_2 (nr__2) ! Central abundance of 14N (mass fraction)
real*4 XO16c_2 (nr__2) ! Central abundance of 16O (mass fraction)
real*4 XO17c_2 (nr__2) ! Central abundance of 17O (mass fraction)
real*4 XO18c_2 (nr__2) ! Central abundance of 18O (mass fraction)
real*4 XF19c_2 (nr__2) ! Central abundance of 19F (mass fraction)
real*4 XNe20c_2 (nr__2) ! Central abundance of 20Ne (mass fraction)
real*4 XNe21c_2 (nr__2) ! Central abundance of 21Ne (mass fraction)
real*4 XNe22c_2 (nr__2) ! Central abundance of 22Ne (mass fraction)
real*4 XNa23c_2 (nr__2) ! Central abundance of 23Na (mass fraction)
real*4 XMg24c_2 (nr__2) ! Central abundance of 24Mg (mass fraction)
real*4 XMg25c_2 (nr__2) ! Central abundance of 25Mg (mass fraction)
real*4 XMg26c_2 (nr__2) ! Central abundance of 26Mg (mass fraction)
real*4 XAl26c_2 (nr__2) ! Central abundance of 26Al (mass fraction)
real*4 XAl27c_2 (nr__2) ! Central abundance of 27Al (mass fraction)
real*4 XSi28c_2 (nr__2) ! Central abundance of 28Si (mass fraction)
real*4 XHs_2 (nr__2) ! Surface abundance of H (mass fraction)
real*4 XH2s_2 (nr__2) ! Surface abundance of 2H (mass fraction)
real*4 XHe3s_2 (nr__2) ! Surface abundance of 3He (mass fraction)
real*4 XHe4s_2 (nr__2) ! Surface abundance of 4He (mass fraction)
real*4 XLi6s_2 (nr__2) ! Surface abundance of 6Li (mass fraction)
real*4 XLi7s_2 (nr__2) ! Surface abundance of 7Li (mass fraction)
real*4 XBe7s_2 (nr__2) ! Surface abundance of 7Be (mass fraction)
real*4 XBe9s_2 (nr__2) ! Surface abundance of 9Be (mass fraction)
real*4 XB10s_2 (nr__2) ! Surface abundance of 10B (mass fraction)
real*4 XB11s_2 (nr__2) ! Surface abundance of 11B (mass fraction)
real*4 XC12s_2 (nr__2) ! Surface abundance of 12C (mass fraction)
real*4 XC13s_2 (nr__2) ! Surface abundance of 13C (mass fraction)
real*4 XC14s_2 (nr__2) ! Surface abundance of 14C (mass fraction)
real*4 XN14s_2 (nr__2) ! Surface abundance of 14N (mass fraction)
real*4 XO16s_2 (nr__2) ! Surface abundance of 16O (mass fraction)
real*4 XO17s_2 (nr__2) ! Surface abundance of 17O (mass fraction)
real*4 XO18s_2 (nr__2) ! Surface abundance of 18O (mass fraction)
real*4 XF19s_2 (nr__2) ! Surface abundance of 19F (mass fraction)
real*4 XNe20s_2 (nr__2) ! Surface abundance of 20Ne (mass fraction)
real*4 XNe21s_2 (nr__2) ! Surface abundance of 21Ne (mass fraction)
real*4 XNe22s_2 (nr__2) ! Surface abundance of 22Ne (mass fraction)
real*4 XNa23s_2 (nr__2) ! Surface abundance of 23Na (mass fraction)
real*4 XMg24s_2 (nr__2) ! Surface abundance of 24Mg (mass fraction)
real*4 XMg25s_2 (nr__2) ! Surface abundance of 25Mg (mass fraction)
real*4 XMg26s_2 (nr__2) ! Surface abundance of 26Mg (mass fraction)
real*4 XAl26s_2 (nr__2) ! Surface abundance of 26Al (mass fraction)
real*4 XAl27s_2 (nr__2) ! Surface abundance of 27Al (mass fraction)
real*4 XSi28s_2 (nr__2) ! Surface abundance of 28Si (mass fraction)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'z014.dat' ! Grid of stellar models for Z=0.014
integer*4 nr__3
parameter (nr__3=6339) ! Number of records
character*900 ar__3 ! Full-size record
real*4 Z_3 (nr__3) ! Initial Z abundance
character*1 T_6 (nr__3) ! [rs] s: standard prescriptions or
* r: th+rot prescriptions
real*4 Mass_3 (nr__3) ! (Msun) [0.85/6] Initial mass
* (0.85, 1.0, 1.5, 2.0, 2.5, 3.0, 4.0, 6.0)
integer*4 Mod_3 (nr__3) ! Model number
real*8 L_3 (nr__3) ! (Lsun) Surface luminosity
real*4 Reff_3 (nr__3) ! (Rsun) Photospheric radius
integer*4 Teff_3 (nr__3) ! (K) Effective temperature
real*4 rhoeff_3 (nr__3) ! (g/cm3) Photospheric density
real*4 log_g_1 (nr__3) ! ([cm/s2]) Photospheric gravity
real*4 Mloss_3 (nr__3) ! (Msun/yr) Mass-loss rate
real*8 M_3 (nr__3) ! (Msun) Stellar mass
real*4 t_7 (nr__3) ! (yr) Age
real*4 tchp_2_3 (nr__3) ! (d) Convective turnover timescale calculated at
* half of the pressure scale height
real*4 tchp_3 (nr__3) ! (d) Convective turnover timescale calculated at
* the pressure scale height
real*4 tcrconv_2_3(nr__3) ! (d) Convective turnover timescale calculated at
* half of the radius of the convective envelope
real*4 tcmconv_2_3(nr__3) ! (d) Convective turnover timescale calculated at
* half of the mass of the convective envelope
real*4 tcmax_3 (nr__3) ! (d) Maximum convective turnover timescale
real*4 tg_3 (nr__3) ! (d) Global convective turnover timescale
real*4 rhp_2_3 (nr__3) ! (Rsun) Radius at which tc_hp/2 was estimated
real*4 rhp_3 (nr__3) ! (Rsun) Radius at which tc_hp was estimated
real*4 rrconv_2_3 (nr__3) ! (Rsun) Radius at which tc_rconv/2 was estimated
real*4 rmconv_2_3 (nr__3) ! (Rsun) Radius at which tc_mconv/2 was estimated
real*4 rmax_3 (nr__3) ! (Rsun) Radius at which tc_max was estimated
real*4 Rohp_2_3 (nr__3) ! Rossby number associated to tc_hp/2
real*4 Rohp_3 (nr__3) ! Rossby number associated to tc_hp
real*4 Rorconv_2_3(nr__3) ! Rossby number associated to tc_rconv/2
real*4 Romconv_2_3(nr__3) ! Rossby number associated to tc_mconv/2
real*4 Romax_3 (nr__3) ! Rossby number associated to tc_max
real*4 Rog_3 (nr__3) ! Rossby number associated to tg
real*4 Tc_3 (nr__3) ! (K) Central temperature
real*4 Tmax_3 (nr__3) ! (K) Maximum of temperature
real*4 MrTmax_3 (nr__3) ! (Msun) Mass coordinate of Tmax
real*8 rhoc_3 (nr__3) ! (g/cm3) Central density
real*4 rhomax_3 (nr__3) ! (g/cm3) Density at the location of Tmax
real*4 Pc_3 (nr__3) ! (dPa) Central pressure (cgs)
real*4 Mbenv_3 (nr__3) ! (Msun) Mass at the base of convective envelope
real*8 XHc_3 (nr__3) ! Central abundance of H (mass fraction)
real*4 XHe3c_3 (nr__3) ! Central abundance of 3He (mass fraction)
real*4 XHe4c_3 (nr__3) ! Central abundance of 4He (mass fraction)
real*4 XC12c_3 (nr__3) ! Central abundance of 12C (mass fraction)
real*4 XC13c_3 (nr__3) ! Central abundance of 13C (mass fraction)
real*4 XC14c_3 (nr__3) ! Central abundance of 14C (mass fraction)
real*4 XN14c_3 (nr__3) ! Central abundance of 14N (mass fraction)
real*4 XO16c_3 (nr__3) ! Central abundance of 16O (mass fraction)
real*4 XO17c_3 (nr__3) ! Central abundance of 17O (mass fraction)
real*4 XO18c_3 (nr__3) ! Central abundance of 18O (mass fraction)
real*4 XF19c_3 (nr__3) ! Central abundance of 19F (mass fraction)
real*4 XNe20c_3 (nr__3) ! Central abundance of 20Ne (mass fraction)
real*4 XNe21c_3 (nr__3) ! Central abundance of 21Ne (mass fraction)
real*4 XNe22c_3 (nr__3) ! Central abundance of 22Ne (mass fraction)
real*4 XNa23c_3 (nr__3) ! Central abundance of 23Na (mass fraction)
real*4 XMg24c_3 (nr__3) ! Central abundance of 24Mg (mass fraction)
real*4 XMg25c_3 (nr__3) ! Central abundance of 25Mg (mass fraction)
real*4 XMg26c_3 (nr__3) ! Central abundance of 26Mg (mass fraction)
real*4 XAl26c_3 (nr__3) ! Central abundance of 26Al (mass fraction)
real*4 XAl27c_3 (nr__3) ! Central abundance of 27Al (mass fraction)
real*4 XSi28c_3 (nr__3) ! Central abundance of 28Si (mass fraction)
real*4 XHs_3 (nr__3) ! Surface abundance of H (mass fraction)
real*4 XH2s_3 (nr__3) ! Surface abundance of 2H (mass fraction)
real*4 XHe3s_3 (nr__3) ! Surface abundance of 3He (mass fraction)
real*4 XHe4s_3 (nr__3) ! Surface abundance of 4He (mass fraction)
real*4 XLi6s_3 (nr__3) ! Surface abundance of 6Li (mass fraction)
real*4 XLi7s_3 (nr__3) ! Surface abundance of 7Li (mass fraction)
real*4 XBe7s_3 (nr__3) ! Surface abundance of 7Be (mass fraction)
real*4 XBe9s_3 (nr__3) ! Surface abundance of 9Be (mass fraction)
real*4 XB10s_3 (nr__3) ! Surface abundance of 10B (mass fraction)
real*4 XB11s_3 (nr__3) ! Surface abundance of 11B (mass fraction)
real*4 XC12s_3 (nr__3) ! Surface abundance of 12C (mass fraction)
real*4 XC13s_3 (nr__3) ! Surface abundance of 13C (mass fraction)
real*4 XC14s_3 (nr__3) ! Surface abundance of 14C (mass fraction)
real*4 XN14s_3 (nr__3) ! Surface abundance of 14N (mass fraction)
real*4 XO16s_3 (nr__3) ! Surface abundance of 16O (mass fraction)
real*4 XO17s_3 (nr__3) ! Surface abundance of 17O (mass fraction)
real*4 XO18s_3 (nr__3) ! Surface abundance of 18O (mass fraction)
real*4 XF19s_3 (nr__3) ! Surface abundance of 19F (mass fraction)
real*4 XNe20s_3 (nr__3) ! Surface abundance of 20Ne (mass fraction)
real*4 XNe21s_3 (nr__3) ! Surface abundance of 21Ne (mass fraction)
real*4 XNe22s_3 (nr__3) ! Surface abundance of 22Ne (mass fraction)
real*4 XNa23s_3 (nr__3) ! Surface abundance of 23Na (mass fraction)
real*4 XMg24s_3 (nr__3) ! Surface abundance of 24Mg (mass fraction)
real*4 XMg25s_3 (nr__3) ! Surface abundance of 25Mg (mass fraction)
real*4 XMg26s_3 (nr__3) ! Surface abundance of 26Mg (mass fraction)
real*4 XAl26s_3 (nr__3) ! Surface abundance of 26Al (mass fraction)
real*4 XAl27s_3 (nr__3) ! Surface abundance of 27Al (mass fraction)
real*4 XSi28s_3 (nr__3) ! Surface abundance of 28Si (mass fraction)
C=============================================================================
C Loading file 'z0001.dat' ! Grid of stellar models for Z=0.0001
C Format for file interpretation
1 format(
+ F6.4,1X,A1,1X,F4.2,1X,I3,1X,F8.2,1X,F6.2,1X,I5,1X,E9.3,1X,
+ E9.3,1X,E10.4,1X,F10.8,1X,E16.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,
+ E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,
+ 1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,
+ E10.4,1X,E10.4,1X,E9.3,1X,E9.3,1X,D12.6,1X,E9.3,1X,E9.3,1X,
+ E9.3,1X,D12.6,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,
+ E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,
+ 1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,
+ E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,
+ 1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,
+ E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,
+ 1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,
+ E10.4,1X,E10.4,1X,E10.4,1X,E10.4)
C Effective file loading
open(unit=1,status='old',file=
+'z0001.dat')
write(6,*) '....Loading file: z0001.dat'
do i__=1,7252
read(1,'(A900)')ar__
read(ar__,1)
+ Z(i__),T(i__),Mass(i__),Mod(i__),L(i__),Reff(i__),Teff(i__),
+ rhoeff(i__),log_g(i__),Mloss(i__),M(i__),t_1(i__),tchp_2(i__),
+ tchp(i__),tcrconv_2(i__),tcmconv_2(i__),tcmax(i__),tg(i__),
+ rhp_2(i__),rhp(i__),rrconv_2(i__),rmconv_2(i__),rmax(i__),
+ Rohp_2(i__),Rohp(i__),Rorconv_2(i__),Romconv_2(i__),
+ Romax(i__),Rog(i__),Tc(i__),Tmax(i__),MrTmax(i__),rhoc(i__),
+ rhomax(i__),Pc(i__),Mbenv(i__),XHc(i__),XHe3c(i__),XHe4c(i__),
+ XC12c(i__),XC13c(i__),XC14c(i__),XN14c(i__),XO16c(i__),
+ XO17c(i__),XO18c(i__),XF19c(i__),XNe20c(i__),XNe21c(i__),
+ XNe22c(i__),XNa23c(i__),XMg24c(i__),XMg25c(i__),XMg26c(i__),
+ XAl26c(i__),XAl27c(i__),XSi28c(i__),XHs(i__),XH2s(i__),
+ XHe3s(i__),XHe4s(i__),XLi6s(i__),XLi7s(i__),XBe7s(i__),
+ XBe9s(i__),XB10s(i__),XB11s(i__),XC12s(i__),XC13s(i__),
+ XC14s(i__),XN14s(i__),XO16s(i__),XO17s(i__),XO18s(i__),
+ XF19s(i__),XNe20s(i__),XNe21s(i__),XNe22s(i__),XNa23s(i__),
+ XMg24s(i__),XMg25s(i__),XMg26s(i__),XAl26s(i__),XAl27s(i__),
+ XSi28s(i__)
c ..............Just test output...........
write(6,1)
+ Z(i__),T(i__),Mass(i__),Mod(i__),L(i__),Reff(i__),Teff(i__),
+ rhoeff(i__),log_g(i__),Mloss(i__),M(i__),t_1(i__),tchp_2(i__),
+ tchp(i__),tcrconv_2(i__),tcmconv_2(i__),tcmax(i__),tg(i__),
+ rhp_2(i__),rhp(i__),rrconv_2(i__),rmconv_2(i__),rmax(i__),
+ Rohp_2(i__),Rohp(i__),Rorconv_2(i__),Romconv_2(i__),
+ Romax(i__),Rog(i__),Tc(i__),Tmax(i__),MrTmax(i__),rhoc(i__),
+ rhomax(i__),Pc(i__),Mbenv(i__),XHc(i__),XHe3c(i__),XHe4c(i__),
+ XC12c(i__),XC13c(i__),XC14c(i__),XN14c(i__),XO16c(i__),
+ XO17c(i__),XO18c(i__),XF19c(i__),XNe20c(i__),XNe21c(i__),
+ XNe22c(i__),XNa23c(i__),XMg24c(i__),XMg25c(i__),XMg26c(i__),
+ XAl26c(i__),XAl27c(i__),XSi28c(i__),XHs(i__),XH2s(i__),
+ XHe3s(i__),XHe4s(i__),XLi6s(i__),XLi7s(i__),XBe7s(i__),
+ XBe9s(i__),XB10s(i__),XB11s(i__),XC12s(i__),XC13s(i__),
+ XC14s(i__),XN14s(i__),XO16s(i__),XO17s(i__),XO18s(i__),
+ XF19s(i__),XNe20s(i__),XNe21s(i__),XNe22s(i__),XNa23s(i__),
+ XMg24s(i__),XMg25s(i__),XMg26s(i__),XAl26s(i__),XAl27s(i__),
+ XSi28s(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'z002.dat' ! Grid of stellar models for Z=0.002
C Format for file interpretation
2 format(
+ F6.4,1X,A1,1X,F4.2,1X,I3,1X,F8.2,1X,F6.2,1X,I5,1X,E9.3,1X,
+ E9.3,1X,E10.4,1X,F10.8,1X,E16.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,
+ E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,
+ 1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,
+ E10.4,1X,E10.4,1X,E9.3,1X,E9.3,1X,D12.6,1X,E9.3,1X,E9.3,1X,
+ E9.3,1X,D12.6,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,
+ E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,
+ 1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,
+ E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,
+ 1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,
+ E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,
+ 1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,
+ E10.4,1X,E10.4,1X,E10.4,1X,E10.4)
C Effective file loading
open(unit=1,status='old',file=
+'z002.dat')
write(6,*) '....Loading file: z002.dat'
do i__=1,7227
read(1,'(A900)')ar__1
read(ar__1,2)
+ Z_1(i__),T_2(i__),Mass_1(i__),Mod_1(i__),L_1(i__),Reff_1(i__),
+ Teff_1(i__),rhoeff_1(i__),log_g_1(i__),Mloss_1(i__),M_1(i__),
+ t_3(i__),tchp_2_1(i__),tchp_1(i__),tcrconv_2_1(i__),
+ tcmconv_2_1(i__),tcmax_1(i__),tg_1(i__),rhp_2_1(i__),
+ rhp_1(i__),rrconv_2_1(i__),rmconv_2_1(i__),rmax_1(i__),
+ Rohp_2_1(i__),Rohp_1(i__),Rorconv_2_1(i__),Romconv_2_1(i__),
+ Romax_1(i__),Rog_1(i__),Tc_1(i__),Tmax_1(i__),MrTmax_1(i__),
+ rhoc_1(i__),rhomax_1(i__),Pc_1(i__),Mbenv_1(i__),XHc_1(i__),
+ XHe3c_1(i__),XHe4c_1(i__),XC12c_1(i__),XC13c_1(i__),
+ XC14c_1(i__),XN14c_1(i__),XO16c_1(i__),XO17c_1(i__),
+ XO18c_1(i__),XF19c_1(i__),XNe20c_1(i__),XNe21c_1(i__),
+ XNe22c_1(i__),XNa23c_1(i__),XMg24c_1(i__),XMg25c_1(i__),
+ XMg26c_1(i__),XAl26c_1(i__),XAl27c_1(i__),XSi28c_1(i__),
+ XHs_1(i__),XH2s_1(i__),XHe3s_1(i__),XHe4s_1(i__),XLi6s_1(i__),
+ XLi7s_1(i__),XBe7s_1(i__),XBe9s_1(i__),XB10s_1(i__),
+ XB11s_1(i__),XC12s_1(i__),XC13s_1(i__),XC14s_1(i__),
+ XN14s_1(i__),XO16s_1(i__),XO17s_1(i__),XO18s_1(i__),
+ XF19s_1(i__),XNe20s_1(i__),XNe21s_1(i__),XNe22s_1(i__),
+ XNa23s_1(i__),XMg24s_1(i__),XMg25s_1(i__),XMg26s_1(i__),
+ XAl26s_1(i__),XAl27s_1(i__),XSi28s_1(i__)
c ..............Just test output...........
write(6,2)
+ Z_1(i__),T_2(i__),Mass_1(i__),Mod_1(i__),L_1(i__),Reff_1(i__),
+ Teff_1(i__),rhoeff_1(i__),log_g_1(i__),Mloss_1(i__),M_1(i__),
+ t_3(i__),tchp_2_1(i__),tchp_1(i__),tcrconv_2_1(i__),
+ tcmconv_2_1(i__),tcmax_1(i__),tg_1(i__),rhp_2_1(i__),
+ rhp_1(i__),rrconv_2_1(i__),rmconv_2_1(i__),rmax_1(i__),
+ Rohp_2_1(i__),Rohp_1(i__),Rorconv_2_1(i__),Romconv_2_1(i__),
+ Romax_1(i__),Rog_1(i__),Tc_1(i__),Tmax_1(i__),MrTmax_1(i__),
+ rhoc_1(i__),rhomax_1(i__),Pc_1(i__),Mbenv_1(i__),XHc_1(i__),
+ XHe3c_1(i__),XHe4c_1(i__),XC12c_1(i__),XC13c_1(i__),
+ XC14c_1(i__),XN14c_1(i__),XO16c_1(i__),XO17c_1(i__),
+ XO18c_1(i__),XF19c_1(i__),XNe20c_1(i__),XNe21c_1(i__),
+ XNe22c_1(i__),XNa23c_1(i__),XMg24c_1(i__),XMg25c_1(i__),
+ XMg26c_1(i__),XAl26c_1(i__),XAl27c_1(i__),XSi28c_1(i__),
+ XHs_1(i__),XH2s_1(i__),XHe3s_1(i__),XHe4s_1(i__),XLi6s_1(i__),
+ XLi7s_1(i__),XBe7s_1(i__),XBe9s_1(i__),XB10s_1(i__),
+ XB11s_1(i__),XC12s_1(i__),XC13s_1(i__),XC14s_1(i__),
+ XN14s_1(i__),XO16s_1(i__),XO17s_1(i__),XO18s_1(i__),
+ XF19s_1(i__),XNe20s_1(i__),XNe21s_1(i__),XNe22s_1(i__),
+ XNa23s_1(i__),XMg24s_1(i__),XMg25s_1(i__),XMg26s_1(i__),
+ XAl26s_1(i__),XAl27s_1(i__),XSi28s_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'z004.dat' ! Grid of stellar models for Z=0.004
C Format for file interpretation
3 format(
+ F6.4,1X,A1,1X,F4.2,1X,I3,1X,F8.2,1X,F6.2,1X,I5,1X,E9.3,1X,
+ E9.3,1X,E10.4,1X,F10.8,1X,E16.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,
+ E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,
+ 1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,
+ E10.4,1X,E10.4,1X,E9.3,1X,E9.3,1X,D12.6,1X,E9.3,1X,E9.3,1X,
+ E9.3,1X,D12.6,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,
+ E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,
+ 1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,
+ E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,
+ 1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,
+ E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,
+ 1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,
+ E10.4,1X,E10.4,1X,E10.4,1X,E10.4)
C Effective file loading
open(unit=1,status='old',file=
+'z004.dat')
write(6,*) '....Loading file: z004.dat'
do i__=1,5811
read(1,'(A900)')ar__2
read(ar__2,3)
+ Z_2(i__),T_4(i__),Mass_2(i__),Mod_2(i__),L_2(i__),Reff_2(i__),
+ Teff_2(i__),rhoeff_2(i__),log_g_1(i__),Mloss_2(i__),M_2(i__),
+ t_5(i__),tchp_2_2(i__),tchp_2(i__),tcrconv_2_2(i__),
+ tcmconv_2_2(i__),tcmax_2(i__),tg_2(i__),rhp_2_2(i__),
+ rhp_2(i__),rrconv_2_2(i__),rmconv_2_2(i__),rmax_2(i__),
+ Rohp_2_2(i__),Rohp_2(i__),Rorconv_2_2(i__),Romconv_2_2(i__),
+ Romax_2(i__),Rog_2(i__),Tc_2(i__),Tmax_2(i__),MrTmax_2(i__),
+ rhoc_2(i__),rhomax_2(i__),Pc_2(i__),Mbenv_2(i__),XHc_2(i__),
+ XHe3c_2(i__),XHe4c_2(i__),XC12c_2(i__),XC13c_2(i__),
+ XC14c_2(i__),XN14c_2(i__),XO16c_2(i__),XO17c_2(i__),
+ XO18c_2(i__),XF19c_2(i__),XNe20c_2(i__),XNe21c_2(i__),
+ XNe22c_2(i__),XNa23c_2(i__),XMg24c_2(i__),XMg25c_2(i__),
+ XMg26c_2(i__),XAl26c_2(i__),XAl27c_2(i__),XSi28c_2(i__),
+ XHs_2(i__),XH2s_2(i__),XHe3s_2(i__),XHe4s_2(i__),XLi6s_2(i__),
+ XLi7s_2(i__),XBe7s_2(i__),XBe9s_2(i__),XB10s_2(i__),
+ XB11s_2(i__),XC12s_2(i__),XC13s_2(i__),XC14s_2(i__),
+ XN14s_2(i__),XO16s_2(i__),XO17s_2(i__),XO18s_2(i__),
+ XF19s_2(i__),XNe20s_2(i__),XNe21s_2(i__),XNe22s_2(i__),
+ XNa23s_2(i__),XMg24s_2(i__),XMg25s_2(i__),XMg26s_2(i__),
+ XAl26s_2(i__),XAl27s_2(i__),XSi28s_2(i__)
c ..............Just test output...........
write(6,3)
+ Z_2(i__),T_4(i__),Mass_2(i__),Mod_2(i__),L_2(i__),Reff_2(i__),
+ Teff_2(i__),rhoeff_2(i__),log_g_1(i__),Mloss_2(i__),M_2(i__),
+ t_5(i__),tchp_2_2(i__),tchp_2(i__),tcrconv_2_2(i__),
+ tcmconv_2_2(i__),tcmax_2(i__),tg_2(i__),rhp_2_2(i__),
+ rhp_2(i__),rrconv_2_2(i__),rmconv_2_2(i__),rmax_2(i__),
+ Rohp_2_2(i__),Rohp_2(i__),Rorconv_2_2(i__),Romconv_2_2(i__),
+ Romax_2(i__),Rog_2(i__),Tc_2(i__),Tmax_2(i__),MrTmax_2(i__),
+ rhoc_2(i__),rhomax_2(i__),Pc_2(i__),Mbenv_2(i__),XHc_2(i__),
+ XHe3c_2(i__),XHe4c_2(i__),XC12c_2(i__),XC13c_2(i__),
+ XC14c_2(i__),XN14c_2(i__),XO16c_2(i__),XO17c_2(i__),
+ XO18c_2(i__),XF19c_2(i__),XNe20c_2(i__),XNe21c_2(i__),
+ XNe22c_2(i__),XNa23c_2(i__),XMg24c_2(i__),XMg25c_2(i__),
+ XMg26c_2(i__),XAl26c_2(i__),XAl27c_2(i__),XSi28c_2(i__),
+ XHs_2(i__),XH2s_2(i__),XHe3s_2(i__),XHe4s_2(i__),XLi6s_2(i__),
+ XLi7s_2(i__),XBe7s_2(i__),XBe9s_2(i__),XB10s_2(i__),
+ XB11s_2(i__),XC12s_2(i__),XC13s_2(i__),XC14s_2(i__),
+ XN14s_2(i__),XO16s_2(i__),XO17s_2(i__),XO18s_2(i__),
+ XF19s_2(i__),XNe20s_2(i__),XNe21s_2(i__),XNe22s_2(i__),
+ XNa23s_2(i__),XMg24s_2(i__),XMg25s_2(i__),XMg26s_2(i__),
+ XAl26s_2(i__),XAl27s_2(i__),XSi28s_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'z014.dat' ! Grid of stellar models for Z=0.014
C Format for file interpretation
4 format(
+ F6.4,1X,A1,1X,F4.2,1X,I3,1X,F8.2,1X,F6.2,1X,I5,1X,E9.3,1X,
+ E9.3,1X,E10.4,1X,F10.8,1X,E16.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,
+ E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,
+ 1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,
+ E10.4,1X,E10.4,1X,E9.3,1X,E9.3,1X,D12.6,1X,E9.3,1X,E9.3,1X,
+ E9.3,1X,D12.6,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,
+ E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,
+ 1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,
+ E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,
+ 1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,
+ E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,
+ 1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,E10.4,1X,
+ E10.4,1X,E10.4,1X,E10.4,1X,E10.4)
C Effective file loading
open(unit=1,status='old',file=
+'z014.dat')
write(6,*) '....Loading file: z014.dat'
do i__=1,6339
read(1,'(A900)')ar__3
read(ar__3,4)
+ Z_3(i__),T_6(i__),Mass_3(i__),Mod_3(i__),L_3(i__),Reff_3(i__),
+ Teff_3(i__),rhoeff_3(i__),log_g_1(i__),Mloss_3(i__),M_3(i__),
+ t_7(i__),tchp_2_3(i__),tchp_3(i__),tcrconv_2_3(i__),
+ tcmconv_2_3(i__),tcmax_3(i__),tg_3(i__),rhp_2_3(i__),
+ rhp_3(i__),rrconv_2_3(i__),rmconv_2_3(i__),rmax_3(i__),
+ Rohp_2_3(i__),Rohp_3(i__),Rorconv_2_3(i__),Romconv_2_3(i__),
+ Romax_3(i__),Rog_3(i__),Tc_3(i__),Tmax_3(i__),MrTmax_3(i__),
+ rhoc_3(i__),rhomax_3(i__),Pc_3(i__),Mbenv_3(i__),XHc_3(i__),
+ XHe3c_3(i__),XHe4c_3(i__),XC12c_3(i__),XC13c_3(i__),
+ XC14c_3(i__),XN14c_3(i__),XO16c_3(i__),XO17c_3(i__),
+ XO18c_3(i__),XF19c_3(i__),XNe20c_3(i__),XNe21c_3(i__),
+ XNe22c_3(i__),XNa23c_3(i__),XMg24c_3(i__),XMg25c_3(i__),
+ XMg26c_3(i__),XAl26c_3(i__),XAl27c_3(i__),XSi28c_3(i__),
+ XHs_3(i__),XH2s_3(i__),XHe3s_3(i__),XHe4s_3(i__),XLi6s_3(i__),
+ XLi7s_3(i__),XBe7s_3(i__),XBe9s_3(i__),XB10s_3(i__),
+ XB11s_3(i__),XC12s_3(i__),XC13s_3(i__),XC14s_3(i__),
+ XN14s_3(i__),XO16s_3(i__),XO17s_3(i__),XO18s_3(i__),
+ XF19s_3(i__),XNe20s_3(i__),XNe21s_3(i__),XNe22s_3(i__),
+ XNa23s_3(i__),XMg24s_3(i__),XMg25s_3(i__),XMg26s_3(i__),
+ XAl26s_3(i__),XAl27s_3(i__),XSi28s_3(i__)
c ..............Just test output...........
write(6,4)
+ Z_3(i__),T_6(i__),Mass_3(i__),Mod_3(i__),L_3(i__),Reff_3(i__),
+ Teff_3(i__),rhoeff_3(i__),log_g_1(i__),Mloss_3(i__),M_3(i__),
+ t_7(i__),tchp_2_3(i__),tchp_3(i__),tcrconv_2_3(i__),
+ tcmconv_2_3(i__),tcmax_3(i__),tg_3(i__),rhp_2_3(i__),
+ rhp_3(i__),rrconv_2_3(i__),rmconv_2_3(i__),rmax_3(i__),
+ Rohp_2_3(i__),Rohp_3(i__),Rorconv_2_3(i__),Romconv_2_3(i__),
+ Romax_3(i__),Rog_3(i__),Tc_3(i__),Tmax_3(i__),MrTmax_3(i__),
+ rhoc_3(i__),rhomax_3(i__),Pc_3(i__),Mbenv_3(i__),XHc_3(i__),
+ XHe3c_3(i__),XHe4c_3(i__),XC12c_3(i__),XC13c_3(i__),
+ XC14c_3(i__),XN14c_3(i__),XO16c_3(i__),XO17c_3(i__),
+ XO18c_3(i__),XF19c_3(i__),XNe20c_3(i__),XNe21c_3(i__),
+ XNe22c_3(i__),XNa23c_3(i__),XMg24c_3(i__),XMg25c_3(i__),
+ XMg26c_3(i__),XAl26c_3(i__),XAl27c_3(i__),XSi28c_3(i__),
+ XHs_3(i__),XH2s_3(i__),XHe3s_3(i__),XHe4s_3(i__),XLi6s_3(i__),
+ XLi7s_3(i__),XBe7s_3(i__),XBe9s_3(i__),XB10s_3(i__),
+ XB11s_3(i__),XC12s_3(i__),XC13s_3(i__),XC14s_3(i__),
+ XN14s_3(i__),XO16s_3(i__),XO17s_3(i__),XO18s_3(i__),
+ XF19s_3(i__),XNe20s_3(i__),XNe21s_3(i__),XNe22s_3(i__),
+ XNa23s_3(i__),XMg24s_3(i__),XMg25s_3(i__),XMg26s_3(i__),
+ XAl26s_3(i__),XAl27s_3(i__),XSi28s_3(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end