Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/559/A3 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-12
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/559/A3 Low-mass stars in Trumpler 37 (Sicilia-Aguilar+, 2013)
*================================================================================
*The low-mass stellar population in the young cluster Tr37:
*Disk evolution, accretion, and environment.
* Sicilia-Aguilar A., Kim J.S., Sobolev A., Getman K., Henning T., Fang M.
* <Astron. Astrophys. 559, A3 (2013)>
* =2013A&A...559A...3S
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea1.dat' ! Observations and spectroscopic information for
all observed objects
integer*4 nr__
parameter (nr__=583) ! Number of records
character*93 ar__ ! Full-size record
character*19 Name (nr__) ! Object name (G1)
integer*4 W_Ha (nr__) ! (0.1nm) ?=99 H{alpha} equivalent width (1)
real*4 W_Li (nr__) ! (0.1nm) ?=9.9 LiI equivalent width (2)
character*1 u_W_Li (nr__) ! [u] Uncertainty label for the W(Li)
character*1 n_SpT (nr__) ! [u/b] Note on spectral type: uncertain/best
character*4 SpT (nr__) ! Spectral type (4)
real*4 AV (nr__) ! (mag) ?=9.9 Extinction in V magnitude
real*4 e_AV (nr__) ! (mag) ?=9.9 Error in the V magnitude extinction
character*6 Disk (nr__) ! Presence of disk/IR excess (5)
character*8 Campgn (nr__) ! Observing campaign and telescope
* (2009-SAO, 2009-MMT, or 2010-MMT)
character*6 Memb (nr__) ! Membership label and references (6)
character*23 Com (nr__) ! Individual comments (7)
*Note (1): in Angstroms; 99 if H{alpha} appears in absorption
*Note (2): in Angstroms;
* 0.0 if Lithium absorption can be excluded (W(Li)<0.1{AA});
* 9.9 if Lithium could not be measured (due to low S/N).
*Note (4): The spectral types are self-explanatory, except!
* E = spectral type earlier than G
* Late = late M-type
* ClsI = Class I object
* X = Uncertain
* The errors in the spectral type are of 0.5 subtypes for objects K5 or later,
* 1 subtype for early K stars, and 2-3 subtypes for G-type objects.
*Note (5): Presence of disk/IR excess code (the number in brackets indicates
* the wavelength at which the excess is observed, especially in case of
* uncertain excess seen at a single wavelength):
* Y = the object has IR excess/disk
* N = the object does not have IR excess/disk
* P = the object has probably IR excess/disk
* pN = the object has probably no IR excess/no disk
*Note (6): Membership label and references as follows:
* Y = member
* N = non-member
* P = probable member
* pN = probable non-member
* X = uncertain membership
* 1 = object previously identified by Sicilia-Aguilar et al.
* (2005, Cat. J/AJ/130/188, 2006, Cat. J/ApJ/638/897, Cat. J/AJ/132/2135)
* 2 = object previously identified by Barentsen et al.
* (2011, Cat. J/MNRAS/415/103)
* 3 = object previously identified by Morales-Calderon et al.
* (2009, Cat. J/ApJ/702/1507)
* 4 = object previously identified in X-ray by Getman et al.
* (2012, Cat. J/MNRAS/426/2917)
* 5 = object previously identified from IR excess by Getman et al.
* (2012, Cat. J/MNRAS/426/2917)
*Note (7): Abbreviations used in comments:
* Neb. = Nebular IR emission
* Nobj. = Nearby object
* Var. = SED suggestive of variability
* oVVI = Object appears old in the V vs V-I diagram
* dMe = dMe star
* UXor = Colors of the object suggestive of UXor variability
* Globu = Object associated to the IC1396A globule
* Mini-cluster = Object associated to one of the Mini-clusters in Tr37
* Any other comments are self-explanatory, other names of the object
* (if any) are also listed in the comment field.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea2.dat' ! Optical and 2MASS data for confirmed members and
probable members
integer*4 nr__1
parameter (nr__1=213) ! Number of records
character*128 ar__1 ! Full-size record
character*19 Name_1 (nr__1) ! Object name (G1)
real*4 Vmag (nr__1) ! (mag) ?=99.999 Optical V magnitude, Johnsons system
real*4 e_Vmag (nr__1) ! (mag) ?=9.999 Error in the optical V magnitude
character*2 n_Vmag (nr__1) ! Label for the V magnitude (1)
real*4 Rmag (nr__1) ! (mag) ?=99.999 Optical R magnitude, Johnsons system
real*4 e_Rmag (nr__1) ! (mag) ?=9.999 Error in the optical R magnitude
character*1 n_Rmag (nr__1) ! Label for the R magnitude (1)
real*4 Imag (nr__1) ! (mag) ?=99.999 Optical I magnitude, Johnsons system
real*4 e_Imag (nr__1) ! (mag) ?=9.999 Error in the optical I magnitude
character*1 n_Imag (nr__1) ! Label for the I magnitude (1)
character*16 v2MASS (nr__1) ! 2MASS identification name
* (Cutri et al., Cat. II/246)
real*4 Jmag (nr__1) ! (mag) ?=99.999 J magnitude from 2MASS
real*4 e_Jmag (nr__1) ! (mag) ?=9.999 Error in the J magnitude from 2MASS
character*1 n_Jmag (nr__1) ! Label for the J magnitude (1)
real*4 Hmag (nr__1) ! (mag) ?=99.999 H magnitude from 2MASS
real*4 e_Hmag (nr__1) ! (mag) ?=9.999 Error in the H magnitude from 2MASS
character*1 n_Hmag (nr__1) ! Label for the H magnitude (1)
real*4 Ksmag (nr__1) ! (mag) ?=99.999 Ks magnitude from 2MASS
real*4 e_Ksmag (nr__1) ! (mag) ?=9.999 Error in the Ks magnitude from 2MASS
character*2 n_Ksmag (nr__1) ! Label for the Ks magnitude (1)
*Note (1): Labels as follows:
* b = Near bright star HD206267
* c = Near 21-840, likely contaminated
* d = Near other bright stars, probably contaminated
* u = Uncertain value
* l = Magnitude lower limit
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea3.dat' ! Spitzer IRAC and MIPS data for confirmed members
and probable members
integer*4 nr__2
parameter (nr__2=205) ! Number of records
character*149 ar__2 ! Full-size record
character*19 Name_2 (nr__2) ! Object name (G1)
integer*4 v_3_6_ (nr__2) ! (uJy) ?=9999 Spitzer/IRAC band 1 (3.6um) flux
integer*4 e__3_6_ (nr__2) ! (uJy) ?=999 Error in the ?=9999 IRAC band 1 flux
integer*4 v_4_5_ (nr__2) ! (uJy) ?=9999 Spitzer/IRAC band 2 (4.5um) flux
integer*4 e__4_5_ (nr__2) ! (uJy) ?=999 Error in the ?=9999 IRAC band 2 flux
integer*4 v_5_8_ (nr__2) ! (uJy) ?=9999 Spitzer/IRAC band 3 (5.8um) flux
integer*4 e__5_8_ (nr__2) ! (uJy) ?=999 Error in the ?=9999 IRAC band 3 flux
integer*4 v_8_0_ (nr__2) ! (uJy) ?=9999 Spitzer/IRAC band 4 (8.0um) flux
integer*4 e__8_0_ (nr__2) ! (uJy) ?=999 Error in the ?=9999 IRAC band 4 flux
character*1 n__8_0_ (nr__2) ! Label for the ?=9999 IRAC band 4 flux (1)
integer*4 v_24_ (nr__2) ! (uJy) ?=9999 Spitzer/MIPS 24um flux
integer*4 e__24_ (nr__2) ! (uJy) ?=999 Error in the MIPS 24um flux
character*1 n__24_ (nr__2) ! Label for the MIPS 24um flux (1)
character*9 SED (nr__2) ! SED type (2)
character*1 n_SED (nr__2) ! Label on the SED type (1)
character*54 Com_1 (nr__2) ! Comments on the sources (3)
*Note (1): Labels in magnitudes as follows:
* u = Uncertain
* l = Upper limit
* e = Edge of field
* o = Out of field
*Note (2): SED type as follows:
* ClsI = Class I object
* Full = Full disk
* PTD = Pre-transitional disk
* TD = Transitional disk
* Dep. = Dust-depleted disk
* ND = No disk
*Note (3): Objects listed as "marginal detection" at 24um have probably IR
* excesses but their photometry is compromised by low S/N or different
* sources of contamination.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea4.dat' ! *WISE data for confirmed members and probable members
integer*4 nr__3
parameter (nr__3=73) ! Number of records
character*77 ar__3 ! Full-size record
character*19 Name_3 (nr__3) ! Object name
real*4 v_3_4_ (nr__3) ! (mJy) Flux in WISE 1 band (3.4um)
real*4 e__3_4_ (nr__3) ! (mJy) Error in the WISE W1 band flux
character*1 n__3_4_ (nr__3) ! [a] Label for the WISE 1 band flux (1)
real*4 v_4_6_ (nr__3) ! (mJy) Flux in WISE 2 band (4.6um)
real*4 e__4_6_ (nr__3) ! (mJy) Error in the WISE W2 band flux
character*1 n__4_6_ (nr__3) ! [a] Label for the WISE 2 band flux (1)
real*8 v_12_ (nr__3) ! (mJy) Flux in WISE 3 band (12um)
real*4 e__12_ (nr__3) ! (mJy) Error in the WISE W3 band flux
character*1 n__12_ (nr__3) ! [a] Label for the WISE 3 band flux (1)
real*4 v_22_ (nr__3) ! (mJy) Flux in WISE 4 band (22um)
real*4 e__22_ (nr__3) ! (mJy) Error in the WISE W4 band flux
character*1 n__22_ (nr__3) ! [a] Label for the WISE 4 band flux (1)
*Note (1): Label as follows:
* a = Uncertain and/or contaminated, value not used to trace the SED
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea5.dat' ! Stars with multiple emission lines and their
identifications
integer*4 nr__4
parameter (nr__4=158) ! Number of records
character*57 ar__4 ! Full-size record
character*19 Name_4 (nr__4) ! Object name (G1)
character*12 Line (nr__4) ! Line name (1)
real*4 W0 (nr__4) ! (0.1nm) ?=9999 Line wavelength (laboratory)
real*4 Wobs (nr__4) ! (0.1nm) Line wavelength (observed)
real*4 EW (nr__4) ! (0.1nm) Equivalent width
character*3 u_EW (nr__4) ! [Unc ] Uncertainty flag (2)
*Note (1): The lines have been identified according to Sicilia-Aguilar et al.
* (2012A&A...544A..93S), the NIST database, Appenzeller et al.
* (1986A&AS...64...65A), and Hamann & Persson (1992ApJS...82..247H).
*Note (2): Unc: Uncertain detection of the line.
C=============================================================================
C Loading file 'tablea1.dat' ! Observations and spectroscopic information for
* all observed objects
C Format for file interpretation
1 format(
+ A19,1X,I4,1X,F3.1,1X,A1,1X,A1,A4,1X,F4.1,1X,F3.1,1X,A6,1X,A8,
+ 1X,A6,2X,A23)
C Effective file loading
open(unit=1,status='old',file=
+'tablea1.dat')
write(6,*) '....Loading file: tablea1.dat'
do i__=1,583
read(1,'(A93)')ar__
read(ar__,1)
+ Name(i__),W_Ha(i__),W_Li(i__),u_W_Li(i__),n_SpT(i__),SpT(i__),
+ AV(i__),e_AV(i__),Disk(i__),Campgn(i__),Memb(i__),Com(i__)
c ..............Just test output...........
write(6,1)
+ Name(i__),W_Ha(i__),W_Li(i__),u_W_Li(i__),n_SpT(i__),SpT(i__),
+ AV(i__),e_AV(i__),Disk(i__),Campgn(i__),Memb(i__),Com(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea2.dat' ! Optical and 2MASS data for confirmed members and
* probable members
C Format for file interpretation
2 format(
+ A19,1X,F6.3,1X,F5.3,1X,A2,1X,F6.3,1X,F5.3,1X,A1,1X,F6.3,1X,
+ F5.3,1X,A1,1X,A16,1X,F6.3,1X,F5.3,1X,A1,1X,F6.3,1X,F5.3,1X,A1,
+ 1X,F6.3,1X,F5.3,1X,A2)
C Effective file loading
open(unit=1,status='old',file=
+'tablea2.dat')
write(6,*) '....Loading file: tablea2.dat'
do i__=1,213
read(1,'(A128)')ar__1
read(ar__1,2)
+ Name_1(i__),Vmag(i__),e_Vmag(i__),n_Vmag(i__),Rmag(i__),
+ e_Rmag(i__),n_Rmag(i__),Imag(i__),e_Imag(i__),n_Imag(i__),
+ v2MASS(i__),Jmag(i__),e_Jmag(i__),n_Jmag(i__),Hmag(i__),
+ e_Hmag(i__),n_Hmag(i__),Ksmag(i__),e_Ksmag(i__),n_Ksmag(i__)
c ..............Just test output...........
write(6,2)
+ Name_1(i__),Vmag(i__),e_Vmag(i__),n_Vmag(i__),Rmag(i__),
+ e_Rmag(i__),n_Rmag(i__),Imag(i__),e_Imag(i__),n_Imag(i__),
+ v2MASS(i__),Jmag(i__),e_Jmag(i__),n_Jmag(i__),Hmag(i__),
+ e_Hmag(i__),n_Hmag(i__),Ksmag(i__),e_Ksmag(i__),n_Ksmag(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea3.dat' ! Spitzer IRAC and MIPS data for confirmed members
* and probable members
C Format for file interpretation
3 format(
+ A19,1X,I6,1X,I3,1X,I5,1X,I3,1X,I6,1X,I3,1X,I5,1X,I3,1X,A1,1X,
+ I7,1X,I5,1X,A1,2X,A9,1X,A1,3X,A54)
C Effective file loading
open(unit=1,status='old',file=
+'tablea3.dat')
write(6,*) '....Loading file: tablea3.dat'
do i__=1,205
read(1,'(A149)')ar__2
read(ar__2,3)
+ Name_2(i__),v_3_6_(i__),e__3_6_(i__),v_4_5_(i__),e__4_5_(i__),
+ v_5_8_(i__),e__5_8_(i__),v_8_0_(i__),e__8_0_(i__),
+ n__8_0_(i__),v_24_(i__),e__24_(i__),n__24_(i__),SED(i__),
+ n_SED(i__),Com_1(i__)
c ..............Just test output...........
write(6,3)
+ Name_2(i__),v_3_6_(i__),e__3_6_(i__),v_4_5_(i__),e__4_5_(i__),
+ v_5_8_(i__),e__5_8_(i__),v_8_0_(i__),e__8_0_(i__),
+ n__8_0_(i__),v_24_(i__),e__24_(i__),n__24_(i__),SED(i__),
+ n_SED(i__),Com_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea4.dat' ! *WISE data for confirmed members and probable members
C Format for file interpretation
4 format(
+ A19,1X,F6.2,1X,F4.2,1X,A1,1X,F6.2,1X,F4.2,1X,A1,1X,F7.2,1X,
+ F4.2,1X,A1,1X,F6.2,1X,F5.2,1X,A1)
C Effective file loading
open(unit=1,status='old',file=
+'tablea4.dat')
write(6,*) '....Loading file: tablea4.dat'
do i__=1,73
read(1,'(A77)')ar__3
read(ar__3,4)
+ Name_3(i__),v_3_4_(i__),e__3_4_(i__),n__3_4_(i__),v_4_6_(i__),
+ e__4_6_(i__),n__4_6_(i__),v_12_(i__),e__12_(i__),n__12_(i__),
+ v_22_(i__),e__22_(i__),n__22_(i__)
c ..............Just test output...........
write(6,4)
+ Name_3(i__),v_3_4_(i__),e__3_4_(i__),n__3_4_(i__),v_4_6_(i__),
+ e__4_6_(i__),n__4_6_(i__),v_12_(i__),e__12_(i__),n__12_(i__),
+ v_22_(i__),e__22_(i__),n__22_(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea5.dat' ! Stars with multiple emission lines and their
* identifications
C Format for file interpretation
5 format(A19,1X,A12,1X,F6.1,1X,F6.1,1X,F6.1,1X,A3)
C Effective file loading
open(unit=1,status='old',file=
+'tablea5.dat')
write(6,*) '....Loading file: tablea5.dat'
do i__=1,158
read(1,'(A57)')ar__4
read(ar__4,5)
+ Name_4(i__),Line(i__),W0(i__),Wobs(i__),EW(i__),u_EW(i__)
c ..............Just test output...........
write(6,5)
+ Name_4(i__),Line(i__),W0(i__),Wobs(i__),EW(i__),u_EW(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end