Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/771/129 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-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/ApJ/771/129 Submillimetric Class II sources of Taurus (Andrews+, 2013)
*================================================================================
*The mass dependence between protoplanetary disks and their stellar hosts.
* Andrews S.M., Rosenfeld K.A., Kraus A.L., Wilner D.J.
* <Astrophys. J., 771, 129 (2013)>
* =2013ApJ...771..129A
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! Millimeter-wave continuum photometry catalog
integer*4 nr__
parameter (nr__=210) ! Number of records
character*84 ar__ ! Full-size record
character*20 Name (nr__) ! Name
character*1 f_Name (nr__) ! [e] Note for V807 Tau B (1)
character*16 SpT (nr__) ! MK spectral type(s) (2)
character*6 r_SpT (nr__) ! Reference(s) or "est"imated flag (3)
character*1 l_F0_89 (nr__) ! 3{sigma} upper limit flag on F0.89
real*4 F0_89 (nr__) ! (Jy) [0.001/1.4] 0.89mm flux (338GHz)
real*4 e_F0_89 (nr__) ! (Jy) [0.001/0.2]? F0.89 uncertainty
character*1 l_F1_3 (nr__) ! 3{sigma} upper limit flag on F1.3
real*4 F1_3 (nr__) ! (Jy) [0.001/0.5] 1.3mm flux (225GHz)
real*4 e_F1_3 (nr__) ! (Jy) [0.0004/0.04]? F1.3 uncertainty
character*6 Notes (nr__) ! Notes (4)
*Note (1):
* e = The companions V807 Tau B (itself a binary) and V710 Tau B have no
* observational disk signatures, and are therefore excluded from this
* Class II sample. The V807 Tau system has no resolved mm-wave measurement,
* but we associate all of the composite emission with the primary
* component. Harris et al. (2012, J/ApJ/751/115) found no mm-wave emission
* associated with V710 Tau B.
*Note (2): Unless otherwise specified, the assumed spectral type uncertainty
* is +/-1 subclass.
*Note (3): The spectral classifications were collated from the references listed
* in refs.dat file. The classifications for FT Tau, ITG 1, and
* IRAS 04370+2559 are unknown: we adopt the wide spectral type bins
* advocated by Luhman et al. (2010ApJS..186..111L), listed here in
* parenthesis.
*Note (4): The "m" and "e" flags denote whether each source has a flux density
* measurement available at (or near) 0.89 and 1.3mm: "m" implies direct
* measurement and "e" implies an extrapolation as described in the text.
* Notes in parenthesis indicate unresolved multiple systems where
* composite measurements were applied to the primaries, and upper
* limits to the secondaries (see Section 3.1).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! Stellar host masses and ages
integer*4 nr__1
parameter (nr__1=212) ! Number of records
character*130 ar__1 ! Full-size record
character*20 Name_1 (nr__1) ! Name
character*1 f_Name_1 (nr__1) ! [ab] Flag on Name (1)
real*4 logM_ (nr__1) ! ([Msun]) [-2.1/0.9] Log of stellar mass from DM97 (2)
real*4 b_logM_ (nr__1) ! ([Msun]) [-2.2/0.8] Lower value of logM* (3)
real*4 B_logM__1 (nr__1) ! ([Msun]) [-1.7/1.1] Upper value of logM* (3)
real*4 logt_ (nr__1) ! ([Myr]) [-1.6/1.7] Log of stellar age from DM97 (2)
real*4 b_logt_ (nr__1) ! ([Myr]) [-2.5/1.1] Lower value of logt* (3)
real*4 B_logt__1 (nr__1) ! ([Myr]) [-1.3/1.9] Upper value of logt* (3)
real*4 logM_2 (nr__1) ! ([Msun]) [-2/1.5] Log of stellar mass from BCAH98 (2)
real*4 b_logM_2 (nr__1) ! ([Msun]) [-2.5/1.1] Lower value of logM*2 (3)
real*4 B_logM_2_1 (nr__1) ! ([Msun]) [-1.7/1.7] Upper value of logM*2 (3)
real*4 logt_2 (nr__1) ! ([Myr]) [-2.6/1.6] Log of stellar age from BCAH98 (2)
real*4 b_logt_2 (nr__1) ! ([Myr]) [-2.8/1.4] Lower value of logt*2 (3)
real*4 B_logt_2_1 (nr__1) ! ([Myr]) [-0.9/1.9] Lower value of logt*2 (3)
real*4 logM_3 (nr__1) ! ([Msun]) [-2.1/0.9] Log of stellar mass from SDF00 (2)
real*4 b_logM_3 (nr__1) ! ([Msun]) [-2.4/0.8] Lower value of logM*3 (3)
real*4 B_logM_3_1 (nr__1) ! ([Msun]) [-1.6/1] Upper value of logM*3 (3)
real*4 logt_3 (nr__1) ! ([Myr]) [-2.1/1.7] Log of stellar age from SDF00 (2)
real*4 b_logt_3 (nr__1) ! ([Myr]) [-2.4/1.1] Lower value of logt*3 (3)
real*4 B_logt_3_1 (nr__1) ! ([Myr]) [-0.6/1.9] Upper value of logt*3 (3)
*Note (1): Flag as follows:
* a = The logM* values (and their confidence intervals) in these composite
* multiple systems were determined by properly combining the marginal
* probability density functions determined for their individual components,
* as described in the text. The corresponding logt* values (and their
* confidence intervals) correspond to the primary components.
* b = These systems are either known or strongly suspected to be associated
* with nearly edge-on disks. Since such a viewing geometry substantially
* obscures the starlight, the derived estimates of L* for these cases are
* not used to determine masses and ages. Instead, we compute the weighted
* mean and standard deviation of the luminosities for the ensemble of
* sources with the same spectral type (+/-1 subclass), and associate them
* with L* and {sigma}_L_ in the Bayesian inference method for
* {M*, t*}: details are available in Appendix B.
*Note (2): References as below:
* DM97 = D'Antona & Mazzitelli 1997MmSAI..68..807D
* BCAH98 = Baraffe et al. 1998, J/A+A/337/403
* SDF00 = Siess et al. 2000A&A...358..593S
*Note (3): Confidence intervals represent the extreme values of the variable that
* span its 68.2% ("1{sigma}") confidence interval, derived from its
* associated marginal probability density function (see Section 3.2.1).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! Stellar spectrum fitting results
integer*4 nr__2
parameter (nr__2=224) ! Number of records
character*62 ar__2 ! Full-size record
character*20 Name_2 (nr__2) ! Name
character*3 f_Name_2 (nr__2) ! [abc, ] Flag on Name (1)
real*4 logT__1 (nr__2) ! ([K]) [3.3/4.1] log of stellar temperature
real*4 e_logT_ (nr__2) ! ([K]) [0.01/0.2] logT* uncertainty
real*4 logL_ (nr__2) ! ([Lsun]) [-3/2.2] log of stellar luminosity
real*4 e_logL_ (nr__2) ! ([Lsun]) [0.03/0.9] logL* uncertainty
real*4 Av (nr__2) ! (mag) [0/28.2] Extinction
real*4 e_Av (nr__2) ! (mag) [0.2/3.3] Av uncertainty
*Note (1): Flag as follows:
* a = The stellar parameters for these close pairs were determined by
* simultaneous fits to both their composite SEDs and any available
* contrast ratios, assuming that both components have the same extinction.
* b = These stars are known or suspected to host edge-on disks. The
* luminosities and extinctions reported here are the results of the formal
* measurements described in Appendix B. However, their L* values are
* modified before use in Section 3.2, as described in the text.
* c = The components V807 Tau Bab and V710 Tau B are included here for
* completeness, but are not considered part of the Class II sample. The
* stellar photosphere parameters listed here were used to subtract off
* their contribution to the composite SEDs shown in Appendix B.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'refs.dat' ! References
integer*4 nr__3
parameter (nr__3=21) ! Number of records
character*263 ar__3 ! Full-size record
character*3 Ref (nr__3) ! Reference code
character*19 BibCode (nr__3) ! Bibcode
character*19 Aut (nr__3) ! Author's name(s)
character*219 Comm (nr__3) ! Comment
C=============================================================================
C Loading file 'table2.dat' ! Millimeter-wave continuum photometry catalog
C Format for file interpretation
1 format(
+ A20,A1,1X,A16,1X,A6,1X,A1,1X,F6.4,1X,F6.4,1X,A1,1X,F6.4,1X,
+ F6.4,1X,A6)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,210
read(1,'(A84)')ar__
read(ar__,1)
+ Name(i__),f_Name(i__),SpT(i__),r_SpT(i__),l_F0_89(i__),
+ F0_89(i__),e_F0_89(i__),l_F1_3(i__),F1_3(i__),e_F1_3(i__),
+ Notes(i__)
if(ar__(56:61) .EQ. '') e_F0_89(i__) = rNULL__
if(ar__(72:77) .EQ. '') e_F1_3(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ Name(i__),f_Name(i__),SpT(i__),r_SpT(i__),l_F0_89(i__),
+ F0_89(i__),e_F0_89(i__),l_F1_3(i__),F1_3(i__),e_F1_3(i__),
+ Notes(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! Stellar host masses and ages
C Format for file interpretation
2 format(
+ A20,1X,A1,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,
+ F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,
+ 1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,212
read(1,'(A130)')ar__1
read(ar__1,2)
+ Name_1(i__),f_Name_1(i__),logM_(i__),b_logM_(i__),
+ B_logM__1(i__),logt_(i__),b_logt_(i__),B_logt__1(i__),
+ logM_2(i__),b_logM_2(i__),B_logM_2_1(i__),logt_2(i__),
+ b_logt_2(i__),B_logt_2_1(i__),logM_3(i__),b_logM_3(i__),
+ B_logM_3_1(i__),logt_3(i__),b_logt_3(i__),B_logt_3_1(i__)
c ..............Just test output...........
write(6,2)
+ Name_1(i__),f_Name_1(i__),logM_(i__),b_logM_(i__),
+ B_logM__1(i__),logt_(i__),b_logt_(i__),B_logt__1(i__),
+ logM_2(i__),b_logM_2(i__),B_logM_2_1(i__),logt_2(i__),
+ b_logt_2(i__),B_logt_2_1(i__),logM_3(i__),b_logM_3(i__),
+ B_logM_3_1(i__),logt_3(i__),b_logt_3(i__),B_logt_3_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! Stellar spectrum fitting results
C Format for file interpretation
3 format(
+ A20,1X,A3,1X,F6.4,1X,F6.4,1X,F6.3,1X,F5.3,1X,F5.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,224
read(1,'(A62)')ar__2
read(ar__2,3)
+ Name_2(i__),f_Name_2(i__),logT__1(i__),e_logT_(i__),
+ logL_(i__),e_logL_(i__),Av(i__),e_Av(i__)
c ..............Just test output...........
write(6,3)
+ Name_2(i__),f_Name_2(i__),logT__1(i__),e_logT_(i__),
+ logL_(i__),e_logL_(i__),Av(i__),e_Av(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'refs.dat' ! References
C Format for file interpretation
4 format(A3,1X,A19,1X,A19,1X,A219)
C Effective file loading
open(unit=1,status='old',file=
+'refs.dat')
write(6,*) '....Loading file: refs.dat'
do i__=1,21
read(1,'(A263)')ar__3
read(ar__3,4)Ref(i__),BibCode(i__),Aut(i__),Comm(i__)
c ..............Just test output...........
write(6,4)Ref(i__),BibCode(i__),Aut(i__),Comm(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end