Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/872/17 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 2026-Aug-16
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/872/17 HAZMAT. V. UV and X-ray evolution of K stars (Richey-Yowell+, 2019)
*================================================================================
*HAZMAT.
*V. The ultraviolet and X-ray evolution of K stars.
* Richey-Yowell T., Shkolnik E.L., Schneider A.C., Osby E., Barman T.,
* Meadows V.S.
* <Astrophys. J., 872, 17 (2019)>
* =2019ApJ...872...17R
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Number of K stars in each group or cluster
integer*4 nr__
parameter (nr__=10) ! Number of records
character*55 ar__ ! Full-size record
character*17 Cluster (nr__) ! Group/Cluster designation as in Table 3
character*1 u_Age (nr__) ! [~] Uncertainty flag on Age
integer*4 Age (nr__) ! (Myr) Age
integer*4 e_Age (nr__) ! (Myr) [3/50]? Negative error on Age
integer*4 E_Age_1 (nr__) ! (Myr) [3/51]? Positive error on Age
character*1 r_Age (nr__) ! Age reference (1)
integer*4 Nin (nr__) ! [3/310] Number of K stars in input sample
integer*4 Nnuv (nr__) ! [1/229] Number of K stars observed in NUV
* by GALEX
integer*4 DetNUV (nr__) ! [1/169] Number of K stars detected in NUV
* by GALEX
integer*4 Nfuv (nr__) ! [1/211] Number of K stars observed in FUV
* by GALEX
integer*4 DetFUV (nr__) ! [1/152] Number of K stars detected in FUV
* by GALEX
integer*4 Nx (nr__) ! [0/107] Number of K stars detected in X-ray
* by ROSAT
*Note (1): References as follows:
* a = Bell et al. (2015, J/MNRAS/454/593),
* b = Shkolnik et al. (2017, J/AJ/154/69),
* c = Jones et al. (2015ApJ...813...58J & 2017AAS...22913105J),
* d = Kraus & Hillenbrand (2007, J/AJ/134/2340),
* e = Perryman et al. (1997, J/A+A/331/81).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Observed flux densities, fluxes, and model
photospheric flux densities of the stars
used in this analysis
integer*4 nr__1
parameter (nr__1=455) ! Number of records
character*180 ar__1 ! Full-size record
C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
real*8 RAdeg (nr__1) ! (deg) Right Ascension J2000
real*8 DEdeg (nr__1) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*17 Cluster_1 (nr__1) ! Group/Cluster
integer*4 Age_1 (nr__1) ! (Myr) [10/5000] Age (1)
character*21 Name (nr__1) ! Simbad Name (2)
integer*4 RAh (nr__1) ! (h) Hour of Right Ascension (J2000)
integer*4 RAm (nr__1) ! (min) Minute of Right Ascension (J2000)
real*8 RAs (nr__1) ! (s) Second of Right Ascension (J2000)
character*1 DE_ (nr__1) ! Sign of the Declination (J2000)
integer*4 DEd (nr__1) ! (deg) Degree of Declination (J2000)
integer*4 DEm (nr__1) ! (arcmin) Arcminute of Declination (J2000)
real*8 DEs (nr__1) ! (arcsec) Arcsecond of Declination (J2000)
real*4 Mass (nr__1) ! (Msun) [0.5/1] Derived stellar mass
real*4 Dist (nr__1) ! (pc) [3.2/355.4] Distance
real*4 SJ (nr__1) ! (uJy) [1870000/568000000] Observed J band flux
* density
real*4 sigJ (nr__1) ! (uJy) [4760000/32400000] J band model photospheric
* flux density ({sigma}_J_)
character*1 l_SNUV (nr__1) ! Limit flag on S-NUV
real*4 SNUV (nr__1) ! (uJy) [5/3460000]? Observed GALEX NUV
* (1771-2831{AA}) flux density
real*4 e_SNUV (nr__1) ! (uJy) [1/13200]? Uncertainty in S-NUV
real*4 sigNUV (nr__1) ! (uJy) [18.7/31100]? NUV model photospheric
* flux density ({sigma}_NUV_)
character*1 l_SFUV (nr__1) ! Limit flag on S-FUV
real*4 SFUV (nr__1) ! (uJy) [2.5/165000]? Observed GALEX FUV
* (1344-1786{AA}) flux density
real*4 e_SFUV (nr__1) ! (uJy) [0.5/1850]? Uncertainty in S-FUV
real*4 sigFUV (nr__1) ! (uJy) [3.2e-5/3.1]? FUV model photospheric
* flux density ({sigma}_FUV_)
real*4 Fx (nr__1) ! (mW/m2) [1e-13/2.4e-10]? Observed X-ray flux
real*4 e_Fx (nr__1) ! (mW/m2) [1.4e-14/2.6e-11]? Uncertainty in Fx
*Note (1): Using the effective stellar temperatures in Pecaut & Mamajek
* (2013, J/ApJS/208/9) and the model isochrones from
* Baraffe+ (2015A&A...577A..42B), we determine mass estimates for the
* spectral types used in this work following Table 1. We used masses
* from 0.6 to 0.9M_{sun}_ but allowed one spectral subtype lower and higher
* to account for spectral type uncertainty. Table 1:
* ------------------------------------------------------------
* SpT TW Hya beta Pic Tuc-Hor AB Dor UMa Hyades Field
* ------------------------------------------------------------
* 10Myr 24Myr 45Myr 149Myr 300Myr 625Myr 5Gyr
* ------------------------------------------------------------
* K0 1.30 0.98 0.90 0.93 0.92 0.92 0.89
* K1 1.26 0.96 0.88 0.89 0.89 0.89 0.86
* K2 1.20 0.93 0.86 0.86 0.86 0.86 0.83
* K3 1.10 0.89 0.83 0.80 0.80 0.80 0.79
* K4 0.98 0.83 0.80 0.75 0.75 0.75 0.74
* K5 0.90 0.78 0.77 0.71 0.71 0.71 0.70
* K6 0.78 0.74 0.72 0.63 0.65 0.66 0.65
* K7 0.75 0.72 0.69 0.59 0.62 0.62 0.61
* K8 0.72 0.71 0.66 0.57 0.60 0.60 0.60
* K9 0.68 0.68 0.63 0.54 0.57 0.57 0.57
* M0 0.59 0.62 0.60 0.51 0.56 0.53 0.53
* M1 0.49 0.52 0.52 0.45 0.50 0.47 0.47
* ------------------------------------------------------------
*Note (2): Names beginning with "J" are 2MASS identifiers; otherwise SIMBAD
* names are used.
C=============================================================================
C Loading file 'table2.dat' ! Number of K stars in each group or cluster
C Format for file interpretation
1 format(
+ A17,1X,A1,I4,1X,I2,1X,I2,1X,A1,1X,I3,1X,I3,1X,I3,1X,I3,1X,I3,
+ 1X,I3)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,10
read(1,'(A55)')ar__
read(ar__,1)
+ Cluster(i__),u_Age(i__),Age(i__),e_Age(i__),E_Age_1(i__),
+ r_Age(i__),Nin(i__),Nnuv(i__),DetNUV(i__),Nfuv(i__),
+ DetFUV(i__),Nx(i__)
if(ar__(25:26) .EQ. '') e_Age(i__) = iNULL__
if(ar__(28:29) .EQ. '') E_Age_1(i__) = iNULL__
c ..............Just test output...........
write(6,1)
+ Cluster(i__),u_Age(i__),Age(i__),e_Age(i__),E_Age_1(i__),
+ r_Age(i__),Nin(i__),Nnuv(i__),DetNUV(i__),Nfuv(i__),
+ DetFUV(i__),Nx(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Observed flux densities, fluxes, and model
* photospheric flux densities of the stars
* used in this analysis
C Format for file interpretation
2 format(
+ A17,1X,I4,1X,A21,1X,I2,1X,I2,1X,F8.5,1X,A1,I2,1X,I2,1X,F7.4,
+ 1X,F4.2,1X,F6.2,1X,E8.2,1X,E8.2,1X,A1,1X,E8.2,1X,E8.2,1X,E8.2,
+ 1X,A1,1X,E8.2,1X,E8.2,1X,E8.2,1X,E8.2,1X,E8.2)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,455
read(1,'(A180)')ar__1
read(ar__1,2)
+ Cluster_1(i__),Age_1(i__),Name(i__),RAh(i__),RAm(i__),
+ RAs(i__),DE_(i__),DEd(i__),DEm(i__),DEs(i__),Mass(i__),
+ Dist(i__),SJ(i__),sigJ(i__),l_SNUV(i__),SNUV(i__),e_SNUV(i__),
+ sigNUV(i__),l_SFUV(i__),SFUV(i__),e_SFUV(i__),sigFUV(i__),
+ Fx(i__),e_Fx(i__)
if(ar__1(108:115) .EQ. '') SNUV(i__) = rNULL__
if(ar__1(117:124) .EQ. '') e_SNUV(i__) = rNULL__
if(ar__1(126:133) .EQ. '') sigNUV(i__) = rNULL__
if(ar__1(137:144) .EQ. '') SFUV(i__) = rNULL__
if(ar__1(146:153) .EQ. '') e_SFUV(i__) = rNULL__
if(ar__1(155:162) .EQ. '') sigFUV(i__) = rNULL__
if(ar__1(164:171) .EQ. '') Fx(i__) = rNULL__
if(ar__1(173:180) .EQ. '') e_Fx(i__) = rNULL__
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,2)
+ Cluster_1(i__),Age_1(i__),Name(i__),RAh(i__),RAm(i__),
+ RAs(i__),DE_(i__),DEd(i__),DEm(i__),DEs(i__),Mass(i__),
+ Dist(i__),SJ(i__),sigJ(i__),l_SNUV(i__),SNUV(i__),e_SNUV(i__),
+ sigNUV(i__),l_SFUV(i__),SFUV(i__),e_SFUV(i__),sigFUV(i__),
+ Fx(i__),e_Fx(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end