FORTRAN Generation
(/./ftp/cats/J/ApJ/715/671)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/715/671 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 2024-Mar-19
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/715/671           X-ray and IR emission from YSOs            (Wolk+, 2010)
*================================================================================
*X-ray and infrared emission from young stellar objects near LkH{alpha} 101.
*    Wolk S.J., Winston E., Bourke T.L., Gutermuth R., Megeath S.T.,
*    Spitzbart B.D., Osten R.
*   <Astrophys. J., 715, 671-695 (2010)>
*   =2010ApJ...715..671W
C=============================================================================

C  Internal variables

      integer*4 i__

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table1.dat'	! X-ray sources detected in the LkH{alpha} 101 field

      integer*4 nr__
      parameter (nr__=213)	! Number of records
      character*41 ar__   	! Full-size record

      integer*4     Seq        (nr__) ! [2/231] Running source number (skips indicate
*                                 the spurious detections)
      character*18  ID         (nr__) ! Source identification; (HHMMSS.ss+DDMMSS.s;
*                                 <LkHa101CXO JHHMMSS.ss+DDMMSS.s> in Simbad)
      real*4        Off        (nr__) ! (arcsec) Offset from axis
      integer*4     RCts       (nr__) ! (ct) Raw counts
      real*4        NCts       (nr__) ! (ct) Net counts

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table2.dat'	! 2MASS photometry of X-ray sources

      integer*4 nr__1
      parameter (nr__1=140)	! Number of records
      character*51 ar__1  	! Full-size record

      integer*4     Seq_1      (nr__1) ! [2/231] Source index number
      real*4        X_2M       (nr__1) ! (arcsec) X-ray/2MASS positional offset
      real*4        Off2M      (nr__1) ! (arcsec) Offset from axis
      real*4        Jmag       (nr__1) ! (mag) 2MASS J band magnitude
      real*4        e_Jmag     (nr__1) ! (mag) ? Uncertainty in Jmag
      real*4        Hmag       (nr__1) ! (mag) 2MASS H band magnitude
      real*4        e_Hmag     (nr__1) ! (mag) ? Uncertainty in Hmag
      real*4        Kmag       (nr__1) ! (mag) 2MASS K band magnitude
      real*4        e_Kmag     (nr__1) ! (mag) ? Uncertainty in Kmag
      character*3   Qual       (nr__1) ! 2MASS photometry quality codes

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table3.dat'	! Spitzer photometry of X-ray sources

      integer*4 nr__2
      parameter (nr__2=158)	! Number of records
      character*89 ar__2  	! Full-size record

      integer*4     Seq_2      (nr__2) ! [2/231] Source index number (1)
      real*4        X_Sptz     (nr__2) ! (arcsec) X-ray/Spitzer positional offset
      real*4        Off_1      (nr__2) ! (arcsec) Offset from axis
      real*4        v_3_6_     (nr__2) ! (mag) ? Spitzer/IRAC 3.6 micron band magnitude
      real*4        v_4_5_     (nr__2) ! (mag) ? Spitzer/IRAC 4.5 micron band magnitude
      real*4        v_5_8_     (nr__2) ! (mag) ? Spitzer/IRAC 5.8 micron band magnitude
      real*4        v_8_0_     (nr__2) ! (mag) ? Spitzer/IRAC 8.0 micron band magnitude
      real*4        v_24_      (nr__2) ! (mag) ? Spitzer/MIPS 24 micron band magnitude
      real*4        e__3_6_    (nr__2) ! (mag) ? Uncertainty in [3.6]
      real*4        e__4_5_    (nr__2) ! (mag) ? Uncertainty in [4.5]
      real*4        e__5_8_    (nr__2) ! (mag) ? Uncertainty in [5.8]
      real*4        e__8_0_    (nr__2) ! (mag) ? Uncertainty in [8.0]
      real*4        e__24_     (nr__2) ! (mag) ? Uncertainty in [24]
      real*4        AK         (nr__2) ! (mag) The K band extinction
      character*14  Class      (nr__2) ! Source classification (1, 2, III, T(transition),
*                                gal (background galaxy), or LkH{alpha} 101)
      character*1   f_Class    (nr__2) ! [b] b: star within 90" of LkH{alpha} 101 (2)
