FORTRAN Generation
(/./ftp/cats/J/ApJ/688/418)

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