Conversion of standardized ReadMe file for
file /./ftp/cats/J/MNRAS/452/173 into FORTRAN code for reading data files line by line.
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-Jun-08
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/MNRAS/452/173 Nearby young stars in Northern hemisphere (Binks+, 2015)
*================================================================================
*A kinematically unbiased search for nearby young stars in the Northern
*hemisphere based on SuperWASP rotation periods.
* Binks A.S., Jeffries R.D., Maxted P.F.L.
* <Mon. Not. R. Astron. Soc. 452, 173 (2015)>
* =2015MNRAS.452..173B
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'mgs.dat' ! Moving Group data
integer*4 nr__
parameter (nr__=10) ! Number of records
character*119 ar__ ! Full-size record
character*7 Name ! Moving group name
integer*4 NObs ! [4/48] Number of observations
integer*4 Age ! (Myr) [3/70] Lower value of Age interval
integer*4 B_Age ! (Myr) [5/120] Upper value of Age interval
character*1 n_Age ! [*] * indicates ages calculated from the
* lithium depletion boundary
character*1 u_Dist ! [~] Uncertainty flag on Dist
integer*4 Dist ! (pc) Lower value of distance interval
integer*4 B_Dist ! (pc) ? Upper value of distance interval
real*4 Uvel ! (km/s) ?=- U velocity
real*4 Vvel ! (km/s) ?=- V velocity
real*4 Wvel ! (km/s) ?=- W velocity
real*4 e_Uvel ! (km/s) rms uncertainty on Uvel
real*4 e_Vvel ! (km/s) rms uncertainty on Vvel
real*4 e_Wvel ! (km/s) rms uncertainty on Wvel
real*4 X ! (pc) ?=- X positions
real*4 Y ! (pc) ?=- Y positions
real*4 Z ! (pc) ?=- Z positions
real*4 e_X ! (pc) rms uncertainty on X
real*4 e_Y ! (pc) rms uncertainty on Y
real*4 e_Z ! (pc) rms uncertainty on Z
real*8 RAconv ! (deg) ?=- Convergent point in right ascension
real*8 DEconv ! (deg) ?=- convergent point in declination
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'age.dat' ! Ages for non-TLSPB young objects
integer*4 nr__1
parameter (nr__1=106) ! Number of records
character*43 ar__1 ! Full-size record
character*11 Name_1 ! Name of Target (SWHHMM+DDMM)
character*1 l_AgeG ! Limit flag on AgeG (1)
integer*4 AgeG ! (Myr) Age from gyrochronology
character*1 l_AgeLi ! Limit flag on AgeLi (1)
integer*4 AgeLi ! (Myr) Age from lithium, or lower value of interval
integer*4 B_AgeLi ! (Myr) ? Upper value of lithium age interval
character*1 l_AgeHa ! Limit flag on AgeHa (1)
integer*4 AgeHa ! (Myr) Age from H{alpha} or lower value of interval
integer*4 B_AgeHa ! (Myr) ? Upper value of H{alpha} age interval
character*1 l_AgeR ! Limit flag on AgeR (1)
integer*4 AgeR ! (Myr) Age from Rsini
*Note (1): Any ages measured as >1000 Myr are lower limits and are at least
* older than the Hyades and for the purpose of this work the precision of
* these ages are unimportant as they are not considered for further analysis.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'age_b.dat' ! Ages for objects in TLSPB systems
integer*4 nr__2
parameter (nr__2=14) ! Number of records
character*43 ar__2 ! Full-size record
character*11 Name_2 ! Name of Target (SWHHMM+DDMM)
character*1 l_AgeG_1 ! Limit flag on AgeG (1)
integer*4 AgeG_1 ! (Myr) Age from gyrochronology
character*1 l_AgeLi_1 ! Limit flag on AgeLi (1)
integer*4 AgeLi_1 ! (Myr) Age from lithium, or lower value of interval
integer*4 B_AgeLi_1 ! (Myr) ? Upper value of lithium age interval
character*1 l_AgeHa_1 ! Limit flag on AgeHa (1)
integer*4 AgeHa_1 ! (Myr) Age from H{alpha} or lower value of interval
integer*4 B_AgeHa_1 ! (Myr) ? Upper value of H{alpha} age interval
character*1 l_AgeR_1 ! Limit flag on AgeR (1)
integer*4 AgeR_1 ! (Myr) Age from Rsini
*Note (1): Any ages measured as >1000 Myr are lower limits and are at least
* older than the Hyades and for the purpose of this work the precision of
* these ages are unimportant as they are not considered for further analysis.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'params.dat' ! EWs, temperatures and abundances
integer*4 nr__3
parameter (nr__3=106) ! Number of records
character*99 ar__3 ! Full-size record
character*11 Name_3 ! Name of Target (SWHHMM+DDMM)
character*1 l_EWLi ! Limit flag on EWLi (2-{sigma} upper limits)
integer*4 EWLi ! (0.1pm) Lithium equivalent width before Fe correction
integer*4 EWLic ! (0.1pm) ?=-999 Lithium equivalent width subsequent to
* deblending the 6707.4 FeI line (m{AA}) (2)
integer*4 e_EWLic ! (0.1pm) ? Error on EWLic calculated using Cayrel de
* Strobel (1988IAUS..132..345C) formulation
real*4 EWHa ! (0.1nm) H{alpha} equivalent width ({AA}) (3)
real*4 Bmag ! (mag) B magnitude
real*4 e_Bmag ! (mag) ? rms uncertainty on Bmag
real*4 Vmag ! (mag) V magnitude
real*4 e_Vmag ! (mag) ? rms uncertainty on Vmag
real*8 Kmag ! (mag) ?=-99 K magnitude
real*4 e_Kmag ! (mag) ? rms uncertainty on Kmag
integer*4 Teff ! (K) Effective temperature (4)
character*2 SpT ! Spectral type
real*4 A_Li ! ([-]) ?=-9.99 Lithium abundance, calculated using
* 12+log(N(Li)/N(H))
real*4 E_A_Li ! ([-]) ?=-9.99 1-{sigma} upper error on A(Li)
real*4 e_A_Li_1 ! ([-]) ?=-9.99 1-{sigma} lower error on A(Li)
*Note (2): using the prescription of Soderblom, 2010ARA&A..48..581S.
*Note (3): Positive/negative values represent absorption/emission.
*Note (4): Effective temperature based on the 5th order V-K calibration in
* Pecaut & Mamajek (2013ApJS..208....9P, Cat. J/ApJS/208/9).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'params_b.dat' ! EWs, temperatures, abundances in TLSPB systems
integer*4 nr__4
parameter (nr__4=14) ! Number of records
character*99 ar__4 ! Full-size record
character*11 Name_4 ! Name of Target (SWHHMM+DDMM)
character*1 l_EWLi_1 ! Limit flag on EWLi (2-{sigma} upper limits)
integer*4 EWLi_1 ! (0.1pm) Lithium equivalent width before Fe correction
integer*4 EWLic_1 ! (0.1pm) ?=-999 Lithium equivalent width subsequent to
* deblending the 6707.4 FeI line (m{AA}) (2)
integer*4 e_EWLic_1 ! (0.1pm) ? Error on EWLic calculated using Cayrel de
* Strobel (1988IAUS..132..345C) formulation
real*4 EWHa_1 ! (0.1nm) H{alpha} equivalent width ({AA}) (3)
real*4 Bmag_1 ! (mag) B magnitude
real*4 e_Bmag_1 ! (mag) ? rms uncertainty on Bmag
real*4 Vmag_1 ! (mag) V magnitude
real*4 e_Vmag_1 ! (mag) ? rms uncertainty on Vmag
real*8 Kmag_1 ! (mag) ?=-99 K magnitude
real*4 e_Kmag_1 ! (mag) ? rms uncertainty on Kmag
integer*4 Teff_1 ! (K) Effective temperature (4)
character*2 SpT_1 ! Spectral type
real*4 A_Li_1 ! ([-]) ?=-9.99 Lithium abundance, calculated using
* 12+log(N(Li)/N(H))
real*4 E_A_Li_1 ! ([-]) ?=-9.99 1-{sigma} upper error on A(Li)
real*4 e_A_Li_1 ! ([-]) ?=-9.99 1-{sigma} lower error on A(Li)
*Note (2): using the prescription of Soderblom, 2010ARA&A..48..581S.
*Note (3): Positive/negative values represent absorption/emission.
*Note (4): Effective temperature based on the 5th order V-K calibration in
* Pecaut & Mamajek (2013ApJS..208....9P, Cat. J/ApJS/208/9).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'kin.dat' ! Kinematics
integer*4 nr__5
parameter (nr__5=106) ! Number of records
character*108 ar__5 ! Full-size record
character*11 Name_5 ! Name of Target (SWHHMM+DDMM)
real*4 pmRA ! (mas/yr) Proper motion in right ascension
real*4 e_pmRA ! (mas/yr) rms uncertainty on pmRA
real*4 pmDE ! (mas/yr) Proper motion in declination
real*4 e_pmDE ! (mas/yr) rms uncertainty on pmDE
real*4 Plx ! (mas) Parallax calculated photometrically
real*4 B_Plx ! (mas) ? Upper value of parallax when interval
real*4 e_Plx ! (mas) ? rms uncertainty on Plx
real*4 Uvel_1 ! (km/s) ?=-999 U velocity
real*4 e_Uvel_1 ! (km/s) ? First error bar on Uvel from the
* uncertainties in proper motion and RV
real*4 s_Uvel ! (km/s) ? Second error bar on Uvel from the
* distance uncertainty
real*4 Vvel_1 ! (km/s) ?=-999 V velocity
real*4 e_Vvel_1 ! (km/s) ? First error bar on Vvel from the
* uncertainties in proper motion and RV
real*4 s_Vvel ! (km/s) ? Second error bar on Vvel from the
* distance uncertainty
real*4 Wvel_1 ! (km/s) ?=-999 W velocity
real*4 e_Wvel_1 ! (km/s) ? First error bar on Wvel from the
* uncertainties in proper motion and RV
real*4 s_Wvel ! (km/s) ? Second error bar on Wvel from the
* distance uncertainty
real*4 Vsini ! (km/s) ?=-999 Projected rotational velocity
real*4 e_Vsini ! (km/s) ? rms uncertainty on Vsini
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'kin_b.dat' ! Kinematics TLSPB systems
integer*4 nr__6
parameter (nr__6=14) ! Number of records
character*60 ar__6 ! Full-size record
character*11 Name_6 ! Name of Target (SWHHMM+DDMM)
real*4 pmRA_1 ! (mas/yr) Proper motion in right ascension
real*4 e_pmRA_1 ! (mas/yr) rms uncertainty on pmRA
real*4 pmDE_1 ! (mas/yr) Proper motion in declination
real*4 e_pmDE_1 ! (mas/yr) rms uncertainty on pmDE
real*4 Plx_1 ! (mas) Parallax calculated photometrically
real*4 B_Plx_1 ! (mas) ? Upper value of parallax when interval
real*4 e_Plx_1 ! (mas) ? rms uncertainty on Plx
real*4 Vsini_1 ! (km/s) ?=-999 Projected rotational velocity
real*4 e_Vsini_1 ! (km/s) ? rms uncertainty on Vsini
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'rvs.dat' ! Radial velocities (for the 106 sources)
integer*4 nr__7
parameter (nr__7=120) ! Number of records
character*49 ar__7 ! Full-size record
character*12 Name_7 ! Name of Target (SWHHMM+DDMM)
real*8 HJD ! (d) Heliocentric Julian Date (HJD-2450000)
real*4 RVi ! (km/s) ?=-999 Individual radial velocity measurement
real*4 e_RVi ! (km/s) ? rms uncertainty on RVi
character*1 n_RVi ! [i] i indicates measurements made using the
* Isaac Newton Telescope
real*4 RVf ! (km/s) ? Final radial velocity measurement
real*4 e_RVf ! (km/s) ? rms uncertainty on RVf
integer*4 Bin ! [1/5] Binary score (1)
*Note (1): Binary score used to distinguish likely single stars (1) from those
* very likely to be binaries (5) as follows:
* 1 = likely single stars (consistent, low-uncertainty RV measurements for 2
* or more spectra)
* 2 = objects which had a single spectrum, an RV uncertainty <5km/s and a
* distinct single peak in the CCF
* 3 = status of the star from the CCF was unclear (presumably as a result of
* poor SNR and/or large vsini), resulting in either an indeterminate RV
* or a RV uncertainty >5km/s
* 4 = only one spectrum which results in (a) a clear, multipeaked CCF and (b)
* an RV error <5km/s
* 5 = very likely to be binaries (detected RV differences >5km/s for a target
* on separate nights or (if there was only one measurement) if there were
* literature sources indicating that either the object is a close binary
* or report a RV measurement >5km/s discrepant from our measurement
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'rvs_b.dat' ! Radial velocities of TLSPB systems
(for the 14 sources)
integer*4 nr__8
parameter (nr__8=26) ! Number of records
character*49 ar__8 ! Full-size record
character*12 Name_8 ! Name of Target (SWHHMM+DDMM)
real*8 HJD_1 ! (d) Heliocentric Julian Date (HJD-2450000)
real*4 RVi_1 ! (km/s) ?=-999 Individual radial velocity measurement
real*4 e_RVi_1 ! (km/s) ? rms uncertainty on RVi
character*1 n_RVi_1 ! [i] i indicates measurements made using the
* Isaac Newton Telescope
real*4 RVf_1 ! (km/s) ? Final radial velocity measurement
real*4 e_RVf_1 ! (km/s) ? rms uncertainty on RVf
integer*4 Bin_1 ! [1/5] Binary score (1)
*Note (1): Binary score used to distinguish likely single stars (1) from those
* very likely to be binaries (5) as follows:
* 1 = likely single stars (consistent, low-uncertainty RV measurements for 2
* or more spectra)
* 2 = objects which had a single spectrum, an RV uncertainty <5km/s and a
* distinct single peak in the CCF
* 3 = status of the star from the CCF was unclear (presumably as a result of
* poor SNR and/or large vsini), resulting in either an indeterminate RV
* or a RV uncertainty >5km/s
* 4 = only one spectrum which results in (a) a clear, multipeaked CCF and (b)
* an RV error <5km/s
* 5 = very likely to be binaries (detected RV differences >5km/s for a target
* on separate nights or (if there was only one measurement) if there were
* literature sources indicating that either the object is a close binary
* or report a RV measurement >5km/s discrepant from our measurement
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'p_xray.dat' ! Periods Xray
integer*4 nr__9
parameter (nr__9=106) ! Number of records
character*89 ar__9 ! Full-size record
character*18 v1SWASP ! Identification in the SuperWASP catalog
* 1SWASP J (HHMMSS.ss+DDMMSS.s)
real*4 Vmag1 ! (mag) V magnitude
real*4 V_K ! (mag) V-K colour index
real*4 Per ! (d) [0.7/20] Rotation period
real*4 e_Per ! (d) [0.002/3.5]? rms uncertainty on Per
integer*4 Chisq ! {Delta}{chi}^2^ value (section 2 of the paper)
integer*4 Ns ! [1/7] Number of seasonal SuperWASP lightcurves
* used in determining the rotation period
character*1 q_Per ! [A-C] Quality of period measurement, A=best
* (section 2 of the paper)
character*3 Date ! Month and year for which spectroscopy was
* first obtained for the object,
* J11 = June 2011, D12 = December 2012
real*4 CR ! (s-1) [0.03/0.4] X-ray count rate
real*4 HR1 ! [-1/1] Hardness ratio 1
real*4 logLx_b ! ([-]) [-4/-2.9] Base 10 logarithm of the fractional
* X-ray to bolometric luminosity
character*11 Name_9 ! Name of Target (SWHHMM+DDMM)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'p_xray_b.dat' ! Periods Xray TLSPB systems
integer*4 nr__10
parameter (nr__10=14) ! Number of records
character*89 ar__10 ! Full-size record
character*18 v1SWASP_1 ! Identification in the SuperWASP catalog
* 1SWASP J (HHMMSS.ss+DDMMSS.s)
real*4 Vmag1_1 ! (mag) V magnitude
real*4 V_K_1 ! (mag) V-K colour index
real*4 Per_1 ! (d) [0.7/20] Rotation period
real*4 e_Per_1 ! (d) [0.002/3.5]? rms uncertainty on Per
integer*4 Chisq_1 ! {Delta}{chi}^2^ value (section 2 of the paper)
integer*4 Ns_1 ! [1/7] Number of seasonal SuperWASP lightcurves
* used in determining the rotation period
character*1 q_Per_1 ! [A-C] Quality of period measurement, A=best
* (section 2 of the paper)
character*3 Date_1 ! Month and year for which spectroscopy was
* first obtained for the object,
* J11 = June 2011, D12 = December 2012
real*4 CR_1 ! (s-1) [0.03/0.4] X-ray count rate
real*4 HR1_1 ! [-1/1] Hardness ratio 1
real*4 logLx_b_1 ! ([-]) [-4/-2.9] Base 10 logarithm of the fractional
* X-ray to bolometric luminosity
character*11 Name_10 ! Name of Target (SWHHMM+DDMM)
C=============================================================================
C Loading file 'mgs.dat' ! Moving Group data
C Format for file interpretation
1 format(
+ A7,1X,I2,1X,I2,1X,I3,A1,A1,I3,1X,I3,1X,F6.2,1X,F6.2,1X,F6.2,
+ 1X,F4.2,1X,F4.2,1X,F4.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F5.2,1X,
+ F5.2,1X,F5.2,1X,F8.4,1X,F8.4)
C Effective file loading
open(unit=1,status='old',file=
+'mgs.dat')
write(6,*) '....Loading file: mgs.dat'
do i__=1,10
read(1,'(A119)')ar__
read(ar__,1)
+ Name,NObs,Age,B_Age,n_Age,u_Dist,Dist,B_Dist,Uvel,Vvel,Wvel,
+ e_Uvel,e_Vvel,e_Wvel,X,Y,Z,e_X,e_Y,e_Z,RAconv,DEconv
if(ar__(24:26) .EQ. '') B_Dist = iNULL__
if (idig(ar__(28:33)).EQ.0) Uvel = rNULL__
if (idig(ar__(35:40)).EQ.0) Vvel = rNULL__
if (idig(ar__(42:47)).EQ.0) Wvel = rNULL__
if (idig(ar__(64:69)).EQ.0) X = rNULL__
if (idig(ar__(71:76)).EQ.0) Y = rNULL__
if (idig(ar__(78:83)).EQ.0) Z = rNULL__
if (idig(ar__(103:110)).EQ.0) RAconv = rNULL__
if (idig(ar__(112:119)).EQ.0) DEconv = rNULL__
c ..............Just test output...........
write(6,1)
+ Name,NObs,Age,B_Age,n_Age,u_Dist,Dist,B_Dist,Uvel,Vvel,Wvel,
+ e_Uvel,e_Vvel,e_Wvel,X,Y,Z,e_X,e_Y,e_Z,RAconv,DEconv
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'age.dat' ! Ages for non-TLSPB young objects
C Format for file interpretation
2 format(A11,1X,A1,I4,1X,A1,I4,1X,I3,2X,A1,I3,1X,I3,1X,A1,I4)
C Effective file loading
open(unit=1,status='old',file=
+'age.dat')
write(6,*) '....Loading file: age.dat'
do i__=1,106
read(1,'(A43)')ar__1
read(ar__1,2)
+ Name_1,l_AgeG,AgeG,l_AgeLi,AgeLi,B_AgeLi,l_AgeHa,AgeHa,
+ B_AgeHa,l_AgeR,AgeR
if(ar__1(25:27) .EQ. '') B_AgeLi = iNULL__
if(ar__1(35:37) .EQ. '') B_AgeHa = iNULL__
c ..............Just test output...........
write(6,2)
+ Name_1,l_AgeG,AgeG,l_AgeLi,AgeLi,B_AgeLi,l_AgeHa,AgeHa,
+ B_AgeHa,l_AgeR,AgeR
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'age_b.dat' ! Ages for objects in TLSPB systems
C Format for file interpretation
3 format(A11,1X,A1,I4,1X,A1,I4,1X,I3,2X,A1,I3,1X,I3,1X,A1,I4)
C Effective file loading
open(unit=1,status='old',file=
+'age_b.dat')
write(6,*) '....Loading file: age_b.dat'
do i__=1,14
read(1,'(A43)')ar__2
read(ar__2,3)
+ Name_2,l_AgeG_1,AgeG_1,l_AgeLi_1,AgeLi_1,B_AgeLi_1,l_AgeHa_1,
+ AgeHa_1,B_AgeHa_1,l_AgeR_1,AgeR_1
if(ar__2(25:27) .EQ. '') B_AgeLi_1 = iNULL__
if(ar__2(35:37) .EQ. '') B_AgeHa_1 = iNULL__
c ..............Just test output...........
write(6,3)
+ Name_2,l_AgeG_1,AgeG_1,l_AgeLi_1,AgeLi_1,B_AgeLi_1,l_AgeHa_1,
+ AgeHa_1,B_AgeHa_1,l_AgeR_1,AgeR_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'params.dat' ! EWs, temperatures and abundances
C Format for file interpretation
4 format(
+ A11,1X,A1,I3,1X,I4,1X,I3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,
+ F5.3,1X,F7.3,1X,F5.3,1X,I4,1X,A2,1X,F5.2,1X,F5.2,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'params.dat')
write(6,*) '....Loading file: params.dat'
do i__=1,106
read(1,'(A99)')ar__3
read(ar__3,4)
+ Name_3,l_EWLi,EWLi,EWLic,e_EWLic,EWHa,Bmag,e_Bmag,Vmag,e_Vmag,
+ Kmag,e_Kmag,Teff,SpT,A_Li,E_A_Li,e_A_Li_1
if(ar__3(23:25) .EQ. '') e_EWLic = iNULL__
if(ar__3(41:46) .EQ. '') e_Bmag = rNULL__
if(ar__3(55:59) .EQ. '') e_Vmag = rNULL__
if(ar__3(69:73) .EQ. '') e_Kmag = rNULL__
c ..............Just test output...........
write(6,4)
+ Name_3,l_EWLi,EWLi,EWLic,e_EWLic,EWHa,Bmag,e_Bmag,Vmag,e_Vmag,
+ Kmag,e_Kmag,Teff,SpT,A_Li,E_A_Li,e_A_Li_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'params_b.dat' ! EWs, temperatures, abundances in TLSPB systems
C Format for file interpretation
5 format(
+ A11,1X,A1,I3,1X,I4,1X,I3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,
+ F5.3,1X,F7.3,1X,F5.3,1X,I4,1X,A2,1X,F5.2,1X,F5.2,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'params_b.dat')
write(6,*) '....Loading file: params_b.dat'
do i__=1,14
read(1,'(A99)')ar__4
read(ar__4,5)
+ Name_4,l_EWLi_1,EWLi_1,EWLic_1,e_EWLic_1,EWHa_1,Bmag_1,
+ e_Bmag_1,Vmag_1,e_Vmag_1,Kmag_1,e_Kmag_1,Teff_1,SpT_1,A_Li_1,
+ E_A_Li_1,e_A_Li_1
if(ar__4(23:25) .EQ. '') e_EWLic_1 = iNULL__
if(ar__4(41:46) .EQ. '') e_Bmag_1 = rNULL__
if(ar__4(55:59) .EQ. '') e_Vmag_1 = rNULL__
if(ar__4(69:73) .EQ. '') e_Kmag_1 = rNULL__
c ..............Just test output...........
write(6,5)
+ Name_4,l_EWLi_1,EWLi_1,EWLic_1,e_EWLic_1,EWHa_1,Bmag_1,
+ e_Bmag_1,Vmag_1,e_Vmag_1,Kmag_1,e_Kmag_1,Teff_1,SpT_1,A_Li_1,
+ E_A_Li_1,e_A_Li_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'kin.dat' ! Kinematics
C Format for file interpretation
6 format(
+ A11,1X,F6.1,1X,F4.1,1X,F6.1,1X,F4.1,1X,F5.2,1X,F5.2,1X,F4.2,
+ 1X,F6.1,1X,F3.1,1X,F3.1,1X,F6.1,1X,F3.1,1X,F3.1,1X,F6.1,1X,
+ F3.1,1X,F3.1,F6.1,1X,F4.1)
C Effective file loading
open(unit=1,status='old',file=
+'kin.dat')
write(6,*) '....Loading file: kin.dat'
do i__=1,106
read(1,'(A108)')ar__5
read(ar__5,6)
+ Name_5,pmRA,e_pmRA,pmDE,e_pmDE,Plx,B_Plx,e_Plx,Uvel_1,
+ e_Uvel_1,s_Uvel,Vvel_1,e_Vvel_1,s_Vvel,Wvel_1,e_Wvel_1,s_Wvel,
+ Vsini,e_Vsini
if(ar__5(43:47) .EQ. '') B_Plx = rNULL__
if(ar__5(49:52) .EQ. '') e_Plx = rNULL__
if(ar__5(61:63) .EQ. '') e_Uvel_1 = rNULL__
if(ar__5(65:67) .EQ. '') s_Uvel = rNULL__
if(ar__5(76:78) .EQ. '') e_Vvel_1 = rNULL__
if(ar__5(80:82) .EQ. '') s_Vvel = rNULL__
if(ar__5(91:93) .EQ. '') e_Wvel_1 = rNULL__
if(ar__5(95:97) .EQ. '') s_Wvel = rNULL__
if(ar__5(105:108) .EQ. '') e_Vsini = rNULL__
c ..............Just test output...........
write(6,6)
+ Name_5,pmRA,e_pmRA,pmDE,e_pmDE,Plx,B_Plx,e_Plx,Uvel_1,
+ e_Uvel_1,s_Uvel,Vvel_1,e_Vvel_1,s_Vvel,Wvel_1,e_Wvel_1,s_Wvel,
+ Vsini,e_Vsini
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'kin_b.dat' ! Kinematics TLSPB systems
C Format for file interpretation
7 format(
+ A11,1X,F5.1,1X,F3.1,1X,F6.1,1X,F3.1,1X,F5.2,1X,F5.2,1X,F4.2,
+ 1X,F6.1,1X,F3.1)
C Effective file loading
open(unit=1,status='old',file=
+'kin_b.dat')
write(6,*) '....Loading file: kin_b.dat'
do i__=1,14
read(1,'(A60)')ar__6
read(ar__6,7)
+ Name_6,pmRA_1,e_pmRA_1,pmDE_1,e_pmDE_1,Plx_1,B_Plx_1,e_Plx_1,
+ Vsini_1,e_Vsini_1
if(ar__6(40:44) .EQ. '') B_Plx_1 = rNULL__
if(ar__6(46:49) .EQ. '') e_Plx_1 = rNULL__
if(ar__6(58:60) .EQ. '') e_Vsini_1 = rNULL__
c ..............Just test output...........
write(6,7)
+ Name_6,pmRA_1,e_pmRA_1,pmDE_1,e_pmDE_1,Plx_1,B_Plx_1,e_Plx_1,
+ Vsini_1,e_Vsini_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'rvs.dat' ! Radial velocities (for the 106 sources)
C Format for file interpretation
8 format(A12,F8.3,2X,F6.1,2X,F4.1,A1,F6.1,3X,F3.1,1X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'rvs.dat')
write(6,*) '....Loading file: rvs.dat'
do i__=1,120
read(1,'(A49)')ar__7
read(ar__7,8)Name_7,HJD,RVi,e_RVi,n_RVi,RVf,e_RVf,Bin
if(ar__7(31:34) .EQ. '') e_RVi = rNULL__
if(ar__7(36:41) .EQ. '') RVf = rNULL__
if(ar__7(45:47) .EQ. '') e_RVf = rNULL__
c ..............Just test output...........
write(6,8)Name_7,HJD,RVi,e_RVi,n_RVi,RVf,e_RVf,Bin
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'rvs_b.dat' ! Radial velocities of TLSPB systems
* (for the 14 sources)
C Format for file interpretation
9 format(A12,F8.3,2X,F6.1,2X,F4.1,A1,F6.1,3X,F3.1,1X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'rvs_b.dat')
write(6,*) '....Loading file: rvs_b.dat'
do i__=1,26
read(1,'(A49)')ar__8
read(ar__8,9)
+ Name_8,HJD_1,RVi_1,e_RVi_1,n_RVi_1,RVf_1,e_RVf_1,Bin_1
if(ar__8(31:34) .EQ. '') e_RVi_1 = rNULL__
if(ar__8(36:41) .EQ. '') RVf_1 = rNULL__
if(ar__8(45:47) .EQ. '') e_RVf_1 = rNULL__
c ..............Just test output...........
write(6,9)
+ Name_8,HJD_1,RVi_1,e_RVi_1,n_RVi_1,RVf_1,e_RVf_1,Bin_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'p_xray.dat' ! Periods Xray
C Format for file interpretation
10 format(
+ A18,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,I5,1X,I2,1X,A1,1X,A3,
+ 1X,F5.3,1X,F5.2,1X,F5.2,1X,A11)
C Effective file loading
open(unit=1,status='old',file=
+'p_xray.dat')
write(6,*) '....Loading file: p_xray.dat'
do i__=1,106
read(1,'(A89)')ar__9
read(ar__9,10)
+ v1SWASP,Vmag1,V_K,Per,e_Per,Chisq,Ns,q_Per,Date,CR,HR1,
+ logLx_b,Name_9
if(ar__9(40:44) .EQ. '') e_Per = rNULL__
c ..............Just test output...........
write(6,10)
+ v1SWASP,Vmag1,V_K,Per,e_Per,Chisq,Ns,q_Per,Date,CR,HR1,
+ logLx_b,Name_9
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'p_xray_b.dat' ! Periods Xray TLSPB systems
C Format for file interpretation
11 format(
+ A18,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,I5,1X,I2,1X,A1,1X,A3,
+ 1X,F5.3,1X,F5.2,1X,F5.2,1X,A11)
C Effective file loading
open(unit=1,status='old',file=
+'p_xray_b.dat')
write(6,*) '....Loading file: p_xray_b.dat'
do i__=1,14
read(1,'(A89)')ar__10
read(ar__10,11)
+ v1SWASP_1,Vmag1_1,V_K_1,Per_1,e_Per_1,Chisq_1,Ns_1,q_Per_1,
+ Date_1,CR_1,HR1_1,logLx_b_1,Name_10
if(ar__10(40:44) .EQ. '') e_Per_1 = rNULL__
c ..............Just test output...........
write(6,11)
+ v1SWASP_1,Vmag1_1,V_K_1,Per_1,e_Per_1,Chisq_1,Ns_1,q_Per_1,
+ Date_1,CR_1,HR1_1,logLx_b_1,Name_10
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end
C Locate position of first digit in string; or return 0
integer function idig(c)
character*(*) c
character*1 c1
integer lc,i
lc=len(c)
idig=0
do i=1,lc
if(c(i:i).ne.' ') go to 1
end do
1 if(i.gt.lc) return
c1=c(i:i)
if(c1.eq.'.'.or.c1.eq.'-'.or.c1.eq.'+') i=i+1
if(i.gt.lc) return
c1=c(i:i)
if(c1.ge.'0'.and.c1.le.'9') idig=i
return
end