*Note (1): In addition to their IAU approved Names:
*          <LkHa101CXO JHHMMSS.ss+DDMMSS.s> given in Table 1, these sources are
*          equally identified as: <LkHa101SST JHHMMSS.ss+DDMMSS.s>.
*Note (2): Stars with a notation "b" are within 90 arcsec of LkH{alpha} 101.
*          Their proximity, led to the conclusion that they were cluster
*          members, despite missing some IRAC data. For these stars 2MASS data
*          played an enhanced role in determining Class identification.

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table4.dat'	! LkH{alpha} 101 cluster spectral classification

      integer*4 nr__3
      parameter (nr__3=87)	! Number of records
      character*57 ar__3  	! Full-size record

      integer*4     Seq_3      (nr__3) ! [2/231]? Chandra identification number
      character*18  ID_1       (nr__3) ! LkHaSST identification (HHMMSS.ss+DDMMSS.s;
*                              <LkHa101SST JHHMMSS.ss+DDMMSS.s> in Simbad)
      character*4   SpT        (nr__3) ! MK spectral type
      real*4        e_SpT      (nr__3) ! SpT uncertainty
      integer*4     Teff       (nr__3) ! (K) Effective temperature
      real*4        logL       (nr__3) ! ([-]) Overall luminosity log(L_bol_/L_0_)
      real*8        EWHa       (nr__3) ! (0.1nm) H{alpha} (6563{AA}) equivalent width
      real*4        EWLi       (nr__3) ! (0.1nm) ? LiI (6707{AA}) equivalent width

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table7.dat'	! Spectral fits to X-ray sources

      integer*4 nr__4
      parameter (nr__4=88)	! Number of records
      character*64 ar__4  	! Full-size record

      integer*4     Seq_4      (nr__4) ! [2/231] Source index number
      real*4        NCts_1     (nr__4) ! (ct) Net counts
      real*4        NH         (nr__4) ! (10+22/cm2) Hydrogen column density
      real*4        e_NH       (nr__4) ! (10+22/cm2) ? The 1{sigma} uncertainty in NH
      character*1   l_kT       (nr__4) ! Limit flag on kT
      real*4        kT         (nr__4) ! (keV) Effective temperature
      real*4        e_kT       (nr__4) ! (keV) ? The 1{sigma} uncertainty in kT
      real*4        logF_A     (nr__4) ! ([mW/m2]) Log of Absorbed flux; erg/s/cm^2^
      real*4        logF_U     (nr__4) ! ([mW/m2]) Log of Unabsorbed flux; erg/s/cm^2^
      real*4        chi2       (nr__4) ! The reduced {chi}^2^ of the fit (1)
      character*9   Model      (nr__4) ! Model used in fit (2)
