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-Jun-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__
ter (rNULL__=-1.e37) ! NULL real number
parameter (iNULL__=-2147483647) ! NULL int number
C=============================================================================
Cat. J/ApJS/152/261 Chromospheric Ca II emission in nearby stars (Wright+, 2004)
*================================================================================
*Chromospheric Ca II emission in nearby F, G, K, and M stars.
* Wright J.T., Marcy G.W., Butler R.P., Vogt S.S.
* <Astrophys. J. Suppl. Ser., 152, 261-295 (2004)>
* =2004ApJS..152..261W
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Sensitive differential S value (S_diff_) for
planet search program stars
integer*4 nr__
parameter (nr__=13808) ! Number of records
character*79 ar__ ! Full-size record
character*8 Name (nr__) ! HD or ADS number
integer*4 HIP (nr__) ! ? The Hipparcos number
character*14 Other (nr__) ! Other name
real*8 Date (nr__) ! (d) Observation date
real*4 Sdiff (nr__) ! Sensitive differential S value (1)
character*28 Remarks (nr__) ! Additional notes (2)
*Note (1): The S-values were constructed as
* S = {alpha}(H+K)/(R+V),
* where H, K, R, and V refer to the flux in the corresponding
* bandpasses and {alpha} was calculated to be 2.4 to make the mean S
* correspond to the mean F determined from the HKP-1 observations.
*Note (2): This column contains alternate names for some stars and J2000
* Hipparcos coordinates for stars for which only one catalog name
* is given.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Measured and derived quantities for all stars
integer*4 nr__1
parameter (nr__1=1205) ! Number of records
character*130 ar__1 ! Full-size record
character*8 Name_1 (nr__1) ! HD or ADS number
integer*4 HIP_1 (nr__1) ! ? The Hipparcos number
character*14 Other_1 (nr__1) ! Other name
real*4 B_V (nr__1) ! (mag) ? The (B-V) color (1)
real*4 VMAG (nr__1) ! (mag) ? Absolute V band magnitude (1)
integer*4 FObs (nr__1) ! (d) First observation date used in calculating
* the grand S value
integer*4 LObs (nr__1) ! (d) Last observation date used in calculating
* the grand S value
integer*4 NObs (nr__1) ! Total number of observations
integer*4 NBins (nr__1) ! Number of month bins used to calculate the
* grand S value as defined in 5.1.1
real*4 GSVal (nr__1) ! The grand S value (2)
real*4 Ratio (nr__1) ! (%) ? The {sigma}_S_diff_/S ratio (3)
real*4 log_R_HK (nr__1) ! ([-]) ? Log of the R'_HK_ activity level
real*4 Per (nr__1) ! (d) ? Rotation period
character*1 l_log_Age (nr__1) ! Limit flag on log(Age)
real*4 log_Age (nr__1) ! ([yr]) ? Log of the age
character*1 Flag (nr__1) ! [1] Flag on log(R'HK), PRot, & log(Age) (4)
character*32 Note (nr__1) ! Additional notes (5)
*Note (1): Values are from Hipparcos (Perryman et al., 1997, Cat. <I/239>),
* where available.
*Note (2): The overall uncertainty in the grand S varies from star to star,
* but we conservatively estimate it to be around 13%, in general
* (see Sec. 5.1).
*Note (3): Ratio is the fractional standard deviation of the differential S
* values for a star (for those stars with Keck data) and measures the
* star's intrinsic variability.
*Note (4): Flag as follows:
* 1: The chromospheric activity relation used to derive this quantity
* is only calibrated for stars with 0.44<B-V<0.9.
*Note (5): The note "d8" refers to an entry based solely on Dewar 8 data, for
* which an the uncertainty of each measurement contributing to the
* grand S is around 10%. All other entries are based on Keck and
* Dewar 6 data which have per-measurement uncertainties of around 6%.
* Alternate names for stars are noted, and J2000 Hipparcos coordinates
* are given for stars with only one catalog name.
C=============================================================================
C Loading file 'table1.dat' ! Sensitive differential S value (S_diff_) for
* planet search program stars
C Format for file interpretation
1 format(A8,1X,I6,1X,A14,1X,F13.5,1X,F5.3,1X,A28)
C Effective file loading
open(unit=1,file='table1.dat', status='old')
write(6,*) '....Loading file: table1.dat'
do i__=1,13808
read(1,'(A79)')ar__
read(ar__,1)
+ Name(i__),HIP(i__),Other(i__),Date(i__),Sdiff(i__),
+ Remarks(i__)
if(ar__(10:15) .EQ. '') HIP(i__) = iNULL__
c ..............Just test output...........
write(6,1)
+ Name(i__),HIP(i__),Other(i__),Date(i__),Sdiff(i__),
+ Remarks(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! Measured and derived quantities for all stars
C Format for file interpretation
2 format(
+ A8,1X,I6,1X,A14,1X,F5.3,1X,F5.2,1X,I7,1X,I7,1X,I3,1X,I2,1X,
+ F5.3,1X,F5.2,1X,F5.2,1X,F4.1,1X,A1,F5.2,1X,A1,1X,A32)
C Effective file loading
open(unit=1,file='table2.dat', status='old')
write(6,*) '....Loading file: table2.dat'
do i__=1,1205
read(1,'(A130)')ar__1
read(ar__1,2)
+ Name_1(i__),HIP_1(i__),Other_1(i__),B_V(i__),VMAG(i__),
+ FObs(i__),LObs(i__),NObs(i__),NBins(i__),GSVal(i__),
+ Ratio(i__),log_R_HK(i__),Per(i__),l_log_Age(i__),log_Age(i__),
+ Flag(i__),Note(i__)
if(ar__1(10:15) .EQ. '') HIP_1(i__) = iNULL__
if(ar__1(32:36) .EQ. '') B_V(i__) = rNULL__
if(ar__1(38:42) .EQ. '') VMAG(i__) = rNULL__
if(ar__1(73:77) .EQ. '') Ratio(i__) = rNULL__
if(ar__1(79:83) .EQ. '') log_R_HK(i__) = rNULL__
if(ar__1(85:88) .EQ. '') Per(i__) = rNULL__
if(ar__1(91:95) .EQ. '') log_Age(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Name_1(i__),HIP_1(i__),Other_1(i__),B_V(i__),VMAG(i__),
+ FObs(i__),LObs(i__),NObs(i__),NBins(i__),GSVal(i__),
+ Ratio(i__),log_R_HK(i__),Per(i__),l_log_Age(i__),log_Age(i__),
+ Flag(i__),Note(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end