Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/157/63 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-May-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__
parameter (rNULL__=--2147483648.) ! NULL real number
parameter (iNULL__=(-2147483647-1)) ! NULL int number
integer idig ! testing NULL number
C=============================================================================
Cat. J/AJ/157/63 Radius relations for low-metallicity M-dwarf stars (Kesseli+, 2019)
*================================================================================
*Radii of 88 M subdwarfs and updated radius relations for low-metallicity
*M-dwarf stars.
* Kesseli A.Y., Kirkpatrick J.D., Fajardo-Acosta S.B., Penny M.T.,
* Gaudi B.S., Veyette M., Boeshaar P.C., Henderson C.B., Cushing M.C.,
* Calchi-Novati S., Shvartzvald Y., Muirhead P.S.
* <Astron. J., 157, 63 (2019)>
* =2019AJ....157...63K (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Photometry
integer*4 nr__
parameter (nr__=94) ! Number of records
character*289 ar__ ! Full-size record
character*14 Star (nr__) ! Star identifier
real*8 umagSD (nr__) ! (mag) [11.0806/24.4154]? Apparent SDSS u band magnitude
real*4 e_umagSD (nr__) ! (mag) [0.0005/1.1093]? Uncertainty in umagSD
real*8 gmagSD (nr__) ! (mag) [8.6298/21.6095]? Apparent SDSS g band magnitude
real*4 e_gmagSD (nr__) ! (mag) [0.0001/1.086]? Uncertainty in gmagSD
real*4 gmagPS (nr__) ! (mag) [14.255/21.15]? Apparent Pan-STARRS g band
* magnitude
real*8 e_gmagPS (nr__) ! (mag) [0.000293/0.034]? Uncertainty in gmagPS
real*8 Bmag (nr__) ! (mag) [7.8828/20.1607]? Apparent Gaia blue band
* magnitude
real*4 e_Bmag (nr__) ! (mag) [0.0012/0.0993]? Uncertainty in Bmag
real*4 rmagPS (nr__) ! (mag) [13.861/19.275]? Apparent Pan-STARRS r band
* magnitude
real*8 e_rmagPS (nr__) ! (mag) [5.1e-05/0.015]? Uncertainty in rmagPS
real*8 rmagSD (nr__) ! (mag) [7.0918/19.2307]? Apparent SDSS r band magnitude
real*4 e_rmagSD (nr__) ! (mag) [0/2.0723]? Uncertainty in rmagSD
real*4 imagPS (nr__) ! (mag) [14.531/17.864]? Apparent Pan-STARRS i band
* magnitude
real*8 e_imagPS (nr__) ! (mag) [0.000976/0.018]? Uncertainty in imagPS
real*8 imagSD (nr__) ! (mag) [6.4159/17.8655]? Apparent SDSS i band magnitude
real*4 e_imagSD (nr__) ! (mag) [0/7.8301]? Uncertainty in imagSD
real*8 Rmag (nr__) ! (mag) [6.0594/17.4027]? Apparent Gaia red band
* magnitude
real*4 e_Rmag (nr__) ! (mag) [0.0006/0.054]? Uncertainty in Rmag
real*4 zmagPS (nr__) ! (mag) [13.534/17.237]? Apparent Pan-STARRS z band
* magnitude
real*8 e_zmagPS (nr__) ! (mag) [1.7e-05/0.02]? Uncertainty in ymagPS
real*8 zmagSD (nr__) ! (mag) [6.7336/17.1567]? Apparent SDSS z band magnitude
real*4 e_zmagSD (nr__) ! (mag) [0.0001/1.0869]? Uncertainty in zmagSD
real*4 ymagPS (nr__) ! (mag) [12.622/16.964]? Apparent Pan-STARRS y band
* magnitude
real*8 e_ymagPS (nr__) ! (mag) [9e-06/0.2]? Uncertainty in ymagPS
real*4 Jmag (nr__) ! (mag) [5.449/15.827]? Apparent 2MASS J band magnitude
real*4 e_Jmag (nr__) ! (mag) [0.018/0.079]? Uncertainty in Jmag
real*4 Hmag (nr__) ! (mag) [5.282/15.539]? Apparent 2MASS H band magnitude
real*4 e_Hmag (nr__) ! (mag) [0.016/0.117]? Uncertainty in Hmag
real*4 Ksmag (nr__) ! (mag) [4.593/15.362]? Apparent 2MASS K_s_ band
* magnitude
real*4 e_Ksmag (nr__) ! (mag) [0.014/0.18]? Uncertainty in Ksmag
real*4 W1mag (nr__) ! (mag) [4.042/14.997]? Apparent WISE 3.4 micron band
* magnitude
real*4 e_W1mag (nr__) ! (mag) [0.022/0.282]? Uncertainty in W1mag
real*4 W2mag (nr__) ! (mag) [3.318/14.754]? Apparent WISE 4.6 micron band
* magnitude
real*4 e_W2mag (nr__) ! (mag) [0.016/0.165]? Uncertainty in W2mag
real*4 W3mag (nr__) ! (mag) [2.495/13.051]? Apparent WISE 12 micron band
* magnitude
real*4 e_W3mag (nr__) ! (mag) [0.013/0.406]? Uncertainty in W3mag
real*4 W4mag (nr__) ! (mag) [2.519/8.816]? Apparent WISE 22 micron band
* magnitude
real*4 e_W4mag (nr__) ! (mag) [0.018/0.423]? Uncertainty in W4mag
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Derived parameters
integer*4 nr__1
parameter (nr__1=87) ! Number of records
character*79 ar__1 ! Full-size record
character*13 Star_1 (nr__1) ! Star identifier
character*6 SpType (nr__1) ! Spectral type (1)
integer*4 Teff (nr__1) ! (K) [2700/4300] Effective temperature
integer*4 e_Teff (nr__1) ! (K) [100] Uncertainty in Teff
real*4 logL (nr__1) ! ([Lsun]) [-3.137/-1.044] Bolometric luminosity
* (log(L_bol_/L_{sun}_))
real*4 e_logL (nr__1) ! ([Lsun]) [0.001/0.333] Uncertainty in logL
real*4 Rad (nr__1) ! (Rsun) [0.104/0.626] Radius
real*4 e_Rad (nr__1) ! (Rsun) [0.007/0.165] Uncertainty in Rad
real*4 v_Fe_H_ (nr__1) ! ([-]) [-1.88/0.25] Metallicity
real*4 e__Fe_H_ (nr__1) ! ([-]) [0.08/0.3] Uncertainty in [Fe/H]
character*11 Meth (nr__1) ! Method for determining the metallicity (2)
character*1 r_Meth (nr__1) ! [abc] Reference for Meth (3)
*Note (1): We restricted our spectral class range to ~K7 through ~M7. The
* Lepine et al. (2007ApJ...669.1235L) subdwarf subclasses - subdwarfs (sd),
* extreme subdwarfs (esd), and ultra subdwarfs (usd), roughly represent objects
* in the metallicity ranges log ([Fe/H])~-0.5, -1.0, and -1.5, respectively.
* Most known late-K through late-M subdwarfs were classified before the
* Lepine et al. system was established, some of which were typed against the
* earlier Gizis (1997, J/AJ/113/806) two-subclass system. Others pre-date both
* of these papers and are classified on a mixture of systems.
*Note (2): Method as follows:
* Lit = Previously measured literature [Fe/H] value;
* Phot = Photometrically determined [Fe/H] value;
* Spec = Spectroscopic metallicity;
* iSHELL_Spec = We used the high-resolution (R~35000) near-infrared iSHELL
* spectra.
*Note (3): Reference as follows:
* a = Mann et al. (2015, J/ApJ/804/64).
* b = Gaidos et al. (2014, J/MNRAS/443/2561).
* c = Schmidt et al. (2016MNRAS.460.2611S).
C=============================================================================
C Loading file 'table2.dat' ! Photometry
C Format for file interpretation
1 format(
+ A14,1X,F7.4,1X,F6.4,1X,F7.4,1X,F6.4,1X,F6.3,1X,F8.6,1X,F7.4,
+ 1X,F6.4,1X,F6.3,1X,F8.6,1X,F7.4,1X,F6.4,1X,F6.3,1X,F8.6,1X,
+ F7.4,1X,F6.4,1X,F7.4,1X,F6.4,1X,F6.3,1X,F8.6,1X,F7.4,1X,F6.4,
+ 1X,F6.3,1X,F8.6,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,
+ F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F5.3,
+ 1X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,94
read(1,'(A289)')ar__
read(ar__,1)
+ Star(i__),umagSD(i__),e_umagSD(i__),gmagSD(i__),e_gmagSD(i__),
+ gmagPS(i__),e_gmagPS(i__),Bmag(i__),e_Bmag(i__),rmagPS(i__),
+ e_rmagPS(i__),rmagSD(i__),e_rmagSD(i__),imagPS(i__),
+ e_imagPS(i__),imagSD(i__),e_imagSD(i__),Rmag(i__),e_Rmag(i__),
+ zmagPS(i__),e_zmagPS(i__),zmagSD(i__),e_zmagSD(i__),
+ ymagPS(i__),e_ymagPS(i__),Jmag(i__),e_Jmag(i__),Hmag(i__),
+ e_Hmag(i__),Ksmag(i__),e_Ksmag(i__),W1mag(i__),e_W1mag(i__),
+ W2mag(i__),e_W2mag(i__),W3mag(i__),e_W3mag(i__),W4mag(i__),
+ e_W4mag(i__)
if(ar__(16:22) .EQ. '') umagSD(i__) = rNULL__
if(ar__(24:29) .EQ. '') e_umagSD(i__) = rNULL__
if(ar__(31:37) .EQ. '') gmagSD(i__) = rNULL__
if(ar__(39:44) .EQ. '') e_gmagSD(i__) = rNULL__
if(ar__(46:51) .EQ. '') gmagPS(i__) = rNULL__
if(ar__(53:60) .EQ. '') e_gmagPS(i__) = rNULL__
if(ar__(62:68) .EQ. '') Bmag(i__) = rNULL__
if(ar__(70:75) .EQ. '') e_Bmag(i__) = rNULL__
if(ar__(77:82) .EQ. '') rmagPS(i__) = rNULL__
if(ar__(84:91) .EQ. '') e_rmagPS(i__) = rNULL__
if(ar__(93:99) .EQ. '') rmagSD(i__) = rNULL__
if(ar__(101:106) .EQ. '') e_rmagSD(i__) = rNULL__
if(ar__(108:113) .EQ. '') imagPS(i__) = rNULL__
if(ar__(115:122) .EQ. '') e_imagPS(i__) = rNULL__
if(ar__(124:130) .EQ. '') imagSD(i__) = rNULL__
if(ar__(132:137) .EQ. '') e_imagSD(i__) = rNULL__
if(ar__(139:145) .EQ. '') Rmag(i__) = rNULL__
if(ar__(147:152) .EQ. '') e_Rmag(i__) = rNULL__
if(ar__(154:159) .EQ. '') zmagPS(i__) = rNULL__
if(ar__(161:168) .EQ. '') e_zmagPS(i__) = rNULL__
if(ar__(170:176) .EQ. '') zmagSD(i__) = rNULL__
if(ar__(178:183) .EQ. '') e_zmagSD(i__) = rNULL__
if(ar__(185:190) .EQ. '') ymagPS(i__) = rNULL__
if(ar__(192:199) .EQ. '') e_ymagPS(i__) = rNULL__
if(ar__(201:206) .EQ. '') Jmag(i__) = rNULL__
if(ar__(208:212) .EQ. '') e_Jmag(i__) = rNULL__
if(ar__(214:219) .EQ. '') Hmag(i__) = rNULL__
if(ar__(221:225) .EQ. '') e_Hmag(i__) = rNULL__
if(ar__(227:232) .EQ. '') Ksmag(i__) = rNULL__
if(ar__(234:238) .EQ. '') e_Ksmag(i__) = rNULL__
if(ar__(240:245) .EQ. '') W1mag(i__) = rNULL__
if(ar__(247:251) .EQ. '') e_W1mag(i__) = rNULL__
if(ar__(253:258) .EQ. '') W2mag(i__) = rNULL__
if(ar__(260:264) .EQ. '') e_W2mag(i__) = rNULL__
if(ar__(266:271) .EQ. '') W3mag(i__) = rNULL__
if(ar__(273:277) .EQ. '') e_W3mag(i__) = rNULL__
if(ar__(279:283) .EQ. '') W4mag(i__) = rNULL__
if(ar__(285:289) .EQ. '') e_W4mag(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Star(i__),umagSD(i__),e_umagSD(i__),gmagSD(i__),e_gmagSD(i__),
+ gmagPS(i__),e_gmagPS(i__),Bmag(i__),e_Bmag(i__),rmagPS(i__),
+ e_rmagPS(i__),rmagSD(i__),e_rmagSD(i__),imagPS(i__),
+ e_imagPS(i__),imagSD(i__),e_imagSD(i__),Rmag(i__),e_Rmag(i__),
+ zmagPS(i__),e_zmagPS(i__),zmagSD(i__),e_zmagSD(i__),
+ ymagPS(i__),e_ymagPS(i__),Jmag(i__),e_Jmag(i__),Hmag(i__),
+ e_Hmag(i__),Ksmag(i__),e_Ksmag(i__),W1mag(i__),e_W1mag(i__),
+ W2mag(i__),e_W2mag(i__),W3mag(i__),e_W3mag(i__),W4mag(i__),
+ e_W4mag(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Derived parameters
C Format for file interpretation
2 format(
+ A13,1X,A6,1X,I4,1X,I3,1X,F6.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.2,
+ 1X,F4.2,1X,A11,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,87
read(1,'(A79)')ar__1
read(ar__1,2)
+ Star_1(i__),SpType(i__),Teff(i__),e_Teff(i__),logL(i__),
+ e_logL(i__),Rad(i__),e_Rad(i__),v_Fe_H_(i__),e__Fe_H_(i__),
+ Meth(i__),r_Meth(i__)
c ..............Just test output...........
write(6,2)
+ Star_1(i__),SpType(i__),Teff(i__),e_Teff(i__),logL(i__),
+ e_logL(i__),Rad(i__),e_Rad(i__),v_Fe_H_(i__),e__Fe_H_(i__),
+ Meth(i__),r_Meth(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end