*Note (1): Included for "cstat" models for completeness. One model was
*          calculated with C statistics, it is unclear how this relates to
*          goodness of fit. The formal errors are unbiased however.
*Note (2): The model used for the fit is coded in two parts. The left hand side
*          indicates the type of model used. The right hand side indicates the
*          data set. The overall model is written: LHS_RHS.
*          LHS as follows:
*  cstat = C-statistics used for the fit with an initial guess of kT=1keV and
*          NH=10^21^cm^-2^.
*  crs   = c-stat used for initial guess with those input parameters and then
*          refit using Raymond-Smith plasma for final fit.
*  Mekal = again c-stat used for initial guess and then Mekal for final fit,
*  APEC  = cstat used for initial guess and then APEC for final fit.
*          RHS as follows:
*  29    = only data from first observation (Obsid = 5429) used in fit,
*  28    = only data from second observation (Obsid = 5428) used in fit,
*  all   = data from both observations used in fit.

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table8.dat'	! Spectral fits to X-ray sources (lower quality)

      integer*4 nr__5
      parameter (nr__5=30)	! Number of records
      character*64 ar__5  	! Full-size record

      integer*4     Seq_5      (nr__5) ! [2/231] Source index number
      real*4        NCts_2     (nr__5) ! (ct) Net counts
      real*4        NH_1       (nr__5) ! (10+22/cm2) Hydrogen column density
      real*4        e_NH_1     (nr__5) ! (10+22/cm2) ? The 1{sigma} uncertainty in NH
      character*1   l_kT_1     (nr__5) ! Limit flag on kT
      real*4        kT_1       (nr__5) ! (keV) Effective temperature
      real*4        e_kT_1     (nr__5) ! (keV) ? The 1{sigma} uncertainty in kT
      real*4        logF_A_1   (nr__5) ! ([mW/m2]) Log of Absorbed flux; erg/s/cm^2^
      real*4        logF_U_1   (nr__5) ! ([mW/m2]) Log of Unabsorbed flux; erg/s/cm^2^
      real*4        chi2_1     (nr__5) ! The reduced {chi}^2^ of the fit (1)
      character*9   Model_1    (nr__5) ! Model used in fit (2)
*Note (1): Included for "cstat" models for completeness. One model was
*          calculated with C statistics, it is unclear how this relates to
*          goodness of fit. The formal errors are unbiased however.
*Note (2): The model used for the fit is coded in two parts. The left hand side
*          indicates the type of model used. The right hand side indicates the
*          data set. The overall model is written: LHS_RHS.
*          LHS as follows:
*  cstat = C-statistics used for the fit with an initial guess of kT=1keV and
*          NH=10^21^cm^-2^.
*  crs   = c-stat used for initial guess with those input parameters and then
*          refit using Raymond-Smith plasma for final fit.
*  Mekal = again c-stat used for initial guess and then Mekal for final fit,
*  APEC  = cstat used for initial guess and then APEC for final fit.
*          RHS as follows:
*  29    = only data from first observation (Obsid = 5429) used in fit,
*  28    = only data from second observation (Obsid = 5428) used in fit,
*  all   = data from both observations used in fit.

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table10.dat'	! PMS stars detected by Spitzer but not detected in
                            X-rays (2MASS colors)

      integer*4 nr__6
      parameter (nr__6=69)	! Number of records
      character*55 ar__6  	! Full-size record

      character*18  ID_2       (nr__6) ! Source identification; (HHMMSS.ss+DDMMSS.s;
*                                 <LkHa101SST JHHMMSS.ss+DDMMSS.s> in Simbad) (1)
      real*4        Jmag_1     (nr__6) ! (mag) ? 2MASS J band magnitude
      real*4        e_Jmag_1   (nr__6) ! (mag) ? Uncertainty in Jmag
      real*4        Hmag_1     (nr__6) ! (mag) ? 2MASS H band magnitude
      real*4        e_Hmag_1   (nr__6) ! (mag) ? Uncertainty in Hmag
      real*4        Ksmag      (nr__6) ! (mag) ? 2MASS Ks band magnitude
      real*4        e_Ksmag    (nr__6) ! (mag) ? Uncertainty in Kmag
      character*3   Qual_1     (nr__6) ! 2MASS photometry quality codes
