FORTRAN Generation
(/./ftp/cats/J/ApJS/193/25)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJS/193/25 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/ApJS/193/25       The NAN complex. II. MIPS observations       (Rebull+, 2011)
*================================================================================
*The North American and Pelican nebulae.
*II. MIPS observations and analysis.
*    Rebull L.M., Guieu S., Stauffer J.R., Hillenbrand L.A., Noriega-Crespo A.,
*    Stapelfeldt K.R., Carey S.J., Carpenter J.M., Cole D.M., Padgett D.L.,
*    Strom S.E., Wolff S.C.
*   <Astrophys. J. Suppl. Ser., 193, 25 (2011)>
*   =2011ApJS..193...25R
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table3.dat'	! Complete list of YSOs I. Spitzer measurements

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

C  J2000 position composed of: RAdeg DEdeg
      character*18  v_RGS2011_ (nr__) ! Source name (HHMMSS.ss+DDMMSS.s on J2000)
      real*8        RAdeg      (nr__) ! (deg) Right Ascension in decimal degrees (J2000)
      real*8        DEdeg      (nr__) ! (deg) Declination in decimal degrees (J2000)
      character*16  v2MASS     (nr__) ! 2MASS name (HHMMSSss+DDMMSSs)
      character*23  AName      (nr__) ! Alternate name
      character*12  SpType     (nr__) ! MK spectral type
      character*22  r_SpType   (nr__) ! Reference for SpType (1)
      character*1   l__3_6_    (nr__) ! Limit flag on 3.6]
      real*4        v_3_6_     (nr__) ! (mag) ? Spitzer/IRAC 3.6 micron band magnitude
      real*4        e__3_6_    (nr__) ! (mag) ? Uncertainty in 3.6mag
      character*1   l__4_5_    (nr__) ! Limit flag on 4.5]
      real*4        v_4_5_     (nr__) ! (mag) ? Spitzer/IRAC 4.5 micron band magnitude
      real*4        e__4_5_    (nr__) ! (mag) ? Uncertainty in 4.5mag
      character*1   l__5_8_    (nr__) ! Limit flag on 5.8]
      real*4        v_5_8_     (nr__) ! (mag) ? Spitzer/IRAC 5.8 micron band magnitude
      real*4        e__5_8_    (nr__) ! (mag) ? Uncertainty in 5.8mag
      character*1   l__8_0_    (nr__) ! Limit flag on 8.0]
      real*4        v_8_0_     (nr__) ! (mag) ? Spitzer/IRAC 8.0 micron band magnitude
      real*4        e__8_0_    (nr__) ! (mag) ? Uncertainty in 8.0mag
      character*1   l__24_     (nr__) ! Limit flag on 24]
      real*4        v_24_      (nr__) ! (mag) ? Spitzer/MIPS 24 micron band magnitude
      real*4        e__24_     (nr__) ! (mag) ? Uncertainty in 24mag
      real*4        v_70_      (nr__) ! (mag) ? Spitzer/MIPS 70 micron band magnitude
      real*4        e__70_     (nr__) ! (mag) ? Uncertainty in 70mag
