Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/797/116 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-07
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/ApJ/797/116 Chemical properties of M31 star clusters (Colucci+, 2014)
*================================================================================
*The detailed chemical properties of M31 star clusters.
*I. Fe, alpha and light elements.
* Colucci J.E., Bernstein R.A., Cohen J.G.
* <Astrophys. J., 797, 116 (2014)>
* =2014ApJ...797..116C (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Observations and cluster properties
integer*4 nr__
parameter (nr__=32) ! Number of records
character*71 ar__ ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg ! (deg) Right Ascension J2000
real*8 DEdeg ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*10 Name ! Cluster name
character*1 f_Name ! [*] Indicates a cluster analyzed in
* Colucci et al. 2009, Cat. J/ApJ/704/385 (1)
integer*4 RAh ! (h) Hour of Right Ascension (J2000)
integer*4 RAm ! (min) Minute of Right Ascension (J2000)
real*4 RAs ! (s) Second of Right Ascension (J2000)
character*1 DE_ ! Sign of the Declination (J2000)
integer*4 DEd ! (deg) Degree of Declination (J2000)
integer*4 DEm ! (arcmin) Arcminute of Declination (J2000)
real*4 DEs ! (arcsec) Arcsecond of Declination (J2000)
real*4 Vmag ! (mag) [14.15/16.91]? V band magnitude
real*4 E_B_V ! (mag) [0.05/0.46]? Reddening
character*1 f_E_B_V ! [a-d] Flag on E(B-V) (2)
real*4 RM31 ! (kpc) [1.77/117.05]? Projected galactocentric
* radius from M31
integer*4 Obs_Y ! (yr) Year of Observation
character*3 Obs_M ! ("month") Month of Observation
real*4 Texp ! (h) [1/6] Exposure time
integer*4 S_N ! [54/114]? Signal-to-Noise ratio
character*1 f_S_N ! [e] Indicates S/N of combined 2009 and 2011
* spectrum
*Note (1): Results from Colucci et al. 2009 (Cat. J/ApJ/704/385) are reproduced
* here for completeness.
*Note (2): Reddening values are taken from Caldwell et al. (2011, J/AJ/141/61),
* with exceptions defined as follows:
* a = Fan et al. (2008MNRAS.385.1973F),
* b = Mackey et al. (2007ApJ...655L..85M),
* c = Barmby et al. (2007, Cat. J/AJ/133/2764),
* d = Mackey et al. (2010MNRAS.401..533M)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Velocity dispersions and radial velocities
integer*4 nr__1
parameter (nr__1=26) ! Number of records
character*61 ar__1 ! Full-size record
character*10 Name_1 ! Cluster name
real*4 sigma ! (km/s) [4.73/25.73] Velocity dispersion
real*4 e_sigma ! (km/s) [0.37/1.24] Error on sigma
real*4 RV ! (km/s) [-555.5/-53.1] Heliocentric radial velocity
real*4 e_RV ! (km/s) [0.1/0.6] Error on RV
real*4 sigmal ! (km/s) [6.8/27.2]? Velocity dispersion
* from the literature
real*4 e_sigmal ! (km/s) [0.3/2.2]? Error on sigmal
real*4 RVl ! (km/s) [-539/-55.8] Radial velocity
* from the literature
real*4 e_RVl ! (km/s) [0.3/48] Error on RVl
character*7 Refs ! References (1)
*Note (1): Where two references are listed, we have calculated the average of
* the two results and an error equal to the standard deviation in the
* mean. The references are defined as follows:
* 1 = Strader et al. 2011AJ....142....8S.
* 2 = Galleti et al. 2004A&A...416..917G.
* 3 = Dubath & Grillmair 1997A&A...321..379D.
* 4 = Alves-Brito et al. 2009MNRAS.395L..34A.
* 5 = Djorgovski et al. 1997ApJ...474L..19D.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! EW versus synthesis analysis
integer*4 nr__2
parameter (nr__2=14) ! Number of records
character*37 ar__2 ! Full-size record
character*10 Name_2 ! Cluster Name
real*4 v_Fe_H_ ! [-1.74/-0.26] [Fe/H] from the equivalent width
real*4 e__Fe_H_ ! [0.15/0.4] Error on [Fe/H]
real*4 v_Fe_H_s ! [-1.81/-0.16] [Fe/H] from line synthesis
real*4 e__Fe_H_s ! [0.1/0.2] Error on [Fe/H]s
real*4 v ! (km/s) [6.5/19.5] Velocity dispersion
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! M31 GC FeI equivalent widths
integer*4 nr__3
parameter (nr__3=164) ! Number of records
character*109 ar__3 ! Full-size record
real*8 lambda ! (0.1nm) [4005.25/8239.13] Wavelength in Angstroms (1)
real*4 EP ! (eV) [0/4.77] Excitation potential
real*4 log_gf ! ([-]) [-5.28/0.65] Log oscillator strength times
* degeneracy
real*4 EWB006 ! (10-13m) [35.7/149.2]? FeI equivalent width in B006 (2)
real*4 EWB048 ! (10-13m) [34.1/148.9]? FeI equivalent width in B048 (2)
real*4 EWB088 ! (10-13m) [45.4/150.1]? FeI equivalent width in B088 (2)
real*4 EWB232 ! (10-13m) [50.5/149.6]? FeI equivalent width in B232 (2)
real*4 EWB235 ! (10-13m) [27.3/144.6]? FeI equivalent width in B235 (2)
real*4 EWB240 ! (10-13m) [32.1/146.8]? FeI equivalent width in B240 (2)
real*4 EWB311 ! (10-13m) [41.5/139.8]? FeI equivalent width in B311 (2)
real*4 EWB312 ! (10-13m) [40.3/199.7]? FeI equivalent width in B312 (2)
real*4 EWB383 ! (10-13m) [22.4/146.4]? FeI equivalent width in B383 (2)
real*4 EWB403 ! (10-13m) [37.2/189.7]? FeI equivalent width in B403 (2)
real*4 EWB514 ! (10-13m) [38.5/143.6]? FeI equivalent width in B514 (2)
real*4 EWG002 ! (10-13m) [35.9/143.5]? FeI equivalent width in G002 (2)
real*4 EWG327 ! (10-13m) [34.1/99.2]? FeI equivalent width in G327 (2)
real*4 EWMCGC5 ! (10-13m) [36.4/144.3]? FeI equivalent width in MCGC5 (2)
real*4 EWMGC1 ! (10-13m) [47.7/97.8]? FeI equivalent width in MGC1 (2)
*Note (1): Lines listed twice were measured in two orders with overlapping
* wavelength coverage.
*Note (2): In units of milli-Angstroms.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! M31 GC FeI synthesis abundances
integer*4 nr__4
parameter (nr__4=83) ! Number of records
character*63 ar__4 ! Full-size record
real*8 lambda_1 ! (0.1nm) [4294.14/7583.79] Wavelength (in Angstroms)
real*4 AB012 ! ([-]) [5.2/6.2]? FeI synthesis log abundance in B012 (1)
real*4 AB029 ! ([-]) [6.7/7.4]? FeI synthesis log abundance in B029 (1)
real*4 AB034 ! ([-]) [6.3/7.1]? FeI synthesis log abundance in B034 (1)
real*4 AB110 ! ([-]) [6.5/7.1]? FeI synthesis log abundance in B110 (1)
real*4 AB163 ! ([-]) [6.5/7.3]? FeI synthesis log abundance in B163 (1)
real*4 AB171 ! ([-]) [6.8/7.9]? FeI synthesis log abundance in B171 (1)
real*4 AB182 ! ([-]) [5.9/7.0]? FeI synthesis log abundance in B182 (1)
real*4 AB193 ! ([-]) [6.8/7.6]? FeI synthesis log abundance in B193 (1)
real*4 AB225 ! ([-]) [6.6/7.3]? FeI synthesis log abundance in B225 (1)
real*4 AB384 ! ([-]) [6.5/7.1]? FeI synthesis log abundance in B384 (1)
real*4 AB457 ! ([-]) [5.8/6.9]? FeI synthesis log abundance in B457 (1)
*Note (1): Abundances for individual lines are measured in steps of 0.1dex from
* the mean abundance using all lines. The typical statistical
* measurement uncertainty of individual lines is 0.05dex. Note that
* these statistical uncertainties do not dominate the total uncertainty.
* See text for further discussion. Abundances for each GC are calculated
* using its CMD solution that has the oldest age.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! M31 GC ages and FeI abundances
integer*4 nr__5
parameter (nr__5=31) ! Number of records
character*50 ar__5 ! Full-size record
character*10 Name_3 ! Cluster name
character*1 f_Name_1 ! [*] Indicates a cluster analyzed in Colucci
* et al. 2009, Cat. J/ApJ/704/385 (1)
real*4 Age ! (Gyr) Age
real*4 DAge ! (Gyr) Delta Age
real*4 v_FeI_H_ ! ([-]) Abundance [FeI/H]
integer*4 N ! Number of FeI lines measured
real*4 e_N ! Statistical error of the mean FeI abundance (2)
real*4 sAge ! Uncertainty due to the dependence of the final
* abundance on the age of the CMD (2)
real*4 sT ! Final total uncertainty in abundance (2)
character*3 Fin ! Final analysis performed with GETJOB (EW)
* or line synthesis (SYN)
character*1 Com ! Denotes whether solar-scaled (S) or
* {alpha}-enhanced (A) isochrones and stellar
* atmospheres were used in the final analysis.
*Note (1): Results from Colucci et al. 2009 (Cat. J/ApJ/704/385) are reproduced
* here for completeness.
*Note (2): Final total uncertainty in abundance:
* {sigma}_T_=(({sigma}_N_/sqrt(N-1))^2^+{sigma}_Age_^2^)(1/2)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table7.dat' ! M31 GC line abundances
integer*4 nr__6
parameter (nr__6=100) ! Number of records
character*182 ar__6 ! Full-size record
character*5 Species ! Species identifier
real*8 lambda_2 ! (0.1nm) [4122.66/7932.35] Wavelength in Angstroms
real*4 EP_1 ! (eV) [0/5.96] Excitation potential
real*4 log_gf_1 ! ([-]) [-5.569/0.504]? Log of oscillator strength
* times degeneracy
real*4 B006 ! ([-]) [4.33/7.48]? 12+log(X/H) abundance in B006
real*4 B012 ! ([-]) [3.51/7.59]? 12+log(X/H) abundance in B012
real*4 B029 ! ([-]) [4.25/7.3]? 12+log(X/H) abundance in B029
real*4 B034 ! ([-]) [3.99/7.57]? 12+log(X/H) abundance in B034
real*4 B045 ! ([-]) [4.14/7.25]? 12+log(X/H) abundance in B045
real*4 B048 ! ([-]) [3.79/7.2]? 12+log(X/H) abundance in B048
real*4 B088 ! ([-]) [2.81/7.15]? 12+log(X/H) abundance in B088
real*4 B110 ! ([-]) [4.04/7.19]? 12+log(X/H) abundance in B110
real*4 B163 ! ([-]) [4.28/7.63]? 12+log(X/H) abundance in B163
real*4 B171 ! ([-]) [3.99/7.73]? 12+log(X/H) abundance in B171
real*4 B182 ! ([-]) [3.85/7.65]? 12+log(X/H) abundance in B182
real*4 B193 ! ([-]) [4.33/7.82]? 12+log(X/H) abundance in B193
real*4 B225 ! ([-]) [4.65/7.93]? 12+log(X/H) abundance in B225
real*4 B232 ! ([-]) [3/7.41]? 12+log(X/H) abundance in B232
real*4 B235 ! ([-]) [3.96/6.81]? 12+log(X/H) abundance in B235
real*4 B240 ! ([-]) [3.34/7.22]? 12+log(X/H) abundance in B240
real*4 B311 ! ([-]) [3.16/6.5]? 12+log(X/H) abundance in B311
real*4 B312 ! ([-]) [3.75/6.99]? 12+log(X/H) abundance in B312
real*4 B358 ! ([-]) [2.64/6.91]? 12+log(X/H) abundance in B358
real*4 B381 ! ([-]) [3.92/6.77]? 12+log(X/H) abundance in B381
real*4 B383 ! ([-]) [4.17/7.22]? 12+log(X/H) abundance in B383
real*4 B384 ! ([-]) [4.17/7.23]? 12+log(X/H) abundance in B384
real*4 B386 ! ([-]) [3.71/7.33]? 12+log(X/H) abundance in B386
real*4 B403 ! ([-]) [3.83/7.08]? 12+log(X/H) abundance in B403
real*4 B405 ! ([-]) [3.49/7.48]? 12+log(X/H) abundance in B405
real*4 B457 ! ([-]) [3.57/6.52]? 12+log(X/H) abundance in B457
real*4 B514 ! ([-]) [3.2/6.51]? 12+log(X/H) abundance in B514
real*4 G002 ! ([-]) [2.77/6.12]? 12+log(X/H) abundance in G002
real*4 G327 ! ([-]) [3.36/6.32]? 12+log(X/H) abundance in G327
real*4 MCGC5 ! ([-]) [3.21/6.67]? 12+log(X/H) abundance in MCGC5
real*4 MGC1 ! ([-]) [3.36/6.42]? 12+log(X/H) abundance in MGC1
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table8.dat' ! M31 GC mean abundances
integer*4 nr__7
parameter (nr__7=31) ! Number of records
character*82 ar__7 ! Full-size record
character*10 Name_4 ! Cluster name
real*4 v_FeII_H_ ! ([-]) Mean [FeII/H] abundance (G1)
real*4 e__FeII_H_ ! ([-]) Statistical error on [FeII/H] (G2)
integer*4 o_FeII ! Number of FeII lines
real*4 eA_FeII_H_ ! ([-]) Uncertainy on [FeII/H] due to age (G3)
real*4 v_NaI_Fe_ ! ([-]) ? Mean abundance ratio [NaI/FeI] (G1)
real*4 e__NaI_Fe_ ! ([-]) ? Statistical error on NaI (G2)
integer*4 o_NaI ! ? Number of NaI lines
real*4 eA_NaI_Fe_ ! ([-]) ? Uncertainy on [NaI/Fe] due to age (G3)
real*4 v_MgI_Fe_ ! ([-]) Mean abundance ratio [MgI/Fe] (G1)
real*4 e__MgI_Fe_ ! ([-]) Statistical error on [MgI/Fe] (G2)
integer*4 o_MgI ! Number of MgI lines
real*4 eA_MgI_Fe_ ! ([-]) Uncertainy on [MgI/Fe] due to age (G3)
real*4 v_AlI_Fe_ ! ([-]) ? Mean abundance ratio [AlI/Fe] (G1)
real*4 e__AlI_Fe_ ! ([-]) ? Statistical error on AlI (G2)
integer*4 o_AlI ! ? Number of AlI lines
real*4 eA_AlI_Fe_ ! ([-]) ? Uncertainy on [AlI/Fe] due to age (G3)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table9.dat' ! M31 GC abundances continued
integer*4 nr__8
parameter (nr__8=31) ! Number of records
character*84 ar__8 ! Full-size record
character*10 Name_5 ! Cluster name
real*4 v_CaI_Fe_ ! ([-]) Mean [CaI/Fe] abundance (G1)
real*4 e__CaI_Fe_ ! ([-]) Statistical error on [CaI/Fe] (G2)
integer*4 o_Ca1 ! Number of CaI lines
real*4 eACa1 ! ([-]) Uncertainy on [CaI/Fe] due to age (G3)
real*4 v_SiI_Fe_ ! ([-]) ? Abundance [SiI/Fe] abundance (G1)
real*4 e__SiI_Fe_ ! ([-]) ? Statistical error on [SiI/Fe] (G2)
integer*4 o_Si1 ! ? Number of Si1 lines
real*4 eA_SiI_Fe_ ! ([-]) ? Uncertainy on [SiI/Fe] due to age (G3)
real*4 v_TiI_FeI_ ! ([-]) ? Abundance [TiI/FeI] abundance (G1)
real*4 e__TiI_FeI_ ! ([-]) ? Statistical error on [TiI/FeI] (G2)
integer*4 o_TiI ! ? Number of TiI lines
real*4 eA_TiI_FeI_ ! ([-]) ? Uncertainy on [TiI/FeI] due to age (G3)
real*4 v_TiII_FeII_ ! ([-]) Abundance [TiII/FeII] abundance (G1)
real*4 e__TiII_FeII_ ! ([-]) Statistical error on [TiII/FeII] (G2)
integer*4 o_TiII ! Number of TiII lines
real*4 eA_TiII_FeII_ ! ([-]) Uncertainy on [TiII/FeII] due to age (G3)
C=============================================================================
C Loading file 'table1.dat' ! Observations and cluster properties
C Format for file interpretation
1 format(
+ A10,A1,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,F4.1,1X,F5.2,1X,
+ F4.2,A1,1X,F6.2,1X,I4,1X,A3,1X,F3.1,1X,I3,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,32
read(1,'(A71)')ar__
read(ar__,1)
+ Name,f_Name,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Vmag,E_B_V,f_E_B_V,
+ RM31,Obs_Y,Obs_M,Texp,S_N,f_S_N
if(ar__(36:40) .EQ. '') Vmag = rNULL__
if(ar__(42:45) .EQ. '') E_B_V = rNULL__
if(ar__(48:53) .EQ. '') RM31 = rNULL__
if(ar__(68:70) .EQ. '') S_N = iNULL__
RAdeg = rNULL__
DEdeg = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh .GT. -180) RAdeg=RAh*15.
if(RAm .GT. -180) RAdeg=RAdeg+RAm/4.
if(RAs .GT. -180) RAdeg=RAdeg+RAs/240.
if(DEd .GE. 0) DEdeg=DEd
if(DEm .GE. 0) DEdeg=DEdeg+DEm/60.
if(DEs .GE. 0) DEdeg=DEdeg+DEs/3600.
if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg
c ..............Just test output...........
write(6,1)
+ Name,f_Name,RAh,RAm,RAs,DE_,DEd,DEm,DEs,Vmag,E_B_V,f_E_B_V,
+ RM31,Obs_Y,Obs_M,Texp,S_N,f_S_N
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Velocity dispersions and radial velocities
C Format for file interpretation
2 format(
+ A10,1X,F5.2,1X,F4.2,1X,F6.1,1X,F3.1,1X,F4.1,1X,F3.1,1X,F6.1,
+ 1X,F4.1,1X,A7)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,26
read(1,'(A61)')ar__1
read(ar__1,2)
+ Name_1,sigma,e_sigma,RV,e_RV,sigmal,e_sigmal,RVl,e_RVl,Refs
if(ar__1(34:37) .EQ. '') sigmal = rNULL__
if(ar__1(39:41) .EQ. '') e_sigmal = rNULL__
c ..............Just test output...........
write(6,2)
+ Name_1,sigma,e_sigma,RV,e_RV,sigmal,e_sigmal,RVl,e_RVl,Refs
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! EW versus synthesis analysis
C Format for file interpretation
3 format(A10,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F4.1)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,14
read(1,'(A37)')ar__2
read(ar__2,3)Name_2,v_Fe_H_,e__Fe_H_,v_Fe_H_s,e__Fe_H_s,v
c ..............Just test output...........
write(6,3)Name_2,v_Fe_H_,e__Fe_H_,v_Fe_H_s,e__Fe_H_s,v
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! M31 GC FeI equivalent widths
C Format for file interpretation
4 format(
+ F8.3,1X,F5.3,1X,F6.3,1X,F5.1,1X,F5.1,1X,F5.1,1X,F5.1,1X,F5.1,
+ 1X,F5.1,1X,F5.1,1X,F5.1,1X,F5.1,1X,F5.1,1X,F5.1,1X,F5.1,1X,
+ F4.1,1X,F5.1,1X,F4.1)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,164
read(1,'(A109)')ar__3
read(ar__3,4)
+ lambda,EP,log_gf,EWB006,EWB048,EWB088,EWB232,EWB235,EWB240,
+ EWB311,EWB312,EWB383,EWB403,EWB514,EWG002,EWG327,EWMCGC5,
+ EWMGC1
if(ar__3(23:27) .EQ. '') EWB006 = rNULL__
if(ar__3(29:33) .EQ. '') EWB048 = rNULL__
if(ar__3(35:39) .EQ. '') EWB088 = rNULL__
if(ar__3(41:45) .EQ. '') EWB232 = rNULL__
if(ar__3(47:51) .EQ. '') EWB235 = rNULL__
if(ar__3(53:57) .EQ. '') EWB240 = rNULL__
if(ar__3(59:63) .EQ. '') EWB311 = rNULL__
if(ar__3(65:69) .EQ. '') EWB312 = rNULL__
if(ar__3(71:75) .EQ. '') EWB383 = rNULL__
if(ar__3(77:81) .EQ. '') EWB403 = rNULL__
if(ar__3(83:87) .EQ. '') EWB514 = rNULL__
if(ar__3(89:93) .EQ. '') EWG002 = rNULL__
if(ar__3(95:98) .EQ. '') EWG327 = rNULL__
if(ar__3(100:104) .EQ. '') EWMCGC5 = rNULL__
if(ar__3(106:109) .EQ. '') EWMGC1 = rNULL__
c ..............Just test output...........
write(6,4)
+ lambda,EP,log_gf,EWB006,EWB048,EWB088,EWB232,EWB235,EWB240,
+ EWB311,EWB312,EWB383,EWB403,EWB514,EWG002,EWG327,EWMCGC5,
+ EWMGC1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! M31 GC FeI synthesis abundances
C Format for file interpretation
5 format(
+ F8.3,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,
+ 1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,83
read(1,'(A63)')ar__4
read(ar__4,5)
+ lambda_1,AB012,AB029,AB034,AB110,AB163,AB171,AB182,AB193,
+ AB225,AB384,AB457
if(ar__4(10:13) .EQ. '') AB012 = rNULL__
if(ar__4(15:18) .EQ. '') AB029 = rNULL__
if(ar__4(20:23) .EQ. '') AB034 = rNULL__
if(ar__4(25:28) .EQ. '') AB110 = rNULL__
if(ar__4(30:33) .EQ. '') AB163 = rNULL__
if(ar__4(35:38) .EQ. '') AB171 = rNULL__
if(ar__4(40:43) .EQ. '') AB182 = rNULL__
if(ar__4(45:48) .EQ. '') AB193 = rNULL__
if(ar__4(50:53) .EQ. '') AB225 = rNULL__
if(ar__4(55:58) .EQ. '') AB384 = rNULL__
if(ar__4(60:63) .EQ. '') AB457 = rNULL__
c ..............Just test output...........
write(6,5)
+ lambda_1,AB012,AB029,AB034,AB110,AB163,AB171,AB182,AB193,
+ AB225,AB384,AB457
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! M31 GC ages and FeI abundances
C Format for file interpretation
6 format(
+ A10,A1,1X,F4.1,1X,F3.1,1X,F5.2,1X,I2,1X,F4.2,1X,F4.2,1X,F4.2,
+ 1X,A3,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,31
read(1,'(A50)')ar__5
read(ar__5,6)
+ Name_3,f_Name_1,Age,DAge,v_FeI_H_,N,e_N,sAge,sT,Fin,Com
c ..............Just test output...........
write(6,6)
+ Name_3,f_Name_1,Age,DAge,v_FeI_H_,N,e_N,sAge,sT,Fin,Com
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table7.dat' ! M31 GC line abundances
C Format for file interpretation
7 format(
+ A5,1X,F8.3,1X,F5.3,1X,F6.3,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,
+ F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,
+ 1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,
+ F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,
+ 1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table7.dat')
write(6,*) '....Loading file: table7.dat'
do i__=1,100
read(1,'(A182)')ar__6
read(ar__6,7)
+ Species,lambda_2,EP_1,log_gf_1,B006,B012,B029,B034,B045,B048,
+ B088,B110,B163,B171,B182,B193,B225,B232,B235,B240,B311,B312,
+ B358,B381,B383,B384,B386,B403,B405,B457,B514,G002,G327,MCGC5,
+ MGC1
if(ar__6(22:27) .EQ. '') log_gf_1 = rNULL__
if(ar__6(29:32) .EQ. '') B006 = rNULL__
if(ar__6(34:37) .EQ. '') B012 = rNULL__
if(ar__6(39:42) .EQ. '') B029 = rNULL__
if(ar__6(44:47) .EQ. '') B034 = rNULL__
if(ar__6(49:52) .EQ. '') B045 = rNULL__
if(ar__6(54:57) .EQ. '') B048 = rNULL__
if(ar__6(59:62) .EQ. '') B088 = rNULL__
if(ar__6(64:67) .EQ. '') B110 = rNULL__
if(ar__6(69:72) .EQ. '') B163 = rNULL__
if(ar__6(74:77) .EQ. '') B171 = rNULL__
if(ar__6(79:82) .EQ. '') B182 = rNULL__
if(ar__6(84:87) .EQ. '') B193 = rNULL__
if(ar__6(89:92) .EQ. '') B225 = rNULL__
if(ar__6(94:97) .EQ. '') B232 = rNULL__
if(ar__6(99:102) .EQ. '') B235 = rNULL__
if(ar__6(104:107) .EQ. '') B240 = rNULL__
if(ar__6(109:112) .EQ. '') B311 = rNULL__
if(ar__6(114:117) .EQ. '') B312 = rNULL__
if(ar__6(119:122) .EQ. '') B358 = rNULL__
if(ar__6(124:127) .EQ. '') B381 = rNULL__
if(ar__6(129:132) .EQ. '') B383 = rNULL__
if(ar__6(134:137) .EQ. '') B384 = rNULL__
if(ar__6(139:142) .EQ. '') B386 = rNULL__
if(ar__6(144:147) .EQ. '') B403 = rNULL__
if(ar__6(149:152) .EQ. '') B405 = rNULL__
if(ar__6(154:157) .EQ. '') B457 = rNULL__
if(ar__6(159:162) .EQ. '') B514 = rNULL__
if(ar__6(164:167) .EQ. '') G002 = rNULL__
if(ar__6(169:172) .EQ. '') G327 = rNULL__
if(ar__6(174:177) .EQ. '') MCGC5 = rNULL__
if(ar__6(179:182) .EQ. '') MGC1 = rNULL__
c ..............Just test output...........
write(6,7)
+ Species,lambda_2,EP_1,log_gf_1,B006,B012,B029,B034,B045,B048,
+ B088,B110,B163,B171,B182,B193,B225,B232,B235,B240,B311,B312,
+ B358,B381,B383,B384,B386,B403,B405,B457,B514,G002,G327,MCGC5,
+ MGC1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table8.dat' ! M31 GC mean abundances
C Format for file interpretation
8 format(
+ A10,1X,F5.2,1X,F4.2,1X,I2,1X,F4.2,1X,F5.2,1X,F4.2,1X,I1,1X,
+ F4.2,1X,F5.2,1X,F4.2,1X,I1,1X,F4.2,1X,F4.2,1X,F4.2,1X,I1,1X,
+ F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table8.dat')
write(6,*) '....Loading file: table8.dat'
do i__=1,31
read(1,'(A82)')ar__7
read(ar__7,8)
+ Name_4,v_FeII_H_,e__FeII_H_,o_FeII,eA_FeII_H_,v_NaI_Fe_,
+ e__NaI_Fe_,o_NaI,eA_NaI_Fe_,v_MgI_Fe_,e__MgI_Fe_,o_MgI,
+ eA_MgI_Fe_,v_AlI_Fe_,e__AlI_Fe_,o_AlI,eA_AlI_Fe_
if(ar__7(31:35) .EQ. '') v_NaI_Fe_ = rNULL__
if(ar__7(37:40) .EQ. '') e__NaI_Fe_ = rNULL__
if(ar__7(42:42) .EQ. '') o_NaI = iNULL__
if(ar__7(44:47) .EQ. '') eA_NaI_Fe_ = rNULL__
if(ar__7(67:70) .EQ. '') v_AlI_Fe_ = rNULL__
if(ar__7(72:75) .EQ. '') e__AlI_Fe_ = rNULL__
if(ar__7(77:77) .EQ. '') o_AlI = iNULL__
if(ar__7(79:82) .EQ. '') eA_AlI_Fe_ = rNULL__
c ..............Just test output...........
write(6,8)
+ Name_4,v_FeII_H_,e__FeII_H_,o_FeII,eA_FeII_H_,v_NaI_Fe_,
+ e__NaI_Fe_,o_NaI,eA_NaI_Fe_,v_MgI_Fe_,e__MgI_Fe_,o_MgI,
+ eA_MgI_Fe_,v_AlI_Fe_,e__AlI_Fe_,o_AlI,eA_AlI_Fe_
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table9.dat' ! M31 GC abundances continued
C Format for file interpretation
9 format(
+ A10,1X,F5.2,1X,F4.2,1X,I2,1X,F4.2,1X,F4.2,1X,F4.2,1X,I1,1X,
+ F4.2,1X,F5.2,1X,F4.2,1X,I2,1X,F4.2,1X,F5.2,1X,F4.2,1X,I2,1X,
+ F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table9.dat')
write(6,*) '....Loading file: table9.dat'
do i__=1,31
read(1,'(A84)')ar__8
read(ar__8,9)
+ Name_5,v_CaI_Fe_,e__CaI_Fe_,o_Ca1,eACa1,v_SiI_Fe_,e__SiI_Fe_,
+ o_Si1,eA_SiI_Fe_,v_TiI_FeI_,e__TiI_FeI_,o_TiI,eA_TiI_FeI_,
+ v_TiII_FeII_,e__TiII_FeII_,o_TiII,eA_TiII_FeII_
if(ar__8(31:34) .EQ. '') v_SiI_Fe_ = rNULL__
if(ar__8(36:39) .EQ. '') e__SiI_Fe_ = rNULL__
if(ar__8(41:41) .EQ. '') o_Si1 = iNULL__
if(ar__8(43:46) .EQ. '') eA_SiI_Fe_ = rNULL__
if(ar__8(48:52) .EQ. '') v_TiI_FeI_ = rNULL__
if(ar__8(54:57) .EQ. '') e__TiI_FeI_ = rNULL__
if(ar__8(59:60) .EQ. '') o_TiI = iNULL__
if(ar__8(62:65) .EQ. '') eA_TiI_FeI_ = rNULL__
c ..............Just test output...........
write(6,9)
+ Name_5,v_CaI_Fe_,e__CaI_Fe_,o_Ca1,eACa1,v_SiI_Fe_,e__SiI_Fe_,
+ o_Si1,eA_SiI_Fe_,v_TiI_FeI_,e__TiI_FeI_,o_TiI,eA_TiI_FeI_,
+ v_TiII_FeII_,e__TiII_FeII_,o_TiII,eA_TiII_FeII_
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end