*Note (1): Source 042955.68+352021.2 has been modified at CDS into
*          042955.68+352021.0 as in table 11 and 2MASS catalog.

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table11.dat'	! PMS stars not detected in X-rays (Spitzer)

      integer*4 nr__7
      parameter (nr__7=69)	! Number of records
      character*81 ar__7  	! Full-size record

      character*18  ID_3       (nr__7) ! Source identification; (HHMMSS.ss+DDMMSS.s;
*                                  <LkHa101SST JHHMMSS.ss+DDMMSS.s> in Simbad)
      character*1   f_ID       (nr__7) ! [*] outside Chandra field of view (1)
      real*4        v_3_6__1   (nr__7) ! (mag) Spitzer/IRAC 3.6 micron band magnitude
      real*4        e__3_6__1  (nr__7) ! (mag) Uncertainty in [3.6]
      real*4        v_4_5__1   (nr__7) ! (mag) ? Spitzer/IRAC 4.5 micron band magnitude
      real*4        e__4_5__1  (nr__7) ! (mag) ? Uncertainty in [4.5]
      real*4        v_5_6_     (nr__7) ! (mag) ? Spitzer/IRAC 5.6 micron band magnitude
      real*4        e__5_6_    (nr__7) ! (mag) ? Uncertainty in [5.6]
      real*4        v_8_0__1   (nr__7) ! (mag) ? Spitzer/IRAC 8.0 micron band magnitude
      real*4        e__8_0__1  (nr__7) ! (mag) ? Uncertainty in [8.0]
      real*4        v_24__1    (nr__7) ! (mag) ? Spitzer/MIPS 24 micron band magnitude
      real*4        e__24__1   (nr__7) ! (mag) ? Uncertainty in [24]
      real*4        AK_1       (nr__7) ! (mag) The K band extinction
      character*1   Cl         (nr__7) ! [012T] Classification (Class 0, I, II or
*                                  T=possible transition objects)
*Note (1):
*    * = These Spitzer identified PMS stars are off the Chandra field of view.

C=============================================================================

C  Loading file 'table1.dat'	! X-ray sources detected in the LkH{alpha} 101 field