*Note (1): References as follows:
*   BallyScoville    = 1980ApJ...239..121B
*   CohenKuhi79      = 1979ApJS...41..743C
*   Comeron05        = 2005A&A...430..541C; <[CP2005] NN> in Simbad
*   Corbally09       = 2009BaltA..18..111C; <[CSL2009] CSL N> in Simbad
*   Laugalys02       = 2002BaltA..11..205L
*   Laugalys06a      = 2006BaltA..15..483L (Cat. J/BaltA/15/483;
*                      <[LSV2006] Area II NNN> in Simbad)
*   Laugalys06b      = 2006BaltA..15..327L
*   MerrillBurwell49 = 1949ApJ...110..387M
*   Ogura2002        = 2002AJ....123.2597O (Cat. J/AJ/123/2597;
*                      <[OSP2002] BRC 31 NN> in Simbad)
*   Shevchenko93     = 1993Ap&SS.202..121S
*   SIMBAD           = http://simbad.u-strasbg.fr/simbad/
*   Straizys08a      = 2008BaltA..17....1S (Cat. J/BaltA/17/1)
*   Straizys08b      = 2008BaltA..17..143S; <[SL2008c] NN> in Simbad
*   Straizys09       = 2009BaltA..18..141S (Cat. J/BaltA/18/141)
*   Straizys89       = 1989VilOB..83....3S
*   Straizys93       = 1993BaltA...2..171S (Cat. II/196)
*   Straizys99       = 1999BaltA...8..355S
*   Terranegra94     = 1994A&AS..104..557T (Cat. J/A+AS/104/557)
*   Vill04           = 2004A&A...419..149V (Cat. J/A+A/419/149;
*                      <Cl* NGC 6996 VBC NNNN> in Simbad)
*   Welin73          = 1973A&AS....9..183W
*   Zdan90           = 1990Ap&SS.173..309Z

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

C  Declarations for 'table4.dat'	! Complete list of YSOs II. YSO criteria

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

      character*18  v_RGS2011__1(nr__1) ! Source name (HHMMSS.ss+DDMMSS.s)
      character*8   Clu        (nr__1) ! Cluster Name (1)
      character*1   Q          (nr__1) ! [ABC] YSO quality (A for good confidence)
      character*4   cl         (nr__1) ! YSO class (2)
      character*3   cr1        (nr__1) ! [yes/no] Selected from [3.6] vs [3.6]-[24] CMD (3)
      character*3   cr2        (nr__1) ! [yes/no] Not too faint in [3.6] vs [3.6]-[24]
*                              CMD (3)
      character*3   cr3        (nr__1) ! [yes/no] Selected from Ks vs Ks-[24] CMD (3)
      character*3   cr4        (nr__1) ! [yes/no] Detected at 70 microns (3)
      character*3   cr5        (nr__1) ! [yes/no] Found in Guieu et al. 2009
*                              Cat. J/ApJ/697/787 ([GRS2009] in Simbad)
      character*3   cr6        (nr__1) ! [yes/no] Close to the middle of the nebula (3)
      character*3   cr7        (nr__1) ! [yes/no] In a cluster (3)
      character*3   cr8        (nr__1) ! [yes/no] Previously identified YSO (3)
      character*3   cr9        (nr__1) ! [yes/no] Selected from Ic vs. V-Ic CMD (3)
      character*3   cr10       (nr__1) ! [yes/no] Selected from IPHAS r' vs. r'-i' CMD (3)
      character*3   cr11       (nr__1) ! [yes/no] A bright IPHAS H{alpha} source (3)
      character*3   cr12       (nr__1) ! [yes/no] Selected from IPHAS r'-H{alpha} vs
*                              r'-i' diagram (3)
*Note (1): Clusters as follows:
*   Pelican = IC 5070. See section 6.3 for further details.
*   PelHat  = Pelican's Hat. North of the Pelican's head, in roughly the same
*             location as would be found a hat on its head, is a previously
*             unidentified dark lane at 24um, similar to but much smaller than
*             the Gulf of Mexico. See section 6.4 for further details.
*   Gulf    = The "Gulf of Mexico" region is the area of highest extinction in
*             the NAN (The North American (NGC 7000) and Pelican (IC 5070)
*             Nebulae) complex. See section 6.2 for further details.
*Note (2): In the spirit of Wilking et al. (2001ApJ...551..357W), we define the
*          near- to mid-IR slope of the SED,
*          {alpha}=dlog{lambda}F_{lambda}_/dlog{lambda}, where {alpha}>0.3 for a
*          Class I, 0.3 to -0.3 for a flat-spectrum source, -0.3 to -1.6 for a
*          Class II, and <-1.6 for a Class III.
*Note (3): For each individual criterion, a "yes" means the object met that
*          given criterion for YSO selection, "no" means it did not, and blank
*          means that it was unable to be assessed for that criterion.
*          For example, column S3624 indicates whether each object was selected
*          as a YSO candidate in the [3.6] vs. [3.6]-[24] diagram ("yes"), or
*          it was not ("no"), or, it was missing a 3.6 or 24 micron
*          measurement, and thus could not be assessed in the [3.6] vs.
*          [3.6]-[24] diagram (blank).
*          See section 4.2 for further details on each criterion.

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

