Conversion of standardized ReadMe file for
file /./ftp/cats/J/AJ/145/66 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-09
C=============================================================================
* This code was generated from the ReadMe file documenting a catalogue
* according to the "Standard for Documentation of Astronomical Catalogues"
* currently in use by the Astronomical Data Centers (CDS, ADC, A&A)
* (see full documentation at URL http://vizier.u-strasbg.fr/doc/catstd.htx)
* Please report problems or questions to
C=============================================================================
implicit none
* Unspecified or NULL values, generally corresponding to blank columns,
* are assigned one of the following special values:
* rNULL__ for unknown or NULL floating-point values
* iNULL__ for unknown or NULL integer values
real*4 rNULL__
integer*4 iNULL__
parameter (rNULL__=--2147483648.) ! NULL real number
parameter (iNULL__=(-2147483647-1)) ! NULL int number
integer idig ! testing NULL number
C=============================================================================
Cat. J/AJ/145/66 Spitzer light curves of YSOs in IC 348 (Flaherty+, 2013)
*================================================================================
*Kinks and dents in protoplanetary disks: rapid infrared variability as evidence
*for large structural perturbations.
* Flaherty K.M., Muzerolle J., Rieke G., Gutermuth R., Balog Z., Herbst W.,
* Megeath S.T.
* <Astron. J., 145, 66 (2013)>
* =2013AJ....145...66F
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table1.dat' ! Cluster member properties
integer*4 nr__
parameter (nr__=377) ! Number of records
character*114 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 LRL (nr__) ! [1/57025] YSO identifier (G1)
character*1 m_LRL (nr__) ! [AB] Multiplicity flag on LRL (A or B)
real*8 RAdeg (nr__) ! (deg) Right Ascension in decimal degrees (J2000)
real*8 DEdeg (nr__) ! (deg) Declination in decimal degrees (J2000)
integer*4 Teff (nr__) ! (K) ?=-100.000 Effective temperature
real*8 Lum (nr__) ! (Lsun) ?=-100.000 Stellar luminosity
real*4 W_Ha (nr__) ! (0.1nm) ?=-10.0 H{alpha} equivalent width; in {AA} (1)
real*8 Per (nr__) ! (d) ?=-10.000 Optical period (2)
real*4 AV (nr__) ! (mag) The V-band extinction
real*4 LX (nr__) ! (ZW) ?=-10.0 X-ray luminosity; in 10^28^erg/s (3)
real*8 E3_6 (nr__) ! (mag) ?=-10.000 Excess emission at 3.6{mu}m
real*8 alpha (nr__) ! ?=-100.00 Slope of the observed infrared SED (4)
real*8 alpha2 (nr__) ! ?=-100.00 Slope of the dereddened infrared
* SED (5)
real*4 Dist (nr__) ! (arcmin) Distance from HH 211 (6)
character*2 f_LRL (nr__) ! Variable? (y, n, y1, n1, p, or b) (7)
character*11 Ref (nr__) ! Reference (8)
*Note (1): Here we actually report the negative of the H{alpha} equivalent width,
* with positive numbers corresponding to emission. H{alpha} equivalent width
* is not reported when it is in absorption.
*Note (2): Optical period is taken from Cohen et al. (2004, cat. J/AJ/127/1602),
* Cieza & Baliber (2006, cat. J/ApJ/649/862) and references therein. Values
* of 0 correspond to stars that were detected in these surveys, but no period
* structure was detected.
*Note (3): From Preibisch & Zinnecker (2002AJ....123.1613P).
*Note (4): Defined as {lambda}F_{lambda}_{propto}{lambda}^{alpha}^ over the
* IRAC (3-8{mu}m) wavelengths.
*Note (5): Unless otherwise stated, this is the value of the infrared SED that
* is used in the text.
*Note (6): IC 348 is composed of an exposed cluster core in the northeast and a
* highly obscured region in the southwest where most of the younger class I
* sources, along with all of the molecular outflows and submillimeter cores,
* reside. The presence of these less-evolved sources indicates that the
* southwestern ridge is a region of very recent star formation. We
* characterize the position as the distance from HH 211, a prominent
* Herbig-Haro object in this region.
*Note (7): Multiplicity flag on Name as follows:
* y = Variable;
* n = Not variable;
* y1 = Variable, but only detected at [3.6];
* n1 = Not variable, but only detected at [3.6];
* p = Periodic;
* b = A nearby source prevents us from accurately assessing variability;
* blank = Source was not detected.
*Note (8): Reference as follows:
* luh03 = Luhman et al. (2003, cat. J/ApJ/593/1093);
* mue07 = Muench et al. (2007, cat. J/AJ/134/411);
* lad06 = Lada et al. (2006, cat. J/AJ/131/1574);
* muz06 = Muzerolle et al. (2006ApJ...643.1003M);
* luh05 = Luhman et al. (2005ApJ...618..810L).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table2.dat' ! [3.6] photometry
integer*4 nr__1
parameter (nr__1=14326) ! Number of records
character*47 ar__1 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 LRL_1 (nr__1) ! [1/57025] YSO identifier (G1)
character*1 m_LRL_1 (nr__1) ! [AB] Multiplicity flag on LRL (A or B)
real*8 RAdeg_1 (nr__1) ! (deg) Right Ascension in decimal degrees (J2000)
real*8 DEdeg_1 (nr__1) ! (deg) Declination in decimal degrees (J2000)
real*8 MJD (nr__1) ! (d) [55106.9/55142.1] Modified Julian Date
real*8 v_3_6_ (nr__1) ! (mag) ?=-100.000 Spitzer/IRAC 3.6{mu}m band magnitude
real*8 e__3_6_ (nr__1) ! (mag) ?=10.0000 Uncertainty in [3.6]
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! [4.5] photometry
integer*4 nr__2
parameter (nr__2=14326) ! Number of records
character*47 ar__2 ! Full-size record
C J2000 position composed of: RAdeg DEdeg
integer*4 LRL_2 (nr__2) ! [1/57025] YSO identifier (G1)
character*1 m_LRL_2 (nr__2) ! [AB] Multiplicity flag on LRL (A or B)
real*8 RAdeg_2 (nr__2) ! (deg) Right Ascension in decimal degrees (J2000)
real*8 DEdeg_2 (nr__2) ! (deg) Declination in decimal degrees (J2000)
real*8 MJD_1 (nr__2) ! (d) [55106.9/55142.1] Modified Julian Date
real*8 v_4_5_ (nr__2) ! (mag) ?=-100.000 Spitzer/IRAC 4.5{mu}m band magnitude
real*8 e__4_5_ (nr__2) ! (mag) ?=10.0000 Uncertainty in [4.5]
C=============================================================================
C Loading file 'table1.dat' ! Cluster member properties
C Format for file interpretation
1 format(
+ I5,A1,1X,F7.4,1X,F7.4,1X,I5,1X,F8.3,1X,F5.1,3X,F7.3,1X,F6.3,
+ 1X,F6.1,3X,F7.3,1X,F7.2,1X,F7.2,1X,F5.2,1X,A2,1X,A11)
C Effective file loading
open(unit=1,status='old',file=
+'table1.dat')
write(6,*) '....Loading file: table1.dat'
do i__=1,377
read(1,'(A114)')ar__
read(ar__,1)
+ LRL(i__),m_LRL(i__),RAdeg(i__),DEdeg(i__),Teff(i__),Lum(i__),
+ W_Ha(i__),Per(i__),AV(i__),LX(i__),E3_6(i__),alpha(i__),
+ alpha2(i__),Dist(i__),f_LRL(i__),Ref(i__)
c ..............Just test output...........
write(6,1)
+ LRL(i__),m_LRL(i__),RAdeg(i__),DEdeg(i__),Teff(i__),Lum(i__),
+ W_Ha(i__),Per(i__),AV(i__),LX(i__),E3_6(i__),alpha(i__),
+ alpha2(i__),Dist(i__),f_LRL(i__),Ref(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table2.dat' ! [3.6] photometry
C Format for file interpretation
2 format(I5,A1,1X,F7.4,1X,F7.4,1X,F7.1,1X,F8.3,1X,F7.4)
C Effective file loading
open(unit=1,status='old',file=
+'table2.dat')
write(6,*) '....Loading file: table2.dat'
do i__=1,14326
read(1,'(A47)')ar__1
read(ar__1,2)
+ LRL_1(i__),m_LRL_1(i__),RAdeg_1(i__),DEdeg_1(i__),MJD(i__),
+ v_3_6_(i__),e__3_6_(i__)
c ..............Just test output...........
write(6,2)
+ LRL_1(i__),m_LRL_1(i__),RAdeg_1(i__),DEdeg_1(i__),MJD(i__),
+ v_3_6_(i__),e__3_6_(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table3.dat' ! [4.5] photometry
C Format for file interpretation
3 format(I5,A1,1X,F7.4,1X,F7.4,1X,F7.1,1X,F8.3,1X,F7.4)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,14326
read(1,'(A47)')ar__2
read(ar__2,3)
+ LRL_2(i__),m_LRL_2(i__),RAdeg_2(i__),DEdeg_2(i__),MJD_1(i__),
+ v_4_5_(i__),e__4_5_(i__)
c ..............Just test output...........
write(6,3)
+ LRL_2(i__),m_LRL_2(i__),RAdeg_2(i__),DEdeg_2(i__),MJD_1(i__),
+ v_4_5_(i__),e__4_5_(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end