Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/530/A115 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/530/A115 Rotating massive MS stars evolutionary models (Brott+, 2011)
*================================================================================
*Rotating massive main-sequence stars.
*I: Grids of evolutionary models and isochrones.
* Brott I., de Mink S.E., Cantiello M., Langer N., de Koter A., Evans C.J.,
* Hunter I., Trundle C., Vink J.S.
* <Astron. Astrophys. 530, A115 (2011)>
* =2011A&A...530A.115B
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'models.dat' ! List of models files
integer*4 nr__
parameter (nr__=663) ! Number of records
character*33 ar__ ! Full-size record
integer*4 Minit (nr__) ! (Msun) [5/60] Inital Mass
integer*4 Vzams (nr__) ! (km/s) [0/600] Inital rotational velocity
character*3 Comp (nr__) ! Chemical mixture (lmc, smc or mw)
character*21 FileName (nr__) ! Name of the file in subdirectory evol
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'isochron.dat' ! List of isochrones files
integer*4 nr__1
parameter (nr__1=528) ! Number of records
character*35 ar__1 ! Full-size record
real*4 Age (nr__1) ! (yr) [0/35e6] Age
character*3 Comp_1 (nr__1) ! Chemical mixture (lmc, smc or mw)
character*21 FileName_1 (nr__1) ! Name of the file in subdirectory iso
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'evol/*' ! Evolutionary models with Galactic, LMC and SMC
initial abundances
integer*4 nr__2
parameter (nr__2=663) ! Number of records
character*1001 ar__2 ! Full-size record
real*4 t (nr__2) ! (yr) Time/age on evolutionary track
real*8 Mass (nr__2) ! (Msun) Mass
real*8 Teff (nr__2) ! (K) Effective temperature
real*8 logL (nr__2) ! ([Lsun]) Luminosity
real*8 R (nr__2) ! (Rsun) Radius
real*8 log_Mdot (nr__2) ! ([Msun/yr]) Mass loss rate
real*8 logg (nr__2) ! ([cm/s2]) Surface gravity
real*8 Vsurf (nr__2) ! (km/s) Rotational velocity at the surface
real*8 Prot (nr__2) ! (d) Rotation period
real*8 Vcrit (nr__2) ! (km/s) ? Beakup velocity (asuming an Eddington
* factor for pure election scattering) (1)
real*8 Ge (nr__2) ! Eddington {Gamma}_e_ factor calculated
* for pure electron scattering
real*8 eps_H (nr__2) ! ([---]) [12] Hydrogen surface abundance (G1)
real*8 eps_He (nr__2) ! ([---]) Helium surface abundance (G1)
real*8 eps_Li (nr__2) ! ([---]) Lithium surface abundance (G1)
real*8 eps_Be (nr__2) ! ([---]) Berilium surface abundance (G1)
real*8 eps_B (nr__2) ! ([---]) Boron surface abundance (G1)
real*8 eps_C (nr__2) ! ([---]) Carbon surface abundance (G1)
real*8 eps_N (nr__2) ! ([---]) Nitrogen surface abundance (G1)
real*8 eps_O (nr__2) ! ([---]) Oxygen surface abundance (G1)
real*8 eps_F (nr__2) ! ([---]) Fluron surface abundance (G1)
real*8 eps_Ne (nr__2) ! ([---]) Neon surface abundance (G1)
real*8 eps_Na (nr__2) ! ([---]) Sodium surface abundance (G1)
real*8 eps_Mg (nr__2) ! ([---]) Magnesium surface abundance (G1)
real*8 eps_Al (nr__2) ! ([---]) Aluminium surface abundance (G1)
real*8 eps_Si (nr__2) ! ([---]) Silicium surface abundance (G1)
real*8 eps_Fe (nr__2) ! ([---]) Iron surface abundance (G1)
real*4 sH1 (nr__2) ! Surface mass fraction of H1 isotope
real*4 sHe3 (nr__2) ! Surface mass fraction of He3 isotope
real*4 sHe4 (nr__2) ! Surface mass fraction of He4 isotope
real*4 sLi7 (nr__2) ! Surface mass fraction of Li7 isotope
real*4 sBe9 (nr__2) ! Surface mass fraction of Be9 isotope
real*4 sB10 (nr__2) ! Surface mass fraction of B10 isotope
real*4 sB11 (nr__2) ! Surface mass fraction of B11 isotope
real*4 sC12 (nr__2) ! Surface mass fraction of C12 isotope
real*4 sC13 (nr__2) ! Surface mass fraction of C13 isotope
real*4 sN14 (nr__2) ! Surface mass fraction of N14 isotope
real*4 sN15 (nr__2) ! Surface mass fraction of N15 isotope
real*4 sO16 (nr__2) ! Surface mass fraction of O16 isotope
real*4 sO17 (nr__2) ! Surface mass fraction of O17 isotope
real*4 sO18 (nr__2) ! Surface mass fraction of O18 isotope
real*4 sF19 (nr__2) ! Surface mass fraction of F19 isotope
real*4 sNe20 (nr__2) ! Surface mass fraction of Ne20 isotope
real*4 sNe21 (nr__2) ! Surface mass fraction of Ne21 isotope
real*4 sNe22 (nr__2) ! Surface mass fraction of Ne22 isotope
real*4 sNa23 (nr__2) ! Surface mass fraction of Na23 isotope
real*4 sMg24 (nr__2) ! Surface mass fraction of Mg24 isotope
real*4 sMg25 (nr__2) ! Surface mass fraction of Mg26 isotope
real*4 sMg26 (nr__2) ! Surface mass fraction of Mg26 isotope
real*4 sAl26 (nr__2) ! Surface mass fraction of Al26 isotope
real*4 sAl27 (nr__2) ! Surface mass fraction of Al27 isotope
real*4 sSi28 (nr__2) ! Surface mass fraction of Si28 isotope
real*4 sSi29 (nr__2) ! Surface mass fraction of Si29 isotope
real*4 sSi30 (nr__2) ! Surface mass fraction of Si30 isotope
real*4 sFe56 (nr__2) ! Surface mass fraction of Fe56 isotope
real*4 cH1 (nr__2) ! Core mass fraction of H1 isotope
real*4 cHe3 (nr__2) ! Core mass fraction of He3 isotope
real*4 cHe4 (nr__2) ! Core mass fraction of He4 isotope
real*4 cLi7 (nr__2) ! Core mass fraction of Li7 isotope
real*4 cBe9 (nr__2) ! Core mass fraction of Be9 isotope
real*4 cB10 (nr__2) ! Core mass fraction of B10 isotope
real*4 cB11 (nr__2) ! Core mass fraction of B11 isotope
real*4 cC12 (nr__2) ! Core mass fraction of C12 isotope
real*4 cC13 (nr__2) ! Core mass fraction of C13 isotope
real*4 cN14 (nr__2) ! Core mass fraction of N14 isotope
real*4 cN15 (nr__2) ! Core mass fraction of N15 isotope
real*4 cO16 (nr__2) ! Core mass fraction of O16 isotope
real*4 cO17 (nr__2) ! Core mass fraction of O17 isotope
real*4 cO18 (nr__2) ! Core mass fraction of O18 isotope
real*4 cF19 (nr__2) ! Core mass fraction of F19 isotope
real*4 cNe20 (nr__2) ! Core mass fraction of Ne20 isotope
real*4 cNe21 (nr__2) ! Core mass fraction of Ne21 isotope
real*4 cNe22 (nr__2) ! Core mass fraction of Ne22 isotope
real*4 cNa23 (nr__2) ! Core mass fraction of Na23 isotope
real*4 cMg24 (nr__2) ! Core mass fraction of Mg24 isotope
real*4 cMg25 (nr__2) ! Core mass fraction of Mg25 isotope
real*4 cMg26 (nr__2) ! Core mass fraction of Mg26 isotope
real*4 cAl26 (nr__2) ! Core mass fraction of Al26 isotope
real*4 cAl27 (nr__2) ! Core mass fraction of Al27 isotope
real*4 cSi28 (nr__2) ! Core mass fraction of Si28 isotope
real*4 cSi29 (nr__2) ! Core mass fraction of Si29 isotope
real*4 cSi30 (nr__2) ! Core mass fraction of Si30 isotope
real*4 cFe56 (nr__2) ! Core mass fraction of Fe56 isotope
*Note (1): Blank for nan (no value)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'iso/*' ! Isochrones generated from evolutionnary models
integer*4 nr__3
parameter (nr__3=528) ! Number of records
character*235 ar__3 ! Full-size record
real*8 t_1 (nr__3) ! (yr) time/age of isochrone
real*8 Mass_1 (nr__3) ! (Msun) Mass
real*8 Teff_1 (nr__3) ! (K) Effective temperature
real*8 logL_1 (nr__3) ! ([Lsun]) Luminosity
real*8 R_1 (nr__3) ! (Rsun) Radius
real*8 log_Mdot_1 (nr__3) ! ([Msun/yr]) Mass loss rate
real*8 logg_1 (nr__3) ! ([s/cm2]) Surface gravity
real*8 Vsurf_1 (nr__3) ! (km/s) Rotational velocity at the surface
real*8 Prot_1 (nr__3) ! (d) Rotation period
real*8 Vcrit_1 (nr__3) ! (km/s) Beakup velocity (asuming an Eddington
* factor for pure election scattering)
real*8 Ge_1 (nr__3) ! Eddington {Gamma}_e_ factor calculated for
* pure electron scattering
real*8 eps_H_1 (nr__3) ! ([---]) [12] Hydrogen surface abundance (G1)
real*8 eps_He_1 (nr__3) ! ([---]) Helium surface abundance (G1)
real*8 eps_B_1 (nr__3) ! ([---]) Boron surface abundance (G1)
real*8 eps_C_1 (nr__3) ! ([---]) Carbon surface abundance (G1)
real*8 eps_N_1 (nr__3) ! ([---]) Nitrogen surface abundance (G1)
real*8 eps_O_1 (nr__3) ! ([---]) Oxygen surface abundance (G1)
real*8 eps_Ne_1 (nr__3) ! ([---]) Neon surface abundance (G1)
real*8 eps_Na_1 (nr__3) ! ([---]) Sodium surface abundance (G1)
real*8 eps_Mg_1 (nr__3) ! ([---]) Magnesium surface abundance (G1)
real*8 eps_Si_1 (nr__3) ! ([---]) Silicium surface abundance (G1)
real*8 H1mf (nr__3) ! Surface mass fraction of H1 isotope
real*8 He4mf (nr__3) ! Surface mass fraction of He4 isotope
real*8 Minit_1 (nr__3) ! (Msun) Inital Mass
real*8 vinit (nr__3) ! (km/s) Inital rotational velocity
C=============================================================================
C Loading file 'models.dat' ! List of models files
C Format for file interpretation
1 format(I2,1X,I3,1X,A3,2X,A21)
C Effective file loading
open(unit=1,status='old',file=
+'models.dat')
write(6,*) '....Loading file: models.dat'
do i__=1,663
read(1,'(A33)')ar__
read(ar__,1)Minit(i__),Vzams(i__),Comp(i__),FileName(i__)
c ..............Just test output...........
write(6,1)Minit(i__),Vzams(i__),Comp(i__),FileName(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'isochron.dat' ! List of isochrones files
C Format for file interpretation
2 format(E8.3,1X,A3,2X,A21)
C Effective file loading
open(unit=1,status='old',file=
+'isochron.dat')
write(6,*) '....Loading file: isochron.dat'
do i__=1,528
read(1,'(A35)')ar__1
read(ar__1,2)Age(i__),Comp_1(i__),FileName_1(i__)
c ..............Just test output...........
write(6,2)Age(i__),Comp_1(i__),FileName_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'evol/*' ! Evolutionary models with Galactic, LMC and SMC
* initial abundances
C Format for file interpretation
3 format(
+ E12.6,1X,F9.5,1X,F12.5,1X,F8.5,1X,F10.5,1X,F9.5,1X,F9.5,1X,
+ F10.5,1X,F15.5,1X,F10.5,1X,F9.5,1X,F9.5,1X,F9.5,1X,F9.5,1X,
+ F9.5,1X,F9.5,1X,F9.5,1X,F9.5,1X,F9.5,1X,F9.5,1X,F9.5,1X,F9.5,
+ 1X,F9.5,1X,F9.5,1X,F9.5,1X,F9.5,1X,E12.6,1X,E12.6,1X,E12.6,1X,
+ E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,
+ 1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,
+ E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,
+ 1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,
+ E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,
+ 1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,
+ E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,
+ 1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6)
C Effective file loading
open(unit=1,status='old',file=
+'evol/*')
write(6,*) '....Loading file: evol/*'
do i__=1,663
read(1,'(A1001)')ar__2
read(ar__2,3)
+ t(i__),Mass(i__),Teff(i__),logL(i__),R(i__),log_Mdot(i__),
+ logg(i__),Vsurf(i__),Prot(i__),Vcrit(i__),Ge(i__),eps_H(i__),
+ eps_He(i__),eps_Li(i__),eps_Be(i__),eps_B(i__),eps_C(i__),
+ eps_N(i__),eps_O(i__),eps_F(i__),eps_Ne(i__),eps_Na(i__),
+ eps_Mg(i__),eps_Al(i__),eps_Si(i__),eps_Fe(i__),sH1(i__),
+ sHe3(i__),sHe4(i__),sLi7(i__),sBe9(i__),sB10(i__),sB11(i__),
+ sC12(i__),sC13(i__),sN14(i__),sN15(i__),sO16(i__),sO17(i__),
+ sO18(i__),sF19(i__),sNe20(i__),sNe21(i__),sNe22(i__),
+ sNa23(i__),sMg24(i__),sMg25(i__),sMg26(i__),sAl26(i__),
+ sAl27(i__),sSi28(i__),sSi29(i__),sSi30(i__),sFe56(i__),
+ cH1(i__),cHe3(i__),cHe4(i__),cLi7(i__),cBe9(i__),cB10(i__),
+ cB11(i__),cC12(i__),cC13(i__),cN14(i__),cN15(i__),cO16(i__),
+ cO17(i__),cO18(i__),cF19(i__),cNe20(i__),cNe21(i__),
+ cNe22(i__),cNa23(i__),cMg24(i__),cMg25(i__),cMg26(i__),
+ cAl26(i__),cAl27(i__),cSi28(i__),cSi29(i__),cSi30(i__),
+ cFe56(i__)
if(ar__2(104:113) .EQ. '') Vcrit(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ t(i__),Mass(i__),Teff(i__),logL(i__),R(i__),log_Mdot(i__),
+ logg(i__),Vsurf(i__),Prot(i__),Vcrit(i__),Ge(i__),eps_H(i__),
+ eps_He(i__),eps_Li(i__),eps_Be(i__),eps_B(i__),eps_C(i__),
+ eps_N(i__),eps_O(i__),eps_F(i__),eps_Ne(i__),eps_Na(i__),
+ eps_Mg(i__),eps_Al(i__),eps_Si(i__),eps_Fe(i__),sH1(i__),
+ sHe3(i__),sHe4(i__),sLi7(i__),sBe9(i__),sB10(i__),sB11(i__),
+ sC12(i__),sC13(i__),sN14(i__),sN15(i__),sO16(i__),sO17(i__),
+ sO18(i__),sF19(i__),sNe20(i__),sNe21(i__),sNe22(i__),
+ sNa23(i__),sMg24(i__),sMg25(i__),sMg26(i__),sAl26(i__),
+ sAl27(i__),sSi28(i__),sSi29(i__),sSi30(i__),sFe56(i__),
+ cH1(i__),cHe3(i__),cHe4(i__),cLi7(i__),cBe9(i__),cB10(i__),
+ cB11(i__),cC12(i__),cC13(i__),cN14(i__),cN15(i__),cO16(i__),
+ cO17(i__),cO18(i__),cF19(i__),cNe20(i__),cNe21(i__),
+ cNe22(i__),cNa23(i__),cMg24(i__),cMg25(i__),cMg26(i__),
+ cAl26(i__),cAl27(i__),cSi28(i__),cSi29(i__),cSi30(i__),
+ cFe56(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'iso/*' ! Isochrones generated from evolutionnary models
C Format for file interpretation
4 format(
+ F13.4,1X,F7.4,1X,F11.4,1X,F7.4,1X,F10.4,1X,F9.4,1X,F9.4,1X,
+ F9.4,1X,F14.4,1X,F9.4,1X,F7.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F7.4,
+ 1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F9.6,1X,
+ F9.6,1X,F7.4,1X,F8.4)
C Effective file loading
open(unit=1,status='old',file=
+'iso/*')
write(6,*) '....Loading file: iso/*'
do i__=1,528
read(1,'(A235)')ar__3
read(ar__3,4)
+ t_1(i__),Mass_1(i__),Teff_1(i__),logL_1(i__),R_1(i__),
+ log_Mdot_1(i__),logg_1(i__),Vsurf_1(i__),Prot_1(i__),
+ Vcrit_1(i__),Ge_1(i__),eps_H_1(i__),eps_He_1(i__),
+ eps_B_1(i__),eps_C_1(i__),eps_N_1(i__),eps_O_1(i__),
+ eps_Ne_1(i__),eps_Na_1(i__),eps_Mg_1(i__),eps_Si_1(i__),
+ H1mf(i__),He4mf(i__),Minit_1(i__),vinit(i__)
c ..............Just test output...........
write(6,4)
+ t_1(i__),Mass_1(i__),Teff_1(i__),logL_1(i__),R_1(i__),
+ log_Mdot_1(i__),logg_1(i__),Vsurf_1(i__),Prot_1(i__),
+ Vcrit_1(i__),Ge_1(i__),eps_H_1(i__),eps_He_1(i__),
+ eps_B_1(i__),eps_C_1(i__),eps_N_1(i__),eps_O_1(i__),
+ eps_Ne_1(i__),eps_Na_1(i__),eps_Mg_1(i__),eps_Si_1(i__),
+ H1mf(i__),He4mf(i__),Minit_1(i__),vinit(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end