C  Declarations for 'table8.dat'	! IRAC-selected YSOs not recovered by our
                             MIPS-based search

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

C  J2000 position composed of: RAdeg DEdeg
      character*18  v_RGS2011__2(nr__2) ! Source name (HHMMSS.ss+DDMMSS.s on J2000)
      real*8        RAdeg_1    (nr__2) ! (deg) Right Ascension in decimal degrees (J2000)
      real*8        DEdeg_1    (nr__2) ! (deg) Declination in decimal degrees (J2000)
      character*16  v2MASS_1   (nr__2) ! 2MASS name (HHMMSSss+DDMMSSs)
      character*23  AName_1    (nr__2) ! Alternate name
      character*12  SpType_1   (nr__2) ! MK spectral type
      character*22  r_SpType_1 (nr__2) ! Reference for SpType (1)
      character*1   l__3_6__1  (nr__2) ! Limit flag on 3.6]
      real*4        v_3_6__1   (nr__2) ! (mag) ? Spitzer/IRAC 3.6 micron band magnitude
      real*4        e__3_6__1  (nr__2) ! (mag) ? Uncertainty in 3.6mag
      character*1   l__4_5__1  (nr__2) ! Limit flag on 4.5]
      real*4        v_4_5__1   (nr__2) ! (mag) ? Spitzer/IRAC 4.5 micron band magnitude
      real*4        e__4_5__1  (nr__2) ! (mag) ? Uncertainty in 4.5mag
      character*1   l__5_8__1  (nr__2) ! Limit flag on 5.8]
      real*4        v_5_8__1   (nr__2) ! (mag) ? Spitzer/IRAC 5.8 micron band magnitude
      real*4        e__5_8__1  (nr__2) ! (mag) ? Uncertainty in 5.8mag
      character*1   l__8_0__1  (nr__2) ! Limit flag on 8.0]
      real*4        v_8_0__1   (nr__2) ! (mag) ? Spitzer/IRAC 8.0 micron band magnitude
      real*4        e__8_0__1  (nr__2) ! (mag) ? Uncertainty in 8.0mag
      character*1   l__24__1   (nr__2) ! Limit flag on 24]
      real*4        v_24__1    (nr__2) ! (mag) ? Spitzer/MIPS 24 micron band magnitude
      real*4        e__24__1   (nr__2) ! (mag) ? Uncertainty in 24mag
      real*4        v_70__1    (nr__2) ! (mag) ? Spitzer/MIPS 70 micron band magnitude
      real*4        e__70__1   (nr__2) ! (mag) ? Uncertainty in 70mag