C  Format for file interpretation

    1 format(I3,1X,A18,1X,F6.1,1X,I4,1X,F6.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,213
        read(1,'(A41)')ar__
        read(ar__,1)Seq(i__),ID(i__),Off(i__),RCts(i__),NCts(i__)
c    ..............Just test output...........
        write(6,1)Seq(i__),ID(i__),Off(i__),RCts(i__),NCts(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table2.dat'	! 2MASS photometry of X-ray sources

C  Format for file interpretation

    2 format(
     +  I3,1X,F4.2,1X,F5.1,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,
     +  F4.2,1X,A3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,140
        read(1,'(A51)')ar__1
        read(ar__1,2)
     +  Seq_1(i__),X_2M(i__),Off2M(i__),Jmag(i__),e_Jmag(i__),
     +  Hmag(i__),e_Hmag(i__),Kmag(i__),e_Kmag(i__),Qual(i__)
        if(ar__1(22:25) .EQ. '') e_Jmag(i__) = rNULL__
        if(ar__1(33:36) .EQ. '') e_Hmag(i__) = rNULL__
        if(ar__1(44:47) .EQ. '') e_Kmag(i__) = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  Seq_1(i__),X_2M(i__),Off2M(i__),Jmag(i__),e_Jmag(i__),
     +  Hmag(i__),e_Hmag(i__),Kmag(i__),e_Kmag(i__),Qual(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table3.dat'	! Spitzer photometry of X-ray sources

C  Format for file interpretation

    3 format(
     +  I3,1X,F4.2,1X,F5.1,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F4.2,1X,
     +  F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,A14,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,158
        read(1,'(A89)')ar__2
        read(ar__2,3)
     +  Seq_2(i__),X_Sptz(i__),Off_1(i__),v_3_6_(i__),v_4_5_(i__),
     +  v_5_8_(i__),v_8_0_(i__),v_24_(i__),e__3_6_(i__),e__4_5_(i__),
     +  e__5_8_(i__),e__8_0_(i__),e__24_(i__),AK(i__),Class(i__),
     +  f_Class(i__)
        if(ar__2(16:20) .EQ. '') v_3_6_(i__) = rNULL__
        if(ar__2(22:26) .EQ. '') v_4_5_(i__) = rNULL__
        if(ar__2(28:32) .EQ. '') v_5_8_(i__) = rNULL__
        if(ar__2(34:38) .EQ. '') v_8_0_(i__) = rNULL__
        if(ar__2(40:43) .EQ. '') v_24_(i__) = rNULL__
        if(ar__2(45:48) .EQ. '') e__3_6_(i__) = rNULL__
        if(ar__2(50:53) .EQ. '') e__4_5_(i__) = rNULL__
        if(ar__2(55:58) .EQ. '') e__5_8_(i__) = rNULL__
        if(ar__2(60:63) .EQ. '') e__8_0_(i__) = rNULL__
        if(ar__2(65:68) .EQ. '') e__24_(i__) = rNULL__
c    ..............Just test output...........
        write(6,3)
     +  Seq_2(i__),X_Sptz(i__),Off_1(i__),v_3_6_(i__),v_4_5_(i__),
     +  v_5_8_(i__),v_8_0_(i__),v_24_(i__),e__3_6_(i__),e__4_5_(i__),
     +  e__5_8_(i__),e__8_0_(i__),e__24_(i__),AK(i__),Class(i__),
     +  f_Class(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table4.dat'	! LkH{alpha} 101 cluster spectral classification

C  Format for file interpretation

    4 format(I3,1X,A18,1X,A4,1X,F3.1,1X,I5,1X,F6.3,1X,F8.2,1X,F3.1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table4.dat')
      write(6,*) '....Loading file: table4.dat'
      do i__=1,87
        read(1,'(A57)')ar__3
        read(ar__3,4)
     +  Seq_3(i__),ID_1(i__),SpT(i__),e_SpT(i__),Teff(i__),logL(i__),
     +  EWHa(i__),EWLi(i__)
        if(ar__3(1:3) .EQ. '') Seq_3(i__) = iNULL__
        if(ar__3(55:57) .EQ. '') EWLi(i__) = rNULL__
c    ..............Just test output...........
        write(6,4)
     +  Seq_3(i__),ID_1(i__),SpT(i__),e_SpT(i__),Teff(i__),logL(i__),
     +  EWHa(i__),EWLi(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table7.dat'	! Spectral fits to X-ray sources

C  Format for file interpretation

    5 format(
     +  I3,1X,F6.1,1X,F5.2,1X,F4.2,1X,A1,1X,F6.2,1X,F4.2,1X,F6.2,1X,
     +  F6.2,1X,F4.2,1X,A9)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table7.dat')
      write(6,*) '....Loading file: table7.dat'
      do i__=1,88
        read(1,'(A64)')ar__4
        read(ar__4,5)
     +  Seq_4(i__),NCts_1(i__),NH(i__),e_NH(i__),l_kT(i__),kT(i__),
     +  e_kT(i__),logF_A(i__),logF_U(i__),chi2(i__),Model(i__)
        if(ar__4(18:21) .EQ. '') e_NH(i__) = rNULL__
        if(ar__4(32:35) .EQ. '') e_kT(i__) = rNULL__
c    ..............Just test output...........
        write(6,5)
     +  Seq_4(i__),NCts_1(i__),NH(i__),e_NH(i__),l_kT(i__),kT(i__),
     +  e_kT(i__),logF_A(i__),logF_U(i__),chi2(i__),Model(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table8.dat'	! Spectral fits to X-ray sources (lower quality)

C  Format for file interpretation

    6 format(
     +  I3,1X,F6.1,1X,F5.2,1X,F4.2,1X,A1,1X,F6.2,1X,F4.2,1X,F6.2,1X,
     +  F6.2,1X,F4.2,1X,A9)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table8.dat')
      write(6,*) '....Loading file: table8.dat'
      do i__=1,30
        read(1,'(A64)')ar__5
        read(ar__5,6)
     +  Seq_5(i__),NCts_2(i__),NH_1(i__),e_NH_1(i__),l_kT_1(i__),
     +  kT_1(i__),e_kT_1(i__),logF_A_1(i__),logF_U_1(i__),chi2_1(i__),
     +  Model_1(i__)
        if(ar__5(18:21) .EQ. '') e_NH_1(i__) = rNULL__
        if(ar__5(32:35) .EQ. '') e_kT_1(i__) = rNULL__
c    ..............Just test output...........
        write(6,6)
     +  Seq_5(i__),NCts_2(i__),NH_1(i__),e_NH_1(i__),l_kT_1(i__),
     +  kT_1(i__),e_kT_1(i__),logF_A_1(i__),logF_U_1(i__),chi2_1(i__),
     +  Model_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table10.dat'	! PMS stars detected by Spitzer but not detected in
*                            X-rays (2MASS colors)

C  Format for file interpretation

    7 format(
     +  A18,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,A3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table10.dat')
      write(6,*) '....Loading file: table10.dat'
      do i__=1,69
        read(1,'(A55)')ar__6
        read(ar__6,7)
     +  ID_2(i__),Jmag_1(i__),e_Jmag_1(i__),Hmag_1(i__),e_Hmag_1(i__),
     +  Ksmag(i__),e_Ksmag(i__),Qual_1(i__)
        if(ar__6(20:24) .EQ. '') Jmag_1(i__) = rNULL__
        if(ar__6(26:29) .EQ. '') e_Jmag_1(i__) = rNULL__
        if(ar__6(31:35) .EQ. '') Hmag_1(i__) = rNULL__
        if(ar__6(37:40) .EQ. '') e_Hmag_1(i__) = rNULL__
        if(ar__6(42:46) .EQ. '') Ksmag(i__) = rNULL__
        if(ar__6(48:51) .EQ. '') e_Ksmag(i__) = rNULL__
c    ..............Just test output...........
        write(6,7)
     +  ID_2(i__),Jmag_1(i__),e_Jmag_1(i__),Hmag_1(i__),e_Hmag_1(i__),
     +  Ksmag(i__),e_Ksmag(i__),Qual_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table11.dat'	! PMS stars not detected in X-rays (Spitzer)

C  Format for file interpretation

    8 format(
     +  A18,1X,A1,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,
     +  F5.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table11.dat')
      write(6,*) '....Loading file: table11.dat'
      do i__=1,69
        read(1,'(A81)')ar__7
        read(ar__7,8)
     +  ID_3(i__),f_ID(i__),v_3_6__1(i__),e__3_6__1(i__),
     +  v_4_5__1(i__),e__4_5__1(i__),v_5_6_(i__),e__5_6_(i__),
     +  v_8_0__1(i__),e__8_0__1(i__),v_24__1(i__),e__24__1(i__),
     +  AK_1(i__),Cl(i__)
        if(ar__7(33:37) .EQ. '') v_4_5__1(i__) = rNULL__
        if(ar__7(39:42) .EQ. '') e__4_5__1(i__) = rNULL__
        if(ar__7(44:48) .EQ. '') v_5_6_(i__) = rNULL__
        if(ar__7(50:53) .EQ. '') e__5_6_(i__) = rNULL__
        if(ar__7(55:59) .EQ. '') v_8_0__1(i__) = rNULL__
        if(ar__7(61:64) .EQ. '') e__8_0__1(i__) = rNULL__
        if(ar__7(66:69) .EQ. '') v_24__1(i__) = rNULL__
        if(ar__7(71:74) .EQ. '') e__24__1(i__) = rNULL__
c    ..............Just test output...........
        write(6,8)
     +  ID_3(i__),f_ID(i__),v_3_6__1(i__),e__3_6__1(i__),
     +  v_4_5__1(i__),e__4_5__1(i__),v_5_6_(i__),e__5_6_(i__),
     +  v_8_0__1(i__),e__8_0__1(i__),v_24__1(i__),e__24__1(i__),
     +  AK_1(i__),Cl(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================
      stop
      end