Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/416/787 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/416/787 Mount Wilson Crv Metallicity index (Soon+ 1993)
*================================================================================
*The Mount Wilson Observatory Metallicity index, Crv: comparison with
*other photometric systems
* SOON W.H., ZHANG Q., BALIUNAS S.L., KURUCZ R.L.
* <Astrophys. J. 416, 787 (1993)>
* =1993ApJ...416..787S
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1' ! Data for dwarfs
integer*4 nr__
parameter (nr__=236) ! Number of records
character*127 ar__ ! Full-size record
integer*4 HD (nr__) ! []? HD designation
character*1 m_HD (nr__) ! Multiplicity index on HD designation
character*1 n_HD (nr__) ! [*] Note for HD 36389 (1)
character*10 Name (nr__) ! Other name
real*4 v_Crv_ (nr__) ! (mag) []? Time averaged spectrophotometric index
* Crv = 2.5 log (N_R/N_V) (2)
integer*4 o__Crv_ (nr__) ! []? Number of observed points
real*4 e__Crv_ (nr__) ! (mag) []? Standard deviation of observed points
real*4 b_y (nr__) ! (mag) Observed b-y index
real*4 m1 (nr__) ! (mag) Observed m1 index
real*4 c1 (nr__) ! (mag) []? Observed c1 index
real*4 V_K (nr__) ! (mag) []? Observed V-K index
real*4 THETAeff (nr__) ! []? Observed stellar effective temperature
* (THETAeff= 5040/Teff, Teff in K)
real*4 log_g (nr__) ! ([cm/s2]) []? Observed surface gravity
real*4 v_Fe_H_1 (nr__) ! []? Scaled iron abundance
real*4 v_Fe_H_2 (nr__) ! []? Highest value of metallicity when
* n_[Fe/H] = ,
real*4 Dm (nr__) ! (mag) []? magnitude difference of the two components
* for double stars
real*4 Sep (nr__) ! (arcsec) []? Maximum observed separation for
* double stars
character*20 Period (nr__) ! Orbital period for double stars
* The symbols SBa and SB identify the primary
* star associated with a known or suspected
* spectroscopic binary system.
*Note (1): Note that the reddest giant on the list HD 36389 had b-y = 1.440
* and is not plotted in all the (b-y)-dependent plots.
*Note (2): N_R and N_V are the counts in the two wider windows (20A) labeled
* as R_4001 and V_3901, centered at 4001 and 3901 A respectively.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2' ! Data for giants
integer*4 nr__1
parameter (nr__1=141) ! Number of records
character*127 ar__1 ! Full-size record
integer*4 HD_1 (nr__1) ! []? HD designation
character*1 m_HD_1 (nr__1) ! Multiplicity index on HD designation
character*1 n_HD_1 (nr__1) ! [*] Note for HD 36389 (1)
character*10 Name_1 (nr__1) ! Other name
real*4 v_Crv__1 (nr__1) ! (mag) []? Time averaged spectrophotometric index
* Crv = 2.5 log (N_R/N_V) (2)
integer*4 o__Crv__1 (nr__1) ! []? Number of observed points
real*4 e__Crv__1 (nr__1) ! (mag) []? Standard deviation of observed points
real*4 b_y_1 (nr__1) ! (mag) Observed b-y index
real*4 m1_1 (nr__1) ! (mag) Observed m1 index
real*4 c1_1 (nr__1) ! (mag) []? Observed c1 index
real*4 V_K_1 (nr__1) ! (mag) []? Observed V-K index
real*4 THETAeff_1 (nr__1) ! []? Observed stellar effective temperature
* (THETAeff= 5040/Teff, Teff in K)
real*4 log_g_1 (nr__1) ! ([cm/s2]) []? Observed surface gravity
real*4 v_Fe_H_1_1 (nr__1) ! []? Scaled iron abundance
real*4 v_Fe_H_2_1 (nr__1) ! []? Highest value of metallicity when
* n_[Fe/H] = ,
real*4 Dm_1 (nr__1) ! (mag) []? magnitude difference of the two components
* for double stars
real*4 Sep_1 (nr__1) ! (arcsec) []? Maximum observed separation for
* double stars
character*20 Period_1 (nr__1) ! Orbital period for double stars
* The symbols SBa and SB identify the primary
* star associated with a known or suspected
* spectroscopic binary system.
*Note (1): Note that the reddest giant on the list HD 36389 had b-y = 1.440
* and is not plotted in all the (b-y)-dependent plots.
*Note (2): N_R and N_V are the counts in the two wider windows (20A) labeled
* as R_4001 and V_3901, centered at 4001 and 3901 A respectively.
C=============================================================================
C Loading file 'table1' ! Data for dwarfs
C Format for file interpretation
1 format(
+ I6,A1,A1,1X,A10,2X,F6.3,1X,I4,1X,F5.3,1X,F5.3,1X,F6.3,1X,F5.3,
+ 1X,F4.2,1X,F4.2,1X,F3.1,2X,F5.2,1X,F5.2,7X,F4.1,6X,F5.1,1X,A20)
C Effective file loading
open(unit=1,status='old',file=
+'table1')
write(6,*) '....Loading file: table1'
do i__=1,236
read(1,'(A127)')ar__
read(ar__,1)
+ HD(i__),m_HD(i__),n_HD(i__),Name(i__),v_Crv_(i__),
+ o__Crv_(i__),e__Crv_(i__),b_y(i__),m1(i__),c1(i__),V_K(i__),
+ THETAeff(i__),log_g(i__),v_Fe_H_1(i__),v_Fe_H_2(i__),Dm(i__),
+ Sep(i__),Period(i__)
if(ar__(1:6) .EQ. '') HD(i__) = iNULL__
if(ar__(22:27) .EQ. '') v_Crv_(i__) = rNULL__
if(ar__(29:32) .EQ. '') o__Crv_(i__) = iNULL__
if(ar__(34:38) .EQ. '') e__Crv_(i__) = rNULL__
if(ar__(53:57) .EQ. '') c1(i__) = rNULL__
if(ar__(59:62) .EQ. '') V_K(i__) = rNULL__
if(ar__(64:67) .EQ. '') THETAeff(i__) = rNULL__
if(ar__(69:71) .EQ. '') log_g(i__) = rNULL__
if(ar__(74:78) .EQ. '') v_Fe_H_1(i__) = rNULL__
if(ar__(80:84) .EQ. '') v_Fe_H_2(i__) = rNULL__
if(ar__(92:95) .EQ. '') Dm(i__) = rNULL__
if(ar__(102:106) .EQ. '') Sep(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ HD(i__),m_HD(i__),n_HD(i__),Name(i__),v_Crv_(i__),
+ o__Crv_(i__),e__Crv_(i__),b_y(i__),m1(i__),c1(i__),V_K(i__),
+ THETAeff(i__),log_g(i__),v_Fe_H_1(i__),v_Fe_H_2(i__),Dm(i__),
+ Sep(i__),Period(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2' ! Data for giants
C Format for file interpretation
2 format(
+ I6,A1,A1,1X,A10,2X,F6.3,1X,I4,1X,F5.3,1X,F5.3,1X,F6.3,1X,F5.3,
+ 1X,F4.2,1X,F4.2,1X,F3.1,2X,F5.2,1X,F5.2,7X,F4.1,6X,F5.1,1X,A20)
C Effective file loading
open(unit=1,status='old',file=
+'table2')
write(6,*) '....Loading file: table2'
do i__=1,141
read(1,'(A127)')ar__1
read(ar__1,2)
+ HD_1(i__),m_HD_1(i__),n_HD_1(i__),Name_1(i__),v_Crv__1(i__),
+ o__Crv__1(i__),e__Crv__1(i__),b_y_1(i__),m1_1(i__),c1_1(i__),
+ V_K_1(i__),THETAeff_1(i__),log_g_1(i__),v_Fe_H_1_1(i__),
+ v_Fe_H_2_1(i__),Dm_1(i__),Sep_1(i__),Period_1(i__)
if(ar__1(1:6) .EQ. '') HD_1(i__) = iNULL__
if(ar__1(22:27) .EQ. '') v_Crv__1(i__) = rNULL__
if(ar__1(29:32) .EQ. '') o__Crv__1(i__) = iNULL__
if(ar__1(34:38) .EQ. '') e__Crv__1(i__) = rNULL__
if(ar__1(53:57) .EQ. '') c1_1(i__) = rNULL__
if(ar__1(59:62) .EQ. '') V_K_1(i__) = rNULL__
if(ar__1(64:67) .EQ. '') THETAeff_1(i__) = rNULL__
if(ar__1(69:71) .EQ. '') log_g_1(i__) = rNULL__
if(ar__1(74:78) .EQ. '') v_Fe_H_1_1(i__) = rNULL__
if(ar__1(80:84) .EQ. '') v_Fe_H_2_1(i__) = rNULL__
if(ar__1(92:95) .EQ. '') Dm_1(i__) = rNULL__
if(ar__1(102:106) .EQ. '') Sep_1(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ HD_1(i__),m_HD_1(i__),n_HD_1(i__),Name_1(i__),v_Crv__1(i__),
+ o__Crv__1(i__),e__Crv__1(i__),b_y_1(i__),m1_1(i__),c1_1(i__),
+ V_K_1(i__),THETAeff_1(i__),log_g_1(i__),v_Fe_H_1_1(i__),
+ v_Fe_H_2_1(i__),Dm_1(i__),Sep_1(i__),Period_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end