*Note (1): References as follows:
*   BallyScoville    = 1980ApJ...239..121B
*   CohenKuhi79      = 1979ApJS...41..743C
*   Comeron05        = 2005A&A...430..541C; <[CP2005] NN> in Simbad
*   Corbally09       = 2009BaltA..18..111C; <[CSL2009] CSL N> in Simbad
*   Laugalys02       = 2002BaltA..11..205L
*   Laugalys06a      = 2006BaltA..15..483L (Cat. J/BaltA/15/483;
*                      <[LSV2006] Area II NNN> in Simbad)
*   Laugalys06b      = 2006BaltA..15..327L
*   MerrillBurwell49 = 1949ApJ...110..387M
*   Ogura2002        = 2002AJ....123.2597O (Cat. J/AJ/123/2597;
*                      <[OSP2002] BRC 31 NN> in Simbad)
*   Shevchenko93     = 1993Ap&SS.202..121S
*   SIMBAD           = http://simbad.u-strasbg.fr/simbad/
*   Straizys08a      = 2008BaltA..17....1S (Cat. J/BaltA/17/1)
*   Straizys08b      = 2008BaltA..17..143S; <[SL2008c] NN> in Simbad
*   Straizys09       = 2009BaltA..18..141S (Cat. J/BaltA/18/141)
*   Straizys89       = 1989VilOB..83....3S
*   Straizys93       = 1993BaltA...2..171S (Cat. II/196)
*   Straizys99       = 1999BaltA...8..355S
*   Terranegra94     = 1994A&AS..104..557T (Cat. J/A+AS/104/557)
*   Vill04           = 2004A&A...419..149V (Cat. J/A+A/419/149;
*                      <Cl* NGC 6996 VBC NNNN> in Simbad)
*   Welin73          = 1973A&AS....9..183W
*   Zdan90           = 1990Ap&SS.173..309Z

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

C  Declarations for 'table9.dat'	! Previously suggested YSOs not recovered by our
                             MIPS or IRAC searches

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

C  J2000 position composed of: RAdeg DEdeg
      character*18  v_RGS2011__3(nr__3) ! Source name (HHMMSS.ss+DDMMSS.s on J2000)
      real*8        RAdeg_2    (nr__3) ! (deg) Right Ascension in decimal degrees (J2000)
      real*8        DEdeg_2    (nr__3) ! (deg) Declination in decimal degrees (J2000)
      character*16  v2MASS_2   (nr__3) ! 2MASS name (HHMMSSss+DDMMSSs)
      character*23  AName_2    (nr__3) ! Alternate name
      character*12  SpType_2   (nr__3) ! MK spectral type
      character*22  r_SpType_2 (nr__3) ! Reference for SpType (1)
      character*1   l__3_6__2  (nr__3) ! Limit flag on 3.6]
      real*4        v_3_6__2   (nr__3) ! (mag) ? Spitzer/IRAC 3.6 micron band magnitude
      real*4        e__3_6__2  (nr__3) ! (mag) ? Uncertainty in 3.6mag
      character*1   l__4_5__2  (nr__3) ! Limit flag on 4.5]
      real*4        v_4_5__2   (nr__3) ! (mag) ? Spitzer/IRAC 4.5 micron band magnitude
      real*4        e__4_5__2  (nr__3) ! (mag) ? Uncertainty in 4.5mag
      character*1   l__5_8__2  (nr__3) ! Limit flag on 5.8]
      real*4        v_5_8__2   (nr__3) ! (mag) ? Spitzer/IRAC 5.8 micron band magnitude
      real*4        e__5_8__2  (nr__3) ! (mag) ? Uncertainty in 5.8mag
      character*1   l__8_0__2  (nr__3) ! Limit flag on 8.0]
      real*4        v_8_0__2   (nr__3) ! (mag) ? Spitzer/IRAC 8.0 micron band magnitude
      real*4        e__8_0__2  (nr__3) ! (mag) ? Uncertainty in 8.0mag
      character*1   l__24__2   (nr__3) ! Limit flag on 24]
      real*4        v_24__2    (nr__3) ! (mag) ? Spitzer/MIPS 24 micron band magnitude
      real*4        e__24__2   (nr__3) ! (mag) ? Uncertainty in 24mag
      real*4        v_70__2    (nr__3) ! (mag) ? Spitzer/MIPS 70 micron band magnitude
      real*4        e__70__2   (nr__3) ! (mag) ? Uncertainty in 70mag
