Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/688/418 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 2024-Mar-29 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/688/418 X-ray flares in Orion young stars. I. (Getman+, 2008) *================================================================================ *X-ray flares in Orion young stars. *I. Flare characteristics. * Getman K.V., Feigelson E.D., Broos P.S., Micela G., Garmire G.P. * <Astrophys. J., 688, 418-436 (2008)> * =2008ApJ...688..418G C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table1.dat' ! Flare properties 1 integer*4 nr__ parameter (nr__=216) ! Number of records character*103 ar__ ! Full-size record integer*4 COUP ! COUP source number integer*4 FNo ! Flare number integer*4 NCf ! (ct) Net flare counts real*4 bgf ! (ct) Background counts in flare real*4 MEf ! (keV) X-ray median energy of flare counts real*4 tchar1 ! (ks) Start time of characteristic state real*4 tchar2 ! (ks) Stop time of characteristic state real*4 tFl1 ! (ks) Flare start time real*4 tFl2 ! (ks) Flare stop time real*4 CRchar ! (ct/ks) Count rate in characteristic state real*4 CRFlpk ! (ct/ks) Peak flare count rate real*4 e_CRFlpk ! (ct/ks) Error in CRFlpk real*4 tCRpk ! (ks) Time of maximum flare count rate real*4 tTpk ! (ks) Time of maximum flare temperature real*4 tEMpk ! (ks) Time of maximum flare emission measure real*4 texp ! (ks) Time point for pure exponential decay (1) real*4 tsep ! (ks) Separation time point (2) *Note (1): Time point until which the lightcurve exhibits a pure exponential * decay. *Note (2): Separation time point between two distinct decay phases. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table2.dat' ! Flare properties 2 integer*4 nr__1 parameter (nr__1=216) ! Number of records character*96 ar__1 ! Full-size record integer*4 COUP_1 ! COUP source number integer*4 FNo_1 ! Flare number real*4 taur ! (ks) Rise e-folding timescale real*4 taud1 ! (ks) Initial decay e-folding timescale real*4 taud2 ! (ks) Full decay e-folding timescale real*4 Lxpk ! ([10-7W]) Flare peak X-ray luminosity (log scale) real*4 Tobspk ! (MK) Flare peak observed plasma temperature real*4 EMpk ! (10+53cm-3) Flare peak emission measure real*4 TEMpk_1 ! (MK) Plasma temperature at maximum emission measure real*4 zeta1 ! Slope over the full decay real*4 zeta2 ! Slope within the initial decay real*4 zeta3 ! Slope within the late decay real*4 L1 ! (10+10cm) Low boundary of the inferred loop length range real*4 L2 ! (10+10cm) High boundary of the inferred loop length range real*8 Llim ! (10+10cm) Limiting maximum loop length using Serio's law integer*4 Flag1 ! Flag for a dominant slope(s) (1) character*3 Type ! Flare morphology (2) *Note (1): Flag defines the best or dominant slope or a combination of slopes * which, in combination with decay timescales taud1 and taud2, are * used to estimate loop size range; as follows: * 01 = {zeta}1; * 02 = {zeta}2; * 03 = {zeta}3; * 12 = {zeta}1, {zeta}2; * 13 = {zeta}1, {zeta}3; * 23 = {zeta}2, {zeta}3; * 44 = unable to estimate loop size. *Note (2): The following morphological classes are defined: * typ = typical; * stp = flares with step decay; * srf = slow-rise-and/or-top-flat flares; * dbl = double flares; * inc = incomplete, interrupted by COUP observational gaps; * oth = other, purely defined shapes generally due to their lower * counting statistics. c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table3.dat' ! COUP star properties integer*4 nr__2 parameter (nr__2=161) ! Number of records character*116 ar__2 ! Full-size record integer*4 COUP_2 ! COUP source number (1) integer*4 NC ! (ct) Source net counts (1) integer*4 bg ! (ct) Background counts in source region (1) real*4 ME ! (keV) Source's X-ray median energy (1) real*4 NH ! ([cm-2]) Source's inferred column density (1) real*4 Jmag ! (mag) ?=-9.99 2MASS J-band magnitude (1) real*4 Hmag ! (mag) 2MASS H-band magnitude (1) real*4 Ksmag ! (mag) 2MASS Ks-band magnitude (1) character*6 Qual ! IR source and quality flag (1) real*4 d_H_K ! (mag) ?=-9.99 IR color excess (2) real*4 v3_6mag ! (mag) ?=-9.999 IR 3.6 micron band magnitude real*4 e_3_6mag ! (mag) ?=-9.999 Error in 3.6mag real*4 v4_5mag ! (mag) ?=-9.999 IR 4.5 micron band magnitude real*4 e_4_5mag ! (mag) ?=-9.999 Error in 4.5mag real*4 Mass ! (Msun) ?=-9.99 Stellar mass (1) real*4 Rad ! (10+8m) ?=-9.99 Stellar radius (1) real*4 Per ! (d) ?=-9.99 Stellar rotation period (3) real*4 Rcor ! (10+8m) ?=-9.99 Keplerian co-rotation radius real*4 EW_CaII ! (0.1nm) ?=-99.99 Equivalent width of Ca II line (1) *Note (1): Values for these quantities are from Getman et al. (2005, Cat. * J/ApJS/160/353): 3 letters for 2MASS quality flags and 3 characters * for the 2MASS contamination and confusion flag. *Note (2): Color excess measured from the rightmost reddening vector on the * J-H vs. H-Ks color-color diagram (originating at 0.1M_{sun}_ assuming * an age of 1Myr and using PMS models of Siess et al. * (2000A&A...358..593S)). *Note (3): Stellar rotation periods with extended rotation data collected in * Flaccomio et al. (2005, Cat. J/ApJS/160/450). c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table4.dat' ! Time-resolved spectroscopy analysis integer*4 nr__3 parameter (nr__3=216) ! Number of records character*55 ar__3 ! Full-size record integer*4 COUP_3 ! COUP source number integer*4 FNo_2 ! Flare number character*3 Type_1 ! Flare morphology (1) integer*4 SNo ! Number of time-resolved spectroscopy (TRS) * segments integer*4 Nseg ! (ct) Average number of counts in TRS segment integer*4 Nchar ! (ct) Number of counts in characteristic segment real*4 Tchar ! (MK) Characteristic plasma temperature real*4 TpknCH ! (MK) Peak flare plasma temperature in "noCH" case (2) real*4 TpkCH ! (MK) Peak flare plasma temperature in "CH" case (2) real*4 znCH ! Full decay slope in "noCH" case (2) real*4 zCH ! Full decay slope in "CH" case (2) *Note (1): The following morphological classes are defined: * typ = typical; * stp = flares with step decay; * srf = slow-rise-and/or-top-flat flares; * dbl = double flares; * inc = incomplete, interrupted by COUP observational gaps; * oth = other, purely defined shapes generally due to their lower * counting statistics. *Note (2): Inferred peak flare plasma temperatures and slopes of temperature- * density trajectories without ("noCH") and with ("CH") account for * the "characteristic" background. C============================================================================= C Loading file 'table1.dat' ! Flare properties 1 C Format for file interpretation 1 format( + I4,1X,I1,1X,I5,1X,F5.1,1X,F4.2,1X,F6.1,1X,F6.1,1X,F6.1,1X, + F6.1,1X,F4.1,1X,F6.1,1X,F4.1,1X,F6.1,1X,F6.1,1X,F6.1,1X,F6.1, + 1X,F6.1) C Effective file loading open(unit=1,status='old',file= +'table1.dat') write(6,*) '....Loading file: table1.dat' do i__=1,216 read(1,'(A103)')ar__ read(ar__,1) + COUP,FNo,NCf,bgf,MEf,tchar1,tchar2,tFl1,tFl2,CRchar,CRFlpk, + e_CRFlpk,tCRpk,tTpk,tEMpk,texp,tsep c ..............Just test output........... write(6,1) + COUP,FNo,NCf,bgf,MEf,tchar1,tchar2,tFl1,tFl2,CRchar,CRFlpk, + e_CRFlpk,tCRpk,tTpk,tEMpk,texp,tsep c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table2.dat' ! Flare properties 2 C Format for file interpretation 2 format( + I4,1X,I1,1X,F5.1,1X,F5.1,1X,F5.1,1X,F5.2,1X,F6.2,1X,F6.2,1X, + F6.2,1X,F4.1,1X,F4.1,1X,F5.1,1X,F6.2,1X,F6.2,1X,F7.2,1X,I2,1X, + A3) C Effective file loading open(unit=1,status='old',file= +'table2.dat') write(6,*) '....Loading file: table2.dat' do i__=1,216 read(1,'(A96)')ar__1 read(ar__1,2) + COUP_1,FNo_1,taur,taud1,taud2,Lxpk,Tobspk,EMpk,TEMpk_1,zeta1, + zeta2,zeta3,L1,L2,Llim,Flag1,Type c ..............Just test output........... write(6,2) + COUP_1,FNo_1,taur,taud1,taud2,Lxpk,Tobspk,EMpk,TEMpk_1,zeta1, + zeta2,zeta3,L1,L2,Llim,Flag1,Type c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table3.dat' ! COUP star properties C Format for file interpretation 3 format( + I4,1X,I5,1X,I3,1X,F4.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,A6, + 1X,F5.2,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F5.2,1X,F5.2,1X, + F5.2,1X,F6.2,1X,F6.2) C Effective file loading open(unit=1,status='old',file= +'table3.dat') write(6,*) '....Loading file: table3.dat' do i__=1,161 read(1,'(A116)')ar__2 read(ar__2,3) + COUP_2,NC,bg,ME,NH,Jmag,Hmag,Ksmag,Qual,d_H_K,v3_6mag, + e_3_6mag,v4_5mag,e_4_5mag,Mass,Rad,Per,Rcor,EW_CaII c ..............Just test output........... write(6,3) + COUP_2,NC,bg,ME,NH,Jmag,Hmag,Ksmag,Qual,d_H_K,v3_6mag, + e_3_6mag,v4_5mag,e_4_5mag,Mass,Rad,Per,Rcor,EW_CaII c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table4.dat' ! Time-resolved spectroscopy analysis C Format for file interpretation 4 format( + I4,1X,I1,1X,A3,1X,I2,1X,I4,1X,I4,1X,F5.2,1X,F6.2,1X,F6.2,1X, + F5.2,1X,F5.2) C Effective file loading open(unit=1,status='old',file= +'table4.dat') write(6,*) '....Loading file: table4.dat' do i__=1,216 read(1,'(A55)')ar__3 read(ar__3,4) + COUP_3,FNo_2,Type_1,SNo,Nseg,Nchar,Tchar,TpknCH,TpkCH,znCH,zCH c ..............Just test output........... write(6,4) + COUP_3,FNo_2,Type_1,SNo,Nseg,Nchar,Tchar,TpknCH,TpkCH,znCH,zCH c .......End.of.Just test output........... end do close(1) C============================================================================= stop end