Conversion of standardized ReadMe file for
catalog 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.5, on 2013-May-18
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__=-1.e37) ! NULL real number
parameter (iNULL__=-2147483647) ! NULL int number
C=============================================================================
Cat. J/MNRAS/400/317 IC 4665 low-mass members (Jeffries+, 2009)
*================================================================================
*Low-mass members of the young cluster IC 4665 and pre-main-sequence
*lithium depletion.
* Jeffries R.D., Jackson R.J., James D.J., Cargile P.A.
* <Mon. Not. R. Astron. Soc., 400, 317-329 (2009)>
* =2009MNRAS.400..317J
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! The spectroscopic measurements of 452 photometric
candidates in IC 4665
integer*4 nr__
parameter (nr__=452) ! Number of records
character*168 ar__ ! Full-size record
character*5 Name (nr__) ! Name (N_NNN)
real*8 RAdeg (nr__) ! (deg) Right ascension in decimal degrees (J2000)
real*8 DEdeg (nr__) ! (deg) Declination in decimal degrees (J2000)
real*4 Vmag (nr__) ! (mag) AF2/WYFFOS V magnitude
real*4 V_I (nr__) ! (mag) AF2/WYFFOS V-I colour index
integer*4 Nobs (nr__) ! Number of observations for radial velocities
real*4 HJD1 (nr__) ! (d) ? Measurement 1, HJD-2454600
real*4 HRV1 (nr__) ! (km/s) ? Measurement 1 heliocentric radial velocity
real*4 e_HRV1 (nr__) ! (km/s) ? rms uncertainty on RV1
real*4 HJD2 (nr__) ! (d) ? Measurement 2, HJD-2454600
real*4 HRV2 (nr__) ! (km/s) ? Measurement 2 heliocentric radial velocity
real*4 e_HRV2 (nr__) ! (km/s) ? rms uncertainty on RV2
real*4 HJD3 (nr__) ! (d) ? Measurement 3, HJD-2454600
real*4 HRV3 (nr__) ! (km/s) ? Measurement 3 heliocentric radial velocity
real*4 e_HRV3 (nr__) ! (km/s) ? rms uncertainty on RV3
real*4 HJD4 (nr__) ! (d) ? Measurement 4, HJD-2454600
real*4 HRV4 (nr__) ! (km/s) ? Measurement 4 heliocentric radial velocity
real*4 e_HRV4 (nr__) ! (km/s) ? rms uncertainty on RV4
real*8 HRV (nr__) ! (km/s) ?=-9999.0 Mean heliocentric radial velocity
real*4 e_HRV (nr__) ! (km/s) ?=99.0 rms uncertainty on RV
integer*4 SNR (nr__) ! Signal-to-noise ratio of the averaged spectrum
integer*4 EWLi (nr__) ! (0.1pm) ? Measured Li equivalent width
integer*4 e_EWLi (nr__) ! (0.1pm) ?=-99 rms uncertainty on EWLi
real*8 vsini (nr__) ! (km/s) ?=-9999.0 Rotational velocity
real*8 e_vsini (nr__) ! (km/s) ?=-9999.0 rms uncertainty on vsini
integer*4 Var (nr__) ! [0/2]? RV variability flag (0=not variable) (1)
integer*4 CCF (nr__) ! [0/2]? Cross-correlation function with the
* template standard star shows evidence for
* more than one star (1=evidence of binarity)
integer*4 P93 (nr__) ! ? Prosser et al., 1993, Cat. J/AJ/105/1441
* <Cl* IC 4665 P NNN> in Simbad
integer*4 M08 (nr__) ! ? Manzi et al., 2008, Cat. J/A+A/479/141
* <Cl* IC 4665 MRD NNN> in Simbad
*Note (1): RV variability flag as follows:
* 1 = variable at >99.9 per cent confidence
* 0 = not variable
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Properties of candidate RV members of IC 4665
integer*4 nr__1
parameter (nr__1=56) ! Number of records
character*86 ar__1 ! Full-size record
character*5 Name_1 (nr__1) ! Name (N_NNN)
real*4 Vmag_1 (nr__1) ! (mag) AF2/WYFFOS V magnitude
real*4 V_I_1 (nr__1) ! (mag) AF2/WYFFOS V-I colour index
integer*4 EWLi_1 (nr__1) ! (0.1pm) Lithium equivalent width
integer*4 e_EWLi_1 (nr__1) ! (0.1pm) ?=-99 rms uncertainty on EWLi (1)
real*4 RV (nr__1) ! (km/s) Radial velocity
real*4 e_RV (nr__1) ! (km/s) rms uncertainty on RV
real*4 vsini_1 (nr__1) ! (km/s) Rotational velocity
real*4 e_vsini_1 (nr__1) ! (km/s) ?=-99.9 rms uncertainty on vsini (1)
real*4 EWHa (nr__1) ! (0.1nm) H{alpha} equivalent width (positive=emission)
real*4 EWCaI1 (nr__1) ! (0.1nm) Equivalent width of the
* gravity-sensitive [CaI] 6439{AA} line
real*4 EWCaI2 (nr__1) ! (0.1nm) Equivalent width of the
* gravity-sensitive [CaI] 6463{AA} line
real*4 pmRA (nr__1) ! (mas/yr) ? Proper motion along right ascension
* (from NOMAD catalog, Cat. I/297)
real*4 e_pmRA (nr__1) ! (mas/yr) ? rms uncertainty on RA
real*4 pmDE (nr__1) ! (mas/yr) ? Proper motion along declination
* (from NOMAD catalog, Cat. I/297)
real*4 e_pmDE (nr__1) ! (mas/yr) ? rms uncertainty on pmDE
character*1 Mm (nr__1) ! [Y/N?] Final assessment of membership flag
*Note (1): An uncertainty on -99 or -99.9 indicates an upper limit.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Lithium abundances in IC 4665
integer*4 nr__2
parameter (nr__2=44) ! Number of records
character*37 ar__2 ! Full-size record
character*5 Name_2 (nr__2) ! Name (N_NNN)
real*4 V_I_2 (nr__2) ! (mag) AF2/WYFFOS V-I colour index
integer*4 Teff (nr__2) ! (K) Effective temperature
character*1 l_A_Li_LTE (nr__2) ! Limit flag on A(Li)LTE
real*4 A_Li_LTE (nr__2) ! (mag) ? LTE Li abundance, A(Li)=12+log[N(Li)/N(H)]
real*4 e_A_Li_LTE (nr__2) ! (mag) ? rms uncertainty on A(Li)LTE
character*1 l_A_Li_NLTE(nr__2) ! Limit flag on A(Li)NLTE
real*4 A_Li_NLTE (nr__2) ! (mag) ? NLTE Li abundance, A(Li)=12+log[N(Li)/N(H)],
* for stars with Teff>4500K
C=============================================================================
C Loading file 'table1.dat' ! The spectroscopic measurements of 452 photometric
* candidates in IC 4665
C Format for file interpretation
1 format(
+ 3X,A5,1X,F9.5,1X,F7.5,1X,F6.3,1X,F5.3,1X,I1,1X,F6.3,1X,F6.1,
+ 1X,F4.1,1X,F6.3,1X,F6.1,1X,F3.1,1X,F6.3,1X,F6.1,1X,F3.1,1X,
+ F6.3,1X,F5.1,1X,F3.1,1X,F7.1,1X,F4.1,1X,I3,1X,I3,1X,I3,1X,
+ F7.1,1X,F7.1,1X,I1,1X,I1,1X,I3,1X,I5)
C Effective file loading
open(unit=1,file='table1.dat', status='old')
write(6,*) '....Loading file: table1.dat'
do i__=1,452
read(1,'(A168)')ar__
read(ar__,1)
+ Name(i__),RAdeg(i__),DEdeg(i__),Vmag(i__),V_I(i__),Nobs(i__),
+ HJD1(i__),HRV1(i__),e_HRV1(i__),HJD2(i__),HRV2(i__),
+ e_HRV2(i__),HJD3(i__),HRV3(i__),e_HRV3(i__),HJD4(i__),
+ HRV4(i__),e_HRV4(i__),HRV(i__),e_HRV(i__),SNR(i__),EWLi(i__),
+ e_EWLi(i__),vsini(i__),e_vsini(i__),Var(i__),CCF(i__),
+ P93(i__),M08(i__)
if(ar__(43:48) .EQ. '') HJD1(i__) = rNULL__
if(ar__(50:55) .EQ. '') HRV1(i__) = rNULL__
if(ar__(57:60) .EQ. '') e_HRV1(i__) = rNULL__
if(ar__(62:67) .EQ. '') HJD2(i__) = rNULL__
if(ar__(69:74) .EQ. '') HRV2(i__) = rNULL__
if(ar__(76:78) .EQ. '') e_HRV2(i__) = rNULL__
if(ar__(80:85) .EQ. '') HJD3(i__) = rNULL__
if(ar__(87:92) .EQ. '') HRV3(i__) = rNULL__
if(ar__(94:96) .EQ. '') e_HRV3(i__) = rNULL__
if(ar__(98:103) .EQ. '') HJD4(i__) = rNULL__
if(ar__(105:109) .EQ. '') HRV4(i__) = rNULL__
if(ar__(111:113) .EQ. '') e_HRV4(i__) = rNULL__
if(ar__(132:134) .EQ. '') EWLi(i__) = iNULL__
if(ar__(156:156) .EQ. '') Var(i__) = iNULL__
if(ar__(158:158) .EQ. '') CCF(i__) = iNULL__
if(ar__(160:162) .EQ. '') P93(i__) = iNULL__
if(ar__(164:168) .EQ. '') M08(i__) = iNULL__
c ..............Just test output...........
write(6,1)
+ Name(i__),RAdeg(i__),DEdeg(i__),Vmag(i__),V_I(i__),Nobs(i__),
+ HJD1(i__),HRV1(i__),e_HRV1(i__),HJD2(i__),HRV2(i__),
+ e_HRV2(i__),HJD3(i__),HRV3(i__),e_HRV3(i__),HJD4(i__),
+ HRV4(i__),e_HRV4(i__),HRV(i__),e_HRV(i__),SNR(i__),EWLi(i__),
+ e_EWLi(i__),vsini(i__),e_vsini(i__),Var(i__),CCF(i__),
+ P93(i__),M08(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Properties of candidate RV members of IC 4665
C Format for file interpretation
2 format(
+ 3X,A5,1X,F5.2,1X,F4.2,1X,I3,1X,I3,1X,F5.1,1X,F3.1,1X,F4.1,1X,
+ F5.1,1X,F5.2,1X,F4.2,1X,F4.2,1X,F5.1,1X,F3.1,1X,F5.1,1X,F3.1,
+ 1X,A1)
C Effective file loading
open(unit=1,file='table2.dat', status='old')
write(6,*) '....Loading file: table2.dat'
do i__=1,56
read(1,'(A86)')ar__1
read(ar__1,2)
+ Name_1(i__),Vmag_1(i__),V_I_1(i__),EWLi_1(i__),e_EWLi_1(i__),
+ RV(i__),e_RV(i__),vsini_1(i__),e_vsini_1(i__),EWHa(i__),
+ EWCaI1(i__),EWCaI2(i__),pmRA(i__),e_pmRA(i__),pmDE(i__),
+ e_pmDE(i__),Mm(i__)
if(ar__1(66:70) .EQ. '') pmRA(i__) = rNULL__
if(ar__1(72:74) .EQ. '') e_pmRA(i__) = rNULL__
if(ar__1(76:80) .EQ. '') pmDE(i__) = rNULL__
if(ar__1(82:84) .EQ. '') e_pmDE(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Name_1(i__),Vmag_1(i__),V_I_1(i__),EWLi_1(i__),e_EWLi_1(i__),
+ RV(i__),e_RV(i__),vsini_1(i__),e_vsini_1(i__),EWHa(i__),
+ EWCaI1(i__),EWCaI2(i__),pmRA(i__),e_pmRA(i__),pmDE(i__),
+ e_pmDE(i__),Mm(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Lithium abundances in IC 4665
C Format for file interpretation
3 format(3X,A5,1X,F5.3,1X,I4,1X,A1,F5.2,1X,F4.2,1X,A1,F4.2)
C Effective file loading
open(unit=1,file='table3.dat', status='old')
write(6,*) '....Loading file: table3.dat'
do i__=1,44
read(1,'(A37)')ar__2
read(ar__2,3)
+ Name_2(i__),V_I_2(i__),Teff(i__),l_A_Li_LTE(i__),
+ A_Li_LTE(i__),e_A_Li_LTE(i__),l_A_Li_NLTE(i__),A_Li_NLTE(i__)
if(ar__2(22:26) .EQ. '') A_Li_LTE(i__) = rNULL__
if(ar__2(28:31) .EQ. '') e_A_Li_LTE(i__) = rNULL__
if(ar__2(34:37) .EQ. '') A_Li_NLTE(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ Name_2(i__),V_I_2(i__),Teff(i__),l_A_Li_LTE(i__),
+ A_Li_LTE(i__),e_A_Li_LTE(i__),l_A_Li_NLTE(i__),A_Li_NLTE(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end