*Note (1): References as follows:
*   BallyScoville    = 1980ApJ...239..121B
*   CohenKuhi79      = 1979ApJS...41..743C
*   Comeron05        = 2005A&A...430..541C; <[CP2005] NN> in Simbad
*   Corbally09       = 2009BaltA..18..111C; <[CSL2009] CSL N> in Simbad
*   Laugalys02       = 2002BaltA..11..205L
*   Laugalys06a      = 2006BaltA..15..483L (Cat. J/BaltA/15/483;
*                      <[LSV2006] Area II NNN> in Simbad)
*   Laugalys06b      = 2006BaltA..15..327L
*   MerrillBurwell49 = 1949ApJ...110..387M
*   Ogura2002        = 2002AJ....123.2597O (Cat. J/AJ/123/2597;
*                      <[OSP2002] BRC 31 NN> in Simbad)
*   Shevchenko93     = 1993Ap&SS.202..121S
*   SIMBAD           = http://simbad.u-strasbg.fr/simbad/
*   Straizys08a      = 2008BaltA..17....1S (Cat. J/BaltA/17/1)
*   Straizys08b      = 2008BaltA..17..143S; <[SL2008c] NN> in Simbad
*   Straizys09       = 2009BaltA..18..141S (Cat. J/BaltA/18/141)
*   Straizys89       = 1989VilOB..83....3S
*   Straizys93       = 1993BaltA...2..171S (Cat. II/196)
*   Straizys99       = 1999BaltA...8..355S
*   Terranegra94     = 1994A&AS..104..557T (Cat. J/A+AS/104/557)
*   Vill04           = 2004A&A...419..149V (Cat. J/A+A/419/149;
*                      <Cl* NGC 6996 VBC NNNN> in Simbad)
*   Welin73          = 1973A&AS....9..183W
*   Zdan90           = 1990Ap&SS.173..309Z

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

C  Loading file 'table3.dat'	! Complete list of YSOs I. Spitzer measurements

