Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/702/A63 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-Oct-14
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/A+A/702/A63 NEMESIS catalogue of YSOs in Orion (Roquette+, 2025)
*================================================================================
*The NEMESIS catalogue of young stellar objects for the Orion star formation
*complex.
* Roquette J., Audard M., Hernandez D., Gezer I., Marton G., Mas C.,
* Madarasz M., Dionatos O.
* <Astron. Astrophys. 702, A63 (2025)>
* =2025A&A...702A..63R (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tabled1.dat' ! Main table (tables D1, D3, D4, D9)
integer*4 nr__
parameter (nr__=27879) ! Number of records
character*260 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 ID (nr__) ! [1/27915] Unique source identifier
real*8 RAdeg (nr__) ! (deg) Adopted Right Ascension (J2000)
real*8 DEdeg (nr__) ! (deg) Adopted Declination (J2000)
integer*4 uv (nr__) ! [0/2] Total number of sed observations
* in the UV range (0.001-0.4um)
integer*4 opt1 (nr__) ! [0/1] Total number of sed observations
* in the optical range (0.4-0.6um)
integer*4 opt2 (nr__) ! [0/1] Total number of sed observations
* in the optical range (0.6-0.8um)
integer*4 opt3 (nr__) ! [0/2] Total number of sed observations
* in the optical range (0.8-1.0um)
integer*4 nir1 (nr__) ! [0/2] Total number of sed observations
* in the NIR range (1.0-2.0um)
integer*4 nir2 (nr__) ! [0/2] Total number of sed observations
* in the NIR range (2.0-3.0um)
integer*4 nir3 (nr__) ! [0/2] Total number of sed observations
* in the NIR range (3.0-4.0um)
integer*4 mir1 (nr__) ! [0/2] Total number of sed observations
* in the MIR range (4.0-5.0um)
integer*4 mir2 (nr__) ! [0/2] Total number of sed observations
* in the MIR range (5.0-6.0um)
integer*4 mir3a (nr__) ! [0/2] Total number of sed observations
* in the MIR range (6.0-8.0um)
integer*4 mir3b (nr__) ! [0/2] Total number of sed observations
* in the MIR range (8.0-10.0um)
integer*4 mir4 (nr__) ! [0/2] Total number of sed observations
* in the MIR range (10.0-20.0um)
integer*4 mir5 (nr__) ! [0/2] Total number of sed observations
* in the MIR range (20.0-30.0um)
integer*4 fir1 (nr__) ! [0/1] Total number of sed observations
* in the FIR range (30.0-60.0um)
integer*4 fir2 (nr__) ! [0/2] Total number of sed observations
* in the FIR range (60.0-100.0um)
integer*4 fir3 (nr__) ! [0/2] Total number of sed observations
* in the FIR range (100.0-1000.0um)
character*4 vizier (nr__) ! If True, an SED file with the name
* vizier_sed_{ID}.csv is
* available in the sed/vizier subdirectory
character*4 nemesis (nr__) ! If True, an SED file with the name
* nemesis_sed_{ID}.csv is
* available in the sed/nemesis subdirectory
real*8 extgalProb (nr__) ! ? Probability of being extragalactic
* (Random forest classifier)
* (extragalactic_prob)
character*5 extgalPred (nr__) ! Predicted class of being extragalactic
* (Random forest classifier)
* (extragalactic_pred)
real*8 giantProb (nr__) ! ? Predicted class of being a giant
* (Random forest classifier) (giant_prob)
character*5 giantPred (nr__) ! Predicted class of being a giant
* (Random forest classifier) (giant_pred)
real*4 giantLabel (nr__) ! [1/1]? Sources labelled as giants
* (giant_label)
real*4 extgalLabel(nr__) ! [-1/1]? Sources labelled as extragalactic
* (extragalactic_label)
real*4 stratLabel (nr__) ! [-4/5]? Label for the stratification of the
* sample based on the origin of labels
* (strat_label) (1)
character*4 massiveStar(nr__) ! Sources identified as massive stars
* (massive_star)
character*4 MScontam (nr__) ! Sources identified as likely main-sequence
* contaminants (ms_contaminant)
real*4 alpha2_24 (nr__) ! ? alpha-IR index estimated using all
* available data in the wavelength
* range 2-24 micron (ms_contaminant)
character*16 Class2_24 (nr__) ! Disk Classification based on alpha_2_24
* (class_2_24) (2)
real*8 lambdaMax2_24(nr__) ! (um) ? Maximum wavelength available for
* alpha_2_24 calculation (lambda_max_2_24)
real*8 lambdaMin2_24(nr__) ! (um) ? Minimum wavelength available for
* alpha_2_24 calculation (lambda_min_2_24)
real*8 AGproxy (nr__) ! ? AG proxy for variability measurement
* estimated as equation (2) in Mowlavi
* et al., 2021A&A...648A..44M
*Note (1): Labels (in append E) as follows:
* 1 = GDR3 Galaxy high-purity sample
* 2 = GDR3 QSO high-purity sample
* 3 = MillionQua catalogue
* 4 = MANGROVE catalogue
* 5 = HyperLEDA database
* 6 = GalaxyZoo catalogue
* -1 = further stellar sample
* -2 = giant sample
* -3 = at Orion (distance) sample
* -4 = YSO sample with spectroscopic support
*Note (2): Disk Classification as:
* 0/I : alpha2-24 > 0.3
* flat spectrum : 0.3 >= alpha2-24 > -0.3
* II : -0.3 >= alpha2-24 > -1.6
* III thin disk : -1.6 >= alpha2-24 > -2.5
* III no disk : -2.5 >= alpha2-24
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tabled2a.dat' ! Spectral types table
integer*4 nr__1
parameter (nr__1=22767) ! Number of records
character*79 ar__1 ! Full-size record
integer*4 spID (nr__1) ! [0/22766] Primary key for SpT_DB (spt_ID)
integer*4 ID_1 (nr__1) ! [2/27878] Unique source identifier
* (Internal_ID)
character*16 SpType (nr__1) ! MK spectral type (1)
character*1 u_SpType (nr__1) ! [:] Uncertainty flag on SpType
character*4 e_SpType (nr__1) ! Uncertainty on MK Spectral type
character*4 photSpType (nr__1) ! [True ] Reference for MK Spectral type
character*19 Refsp (nr__1) ! Reference for MK Spectral type
character*16 n_SpType (nr__1) ! Additional Spectral information reported
*Note (1): Spectral types may be reported as a range with two MK SpTs separated
* by a symbol '-'. Unresolved multiple systems may have each components'
* spectral types reported separated by a symbol '+'.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tabled2b.dat' ! Temperatures table
integer*4 nr__2
parameter (nr__2=74508) ! Number of records
character*122 ar__2 ! Full-size record
integer*4 TempID (nr__2) ! [0/74507] Primary key for temperature_DB
* (temperature_ID)
integer*4 ID_2 (nr__2) ! [1/27913] Unique source identifier
* (Internal_ID)
character*19 RefTemp (nr__2) ! Reference for temperature measurement
real*8 Teff (nr__2) ! (K) ? Effective temperature
real*8 e_Teff (nr__2) ! (K) ? Uncertainty on effective temperature or
* lower uncertainty if E_Teff is also reported
character*8 n_Teff (nr__2) ! Note on Teff (1)
real*8 E_Teff_1 (nr__2) ! (K) ? Upper uncertainty on effective temperature
real*4 Tbol (nr__2) ! (K) ? Bolometric temperature
real*4 e_Tbol (nr__2) ! (K) ? Uncertainty on bolometric temperature
*Note (1): Note as follows:
* SpT-PHOT = combined spectral typing and photometric data
* SED = based on SED model-fit
* SPEC = spectra-library fit to spectra
* PHOT = fit of photometric data to models using MCMC or Bayesian methods
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tabled4.dat' ! Disk literature table
integer*4 nr__3
parameter (nr__3=65647) ! Number of records
character*110 ar__3 ! Full-size record
integer*4 DiskID (nr__3) ! [0/65646] Primary key for disk_DB (disk_ID)
integer*4 ID_3 (nr__3) ! [1/27471] Unique source identifier
* (Internal_ID)
character*19 RefDisk (nr__3) ! Reference for disk measurement
character*10 Class (nr__3) ! Disk IR class compiled from the literature
character*15 Disk (nr__3) ! Disk description keyword
character*9 n_Class (nr__3) ! Method used to classify the disk (1)
character*4 IRext (nr__3) ! Furhter infrared information available
real*8 alpha (nr__3) ! ? alpha IR index
real*4 e_alpha (nr__3) ! ? Uncertainty on alpha IR index
character*11 IRrange (nr__3) ! Wavelengths range used to compute the
* alpha IR index (2)
*Note (1): Method as follows:
* standard = standard classification scheme (Lada, 1987IAUS..115....1L,
* Andre et al., 1993ApJ...406..122A)
* Hernandez = Hernandez et al., 2007ApJ...662.1067H, Cat. J/ApJ/662/1067
* MGM = classification as in MEGEATH et al., 2012AJ....144..192M,
* Cat. J/AJ/144/192
* TD = classification schemes focused on Transition discs
* Robitaille = evolution classes derived with Robitaille et al.
* (2006ApJS..167..256R) models
* misc = miscellaneous classification schemes
*Note (2): >2 indicates a generic IR-range with available photometry from 2um,
* and 3-8 indicates that photometry in the range 3-8um was used
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tabled5.dat' ! Lithum data table
integer*4 nr__4
parameter (nr__4=9185) ! Number of records
character*153 ar__4 ! Full-size record
integer*4 LiID (nr__4) ! [0/9184] Primary key for Li_DB
integer*4 ID_4 (nr__4) ! [1/27887] Unique source identifier (Internal_ID)
character*19 RefLi (nr__4) ! Reference for Li measurement
real*4 EWLi (nr__4) ! (0.1nm) ? Equivalent width at {lambda}6708{AA}
real*8 e_EWLi (nr__4) ! (0.1nm) ? Uncertainty on Equivalent width
* at {lambda}6708{AA}
character*3 f_EWLi (nr__4) ! Flag for EW measurement as reported
* in source paper
real*8 E_EWLi_1 (nr__4) ! (0.1nm) ? Upper uncertainty for Equivalent width
* at {lambda}6708{AA}
real*8 ALi (nr__4) ! ? Li abundance
real*8 s_ALi (nr__4) ! ? Uncertainty on Li abundance
character*3 n_EWLi (nr__4) ! Limit flag for Equivalent width
* at {lambda}6708{AA}
real*8 e_ALi (nr__4) ! ? Lower uncertainty for Li abundance
real*8 E_ALi_1 (nr__4) ! ? Upper uncertainty for Li abundance
real*4 l_ALi (nr__4) ! [0/1]? Limit flag for Li abundance
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tabled6a.dat' ! Equivalent widths, gravity, table
integer*4 nr__5
parameter (nr__5=1880) ! Number of records
character*232 ar__5 ! Full-size record
integer*4 EWID (nr__5) ! [0/1879] Primary key for EW_gravity_DB
* (ewg_ID)
integer*4 ID_5 (nr__5) ! [6/16317] Unique source identifier
* (Internal_ID)
real*8 EWNaI8200 (nr__5) ! (0.1nm) ? Equivalent width at {lambda}8200{AA}
* (unresolved doublet)
real*8 e_EWNaI8200(nr__5) ! (0.1nm) ? Uncertainty on Equivalent width
* at {lambda}8200{AA} (unresolved doublet)
character*2 f_EWNaI8200(nr__5) ! [OO ] Comment from original source on
* Equivalent width at {lambda}8200{AA}
character*19 RefEW (nr__5) ! Reference for EW measurement
real*4 EWNaI8195 (nr__5) ! (0.1nm) ? Equivalent width at {lambda}8195{AA}
real*4 e_EWNaI8195(nr__5) ! (0.1nm) ? Uncertainty on Equivalent width
* at {lambda}8195{AA}
character*1 f_EWNaI8195(nr__5) ! [Y/N] Comment from original source on
* Equivalent width at {lambda}8195{AA}
real*4 EWCaH3_6965(nr__5) ! (0.1nm) ? Equivalent width at {lambda}6965{AA}
real*4 e_EWCaH3_6965(nr__5) ! (0.1nm) ? Uncertainty on Equivalent width
* at {lambda}6965{AA}
real*4 EWKI7665_7699(nr__5) ! (0.1nm) ? Equivalent width at
* {lambda}7665-{lambda}7699{AA}
real*4 e_EWKI7665_7699(nr__5) ! (0.1nm) ? Uncertainty on Equivalent width
* at {lambda}7665-{lambda}7699{AA}
real*8 EWKI7665 (nr__5) ! (0.1nm) ? Equivalent width at {lambda}7665{AA}
real*8 e_EWKI7665 (nr__5) ! (0.1nm) ? Uncertainty on Equivalent width
* at {lambda}7665{AA}
real*8 EWKI7699 (nr__5) ! (0.1nm) ? Equivalent width at {lambda}7699{AA}
real*8 e_EWKI7699 (nr__5) ! (0.1nm) ? Uncertainty on Equivalent width
* at {lambda}7699{AA}
real*8 EWTiO8442 (nr__5) ! (0.1nm) ? Equivalent width at {lambda}8442{AA}
real*8 e_EWTiO8442(nr__5) ! (0.1nm) ? Uncertainty on Equivalent width
* at {lambda}8442{AA}
real*4 EWTiO57130 (nr__5) ! (0.1nm) ? Equivalent width at {lambda}7130{AA}
real*4 e_EWTiO57130(nr__5) ! (0.1nm) ? Uncertainty on Equivalent width
* at {lambda}7130{AA}
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tabled6b.dat' ! Surface gravity from spectra fit table
integer*4 nr__6
parameter (nr__6=57649) ! Number of records
character*61 ar__6 ! Full-size record
integer*4 logID (nr__6) ! [0/57648] Primary key for logg table
* (logg_ID)
integer*4 ID_6 (nr__6) ! [1/27913] Unique source identifier
* (Internal_ID)
real*8 logg (nr__6) ! ([cm/s2]) Surface gravity derived from spectra fit
real*8 e_logg (nr__6) ! ([cm/s2]) ? Uncertainty for logg derived from spectra
* fit or lower uncertainty if E_logg is also
* reported
character*19 Reflogg (nr__6) ! Reference for logg measurement
real*4 E_logg_1 (nr__6) ! ([cm/s2]) ? Upper uncertainty for logg derived
* from spectra fit
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tabled7a.dat' ! Equivalent width accretion, other, table
integer*4 nr__7
parameter (nr__7=5486) ! Number of records
character*469 ar__7 ! Full-size record
integer*4 accID (nr__7) ! [0/5485] Primary key for
* EWaccretion_other_DB (acc_ID)
integer*4 ID_7 (nr__7) ! [4/27887] Unique source identifier
* (Internal_ID)
character*19 RefEW_1 (nr__7) ! Reference for EW measurement
real*8 EWHeI6678 (nr__7) ! (0.1nm) ? Equivalent width for HeI line
* at {lambda}6678{AA}
real*4 e_EWHeI6678(nr__7) ! (0.1nm) ? Uncertainty on Equivalent width for HeI
* line at {lambda}6678{AA}
character*1 f_EWHeI6678(nr__7) ! Flag for Equivalent width for HeI line
* at {lambda}6678{AA}
real*4 EWHeI5876 (nr__7) ! (0.1nm) ? Equivalent width for HeI line
* at {lambda}5876{AA}
real*4 e_EWHeI5876(nr__7) ! (0.1nm) ? Lower uncertainty for Equivalent width
* for HeI line at {lambda}5876{AA} or lower
* uncertainty if e_EWHeI5876 is also
* reported
real*4 E_EWHeI5876_1(nr__7) ! (0.1nm) ? Upper uncertainty for Equivalent width
* for HeI line at {lambda}5876{AA}
real*8 EWHeI10830 (nr__7) ! (0.1nm) ? Equivalent width for HeI line
* at {lambda}10830{AA}
real*4 e_EWHeI10830(nr__7) ! (0.1nm) []? Uncertainty on Equivalent width
* for HeI line at {lambda}10830{AA}
real*4 EWCaII8498 (nr__7) ! (0.1nm) ? Equivalent width for CaII line triplet
* at {lambda}8498{AA}
real*4 e_EWCaII8498(nr__7) ! (0.1nm) ? Uncertainty on Equivalent width for
* CaII triplet at {lambda}8498{AA}
real*8 EWCaII8542 (nr__7) ! (0.1nm) ? Equivalent width for CaII triplet
* at {lambda}8542{AA}
real*4 e_EWCaII8542(nr__7) ! (0.1nm) ? Uncertainty on Equivalent width for
* CaII triplet at {lambda}8542{AA}
character*1 f_EWCaII8542(nr__7) ! Flag for Equivalent width for CaII triplet
* at {lambda}8542{AA}
real*8 EWCaII8662 (nr__7) ! (0.1nm) ? Equivalent width for CaII triplet
* at {lambda}8662{AA}
real*4 e_EWCaII8662(nr__7) ! (0.1nm) ? Uncertainty on Equivalent width for
* CaII triplet at {lambda}8662{AA}
character*1 f_EWCaII8662(nr__7) ! Flag for Equivalent width for CaII triplet
* at {lambda}8662{AA}
real*4 EWCaH3934 (nr__7) ! (0.1nm) ? Equivalent width for CaII H&K doublet
* at {lambda}3934{AA}
real*4 e_EWCaH3934(nr__7) ! (0.1nm) ? Uncertainty on Equivalent width for
* CaII H&K doublet at {lambda}3934{AA} or
* lower uncertainty if e_EWCaH3934 is also
* reported
character*1 f_EWCaH3934(nr__7) ! Flag for Equivalent width for CaII H&K
* doublet at {lambda}3934{AA}
real*4 E_EWCaH3934_1(nr__7) ! (0.1nm) ? Upper uncertainty for Equivalent width
* for CaII H&K doublet at {lambda}3934{AA}
real*4 EWCaK3969 (nr__7) ! (0.1nm) ? Equivalent width for CaII H&K doublet
* at {lambda}3969{AA}
real*4 e_EWCaK3969(nr__7) ! (0.1nm) ? Uncertainty on Equivalent width for
* CaII H&K doublet at {lambda}3969{AA} or
* lower uncertainty if e_EWCaK3969 is also
* reported
character*1 f_EWCaK3969(nr__7) ! Flag for Equivalent width for CaII H&K
* doublet at {lambda}3969{AA}
real*4 E_EWCaK3969_1(nr__7) ! (0.1nm) ? Upper uncertainty for Equivalent width
* for CaII H&K doublet at {lambda}3969{AA}
real*8 EWOI5577 (nr__7) ! (0.1nm) ? Equivalent width for [OI] line
* at {lambda}5577{AA}
real*4 e_EWOI5577 (nr__7) ! (0.1nm) ? Uncertainty on Equivalent width for [OI]
* line at {lambda}5577{AA}
real*8 EWOI6300 (nr__7) ! (0.1nm) ? Equivalent width for [OI] line
* at {lambda}6300{AA}
real*4 e_EWOI6300 (nr__7) ! (0.1nm) ? Uncertainty on Equivalent width for [OI]
* line at {lambda}6300{AA}
real*8 EWOI6364 (nr__7) ! (0.1nm) ? Equivalent width for [OI] line
* at {lambda}6364{AA}
real*4 e_EWOI6364 (nr__7) ! (0.1nm) ? Uncertainty on Equivalent width for [OI]
* line at {lambda}6364{AA}
real*4 EWOI7773 (nr__7) ! (0.1nm) ? Equivalent width for [OI] line
* at {lambda}7773{AA}
real*4 EWOI8446 (nr__7) ! (0.1nm) ? Equivalent width for [OI] line
* at {lambda}8446{AA}
real*8 EWSII6716 (nr__7) ! (0.1nm) ? Equivalent width for [SII] line
* at {lambda}6716{AA}
real*4 e_EWSII6716(nr__7) ! (0.1nm) ? Uncertainty on Equivalent width for [SII]
* line at {lambda}6716{AA}
character*1 f_EWSII6716(nr__7) ! Flag for Equivalent width for [SII] line
* at {lambda}6716{AA}
real*8 EWSII6731 (nr__7) ! (0.1nm) ? Equivalent width for [SII] line
* at {lambda}6731{AA}
real*4 e_EWSII6731(nr__7) ! (0.1nm) ? Uncertainty on Equivalent width for [SII]
* line at {lambda}6731{AA}
character*1 f_EWSII6731(nr__7) ! Flag for Equivalent width for [SII] line
* at {lambda}6731{AA}
real*4 EWFeII16440(nr__7) ! (0.1nm) ? Equivalent width for [FeII] line
* at {lambda}16440{AA}
real*4 e_EWFeII16440(nr__7) ! (0.1nm) ? Uncertainty on Equivalent width for
* [FeII] line at {lambda}16440{AA}
real*4 EWNII6548 (nr__7) ! (0.1nm) ? Equivalent width for [NII] line
* at {lambda}6548{AA}
real*4 e_EWNII6548(nr__7) ! (0.1nm) ? Uncertainty on Equivalent width for
* [NII] line at {lambda}6548{AA}
real*4 EWNII6581 (nr__7) ! (0.1nm) ? Equivalent width for [NII] line
* at {lambda}6581{AA}
real*4 e_EWNII6581(nr__7) ! (0.1nm) ? Uncertainty on Equivalent width
* for [NII] line at {lambda}6581{AA}
real*8 EWNII6583 (nr__7) ! (0.1nm) ? Equivalent width for [NII] line
* at {lambda}6583{AA}
character*1 f_EWNII6583(nr__7) ! Flag for Equivalent width for [NII] line
* at {lambda}6583{AA}
real*8 EWPab1282 (nr__7) ! (0.1nm) ? Equivalent width for Pab line
* at {lambda}1282{AA}
real*4 e_EWPab1282(nr__7) ! (0.1nm) ? Uncertainty on Equivalent width for
* Pab line at {lambda}1282{AA}
real*4 EWPag1090 (nr__7) ! (0.1nm) ? Equivalent width for Pag line
* at {lambda}1090{AA}
real*4 e_EWPag1090(nr__7) ! (0.1nm) ? Uncertainty on Equivalent width for
* Pag line at {lambda}1090{AA}
real*4 EWBrg21661 (nr__7) ! (0.1nm) ? Equivalent width for Brg line
* at {lambda}21661{AA}
real*4 e_EWBrg21661(nr__7) ! (0.1nm) ? Uncertainty on Equivalent width for
* Brg line at {lambda}21661{AA}
real*4 EWBr11_16810(nr__7) ! (0.1nm) ? Equivalent width for Br11 line
* at {lambda}16810{AA}
real*4 e_EWBr11_16810(nr__7) ! (0.1nm) ? Uncertainty on Equivalent width for
* Br11 line at {lambda}16810{AA}
real*4 EWHb4861 (nr__7) ! (0.1nm) ? Equivalent width for Hb line
* at {lambda}4861{AA}
real*4 e_EWHb4861 (nr__7) ! (0.1nm) ? Uncertainty on Equivalent width for
* Hb line at {lambda}4861{AA} or lower
* uncertainty if e_EWHb4861 is also reported
character*1 f_EWHb4861 (nr__7) ! Flag for Equivalent width for Hb line
* at {lambda}4861{AA}
real*4 E_EWHb4861_1(nr__7) ! (0.1nm) ? Upper uncertainty for Equivalent width
* for Hb line at {lambda}4861{AA}
real*4 EWHg4340 (nr__7) ! (0.1nm) ? Equivalent width for Hg line
* at {lambda}4340{AA}
real*4 e_EWHg4340 (nr__7) ! (0.1nm) ? Uncertainty on Equivalent width for
* Hg line at {lambda}4340{AA} or lower
* uncertainty if e_EWHg4340 is also reported
character*1 f_EWHg4340 (nr__7) ! Flag for Equivalent width for Hg line
* at {lambda}4340{AA}
real*4 E_EWHg4340_1(nr__7) ! (0.1nm) ? Upper uncertainty for Equivalent width
* for Hg line at {lambda}4340{AA}
real*4 EWH113771 (nr__7) ! (0.1nm) ? Equivalent width for H11 line
* at {lambda}3771{AA}
real*4 e_EWH113771(nr__7) ! (0.1nm) ? Uncertainty on Equivalent width for
* H11 line at {lambda}3771{AA}
real*8 EWH2_1_0S_1(nr__7) ! (0.1nm) ? Equivalent width for H2 band
* at {lambda}2.118um
real*4 e_EWH2_1_0S_1(nr__7) ! (0.1nm) ? Uncertainty on Equivalent width
* for H2 band at {lambda}2.118um
real*4 EWCO_2_0 (nr__7) ! (0.1nm) ? Equivalent width for CO band
* at {lambda}2.293um
real*4 e_EWCO_2_0 (nr__7) ! (0.1nm) ? Uncertainty on Equivalent width for
* CO band at {lambda}2.293um
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tabled7b.dat' ! H{alpha} data table
integer*4 nr__8
parameter (nr__8=37449) ! Number of records
character*122 ar__8 ! Full-size record
integer*4 HaID (nr__8) ! [0/37448] Primary key for Ha_DB (ha_ID)
integer*4 ID_8 (nr__8) ! [1/27915] Unique source identifier
* (Internal_ID)
character*19 RefEWHa (nr__8) ! Reference for H{alpha} measurement
real*4 EWHa (nr__8) ! (0.1nm) ? Equivalent width for H{alpha} line
* at {lambda}6563{AA}
real*4 e_EWHa (nr__8) ! (0.1nm) []? Uncertainty on Equivalent width for
* H{alpha} line at {lambda}6563{AA}
character*1 f_EWHa (nr__8) ! Flag for Equivalent width for H{alpha} line
* at {lambda}6563{AA}
real*8 EWHac (nr__8) ! (0.1nm) ? Equivalent width for H{alpha} line
* at {lambda}6563{AA} with possible background
* contamination
real*4 e_EWHac (nr__8) ! (0.1nm) ? Uncertainty on Equivalent width for
* H{alpha} line at {lambda}6563{AA} with
* possible background contamination
character*3 Halineflag (nr__8) ! Flag for H{alpha} line measurement (1)
real*8 W10Ha (nr__8) ! (km/s) ? Full-width at 10% of the H{alpha} profile
* peak at {lambda}6563{AA}
real*8 e_W10Ha (nr__8) ! (km/s) ? Uncertainty on full-width at 10% of the
* H{alpha} profile peak at {lambda}6563{AA}
real*8 EWHaPhot (nr__8) ! (0.1nm) ? Equivalent width for H{alpha} line at
* {lambda}6563{AA} from photometry
*Note (1): Qualitative flags describing the H_alpha line profile as reported by
* Pettersson et al., 2014A&A...570A..30P, Cat. J/A+A/570/A30.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tabled7c.dat' ! Veiling table
integer*4 nr__9
parameter (nr__9=34358) ! Number of records
character*83 ar__9 ! Full-size record
integer*4 veilingID (nr__9) ! [0/34357] Primary key for veiling_DB
* (veiling_ID)
integer*4 ID_9 (nr__9) ! [1/27507] Unique source identifier
* (Internal_ID)
character*19 Refveiling (nr__9) ! Reference for veiling measurement
real*4 veiling_opt(nr__9) ! ? Veiling at optical wavelengths
* (6390-6710{AA})
real*4 e_veiling_opt(nr__9) ! ? Uncertainty on veiling at optical
* wavelengths (6390-6710{AA})
real*4 veiling_K (nr__9) ! ? Veiling at K band (~2.2um)
real*4 E_veiling_K(nr__9) ! ? Upper uncertainty on veiling
* at K band (~2.2um)
real*4 e_veiling_K_1(nr__9) ! ? Lower uncertainty on veiling
* at K band (~2.2um)
real*4 veiling_H (nr__9) ! ? Veiling at H band (~1.5um)
real*8 e_veiling_H(nr__9) ! ? Uncertainty on veiling
* at H band (~1.5 micron)
real*4 veiling_7465(nr__9) ! ? Veiling around 7465{AA}
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tabled8.dat' ! Radial velocities table
integer*4 nr__10
parameter (nr__10=60294) ! Number of records
character*81 ar__10 ! Full-size record
integer*4 rvID (nr__10) ! [0/60293] Primary key for RV_DB (rv_ID)
integer*4 ID_10 (nr__10) ! [1/27913] Unique source identifier
* (Internal_ID)
real*4 RV (nr__10) ! (km/s) Radial velocity measurement
real*8 e_RV (nr__10) ! (km/s) ? Uncertainty on radial velocity measurement
character*19 RefRV (nr__10) ! Reference for radial velocity measurement
character*2 f_RV (nr__10) ! Flag for radial velocity measurement
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tabled10.dat' ! Rotation table
integer*4 nr__11
parameter (nr__11=42132) ! Number of records
character*237 ar__11 ! Full-size record
integer*4 rotID (nr__11) ! [0/42132] Primary key for rotation_DB
* (rot_ID)
integer*4 ID_11 (nr__11) ! [1/27913] Unique source identifier
* (Internal_ID)
real*8 Per (nr__11) ! (d) ? Variability period
real*4 Per2 (nr__11) ! (d) ? Other variability period
real*4 Per3 (nr__11) ! (d) ? Other variability period
real*4 Per4 (nr__11) ! (d) ? Other variability period
real*4 Per5 (nr__11) ! (d) ? Other variability period
real*4 Per6 (nr__11) ! (d) ? Other variability period
real*4 Per7 (nr__11) ! (d) ? Other variability period
real*4 e_Per (nr__11) ! (d) ? Uncertainty on variability period or
* lower uncertainty if E_Per is also reported
character*10 n_Per (nr__11) ! Photometric band used for variability
* period measurement
real*8 vsini (nr__11) ! (km/s) ? vsini measurement
character*4 f_vsini (nr__11) ! Flag for vsini measurement
real*8 e_vsini (nr__11) ! (km/s) ? Uncertainty on vsini measurement or lower
* uncertainty if E_vsini is also reported
character*21 Refrot (nr__11) ! Reference for rotation measurement
real*8 RotSep (nr__11) ! (arcsec) ? Separation between source measured and
* source in this catalogue
character*22 E_Per_1 (nr__11) ! (d) Upper uncertainty on variability period
real*8 E_vsini_1 (nr__11) ! (km/s) ? Upper uncertainty on vsini measurement
character*1 f_Per (nr__11) ! Flag for variability period measurement
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tabled11.dat' ! X-ray table
integer*4 nr__12
parameter (nr__12=7926) ! Number of records
character*420 ar__12 ! Full-size record
integer*4 xrayID (nr__12) ! [0/7925] Primary key for Xray_DB (xray_ID)
integer*4 ID_12 (nr__12) ! [11/27915] Unique source identifier
* (Internal_ID)
character*19 Refxray (nr__12) ! Reference for Xray measurement
real*4 FX (nr__12) ! (mW/m2) ? X-ray flux
character*11 FXband (nr__12) ! X-ray band for Fx measurement (1)
real*4 e_FX (nr__12) ! (mW/m2) ? Uncertainty on X-ray flux
character*1 l_FX (nr__12) ! limit flag for X-ray flux
real*8 logLX (nr__12) ! (0.1uW) ? Total X-ray luminosity
real*4 e_logLX (nr__12) ! (0.1uW) ? Uncertainty on total X-ray luminosity
character*11 logLXband (nr__12) ! X-ray band for logLx measurement (2)
character*1 l_logLX (nr__12) ! limit flag for total X-ray luminosity
real*8 logNH (nr__12) ! ([cm-2]) ? Hydrogen column density
real*8 e_logNH (nr__12) ! ([cm-2]) ? Uncertainty on hydrogen column density
real*8 E_logNH_1 (nr__12) ! ([cm-2]) ? upper uncertainty on hydrogen column
* density
character*1 l_logNH (nr__12) ! Limit flag for hydrogen column density
real*8 HR1 (nr__12) ! ? Hardness ratio 1
real*4 e_HR1 (nr__12) ! ? Uncertainty on hardness ratio 1
real*4 E_HR1_1 (nr__12) ! ? Upper uncertainty on hardness ratio 1
character*55 HR1band (nr__12) ! Hardness ratio 1 (3)
character*2 l_HR1 (nr__12) ! [< ] Limit flag for hardness ratio 1
real*4 HR2 (nr__12) ! ? Hardness ratio 2
real*4 e_HR2 (nr__12) ! ? Uncertainty on hardness ratio 2
real*4 E_HR2_1 (nr__12) ! ? Upper uncertainty on hardness ratio
character*55 HR2band (nr__12) ! Hardness ratio 2 (4)
real*4 HR3 (nr__12) ! ? Hardness ratio 3
real*4 e_HR3 (nr__12) ! ? Lower uncertainty on hardness ratio 3
real*4 E_HR3_1 (nr__12) ! ? Upper uncertainty on hardness ratio 3
character*18 HR3band (nr__12) ! Hardness ratio 3 (5)
character*10 Survey (nr__12) ! Survey name for X-ray measurement (6)
real*4 sepX (nr__12) ! (arcsec) ? Separation between X-ray and optical
* position
*Note (1): X-ray band possible values: 0.3-10keV, 0.3-7.5keV, 0.3-8.0keV,
* 0.3-8keV, 0.5-19.5keV, 0.5-7keV and 0.5-8keV.
*Note (2): X-ray band for logLx measurement possible values: 0.1-4keV,
* 0.2-10keV, 0.2-12keV, 0.3-10keV, 0.3-7.6keV, 0.3-7.8keV, 0.5-19.5keV,
* 0.5-7.3keV, 0.5-7keV and 0.5-8keV.
*Note (3): HR1 band possible values: (1.0-2.0)-(0.5-1.0)/(0.5-2.0)keV,
* (2.0-4.5 keV)-(0.3-2.0 keV)/(2.0-4.5 keV)+(0.3-2.0 keV) and
* 0.5-2.0/2.0-8.0keV.
*Note (4): HR2 band possible values: (1.0-7.3)-(0.5-1.0)/(0.5-7.3)keV,
* (4.5-7.5 keV)-(2.0-4.5 keV)/(4.5-7.5 keV)+(2.0-4.5 keV) and
* 0.5-1.7/1.7-2.8keV.
*Note (5): HR3 band possible values: 1.7-2.8/2.8-8.0keV.
*Note (6): Surveys are COUP, Chandra, Einstein, ROSAT and XMM-Newton.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tabled12.dat' ! Multiplicity table
integer*4 nr__13
parameter (nr__13=18929) ! Number of records
character*158 ar__13 ! Full-size record
integer*4 binID (nr__13) ! [0/18928] Primary key for multiplicity_DB
integer*4 ID_13 (nr__13) ! [3/27913] Unique source identifier
* (Internal_ID)
character*2 binType (nr__13) ! Multiplicity flag (1)
character*120 binExtType (nr__13) ! Multiplicity type in source paper
character*20 Refbin (nr__13) ! Reference for multiplicity measurement
*Note (1): Multiplicity flag as follows:
* S = Binary identified from spectroscopy
* S? = Binary candidate identified from spectroscopy
* V = Visual pair
* V? = Visual pair candidate
* B = Generic binary reported from the literature
* B? = Generic binary candidate reported from the literature
* E = Eclipsing binary
* E? = Eclipsing binary candidate
* U = Unresolved pair
* U? = Unresolved pair candidate
* R = RUWE binary candidate
* A = Astrometric binary
* Bl = Likely blended source
C=============================================================================
C Loading file 'tabled1.dat' ! Main table (tables D1, D3, D4, D9)
C Format for file interpretation
1 format(
+ I5,1X,F17.14,1X,F20.16,1X,I1,1X,I1,1X,I1,1X,I1,1X,I1,1X,I1,1X,
+ I1,1X,I1,1X,I1,1X,I1,1X,I1,1X,I1,1X,I1,1X,I1,1X,I1,1X,I1,1X,
+ A4,1X,A4,1X,F18.16,1X,A5,1X,F18.16,1X,A5,1X,F3.1,1X,F4.1,1X,
+ F4.1,1X,A4,1X,A4,1X,E23.17,1X,A16,1X,F19.16,1X,F19.16,1X,
+ F18.16)
C Effective file loading
open(unit=1,status='old',file=
+'tabled1.dat')
write(6,*) '....Loading file: tabled1.dat'
do i__=1,27879
read(1,'(A260)')ar__
read(ar__,1)
+ ID(i__),RAdeg(i__),DEdeg(i__),uv(i__),opt1(i__),opt2(i__),
+ opt3(i__),nir1(i__),nir2(i__),nir3(i__),mir1(i__),mir2(i__),
+ mir3a(i__),mir3b(i__),mir4(i__),mir5(i__),fir1(i__),fir2(i__),
+ fir3(i__),vizier(i__),nemesis(i__),extgalProb(i__),
+ extgalPred(i__),giantProb(i__),giantPred(i__),giantLabel(i__),
+ extgalLabel(i__),stratLabel(i__),massiveStar(i__),
+ MScontam(i__),alpha2_24(i__),Class2_24(i__),
+ lambdaMax2_24(i__),lambdaMin2_24(i__),AGproxy(i__)
if(ar__(88:105) .EQ. '') extgalProb(i__) = rNULL__
if(ar__(113:130) .EQ. '') giantProb(i__) = rNULL__
if(ar__(138:140) .EQ. '') giantLabel(i__) = rNULL__
if(ar__(142:145) .EQ. '') extgalLabel(i__) = rNULL__
if(ar__(147:150) .EQ. '') stratLabel(i__) = rNULL__
if(ar__(162:184) .EQ. '') alpha2_24(i__) = rNULL__
if(ar__(203:221) .EQ. '') lambdaMax2_24(i__) = rNULL__
if(ar__(223:241) .EQ. '') lambdaMin2_24(i__) = rNULL__
if(ar__(243:260) .EQ. '') AGproxy(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ ID(i__),RAdeg(i__),DEdeg(i__),uv(i__),opt1(i__),opt2(i__),
+ opt3(i__),nir1(i__),nir2(i__),nir3(i__),mir1(i__),mir2(i__),
+ mir3a(i__),mir3b(i__),mir4(i__),mir5(i__),fir1(i__),fir2(i__),
+ fir3(i__),vizier(i__),nemesis(i__),extgalProb(i__),
+ extgalPred(i__),giantProb(i__),giantPred(i__),giantLabel(i__),
+ extgalLabel(i__),stratLabel(i__),massiveStar(i__),
+ MScontam(i__),alpha2_24(i__),Class2_24(i__),
+ lambdaMax2_24(i__),lambdaMin2_24(i__),AGproxy(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tabled2a.dat' ! Spectral types table
C Format for file interpretation
2 format(I5,1X,I5,1X,A16,3X,A1,1X,A4,1X,A4,1X,A19,1X,A16)
C Effective file loading
open(unit=1,status='old',file=
+'tabled2a.dat')
write(6,*) '....Loading file: tabled2a.dat'
do i__=1,22767
read(1,'(A79)')ar__1
read(ar__1,2)
+ spID(i__),ID_1(i__),SpType(i__),u_SpType(i__),e_SpType(i__),
+ photSpType(i__),Refsp(i__),n_SpType(i__)
c ..............Just test output...........
write(6,2)
+ spID(i__),ID_1(i__),SpType(i__),u_SpType(i__),e_SpType(i__),
+ photSpType(i__),Refsp(i__),n_SpType(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tabled2b.dat' ! Temperatures table
C Format for file interpretation
3 format(
+ I5,1X,I5,1X,A19,3X,F19.13,1X,F24.16,1X,A8,1X,F21.16,1X,F6.1,
+ 1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'tabled2b.dat')
write(6,*) '....Loading file: tabled2b.dat'
do i__=1,74508
read(1,'(A122)')ar__2
read(ar__2,3)
+ TempID(i__),ID_2(i__),RefTemp(i__),Teff(i__),e_Teff(i__),
+ n_Teff(i__),E_Teff_1(i__),Tbol(i__),e_Tbol(i__)
if(ar__2(35:53) .EQ. '') Teff(i__) = rNULL__
if(ar__2(55:78) .EQ. '') e_Teff(i__) = rNULL__
if(ar__2(89:109) .EQ. '') E_Teff_1(i__) = rNULL__
if(ar__2(111:116) .EQ. '') Tbol(i__) = rNULL__
if(ar__2(118:122) .EQ. '') e_Tbol(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ TempID(i__),ID_2(i__),RefTemp(i__),Teff(i__),e_Teff(i__),
+ n_Teff(i__),E_Teff_1(i__),Tbol(i__),e_Tbol(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tabled4.dat' ! Disk literature table
C Format for file interpretation
4 format(
+ I5,1X,I5,1X,A19,1X,A10,1X,A15,1X,A9,1X,A4,1X,F19.16,1X,F4.2,
+ 1X,A11)
C Effective file loading
open(unit=1,status='old',file=
+'tabled4.dat')
write(6,*) '....Loading file: tabled4.dat'
do i__=1,65647
read(1,'(A110)')ar__3
read(ar__3,4)
+ DiskID(i__),ID_3(i__),RefDisk(i__),Class(i__),Disk(i__),
+ n_Class(i__),IRext(i__),alpha(i__),e_alpha(i__),IRrange(i__)
if(ar__3(75:93) .EQ. '') alpha(i__) = rNULL__
if(ar__3(95:98) .EQ. '') e_alpha(i__) = rNULL__
c ..............Just test output...........
write(6,4)
+ DiskID(i__),ID_3(i__),RefDisk(i__),Class(i__),Disk(i__),
+ n_Class(i__),IRext(i__),alpha(i__),e_alpha(i__),IRrange(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tabled5.dat' ! Lithum data table
C Format for file interpretation
5 format(
+ I4,1X,I5,1X,A19,1X,E23.17,1X,F20.16,1X,A3,1X,F12.8,1X,F13.10,
+ 1X,F12.8,1X,A3,1X,F12.10,1X,F12.10,1X,F3.1)
C Effective file loading
open(unit=1,status='old',file=
+'tabled5.dat')
write(6,*) '....Loading file: tabled5.dat'
do i__=1,9185
read(1,'(A153)')ar__4
read(ar__4,5)
+ LiID(i__),ID_4(i__),RefLi(i__),EWLi(i__),e_EWLi(i__),
+ f_EWLi(i__),E_EWLi_1(i__),ALi(i__),s_ALi(i__),n_EWLi(i__),
+ e_ALi(i__),E_ALi_1(i__),l_ALi(i__)
if(ar__4(32:54) .EQ. '') EWLi(i__) = rNULL__
if(ar__4(56:75) .EQ. '') e_EWLi(i__) = rNULL__
if(ar__4(81:92) .EQ. '') E_EWLi_1(i__) = rNULL__
if(ar__4(94:106) .EQ. '') ALi(i__) = rNULL__
if(ar__4(108:119) .EQ. '') s_ALi(i__) = rNULL__
if(ar__4(125:136) .EQ. '') e_ALi(i__) = rNULL__
if(ar__4(138:149) .EQ. '') E_ALi_1(i__) = rNULL__
if(ar__4(151:153) .EQ. '') l_ALi(i__) = rNULL__
c ..............Just test output...........
write(6,5)
+ LiID(i__),ID_4(i__),RefLi(i__),EWLi(i__),e_EWLi(i__),
+ f_EWLi(i__),E_EWLi_1(i__),ALi(i__),s_ALi(i__),n_EWLi(i__),
+ e_ALi(i__),E_ALi_1(i__),l_ALi(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tabled6a.dat' ! Equivalent widths, gravity, table
C Format for file interpretation
6 format(
+ I4,1X,I5,1X,F19.16,1X,F18.16,1X,A2,1X,A19,1X,F4.2,1X,F4.2,1X,
+ A1,1X,F5.3,1X,F5.3,1X,F4.2,1X,F4.2,1X,F18.16,1X,F18.16,1X,
+ F18.16,1X,F18.16,1X,F18.16,1X,F18.16,1X,F5.3,1X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'tabled6a.dat')
write(6,*) '....Loading file: tabled6a.dat'
do i__=1,1880
read(1,'(A232)')ar__5
read(ar__5,6)
+ EWID(i__),ID_5(i__),EWNaI8200(i__),e_EWNaI8200(i__),
+ f_EWNaI8200(i__),RefEW(i__),EWNaI8195(i__),e_EWNaI8195(i__),
+ f_EWNaI8195(i__),EWCaH3_6965(i__),e_EWCaH3_6965(i__),
+ EWKI7665_7699(i__),e_EWKI7665_7699(i__),EWKI7665(i__),
+ e_EWKI7665(i__),EWKI7699(i__),e_EWKI7699(i__),EWTiO8442(i__),
+ e_EWTiO8442(i__),EWTiO57130(i__),e_EWTiO57130(i__)
if(ar__5(12:30) .EQ. '') EWNaI8200(i__) = rNULL__
if(ar__5(32:49) .EQ. '') e_EWNaI8200(i__) = rNULL__
if(ar__5(74:77) .EQ. '') EWNaI8195(i__) = rNULL__
if(ar__5(79:82) .EQ. '') e_EWNaI8195(i__) = rNULL__
if(ar__5(86:90) .EQ. '') EWCaH3_6965(i__) = rNULL__
if(ar__5(92:96) .EQ. '') e_EWCaH3_6965(i__) = rNULL__
if(ar__5(98:101) .EQ. '') EWKI7665_7699(i__) = rNULL__
if(ar__5(103:106) .EQ. '') e_EWKI7665_7699(i__) = rNULL__
if(ar__5(108:125) .EQ. '') EWKI7665(i__) = rNULL__
if(ar__5(127:144) .EQ. '') e_EWKI7665(i__) = rNULL__
if(ar__5(146:163) .EQ. '') EWKI7699(i__) = rNULL__
if(ar__5(165:182) .EQ. '') e_EWKI7699(i__) = rNULL__
if(ar__5(184:201) .EQ. '') EWTiO8442(i__) = rNULL__
if(ar__5(203:220) .EQ. '') e_EWTiO8442(i__) = rNULL__
if(ar__5(222:226) .EQ. '') EWTiO57130(i__) = rNULL__
if(ar__5(228:232) .EQ. '') e_EWTiO57130(i__) = rNULL__
c ..............Just test output...........
write(6,6)
+ EWID(i__),ID_5(i__),EWNaI8200(i__),e_EWNaI8200(i__),
+ f_EWNaI8200(i__),RefEW(i__),EWNaI8195(i__),e_EWNaI8195(i__),
+ f_EWNaI8195(i__),EWCaH3_6965(i__),e_EWCaH3_6965(i__),
+ EWKI7665_7699(i__),e_EWKI7665_7699(i__),EWKI7665(i__),
+ e_EWKI7665(i__),EWKI7699(i__),e_EWKI7699(i__),EWTiO8442(i__),
+ e_EWTiO8442(i__),EWTiO57130(i__),e_EWTiO57130(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tabled6b.dat' ! Surface gravity from spectra fit table
C Format for file interpretation
7 format(I5,1X,I5,1X,F11.8,1X,F12.10,1X,A19,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'tabled6b.dat')
write(6,*) '....Loading file: tabled6b.dat'
do i__=1,57649
read(1,'(A61)')ar__6
read(ar__6,7)
+ logID(i__),ID_6(i__),logg(i__),e_logg(i__),Reflogg(i__),
+ E_logg_1(i__)
if(ar__6(25:36) .EQ. '') e_logg(i__) = rNULL__
if(ar__6(58:61) .EQ. '') E_logg_1(i__) = rNULL__
c ..............Just test output...........
write(6,7)
+ logID(i__),ID_6(i__),logg(i__),e_logg(i__),Reflogg(i__),
+ E_logg_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tabled7a.dat' ! Equivalent width accretion, other, table
C Format for file interpretation
8 format(
+ I4,1X,I5,1X,A19,1X,F19.16,1X,F4.2,1X,A1,1X,F6.2,1X,F4.2,1X,
+ F4.2,1X,F7.3,1X,F5.2,1X,F6.2,1X,F4.2,1X,F7.3,1X,F4.2,1X,A1,1X,
+ F7.3,1X,F4.2,1X,A1,1X,F5.1,1X,F4.1,1X,A1,1X,F4.1,1X,F6.2,1X,
+ F5.2,1X,A1,1X,F5.2,1X,F7.2,1X,F5.2,1X,F7.2,1X,F5.2,1X,F7.2,1X,
+ F5.2,1X,F5.2,1X,F6.2,1X,F20.16,1X,F4.2,1X,A1,1X,F20.16,1X,
+ F4.2,1X,A1,1X,F6.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F6.2,1X,F4.2,1X,
+ F20.16,1X,A1,1X,F7.3,1X,F4.2,1X,F6.3,1X,F4.2,1X,F6.2,1X,F4.2,
+ 1X,F5.2,1X,F4.2,1X,E14.11,1X,F5.3,1X,A1,1X,F3.1,1X,F6.2,1X,
+ F5.2,1X,A1,1X,F4.1,1X,F5.2,1X,F4.2,1X,F7.2,1X,F4.2,1X,F6.2,1X,
+ F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'tabled7a.dat')
write(6,*) '....Loading file: tabled7a.dat'
do i__=1,5486
read(1,'(A469)')ar__7
read(ar__7,8)
+ accID(i__),ID_7(i__),RefEW_1(i__),EWHeI6678(i__),
+ e_EWHeI6678(i__),f_EWHeI6678(i__),EWHeI5876(i__),
+ e_EWHeI5876(i__),E_EWHeI5876_1(i__),EWHeI10830(i__),
+ e_EWHeI10830(i__),EWCaII8498(i__),e_EWCaII8498(i__),
+ EWCaII8542(i__),e_EWCaII8542(i__),f_EWCaII8542(i__),
+ EWCaII8662(i__),e_EWCaII8662(i__),f_EWCaII8662(i__),
+ EWCaH3934(i__),e_EWCaH3934(i__),f_EWCaH3934(i__),
+ E_EWCaH3934_1(i__),EWCaK3969(i__),e_EWCaK3969(i__),
+ f_EWCaK3969(i__),E_EWCaK3969_1(i__),EWOI5577(i__),
+ e_EWOI5577(i__),EWOI6300(i__),e_EWOI6300(i__),EWOI6364(i__),
+ e_EWOI6364(i__),EWOI7773(i__),EWOI8446(i__),EWSII6716(i__),
+ e_EWSII6716(i__),f_EWSII6716(i__),EWSII6731(i__),
+ e_EWSII6731(i__),f_EWSII6731(i__),EWFeII16440(i__),
+ e_EWFeII16440(i__),EWNII6548(i__),e_EWNII6548(i__),
+ EWNII6581(i__),e_EWNII6581(i__),EWNII6583(i__),
+ f_EWNII6583(i__),EWPab1282(i__),e_EWPab1282(i__),
+ EWPag1090(i__),e_EWPag1090(i__),EWBrg21661(i__),
+ e_EWBrg21661(i__),EWBr11_16810(i__),e_EWBr11_16810(i__),
+ EWHb4861(i__),e_EWHb4861(i__),f_EWHb4861(i__),
+ E_EWHb4861_1(i__),EWHg4340(i__),e_EWHg4340(i__),
+ f_EWHg4340(i__),E_EWHg4340_1(i__),EWH113771(i__),
+ e_EWH113771(i__),EWH2_1_0S_1(i__),e_EWH2_1_0S_1(i__),
+ EWCO_2_0(i__),e_EWCO_2_0(i__)
if(ar__7(32:50) .EQ. '') EWHeI6678(i__) = rNULL__
if(ar__7(52:55) .EQ. '') e_EWHeI6678(i__) = rNULL__
if(ar__7(59:64) .EQ. '') EWHeI5876(i__) = rNULL__
if(ar__7(66:69) .EQ. '') e_EWHeI5876(i__) = rNULL__
if(ar__7(71:74) .EQ. '') E_EWHeI5876_1(i__) = rNULL__
if(ar__7(76:82) .EQ. '') EWHeI10830(i__) = rNULL__
if(ar__7(84:88) .EQ. '') e_EWHeI10830(i__) = rNULL__
if(ar__7(90:95) .EQ. '') EWCaII8498(i__) = rNULL__
if(ar__7(97:100) .EQ. '') e_EWCaII8498(i__) = rNULL__
if(ar__7(102:108) .EQ. '') EWCaII8542(i__) = rNULL__
if(ar__7(110:113) .EQ. '') e_EWCaII8542(i__) = rNULL__
if(ar__7(117:123) .EQ. '') EWCaII8662(i__) = rNULL__
if(ar__7(125:128) .EQ. '') e_EWCaII8662(i__) = rNULL__
if(ar__7(132:136) .EQ. '') EWCaH3934(i__) = rNULL__
if(ar__7(138:141) .EQ. '') e_EWCaH3934(i__) = rNULL__
if(ar__7(145:148) .EQ. '') E_EWCaH3934_1(i__) = rNULL__
if(ar__7(150:155) .EQ. '') EWCaK3969(i__) = rNULL__
if(ar__7(157:161) .EQ. '') e_EWCaK3969(i__) = rNULL__
if(ar__7(165:169) .EQ. '') E_EWCaK3969_1(i__) = rNULL__
if(ar__7(171:177) .EQ. '') EWOI5577(i__) = rNULL__
if(ar__7(179:183) .EQ. '') e_EWOI5577(i__) = rNULL__
if(ar__7(185:191) .EQ. '') EWOI6300(i__) = rNULL__
if(ar__7(193:197) .EQ. '') e_EWOI6300(i__) = rNULL__
if(ar__7(199:205) .EQ. '') EWOI6364(i__) = rNULL__
if(ar__7(207:211) .EQ. '') e_EWOI6364(i__) = rNULL__
if(ar__7(213:217) .EQ. '') EWOI7773(i__) = rNULL__
if(ar__7(219:224) .EQ. '') EWOI8446(i__) = rNULL__
if(ar__7(226:245) .EQ. '') EWSII6716(i__) = rNULL__
if(ar__7(247:250) .EQ. '') e_EWSII6716(i__) = rNULL__
if(ar__7(254:273) .EQ. '') EWSII6731(i__) = rNULL__
if(ar__7(275:278) .EQ. '') e_EWSII6731(i__) = rNULL__
if(ar__7(282:287) .EQ. '') EWFeII16440(i__) = rNULL__
if(ar__7(289:292) .EQ. '') e_EWFeII16440(i__) = rNULL__
if(ar__7(294:298) .EQ. '') EWNII6548(i__) = rNULL__
if(ar__7(300:303) .EQ. '') e_EWNII6548(i__) = rNULL__
if(ar__7(305:310) .EQ. '') EWNII6581(i__) = rNULL__
if(ar__7(312:315) .EQ. '') e_EWNII6581(i__) = rNULL__
if(ar__7(317:336) .EQ. '') EWNII6583(i__) = rNULL__
if(ar__7(340:346) .EQ. '') EWPab1282(i__) = rNULL__
if(ar__7(348:351) .EQ. '') e_EWPab1282(i__) = rNULL__
if(ar__7(353:358) .EQ. '') EWPag1090(i__) = rNULL__
if(ar__7(360:363) .EQ. '') e_EWPag1090(i__) = rNULL__
if(ar__7(365:370) .EQ. '') EWBrg21661(i__) = rNULL__
if(ar__7(372:375) .EQ. '') e_EWBrg21661(i__) = rNULL__
if(ar__7(377:381) .EQ. '') EWBr11_16810(i__) = rNULL__
if(ar__7(383:386) .EQ. '') e_EWBr11_16810(i__) = rNULL__
if(ar__7(388:401) .EQ. '') EWHb4861(i__) = rNULL__
if(ar__7(403:407) .EQ. '') e_EWHb4861(i__) = rNULL__
if(ar__7(411:413) .EQ. '') E_EWHb4861_1(i__) = rNULL__
if(ar__7(415:420) .EQ. '') EWHg4340(i__) = rNULL__
if(ar__7(422:426) .EQ. '') e_EWHg4340(i__) = rNULL__
if(ar__7(430:433) .EQ. '') E_EWHg4340_1(i__) = rNULL__
if(ar__7(435:439) .EQ. '') EWH113771(i__) = rNULL__
if(ar__7(441:444) .EQ. '') e_EWH113771(i__) = rNULL__
if(ar__7(446:452) .EQ. '') EWH2_1_0S_1(i__) = rNULL__
if(ar__7(454:457) .EQ. '') e_EWH2_1_0S_1(i__) = rNULL__
if(ar__7(459:464) .EQ. '') EWCO_2_0(i__) = rNULL__
if(ar__7(466:469) .EQ. '') e_EWCO_2_0(i__) = rNULL__
c ..............Just test output...........
write(6,8)
+ accID(i__),ID_7(i__),RefEW_1(i__),EWHeI6678(i__),
+ e_EWHeI6678(i__),f_EWHeI6678(i__),EWHeI5876(i__),
+ e_EWHeI5876(i__),E_EWHeI5876_1(i__),EWHeI10830(i__),
+ e_EWHeI10830(i__),EWCaII8498(i__),e_EWCaII8498(i__),
+ EWCaII8542(i__),e_EWCaII8542(i__),f_EWCaII8542(i__),
+ EWCaII8662(i__),e_EWCaII8662(i__),f_EWCaII8662(i__),
+ EWCaH3934(i__),e_EWCaH3934(i__),f_EWCaH3934(i__),
+ E_EWCaH3934_1(i__),EWCaK3969(i__),e_EWCaK3969(i__),
+ f_EWCaK3969(i__),E_EWCaK3969_1(i__),EWOI5577(i__),
+ e_EWOI5577(i__),EWOI6300(i__),e_EWOI6300(i__),EWOI6364(i__),
+ e_EWOI6364(i__),EWOI7773(i__),EWOI8446(i__),EWSII6716(i__),
+ e_EWSII6716(i__),f_EWSII6716(i__),EWSII6731(i__),
+ e_EWSII6731(i__),f_EWSII6731(i__),EWFeII16440(i__),
+ e_EWFeII16440(i__),EWNII6548(i__),e_EWNII6548(i__),
+ EWNII6581(i__),e_EWNII6581(i__),EWNII6583(i__),
+ f_EWNII6583(i__),EWPab1282(i__),e_EWPab1282(i__),
+ EWPag1090(i__),e_EWPag1090(i__),EWBrg21661(i__),
+ e_EWBrg21661(i__),EWBr11_16810(i__),e_EWBr11_16810(i__),
+ EWHb4861(i__),e_EWHb4861(i__),f_EWHb4861(i__),
+ E_EWHb4861_1(i__),EWHg4340(i__),e_EWHg4340(i__),
+ f_EWHg4340(i__),E_EWHg4340_1(i__),EWH113771(i__),
+ e_EWH113771(i__),EWH2_1_0S_1(i__),e_EWH2_1_0S_1(i__),
+ EWCO_2_0(i__),e_EWCO_2_0(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tabled7b.dat' ! H{alpha} data table
C Format for file interpretation
9 format(
+ I5,1X,I5,1X,A19,3X,E19.16,1X,E19.16,1X,A1,1X,F8.3,1X,F5.3,1X,
+ A3,1X,F7.3,1X,F7.3,1X,F11.4)
C Effective file loading
open(unit=1,status='old',file=
+'tabled7b.dat')
write(6,*) '....Loading file: tabled7b.dat'
do i__=1,37449
read(1,'(A122)')ar__8
read(ar__8,9)
+ HaID(i__),ID_8(i__),RefEWHa(i__),EWHa(i__),e_EWHa(i__),
+ f_EWHa(i__),EWHac(i__),e_EWHac(i__),Halineflag(i__),
+ W10Ha(i__),e_W10Ha(i__),EWHaPhot(i__)
if(ar__8(35:53) .EQ. '') EWHa(i__) = rNULL__
if(ar__8(55:73) .EQ. '') e_EWHa(i__) = rNULL__
if(ar__8(77:84) .EQ. '') EWHac(i__) = rNULL__
if(ar__8(86:90) .EQ. '') e_EWHac(i__) = rNULL__
if(ar__8(96:102) .EQ. '') W10Ha(i__) = rNULL__
if(ar__8(104:110) .EQ. '') e_W10Ha(i__) = rNULL__
if(ar__8(112:122) .EQ. '') EWHaPhot(i__) = rNULL__
c ..............Just test output...........
write(6,9)
+ HaID(i__),ID_8(i__),RefEWHa(i__),EWHa(i__),e_EWHa(i__),
+ f_EWHa(i__),EWHac(i__),e_EWHac(i__),Halineflag(i__),
+ W10Ha(i__),e_W10Ha(i__),EWHaPhot(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tabled7c.dat' ! Veiling table
C Format for file interpretation
10 format(
+ I5,1X,I5,1X,A19,1X,F6.4,1X,F6.4,1X,F3.1,1X,F3.1,1X,F3.1,1X,
+ E8.6,1X,F11.6,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'tabled7c.dat')
write(6,*) '....Loading file: tabled7c.dat'
do i__=1,34358
read(1,'(A83)')ar__9
read(ar__9,10)
+ veilingID(i__),ID_9(i__),Refveiling(i__),veiling_opt(i__),
+ e_veiling_opt(i__),veiling_K(i__),E_veiling_K(i__),
+ e_veiling_K_1(i__),veiling_H(i__),e_veiling_H(i__),
+ veiling_7465(i__)
if(ar__9(33:38) .EQ. '') veiling_opt(i__) = rNULL__
if(ar__9(40:45) .EQ. '') e_veiling_opt(i__) = rNULL__
if(ar__9(47:49) .EQ. '') veiling_K(i__) = rNULL__
if(ar__9(51:53) .EQ. '') E_veiling_K(i__) = rNULL__
if(ar__9(55:57) .EQ. '') e_veiling_K_1(i__) = rNULL__
if(ar__9(59:66) .EQ. '') veiling_H(i__) = rNULL__
if(ar__9(68:78) .EQ. '') e_veiling_H(i__) = rNULL__
if(ar__9(80:83) .EQ. '') veiling_7465(i__) = rNULL__
c ..............Just test output...........
write(6,10)
+ veilingID(i__),ID_9(i__),Refveiling(i__),veiling_opt(i__),
+ e_veiling_opt(i__),veiling_K(i__),E_veiling_K(i__),
+ e_veiling_K_1(i__),veiling_H(i__),e_veiling_H(i__),
+ veiling_7465(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tabled8.dat' ! Radial velocities table
C Format for file interpretation
11 format(I5,1X,I5,1X,E19.16,1X,F26.16,1X,A19,1X,A2)
C Effective file loading
open(unit=1,status='old',file=
+'tabled8.dat')
write(6,*) '....Loading file: tabled8.dat'
do i__=1,60294
read(1,'(A81)')ar__10
read(ar__10,11)
+ rvID(i__),ID_10(i__),RV(i__),e_RV(i__),RefRV(i__),f_RV(i__)
if(ar__10(33:58) .EQ. '') e_RV(i__) = rNULL__
c ..............Just test output...........
write(6,11)
+ rvID(i__),ID_10(i__),RV(i__),e_RV(i__),RefRV(i__),f_RV(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tabled10.dat' ! Rotation table
C Format for file interpretation
12 format(
+ I5,1X,I5,2X,F20.16,1X,F6.2,1X,F6.2,1X,F6.2,1X,F5.2,1X,F4.2,1X,
+ F4.2,1X,E23.17,1X,A10,1X,F20.15,1X,A4,1X,F20.16,1X,A21,1X,
+ F18.16,1X,A22,1X,F18.15,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'tabled10.dat')
write(6,*) '....Loading file: tabled10.dat'
do i__=1,42132
read(1,'(A237)')ar__11
read(ar__11,12)
+ rotID(i__),ID_11(i__),Per(i__),Per2(i__),Per3(i__),Per4(i__),
+ Per5(i__),Per6(i__),Per7(i__),e_Per(i__),n_Per(i__),
+ vsini(i__),f_vsini(i__),e_vsini(i__),Refrot(i__),RotSep(i__),
+ E_Per_1(i__),E_vsini_1(i__),f_Per(i__)
if(ar__11(14:33) .EQ. '') Per(i__) = rNULL__
if(ar__11(35:40) .EQ. '') Per2(i__) = rNULL__
if(ar__11(42:47) .EQ. '') Per3(i__) = rNULL__
if(ar__11(49:54) .EQ. '') Per4(i__) = rNULL__
if(ar__11(56:60) .EQ. '') Per5(i__) = rNULL__
if(ar__11(62:65) .EQ. '') Per6(i__) = rNULL__
if(ar__11(67:70) .EQ. '') Per7(i__) = rNULL__
if(ar__11(72:94) .EQ. '') e_Per(i__) = rNULL__
if(ar__11(107:126) .EQ. '') vsini(i__) = rNULL__
if(ar__11(133:152) .EQ. '') e_vsini(i__) = rNULL__
if(ar__11(176:193) .EQ. '') RotSep(i__) = rNULL__
if(ar__11(218:235) .EQ. '') E_vsini_1(i__) = rNULL__
c ..............Just test output...........
write(6,12)
+ rotID(i__),ID_11(i__),Per(i__),Per2(i__),Per3(i__),Per4(i__),
+ Per5(i__),Per6(i__),Per7(i__),e_Per(i__),n_Per(i__),
+ vsini(i__),f_vsini(i__),e_vsini(i__),Refrot(i__),RotSep(i__),
+ E_Per_1(i__),E_vsini_1(i__),f_Per(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tabled11.dat' ! X-ray table
C Format for file interpretation
13 format(
+ I4,1X,I5,3X,A19,1X,E22.17,1X,A11,1X,E22.17,1X,A1,1X,F18.15,1X,
+ F4.2,1X,A11,1X,A1,1X,F18.15,1X,F19.16,1X,F18.15,1X,A1,1X,
+ F19.16,1X,F4.2,1X,F4.2,1X,A55,1X,A2,1X,F5.2,1X,F4.2,1X,F4.2,
+ 1X,A55,1X,F5.2,1X,F4.2,1X,F4.2,1X,A18,1X,A10,1X,E22.17)
C Effective file loading
open(unit=1,status='old',file=
+'tabled11.dat')
write(6,*) '....Loading file: tabled11.dat'
do i__=1,7926
read(1,'(A420)')ar__12
read(ar__12,13)
+ xrayID(i__),ID_12(i__),Refxray(i__),FX(i__),FXband(i__),
+ e_FX(i__),l_FX(i__),logLX(i__),e_logLX(i__),logLXband(i__),
+ l_logLX(i__),logNH(i__),e_logNH(i__),E_logNH_1(i__),
+ l_logNH(i__),HR1(i__),e_HR1(i__),E_HR1_1(i__),HR1band(i__),
+ l_HR1(i__),HR2(i__),e_HR2(i__),E_HR2_1(i__),HR2band(i__),
+ HR3(i__),e_HR3(i__),E_HR3_1(i__),HR3band(i__),Survey(i__),
+ sepX(i__)
if(ar__12(34:55) .EQ. '') FX(i__) = rNULL__
if(ar__12(69:90) .EQ. '') e_FX(i__) = rNULL__
if(ar__12(94:111) .EQ. '') logLX(i__) = rNULL__
if(ar__12(113:116) .EQ. '') e_logLX(i__) = rNULL__
if(ar__12(132:149) .EQ. '') logNH(i__) = rNULL__
if(ar__12(151:169) .EQ. '') e_logNH(i__) = rNULL__
if(ar__12(171:188) .EQ. '') E_logNH_1(i__) = rNULL__
if(ar__12(192:210) .EQ. '') HR1(i__) = rNULL__
if(ar__12(212:215) .EQ. '') e_HR1(i__) = rNULL__
if(ar__12(217:220) .EQ. '') E_HR1_1(i__) = rNULL__
if(ar__12(281:285) .EQ. '') HR2(i__) = rNULL__
if(ar__12(287:290) .EQ. '') e_HR2(i__) = rNULL__
if(ar__12(292:295) .EQ. '') E_HR2_1(i__) = rNULL__
if(ar__12(353:357) .EQ. '') HR3(i__) = rNULL__
if(ar__12(359:362) .EQ. '') e_HR3(i__) = rNULL__
if(ar__12(364:367) .EQ. '') E_HR3_1(i__) = rNULL__
if(ar__12(399:420) .EQ. '') sepX(i__) = rNULL__
c ..............Just test output...........
write(6,13)
+ xrayID(i__),ID_12(i__),Refxray(i__),FX(i__),FXband(i__),
+ e_FX(i__),l_FX(i__),logLX(i__),e_logLX(i__),logLXband(i__),
+ l_logLX(i__),logNH(i__),e_logNH(i__),E_logNH_1(i__),
+ l_logNH(i__),HR1(i__),e_HR1(i__),E_HR1_1(i__),HR1band(i__),
+ l_HR1(i__),HR2(i__),e_HR2(i__),E_HR2_1(i__),HR2band(i__),
+ HR3(i__),e_HR3(i__),E_HR3_1(i__),HR3band(i__),Survey(i__),
+ sepX(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tabled12.dat' ! Multiplicity table
C Format for file interpretation
14 format(I5,1X,I5,1X,A2,3X,A120,1X,A20)
C Effective file loading
open(unit=1,status='old',file=
+'tabled12.dat')
write(6,*) '....Loading file: tabled12.dat'
do i__=1,18929
read(1,'(A158)')ar__13
read(ar__13,14)
+ binID(i__),ID_13(i__),binType(i__),binExtType(i__),Refbin(i__)
c ..............Just test output...........
write(6,14)
+ binID(i__),ID_13(i__),binType(i__),binExtType(i__),Refbin(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end