Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/797/116 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-19
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 (nr__) ! (deg) Right Ascension J2000
real*8 DEdeg (nr__) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*10 Name (nr__) ! Cluster name
character*1 f_Name (nr__) ! [*] Indicates a cluster analyzed in
* Colucci et al. 2009, Cat. J/ApJ/704/385 (1)
integer*4 RAh (nr__) ! (h) Hour of Right Ascension (J2000)
integer*4 RAm (nr__) ! (min) Minute of Right Ascension (J2000)
real*4 RAs (nr__) ! (s) Second of Right Ascension (J2000)
character*1 DE_ (nr__) ! Sign of the Declination (J2000)
integer*4 DEd (nr__) ! (deg) Degree of Declination (J2000)
integer*4 DEm (nr__) ! (arcmin) Arcminute of Declination (J2000)
real*4 DEs (nr__) ! (arcsec) Arcsecond of Declination (J2000)
real*4 Vmag (nr__) ! (mag) [14.15/16.91]? V band magnitude
real*4 E_B_V (nr__) ! (mag) [0.05/0.46]? Reddening
character*1 f_E_B_V (nr__) ! [a-d] Flag on E(B-V) (2)
real*4 RM31 (nr__) ! (kpc) [1.77/117.05]? Projected galactocentric
* radius from M31
integer*4 Obs_Y (nr__) ! (yr) Year of Observation
character*3 Obs_M (nr__) ! ("month") Month of Observation
real*4 Texp (nr__) ! (h) [1/6] Exposure time
integer*4 S_N (nr__) ! [54/114]? Signal-to-Noise ratio
character*1 f_S_N (nr__) ! [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 (nr__1) ! Cluster name
real*4 sigma (nr__1) ! (km/s) [4.73/25.73] Velocity dispersion
real*4 e_sigma (nr__1) ! (km/s) [0.37/1.24] Error on sigma
real*4 RV (nr__1) ! (km/s) [-555.5/-53.1] Heliocentric radial velocity
real*4 e_RV (nr__1) ! (km/s) [0.1/0.6] Error on RV
real*4 sigmal (nr__1) ! (km/s) [6.8/27.2]? Velocity dispersion
* from the literature
real*4 e_sigmal (nr__1) ! (km/s) [0.3/2.2]? Error on sigmal
real*4 RVl (nr__1) ! (km/s) [-539/-55.8] Radial velocity
* from the literature
real*4 e_RVl (nr__1) ! (km/s) [0.3/48] Error on RVl
character*7 Refs (nr__1) ! 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 (nr__2) ! Cluster Name
real*4 v_Fe_H_ (nr__2) ! [-1.74/-0.26] [Fe/H] from the equivalent width
real*4 e__Fe_H_ (nr__2) ! [0.15/0.4] Error on [Fe/H]
real*4 v_Fe_H_s (nr__2) ! [-1.81/-0.16] [Fe/H] from line synthesis
real*4 e__Fe_H_s (nr__2) ! [0.1/0.2] Error on [Fe/H]s
real*4 v (nr__2) ! (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 (nr__3) ! (0.1nm) [4005.25/8239.13] Wavelength in Angstroms (1)
real*4 EP (nr__3) ! (eV) [0/4.77] Excitation potential
real*4 log_gf (nr__3) ! ([-]) [-5.28/0.65] Log oscillator strength times
* degeneracy
real*4 EWB006 (nr__3) ! (10-13m) [35.7/149.2]? FeI equivalent width in B006 (2)
real*4 EWB048 (nr__3) ! (10-13m) [34.1/148.9]? FeI equivalent width in B048 (2)
real*4 EWB088 (nr__3) ! (10-13m) [45.4/150.1]? FeI equivalent width in B088 (2)
real*4 EWB232 (nr__3) ! (10-13m) [50.5/149.6]? FeI equivalent width in B232 (2)
real*4 EWB235 (nr__3) ! (10-13m) [27.3/144.6]? FeI equivalent width in B235 (2)
real*4 EWB240 (nr__3) ! (10-13m) [32.1/146.8]? FeI equivalent width in B240 (2)
real*4 EWB311 (nr__3) ! (10-13m) [41.5/139.8]? FeI equivalent width in B311 (2)
real*4 EWB312 (nr__3) ! (10-13m) [40.3/199.7]? FeI equivalent width in B312 (2)
real*4 EWB383 (nr__3) ! (10-13m) [22.4/146.4]? FeI equivalent width in B383 (2)
real*4 EWB403 (nr__3) ! (10-13m) [37.2/189.7]? FeI equivalent width in B403 (2)
real*4 EWB514 (nr__3) ! (10-13m) [38.5/143.6]? FeI equivalent width in B514 (2)
real*4 EWG002 (nr__3) ! (10-13m) [35.9/143.5]? FeI equivalent width in G002 (2)
real*4 EWG327 (nr__3) ! (10-13m) [34.1/99.2]? FeI equivalent width in G327 (2)
real*4 EWMCGC5 (nr__3) ! (10-13m) [36.4/144.3]? FeI equivalent width in MCGC5 (2)
real*4 EWMGC1 (nr__3) ! (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 (nr__4) ! (0.1nm) [4294.14/7583.79] Wavelength (in Angstroms)
real*4 AB012 (nr__4) ! ([-]) [5.2/6.2]? FeI synthesis log abundance in B012 (1)
real*4 AB029 (nr__4) ! ([-]) [6.7/7.4]? FeI synthesis log abundance in B029 (1)
real*4 AB034 (nr__4) ! ([-]) [6.3/7.1]? FeI synthesis log abundance in B034 (1)
real*4 AB110 (nr__4) ! ([-]) [6.5/7.1]? FeI synthesis log abundance in B110 (1)
real*4 AB163 (nr__4) ! ([-]) [6.5/7.3]? FeI synthesis log abundance in B163 (1)
real*4 AB171 (nr__4) ! ([-]) [6.8/7.9]? FeI synthesis log abundance in B171 (1)
real*4 AB182 (nr__4) ! ([-]) [5.9/7.0]? FeI synthesis log abundance in B182 (1)
real*4 AB193 (nr__4) ! ([-]) [6.8/7.6]? FeI synthesis log abundance in B193 (1)
real*4 AB225 (nr__4) ! ([-]) [6.6/7.3]? FeI synthesis log abundance in B225 (1)
real*4 AB384 (nr__4) ! ([-]) [6.5/7.1]? FeI synthesis log abundance in B384 (1)
real*4 AB457 (nr__4) ! ([-]) [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 (nr__5) ! Cluster name
character*1 f_Name_1 (nr__5) ! [*] Indicates a cluster analyzed in Colucci
* et al. 2009, Cat. J/ApJ/704/385 (1)
real*4 Age (nr__5) ! (Gyr) Age
real*4 DAge (nr__5) ! (Gyr) Delta Age
real*4 v_FeI_H_ (nr__5) ! ([-]) Abundance [FeI/H]
integer*4 N (nr__5) ! Number of FeI lines measured
real*4 e_N (nr__5) ! Statistical error of the mean FeI abundance (2)
real*4 sAge (nr__5) ! Uncertainty due to the dependence of the final
* abundance on the age of the CMD (2)
real*4 sT (nr__5) ! Final total uncertainty in abundance (2)
character*3 Fin (nr__5) ! Final analysis performed with GETJOB (EW)
* or line synthesis (SYN)
character*1 Com (nr__5) ! 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 (nr__6) ! Species identifier
real*8 lambda_2 (nr__6) ! (0.1nm) [4122.66/7932.35] Wavelength in Angstroms
real*4 EP_1 (nr__6) ! (eV) [0/5.96] Excitation potential
real*4 log_gf_1 (nr__6) ! ([-]) [-5.569/0.504]? Log of oscillator strength
* times degeneracy
real*4 B006 (nr__6) ! ([-]) [4.33/7.48]? 12+log(X/H) abundance in B006
real*4 B012 (nr__6) ! ([-]) [3.51/7.59]? 12+log(X/H) abundance in B012
real*4 B029 (nr__6) ! ([-]) [4.25/7.3]? 12+log(X/H) abundance in B029
real*4 B034 (nr__6) ! ([-]) [3.99/7.57]? 12+log(X/H) abundance in B034
real*4 B045 (nr__6) ! ([-]) [4.14/7.25]? 12+log(X/H) abundance in B045
real*4 B048 (nr__6) ! ([-]) [3.79/7.2]? 12+log(X/H) abundance in B048
real*4 B088 (nr__6) ! ([-]) [2.81/7.15]? 12+log(X/H) abundance in B088
real*4 B110 (nr__6) ! ([-]) [4.04/7.19]? 12+log(X/H) abundance in B110
real*4 B163 (nr__6) ! ([-]) [4.28/7.63]? 12+log(X/H) abundance in B163
real*4 B171 (nr__6) ! ([-]) [3.99/7.73]? 12+log(X/H) abundance in B171
real*4 B182 (nr__6) ! ([-]) [3.85/7.65]? 12+log(X/H) abundance in B182
real*4 B193 (nr__6) ! ([-]) [4.33/7.82]? 12+log(X/H) abundance in B193
real*4 B225 (nr__6) ! ([-]) [4.65/7.93]? 12+log(X/H) abundance in B225
real*4 B232 (nr__6) ! ([-]) [3/7.41]? 12+log(X/H) abundance in B232
real*4 B235 (nr__6) ! ([-]) [3.96/6.81]? 12+log(X/H) abundance in B235
real*4 B240 (nr__6) ! ([-]) [3.34/7.22]? 12+log(X/H) abundance in B240
real*4 B311 (nr__6) ! ([-]) [3.16/6.5]? 12+log(X/H) abundance in B311
real*4 B312 (nr__6) ! ([-]) [3.75/6.99]? 12+log(X/H) abundance in B312
real*4 B358 (nr__6) ! ([-]) [2.64/6.91]? 12+log(X/H) abundance in B358
real*4 B381 (nr__6) ! ([-]) [3.92/6.77]? 12+log(X/H) abundance in B381
real*4 B383 (nr__6) ! ([-]) [4.17/7.22]? 12+log(X/H) abundance in B383
real*4 B384 (nr__6) ! ([-]) [4.17/7.23]? 12+log(X/H) abundance in B384
real*4 B386 (nr__6) ! ([-]) [3.71/7.33]? 12+log(X/H) abundance in B386
real*4 B403 (nr__6) ! ([-]) [3.83/7.08]? 12+log(X/H) abundance in B403
real*4 B405 (nr__6) ! ([-]) [3.49/7.48]? 12+log(X/H) abundance in B405
real*4 B457 (nr__6) ! ([-]) [3.57/6.52]? 12+log(X/H) abundance in B457
real*4 B514 (nr__6) ! ([-]) [3.2/6.51]? 12+log(X/H) abundance in B514
real*4 G002 (nr__6) ! ([-]) [2.77/6.12]? 12+log(X/H) abundance in G002
real*4 G327 (nr__6) ! ([-]) [3.36/6.32]? 12+log(X/H) abundance in G327
real*4 MCGC5 (nr__6) ! ([-]) [3.21/6.67]? 12+log(X/H) abundance in MCGC5
real*4 MGC1 (nr__6) ! ([-]) [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 (nr__7) ! Cluster name
real*4 v_FeII_H_ (nr__7) ! ([-]) Mean [FeII/H] abundance (G1)
real*4 e__FeII_H_ (nr__7) ! ([-]) Statistical error on [FeII/H] (G2)
integer*4 o_FeII (nr__7) ! Number of FeII lines
real*4 eA_FeII_H_ (nr__7) ! ([-]) Uncertainy on [FeII/H] due to age (G3)
real*4 v_NaI_Fe_ (nr__7) ! ([-]) ? Mean abundance ratio [NaI/FeI] (G1)
real*4 e__NaI_Fe_ (nr__7) ! ([-]) ? Statistical error on NaI (G2)
integer*4 o_NaI (nr__7) ! ? Number of NaI lines
real*4 eA_NaI_Fe_ (nr__7) ! ([-]) ? Uncertainy on [NaI/Fe] due to age (G3)
real*4 v_MgI_Fe_ (nr__7) ! ([-]) Mean abundance ratio [MgI/Fe] (G1)
real*4 e__MgI_Fe_ (nr__7) ! ([-]) Statistical error on [MgI/Fe] (G2)
integer*4 o_MgI (nr__7) ! Number of MgI lines
real*4 eA_MgI_Fe_ (nr__7) ! ([-]) Uncertainy on [MgI/Fe] due to age (G3)
real*4 v_AlI_Fe_ (nr__7) ! ([-]) ? Mean abundance ratio [AlI/Fe] (G1)
real*4 e__AlI_Fe_ (nr__7) ! ([-]) ? Statistical error on AlI (G2)
integer*4 o_AlI (nr__7) ! ? Number of AlI lines
real*4 eA_AlI_Fe_ (nr__7) ! ([-]) ? 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 (nr__8) ! Cluster name
real*4 v_CaI_Fe_ (nr__8) ! ([-]) Mean [CaI/Fe] abundance (G1)
real*4 e__CaI_Fe_ (nr__8) ! ([-]) Statistical error on [CaI/Fe] (G2)
integer*4 o_Ca1 (nr__8) ! Number of CaI lines
real*4 eACa1 (nr__8) ! ([-]) Uncertainy on [CaI/Fe] due to age (G3)
real*4 v_SiI_Fe_ (nr__8) ! ([-]) ? Abundance [SiI/Fe] abundance (G1)
real*4 e__SiI_Fe_ (nr__8) ! ([-]) ? Statistical error on [SiI/Fe] (G2)
integer*4 o_Si1 (nr__8) ! ? Number of Si1 lines
real*4 eA_SiI_Fe_ (nr__8) ! ([-]) ? Uncertainy on [SiI/Fe] due to age (G3)
real*4 v_TiI_FeI_ (nr__8) ! ([-]) ? Abundance [TiI/FeI] abundance (G1)
real*4 e__TiI_FeI_(nr__8) ! ([-]) ? Statistical error on [TiI/FeI] (G2)
integer*4 o_TiI (nr__8) ! ? Number of TiI lines
real*4 eA_TiI_FeI_(nr__8) ! ([-]) ? Uncertainy on [TiI/FeI] due to age (G3)
real*4 v_TiII_FeII_(nr__8) ! ([-]) Abundance [TiII/FeII] abundance (G1)
real*4 e__TiII_FeII_(nr__8) ! ([-]) Statistical error on [TiII/FeII] (G2)
integer*4 o_TiII (nr__8) ! Number of TiII lines
real*4 eA_TiII_FeII_(nr__8) ! ([-]) 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(i__),f_Name(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),
+ DEd(i__),DEm(i__),DEs(i__),Vmag(i__),E_B_V(i__),f_E_B_V(i__),
+ RM31(i__),Obs_Y(i__),Obs_M(i__),Texp(i__),S_N(i__),f_S_N(i__)
if(ar__(36:40) .EQ. '') Vmag(i__) = rNULL__
if(ar__(42:45) .EQ. '') E_B_V(i__) = rNULL__
if(ar__(48:53) .EQ. '') RM31(i__) = rNULL__
if(ar__(68:70) .EQ. '') S_N(i__) = iNULL__
RAdeg(i__) = rNULL__
DEdeg(i__) = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15.
if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4.
if(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240.
if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__)
if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60.
if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600.
if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__)
c ..............Just test output...........
write(6,1)
+ Name(i__),f_Name(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__),
+ DEd(i__),DEm(i__),DEs(i__),Vmag(i__),E_B_V(i__),f_E_B_V(i__),
+ RM31(i__),Obs_Y(i__),Obs_M(i__),Texp(i__),S_N(i__),f_S_N(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__)
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(i__),sigma(i__),e_sigma(i__),RV(i__),e_RV(i__),
+ sigmal(i__),e_sigmal(i__),RVl(i__),e_RVl(i__),Refs(i__)
if(ar__1(34:37) .EQ. '') sigmal(i__) = rNULL__
if(ar__1(39:41) .EQ. '') e_sigmal(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Name_1(i__),sigma(i__),e_sigma(i__),RV(i__),e_RV(i__),
+ sigmal(i__),e_sigmal(i__),RVl(i__),e_RVl(i__),Refs(i__)
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(i__),v_Fe_H_(i__),e__Fe_H_(i__),v_Fe_H_s(i__),
+ e__Fe_H_s(i__),v(i__)
c ..............Just test output...........
write(6,3)
+ Name_2(i__),v_Fe_H_(i__),e__Fe_H_(i__),v_Fe_H_s(i__),
+ e__Fe_H_s(i__),v(i__)
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(i__),EP(i__),log_gf(i__),EWB006(i__),EWB048(i__),
+ EWB088(i__),EWB232(i__),EWB235(i__),EWB240(i__),EWB311(i__),
+ EWB312(i__),EWB383(i__),EWB403(i__),EWB514(i__),EWG002(i__),
+ EWG327(i__),EWMCGC5(i__),EWMGC1(i__)
if(ar__3(23:27) .EQ. '') EWB006(i__) = rNULL__
if(ar__3(29:33) .EQ. '') EWB048(i__) = rNULL__
if(ar__3(35:39) .EQ. '') EWB088(i__) = rNULL__
if(ar__3(41:45) .EQ. '') EWB232(i__) = rNULL__
if(ar__3(47:51) .EQ. '') EWB235(i__) = rNULL__
if(ar__3(53:57) .EQ. '') EWB240(i__) = rNULL__
if(ar__3(59:63) .EQ. '') EWB311(i__) = rNULL__
if(ar__3(65:69) .EQ. '') EWB312(i__) = rNULL__
if(ar__3(71:75) .EQ. '') EWB383(i__) = rNULL__
if(ar__3(77:81) .EQ. '') EWB403(i__) = rNULL__
if(ar__3(83:87) .EQ. '') EWB514(i__) = rNULL__
if(ar__3(89:93) .EQ. '') EWG002(i__) = rNULL__
if(ar__3(95:98) .EQ. '') EWG327(i__) = rNULL__
if(ar__3(100:104) .EQ. '') EWMCGC5(i__) = rNULL__
if(ar__3(106:109) .EQ. '') EWMGC1(i__) = rNULL__
c ..............Just test output...........
write(6,4)
+ lambda(i__),EP(i__),log_gf(i__),EWB006(i__),EWB048(i__),
+ EWB088(i__),EWB232(i__),EWB235(i__),EWB240(i__),EWB311(i__),
+ EWB312(i__),EWB383(i__),EWB403(i__),EWB514(i__),EWG002(i__),
+ EWG327(i__),EWMCGC5(i__),EWMGC1(i__)
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(i__),AB012(i__),AB029(i__),AB034(i__),AB110(i__),
+ AB163(i__),AB171(i__),AB182(i__),AB193(i__),AB225(i__),
+ AB384(i__),AB457(i__)
if(ar__4(10:13) .EQ. '') AB012(i__) = rNULL__
if(ar__4(15:18) .EQ. '') AB029(i__) = rNULL__
if(ar__4(20:23) .EQ. '') AB034(i__) = rNULL__
if(ar__4(25:28) .EQ. '') AB110(i__) = rNULL__
if(ar__4(30:33) .EQ. '') AB163(i__) = rNULL__
if(ar__4(35:38) .EQ. '') AB171(i__) = rNULL__
if(ar__4(40:43) .EQ. '') AB182(i__) = rNULL__
if(ar__4(45:48) .EQ. '') AB193(i__) = rNULL__
if(ar__4(50:53) .EQ. '') AB225(i__) = rNULL__
if(ar__4(55:58) .EQ. '') AB384(i__) = rNULL__
if(ar__4(60:63) .EQ. '') AB457(i__) = rNULL__
c ..............Just test output...........
write(6,5)
+ lambda_1(i__),AB012(i__),AB029(i__),AB034(i__),AB110(i__),
+ AB163(i__),AB171(i__),AB182(i__),AB193(i__),AB225(i__),
+ AB384(i__),AB457(i__)
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(i__),f_Name_1(i__),Age(i__),DAge(i__),v_FeI_H_(i__),
+ N(i__),e_N(i__),sAge(i__),sT(i__),Fin(i__),Com(i__)
c ..............Just test output...........
write(6,6)
+ Name_3(i__),f_Name_1(i__),Age(i__),DAge(i__),v_FeI_H_(i__),
+ N(i__),e_N(i__),sAge(i__),sT(i__),Fin(i__),Com(i__)
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(i__),lambda_2(i__),EP_1(i__),log_gf_1(i__),B006(i__),
+ B012(i__),B029(i__),B034(i__),B045(i__),B048(i__),B088(i__),
+ B110(i__),B163(i__),B171(i__),B182(i__),B193(i__),B225(i__),
+ B232(i__),B235(i__),B240(i__),B311(i__),B312(i__),B358(i__),
+ B381(i__),B383(i__),B384(i__),B386(i__),B403(i__),B405(i__),
+ B457(i__),B514(i__),G002(i__),G327(i__),MCGC5(i__),MGC1(i__)
if(ar__6(22:27) .EQ. '') log_gf_1(i__) = rNULL__
if(ar__6(29:32) .EQ. '') B006(i__) = rNULL__
if(ar__6(34:37) .EQ. '') B012(i__) = rNULL__
if(ar__6(39:42) .EQ. '') B029(i__) = rNULL__
if(ar__6(44:47) .EQ. '') B034(i__) = rNULL__
if(ar__6(49:52) .EQ. '') B045(i__) = rNULL__
if(ar__6(54:57) .EQ. '') B048(i__) = rNULL__
if(ar__6(59:62) .EQ. '') B088(i__) = rNULL__
if(ar__6(64:67) .EQ. '') B110(i__) = rNULL__
if(ar__6(69:72) .EQ. '') B163(i__) = rNULL__
if(ar__6(74:77) .EQ. '') B171(i__) = rNULL__
if(ar__6(79:82) .EQ. '') B182(i__) = rNULL__
if(ar__6(84:87) .EQ. '') B193(i__) = rNULL__
if(ar__6(89:92) .EQ. '') B225(i__) = rNULL__
if(ar__6(94:97) .EQ. '') B232(i__) = rNULL__
if(ar__6(99:102) .EQ. '') B235(i__) = rNULL__
if(ar__6(104:107) .EQ. '') B240(i__) = rNULL__
if(ar__6(109:112) .EQ. '') B311(i__) = rNULL__
if(ar__6(114:117) .EQ. '') B312(i__) = rNULL__
if(ar__6(119:122) .EQ. '') B358(i__) = rNULL__
if(ar__6(124:127) .EQ. '') B381(i__) = rNULL__
if(ar__6(129:132) .EQ. '') B383(i__) = rNULL__
if(ar__6(134:137) .EQ. '') B384(i__) = rNULL__
if(ar__6(139:142) .EQ. '') B386(i__) = rNULL__
if(ar__6(144:147) .EQ. '') B403(i__) = rNULL__
if(ar__6(149:152) .EQ. '') B405(i__) = rNULL__
if(ar__6(154:157) .EQ. '') B457(i__) = rNULL__
if(ar__6(159:162) .EQ. '') B514(i__) = rNULL__
if(ar__6(164:167) .EQ. '') G002(i__) = rNULL__
if(ar__6(169:172) .EQ. '') G327(i__) = rNULL__
if(ar__6(174:177) .EQ. '') MCGC5(i__) = rNULL__
if(ar__6(179:182) .EQ. '') MGC1(i__) = rNULL__
c ..............Just test output...........
write(6,7)
+ Species(i__),lambda_2(i__),EP_1(i__),log_gf_1(i__),B006(i__),
+ B012(i__),B029(i__),B034(i__),B045(i__),B048(i__),B088(i__),
+ B110(i__),B163(i__),B171(i__),B182(i__),B193(i__),B225(i__),
+ B232(i__),B235(i__),B240(i__),B311(i__),B312(i__),B358(i__),
+ B381(i__),B383(i__),B384(i__),B386(i__),B403(i__),B405(i__),
+ B457(i__),B514(i__),G002(i__),G327(i__),MCGC5(i__),MGC1(i__)
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(i__),v_FeII_H_(i__),e__FeII_H_(i__),o_FeII(i__),
+ eA_FeII_H_(i__),v_NaI_Fe_(i__),e__NaI_Fe_(i__),o_NaI(i__),
+ eA_NaI_Fe_(i__),v_MgI_Fe_(i__),e__MgI_Fe_(i__),o_MgI(i__),
+ eA_MgI_Fe_(i__),v_AlI_Fe_(i__),e__AlI_Fe_(i__),o_AlI(i__),
+ eA_AlI_Fe_(i__)
if(ar__7(31:35) .EQ. '') v_NaI_Fe_(i__) = rNULL__
if(ar__7(37:40) .EQ. '') e__NaI_Fe_(i__) = rNULL__
if(ar__7(42:42) .EQ. '') o_NaI(i__) = iNULL__
if(ar__7(44:47) .EQ. '') eA_NaI_Fe_(i__) = rNULL__
if(ar__7(67:70) .EQ. '') v_AlI_Fe_(i__) = rNULL__
if(ar__7(72:75) .EQ. '') e__AlI_Fe_(i__) = rNULL__
if(ar__7(77:77) .EQ. '') o_AlI(i__) = iNULL__
if(ar__7(79:82) .EQ. '') eA_AlI_Fe_(i__) = rNULL__
c ..............Just test output...........
write(6,8)
+ Name_4(i__),v_FeII_H_(i__),e__FeII_H_(i__),o_FeII(i__),
+ eA_FeII_H_(i__),v_NaI_Fe_(i__),e__NaI_Fe_(i__),o_NaI(i__),
+ eA_NaI_Fe_(i__),v_MgI_Fe_(i__),e__MgI_Fe_(i__),o_MgI(i__),
+ eA_MgI_Fe_(i__),v_AlI_Fe_(i__),e__AlI_Fe_(i__),o_AlI(i__),
+ eA_AlI_Fe_(i__)
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(i__),v_CaI_Fe_(i__),e__CaI_Fe_(i__),o_Ca1(i__),
+ eACa1(i__),v_SiI_Fe_(i__),e__SiI_Fe_(i__),o_Si1(i__),
+ eA_SiI_Fe_(i__),v_TiI_FeI_(i__),e__TiI_FeI_(i__),o_TiI(i__),
+ eA_TiI_FeI_(i__),v_TiII_FeII_(i__),e__TiII_FeII_(i__),
+ o_TiII(i__),eA_TiII_FeII_(i__)
if(ar__8(31:34) .EQ. '') v_SiI_Fe_(i__) = rNULL__
if(ar__8(36:39) .EQ. '') e__SiI_Fe_(i__) = rNULL__
if(ar__8(41:41) .EQ. '') o_Si1(i__) = iNULL__
if(ar__8(43:46) .EQ. '') eA_SiI_Fe_(i__) = rNULL__
if(ar__8(48:52) .EQ. '') v_TiI_FeI_(i__) = rNULL__
if(ar__8(54:57) .EQ. '') e__TiI_FeI_(i__) = rNULL__
if(ar__8(59:60) .EQ. '') o_TiI(i__) = iNULL__
if(ar__8(62:65) .EQ. '') eA_TiI_FeI_(i__) = rNULL__
c ..............Just test output...........
write(6,9)
+ Name_5(i__),v_CaI_Fe_(i__),e__CaI_Fe_(i__),o_Ca1(i__),
+ eACa1(i__),v_SiI_Fe_(i__),e__SiI_Fe_(i__),o_Si1(i__),
+ eA_SiI_Fe_(i__),v_TiI_FeI_(i__),e__TiI_FeI_(i__),o_TiI(i__),
+ eA_TiI_FeI_(i__),v_TiII_FeII_(i__),e__TiII_FeII_(i__),
+ o_TiII(i__),eA_TiII_FeII_(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end