C  Format for file interpretation

    1 format(
     +  A18,1X,F10.6,1X,F9.6,1X,A16,1X,A23,1X,A12,1X,A22,1X,A1,F5.2,
     +  1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,A1,F5.2,1X,
     +  F4.2,1X,A1,F5.2,1X,F4.2,1X,F5.2,1X,F4.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,1286
        read(1,'(A187)')ar__
        read(ar__,1)
     +  v_RGS2011_(i__),RAdeg(i__),DEdeg(i__),v2MASS(i__),AName(i__),
     +  SpType(i__),r_SpType(i__),l__3_6_(i__),v_3_6_(i__),
     +  e__3_6_(i__),l__4_5_(i__),v_4_5_(i__),e__4_5_(i__),
     +  l__5_8_(i__),v_5_8_(i__),e__5_8_(i__),l__8_0_(i__),
     +  v_8_0_(i__),e__8_0_(i__),l__24_(i__),v_24_(i__),e__24_(i__),
     +  v_70_(i__),e__70_(i__)
        if(ar__(119:123) .EQ. '') v_3_6_(i__) = rNULL__
        if(ar__(125:128) .EQ. '') e__3_6_(i__) = rNULL__
        if(ar__(131:135) .EQ. '') v_4_5_(i__) = rNULL__
        if(ar__(137:140) .EQ. '') e__4_5_(i__) = rNULL__
        if(ar__(143:147) .EQ. '') v_5_8_(i__) = rNULL__
        if(ar__(149:152) .EQ. '') e__5_8_(i__) = rNULL__
        if(ar__(155:159) .EQ. '') v_8_0_(i__) = rNULL__
        if(ar__(161:164) .EQ. '') e__8_0_(i__) = rNULL__
        if(ar__(167:171) .EQ. '') v_24_(i__) = rNULL__
        if(ar__(173:176) .EQ. '') e__24_(i__) = rNULL__
        if(ar__(178:182) .EQ. '') v_70_(i__) = rNULL__
        if(ar__(184:187) .EQ. '') e__70_(i__) = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  v_RGS2011_(i__),RAdeg(i__),DEdeg(i__),v2MASS(i__),AName(i__),
     +  SpType(i__),r_SpType(i__),l__3_6_(i__),v_3_6_(i__),
     +  e__3_6_(i__),l__4_5_(i__),v_4_5_(i__),e__4_5_(i__),
     +  l__5_8_(i__),v_5_8_(i__),e__5_8_(i__),l__8_0_(i__),
     +  v_8_0_(i__),e__8_0_(i__),l__24_(i__),v_24_(i__),e__24_(i__),
     +  v_70_(i__),e__70_(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table4.dat'	! Complete list of YSOs II. YSO criteria

C  Format for file interpretation

    2 format(
     +  A18,1X,A8,1X,A1,1X,A4,1X,A3,1X,A3,1X,A3,1X,A3,1X,A3,1X,A3,1X,
     +  A3,1X,A3,1X,A3,1X,A3,1X,A3,1X,A3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table4.dat')
      write(6,*) '....Loading file: table4.dat'
      do i__=1,1286
        read(1,'(A82)')ar__1
        read(ar__1,2)
     +  v_RGS2011__1(i__),Clu(i__),Q(i__),cl(i__),cr1(i__),cr2(i__),
     +  cr3(i__),cr4(i__),cr5(i__),cr6(i__),cr7(i__),cr8(i__),
     +  cr9(i__),cr10(i__),cr11(i__),cr12(i__)
c    ..............Just test output...........
        write(6,2)
     +  v_RGS2011__1(i__),Clu(i__),Q(i__),cl(i__),cr1(i__),cr2(i__),
     +  cr3(i__),cr4(i__),cr5(i__),cr6(i__),cr7(i__),cr8(i__),
     +  cr9(i__),cr10(i__),cr11(i__),cr12(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table8.dat'	! IRAC-selected YSOs not recovered by our
*                             MIPS-based search

C  Format for file interpretation

    3 format(
     +  A18,1X,F10.6,1X,F9.6,1X,A16,1X,A23,1X,A12,1X,A22,1X,A1,F5.2,
     +  1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,A1,F5.2,1X,
     +  F4.2,1X,A1,F5.2,1X,F4.2,1X,F5.2,1X,F4.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table8.dat')
      write(6,*) '....Loading file: table8.dat'
      do i__=1,796
        read(1,'(A187)')ar__2
        read(ar__2,3)
     +  v_RGS2011__2(i__),RAdeg_1(i__),DEdeg_1(i__),v2MASS_1(i__),
     +  AName_1(i__),SpType_1(i__),r_SpType_1(i__),l__3_6__1(i__),
     +  v_3_6__1(i__),e__3_6__1(i__),l__4_5__1(i__),v_4_5__1(i__),
     +  e__4_5__1(i__),l__5_8__1(i__),v_5_8__1(i__),e__5_8__1(i__),
     +  l__8_0__1(i__),v_8_0__1(i__),e__8_0__1(i__),l__24__1(i__),
     +  v_24__1(i__),e__24__1(i__),v_70__1(i__),e__70__1(i__)
        if(ar__2(119:123) .EQ. '') v_3_6__1(i__) = rNULL__
        if(ar__2(125:128) .EQ. '') e__3_6__1(i__) = rNULL__
        if(ar__2(131:135) .EQ. '') v_4_5__1(i__) = rNULL__
        if(ar__2(137:140) .EQ. '') e__4_5__1(i__) = rNULL__
        if(ar__2(143:147) .EQ. '') v_5_8__1(i__) = rNULL__
        if(ar__2(149:152) .EQ. '') e__5_8__1(i__) = rNULL__
        if(ar__2(155:159) .EQ. '') v_8_0__1(i__) = rNULL__
        if(ar__2(161:164) .EQ. '') e__8_0__1(i__) = rNULL__
        if(ar__2(167:171) .EQ. '') v_24__1(i__) = rNULL__
        if(ar__2(173:176) .EQ. '') e__24__1(i__) = rNULL__
        if(ar__2(178:182) .EQ. '') v_70__1(i__) = rNULL__
        if(ar__2(184:187) .EQ. '') e__70__1(i__) = rNULL__
c    ..............Just test output...........
        write(6,3)
     +  v_RGS2011__2(i__),RAdeg_1(i__),DEdeg_1(i__),v2MASS_1(i__),
     +  AName_1(i__),SpType_1(i__),r_SpType_1(i__),l__3_6__1(i__),
     +  v_3_6__1(i__),e__3_6__1(i__),l__4_5__1(i__),v_4_5__1(i__),
     +  e__4_5__1(i__),l__5_8__1(i__),v_5_8__1(i__),e__5_8__1(i__),
     +  l__8_0__1(i__),v_8_0__1(i__),e__8_0__1(i__),l__24__1(i__),
     +  v_24__1(i__),e__24__1(i__),v_70__1(i__),e__70__1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table9.dat'	! Previously suggested YSOs not recovered by our
*                             MIPS or IRAC searches

C  Format for file interpretation

    4 format(
     +  A18,1X,F10.6,1X,F9.6,1X,A16,1X,A23,1X,A12,1X,A22,1X,A1,F5.2,
     +  1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,A1,F5.2,1X,
     +  F4.2,1X,A1,F5.2,1X,F4.2,1X,F5.2,1X,F4.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table9.dat')
      write(6,*) '....Loading file: table9.dat'
      do i__=1,114
        read(1,'(A187)')ar__3
        read(ar__3,4)
     +  v_RGS2011__3(i__),RAdeg_2(i__),DEdeg_2(i__),v2MASS_2(i__),
     +  AName_2(i__),SpType_2(i__),r_SpType_2(i__),l__3_6__2(i__),
     +  v_3_6__2(i__),e__3_6__2(i__),l__4_5__2(i__),v_4_5__2(i__),
     +  e__4_5__2(i__),l__5_8__2(i__),v_5_8__2(i__),e__5_8__2(i__),
     +  l__8_0__2(i__),v_8_0__2(i__),e__8_0__2(i__),l__24__2(i__),
     +  v_24__2(i__),e__24__2(i__),v_70__2(i__),e__70__2(i__)
        if(ar__3(119:123) .EQ. '') v_3_6__2(i__) = rNULL__
        if(ar__3(125:128) .EQ. '') e__3_6__2(i__) = rNULL__
        if(ar__3(131:135) .EQ. '') v_4_5__2(i__) = rNULL__
        if(ar__3(137:140) .EQ. '') e__4_5__2(i__) = rNULL__
        if(ar__3(143:147) .EQ. '') v_5_8__2(i__) = rNULL__
        if(ar__3(149:152) .EQ. '') e__5_8__2(i__) = rNULL__
        if(ar__3(155:159) .EQ. '') v_8_0__2(i__) = rNULL__
        if(ar__3(161:164) .EQ. '') e__8_0__2(i__) = rNULL__
        if(ar__3(167:171) .EQ. '') v_24__2(i__) = rNULL__
        if(ar__3(173:176) .EQ. '') e__24__2(i__) = rNULL__
        if(ar__3(178:182) .EQ. '') v_70__2(i__) = rNULL__
        if(ar__3(184:187) .EQ. '') e__70__2(i__) = rNULL__
c    ..............Just test output...........
        write(6,4)
     +  v_RGS2011__3(i__),RAdeg_2(i__),DEdeg_2(i__),v2MASS_2(i__),
     +  AName_2(i__),SpType_2(i__),r_SpType_2(i__),l__3_6__2(i__),
     +  v_3_6__2(i__),e__3_6__2(i__),l__4_5__2(i__),v_4_5__2(i__),
     +  e__4_5__2(i__),l__5_8__2(i__),v_5_8__2(i__),e__5_8__2(i__),
     +  l__8_0__2(i__),v_8_0__2(i__),e__8_0__2(i__),l__24__2(i__),
     +  v_24__2(i__),e__24__2(i__),v_70__2(i__),e__70__2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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