FORTRAN Generation
(/./ftp/cats/J/AJ/160/44)

Conversion of standardized ReadMe file for file /./ftp/cats/J/AJ/160/44 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.



++++Only one proper motion (pmDec) found for: table1.dat
      program load_ReadMe
C=============================================================================
C  F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2026-Sep-07
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/160/44      Upper Scorpius spectroscopy and photometry      (Luhman+, 2020)
*================================================================================
*Refining the Census of the Upper Scorpius Association with Gaia.
*    Luhman K.L., Esplin T.L.
*   <Astron. J., 160, 44 (2020)>
*   =2020AJ....160...44L
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! Candidate members of Upper Sco

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

C  J2000 position composed of: RAdeg DEdeg
      character*17  v2MASS      ! 2MASS Point Source Catalog source name
      character*19  WISEA       ! AllWISE Source Catalog source name
      character*19  UGCS        ! UKIDSS Galactic Clusters Survey source name (1)
      integer*8     Gaia        ! ? Gaia DR2 source name
      real*8        RAdeg       ! (deg) [233/252] Right Ascension (J2000)
      real*8        DEdeg       ! (deg) [-30/-16] Declination (J2000)
      character*6   r_RAdeg     ! Reference for position (2)
      character*5   SpType      ! Spectral type
      integer*4     r_SpType    ! ? Spectral type reference
      real*8        pmRA        ! (mas/yr) [-33/3]? Proper motion in RA from Gaia DR2
      real*4        e_pmRA      ! (mas/yr) [0.07/4.1]? Error in pmRA
      real*8        pmDec       ! (mas/yr) [-40/-3.1]? Proper motion in Dec from Gaia DR2
      real*4        e_pmDec     ! (mas/yr) [0.03/3.2]? Error in pmDec
      real*4        plx         ! (mas) [0.25/16]? Parallax from Gaia DR2
      real*4        e_plx       ! (mas) [0.03/1.8]? Error in plx
      real*4        Gmag        ! (mag) [9.02/21.2]? G magnitude from Gaia DR2
      real*4        e_Gmag      ! (mag) [0.001/0.09]? Error in Gmag
      real*4        GBPmag      ! (mag) [9.25/22.4]? GBP magnitude from Gaia DR2
      real*4        e_GBPmag    ! (mag) [0.001/0.42]? Error in GBPmag
      real*4        GRPmag      ! (mag) [8.66/19.3]? GRP magnitude from Gaia DR2
      real*4        e_GRPmag    ! (mag) [0.001/0.21]? Error in GRPmag
      real*4        RUWE        ! [0.82/58.5]? Renormalized unit weight error (3)
      real*4        Jmag        ! (mag) [8.29/19.5]? J band magnitude
      real*4        e_Jmag      ! (mag) [0.02/0.15]? Error in Jmag
      real*4        Hmag        ! (mag) [7.9/18.5]? H band magnitude
      real*4        e_Hmag      ! (mag) [0.02/0.15]? Error in Hmag
      real*4        Ksmag       ! (mag) [7.48/17.3]? Ks band magnitude
      real*4        e_Ksmag     ! (mag) [0.01/0.1]? Error in Ksmag
      character*3   r_JHKmag    ! ? JHK reference (4)
      character*49  Sel         ! Selection criteria satisfied by candidate (5)
      real*4        Sep         ! (arcsec) [0.66/4.29]? Angular separation (6)
      integer*8     CompGaia    ! ? Gaia DR2 source name (6)
*Note (1): Based on coordinates from DR10 of the UKIDSS Galactic Clusters
*          Survey for stars with Ks>10mag from 2MASS.
*Note (2): References as follows:
*          2MASS = 2MASS Point Source Catalog; II/246 (2 occurrences)
*          Gaia = Gaia DR2; I/345 (481 occurrences)
*          UKIDSS = UKIDSS Galactic Clusters Survey DR10 (76 occurrences)
*          VHS = VISTA VHS DR6 (8 occurrences)
*Note (3): From Lindegren, 2018, http://www.rssd.esa.int/doc_fetch.php?id=3757412
*Note (4):
*    2 = 2MASS Point Source Catalog, II/246
*    u = DR10 of the UKIDSS Galactic Clusters Survey;
*    v = DR6 of VISTA VHS.
*Note (5):
*                    G/W/i/Y/Z/ip/zp/yp = CMDs in Luhman et al. (2018);
*                                    pi = parallax from Gaia DR2;
*    gaia/gps/ucac/2m-gaia/2m-ps/ukidss = proper motions in
*                                         Luhman et al. (2018) and this work.
*Note (6): Of nearest known young star within 5".

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

C  Declarations for 'table2.dat'	! Observing log

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

      character*14  Tel         ! Telescope used
      character*11  Inst        ! Instrument used
      character*9   Disp        ! Disperser
      character*10  Aper        ! Aperture
      character*9   lambda      ! (um) Wavelength coverage
      character*6   R           ! Resolution
      integer*4     Ntarget     ! [3/179] Number of observed targets

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

C  Declarations for 'table3.dat'	! Spectroscopic data for candidate members of
                              Upper Sco

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

C  J2000 position composed of: RAdeg DEdeg
      character*24  Name        ! Source Name (1)
      real*8        RAdeg_1     ! (deg) [235/251] Right Ascension (J2000)
      real*8        DEdeg_1     ! (deg) [-30/-16] Declination (J2000)
      character*11  SpType_1    ! Spectral type (2)
      character*1   l_EWLi      ! Limit flag on EWLi
      real*4        EWLi        ! (0.1nm) [0.1/0.8]? Equivalent Width of Li (3)
      character*1   f_EWLi      ! [:] Uncertainty flag on EWLi (1 occurrence)
      character*11  Inst_1      ! Instrument used (4)
      integer*4     Obs_Y       ! (yr) [2018/2020] Year of observation
      character*3   Obs_M       ! ("month") Month of observation
      integer*4     Obs_D       ! (d) Day of observation
      integer*4     Obs_Y_2     ! (yr) [2019]? Year of second observation
      character*3   Obs_M_2     ! ("month") Month of second observation
      integer*4     Obs_D_2     ! (d) ? Day of second observation
      character*2   Young       ! Young? Y=yes, N=No or ?=unknown
*Note (1): Source names are from DR2 of Gaia, DR10 of the UKIDSS Galactic
*          Clusters Survey, and the 2MASS Point Source Catalog.
*Note (2): Uncertainties are 0.25 and 0.5 subclass for optical and IR spectral
*          types, respectively, unless indicated otherwise.
*Note (3): Typical uncertainties are 0.05{AA}.
*Note (4): Instruments as follows:
*          BC = Bok Reflector /B&C at the Steward Observatory (24 occurrences)
*          COSMOS = 4m Blanco telescope at the Cerro Tololo Inter-American
*                   Observatory (177 occurrences)
*          SpeX = The NASA Infrared Telescope Facility (179 occurrences)
*          FIRE = The Magellan Baade telescope at Las Campanas Observatory
*                 (3 occurrences)
*          FLAMINGOS-2 = Gemini South observatory (3 occurrences)
*          GNIRS = Gemini North observatory (47 occurrences)

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

C  Declarations for 'table4.dat'	! Young stars toward Upper Sco

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

C  J2000 position composed of: RAdeg DEdeg
      character*18  v2MASS_1    ! 2MASS Point Source Catalog source name
      character*20  WISEA_1     ! AllWISE Source Catalog source name (1)
      character*19  UGCS_1      ! UKIDSS Galactic Clusters Survey source name (2)
      integer*8     Gaia_1      ! ? Gaia DR2 source name
      character*54  Name_1      ! Other source name
      real*8        RAdeg_2     ! (deg) [233/254] Right Ascension (J2000)
      real*8        DEdeg_2     ! (deg) [-34/-14] Declination  (J2000)
      character*6   r_RAdeg_1   ! Reference for RA & Dec (3)
      character*23  SpType_2    ! Spectral type
      character*13  r_SpType_1  ! Spectral type reference
      character*6   SpT         ! Adopted spectral type
      real*4        Ak          ! (mag) [0/1.469]? Extinction in Ks
      character*5   f_Ak        ! Method of extinction estimation (4)
      real*8        pmRA_1      ! (mas/yr) [-39.768/2.768]? Proper motion in RA from
*                                Gaia DR2
      real*4        e_pmRA_1    ! (mas/yr) [0.056/4.184]? Error in pmRA
      real*8        pmDE        ! (mas/yr) [-48.798/-9.414]? Proper motion in DE from
*                                Gaia DR2
      real*4        e_pmDE      ! (mas/yr) [0.039/2.619]? Error in pmDE
      real*4        plx_1       ! (mas) [0.764/14.907]? Parallax from Gaia DR2
      real*4        e_plx_1     ! (mas) [0.027/2.647]? Error in plx
      real*8        RVel        ! (km/s) [-40/114]? Radial velocity
      real*4        e_RVel      ! (km/s) [0.04/4]? Error in RVel
      character*1   f_RVel      ! [r] Indicates RVel ranges between 1 and 2
      integer*4     r_RVel      ! [1/21]? Reference in RVel (5)
      real*4        U           ! (km/s) [-40/12.2]? U component of space velocity
      real*4        e_U         ! (km/s) [0.07/3.8]? Error in U
      real*4        V           ! (km/s) [-23.1/-9.03]? V component of space velocity
      real*4        e_V         ! (km/s) [0.09/2.1]? Error in V
      real*4        W           ! (km/s) [-16.9/3.52]? W component of space velocity
      real*4        e_W         ! (km/s) [0.07/1.8]? Error in W
      real*4        Gmag_1      ! (mag) [2.26/21.5]? G magnitude from Gaia DR2
      real*4        e_Gmag_1    ! (mag) [0/0.08]? Error in Gmag
      real*4        GBPmag_1    ! (mag) [2.98/22.2]? GBP magnitude from Gaia DR2
      real*4        e_GBPmag_1  ! (mag) [0.001/0.97]? Error in GBPmag
      real*4        GRPmag_1    ! (mag) [2.64/19.3]? GRP magnitude from Gaia DR2
      real*4        e_GRPmag_1  ! (mag) [0.001/0.33]? Error in GRPmag
      real*4        RUWE_1      ! [0.69/98]? Renormalized unit weight error (6)
      real*4        Jmag_1      ! (mag) [-2.73/19.5]? J band magnitude
      real*4        eJmag       ! (mag) [0.01/0.24]? Error in Jmag
      real*4        Hmag_1      ! (mag) [-3.49/18.2]? H band magnitude
      real*4        eHmag       ! (mag) [0.01/0.23]? Error in Hmag
      real*4        Ksmag_1     ! (mag) [-3.79/17.2]? Ks band magnitude
      real*4        eKsmag      ! (mag) [0.01/0.29]? Error in Ksmag
      character*3   r_JHKmag_1  ! ? JHKs reference (7)
      character*5   Pops        ! Gaia kinematic populations (8)
      character*1   memb_       ! Adopted member of Upper Sco? Y=Yes
*                                (1761 occurrences) or N=No (259 occurrences)
*Note (1): The following names are from the WISE All-Sky Source Catalog:
*          J160027.15-223850.5, J160414.16-212915.5, J161320.78-175752.3,
*          J161317.38-292220.0, J161837.22-240522.8, J162210.14-240905.4,
*          J162620.15-223312.8.
*Note (2): Based on coordinates from DR10 of the UKIDSS Galactic Clusters
*          Survey for stars with Ks>10 from 2MASS.
*Note (3): References as follows :
*          2MASS = 2MASS Point Source Catalog; II/246 (6 occurrences)
*          Gaia = Gaia = Gaia DR2; I/345 (1940 occurrences)
*          UKIDSS = UKIDSS Galactic Clusters Survey DR10 (70 occurrences)
*          bry = Bryan+, 2016ApJ...827..100B (1 occurrence)
*          ire = Ireland+, 2011ApJ...726..113I (1 occurrence)
*          lac = Lachapelle+, 2015ApJ...802...61L (2 occurrences)
*Note (4): Extinction estimated from a near-IR spectrum or the indicated color
*          assuming the intrinsic spectrum from Luhman+, 2017, J/AJ/153/46 or
*          the intrinsic color from the Appendix for the spectral type in
*          question.
*Note (5): References as follows
*    1 = Dahm+, 2012, J/ApJ/745/56
*    2 = Gaia DR2, 2018, I/345
*    3 = Jilinski+, 2006, J/A+A/448/1001
*    4 = Gontcharov, 2006, III/252
*    5 = Torres+, 2006, J/A+A/460/695
*    6 = Chen+, 2011AJ....738..122C
*    7 = Kurosawa+, 2006MNRAS.372.1879K
*    8 = Kunder+, 2017, III/279
*    9 = Durkan+, 2018A&A...618....5D
*   10 = Andersen & Nordstrom, 1983 III/127
*   11 = Muzerolle+, 2003ApJ...592..266M
*   12 = David+, 2016ApJ...827..161D
*   13 = White+, 2007, J/AJ/133/2524
*   14 = Mace+, 2012, J/AJ/144/55
*   15 = Prato+, 2002ApJ...569..863P
*   16 = Rice+, 2010ApJS..186...63R
*   17 = Alonso+, 2015A&A...584L...8A
*   18 = David+, 2016Natur.534..658D
*   19 = Eisner+, 2005ApJ...623..952E
*   20 = Wang+, 2018ApJ...865..141W
*   21 = David+, 2017ApJ...835..168D
*Note (6): From Lindegren, 2018, http://www.rssd.esa.int/doc_fetch.php?id=3757412
*Note (7): References as follows:
*    2 = 2MASS Point Source Catalog, II/246
*    u = UKIDSS Galactic Clusters Survey DR10
*    v = VISTA VHS DR6;
*    van = van der Bliek+, 1996, J/A+AS/119/547
*    duc = Ducati, 2002, II/237
*    luh = Luhman+, 2007ApJ...666.1219L
*    bow = Bowler+, 2011ApJ...743..148B
*    kra = Kraus+, 2014ApJ...781...20K
*    lac = Lachapelle+, 2015ApJ...802...61L
*Note (8): The Gaia parallax and proper motion are consistent with
*          populations as follows:
*    O/o = Ophiuchus;
*    U/u = Upper Sco;
*    S/s = remainder of Sco-Cen;
*      n = none of those populations. Upper case letters
*          Upper case letters indicate that both proper motion and parallax
*          support membership in that population. Lower case letters indicate
*          that either proper motion or parallax supports membership while
*          the other parameter is inaccurate, unreliable, or unavailable.

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

C  Declarations for 'table5.dat'	! Mid-IR photometry for adopted Members of Upper Sco

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

C  J2000 position composed of: RAdeg DEdeg
      character*18  v2MASS_2    ! 2MASS Point Source Catalog source name
      character*20  WISEA_2     ! AllWISE Source Catalog source name
      character*19  UGCS_2      ! UKIDSS Galactic Clusters Survey source name
      integer*8     Gaia_2      ! ? Gaia DR2 source name
      character*54  Name_2      ! Other source name
      real*8        RAdeg_3     ! (deg) [233/253] Right Ascension (J2000)
      real*8        DEdeg_3     ! (deg) [-32/-16] Declination (J2000)
      character*6   r_RAdeg_2   ! Reference for RA & Dec (1)
      real*4        v3_6mag     ! (mag) [6.03/15.3]? Spitzer/IRAC 3.6um magnitude
      real*4        e_3_6mag    ! (mag) [0.02/0.21]? Error in 3.6mag
      character*5   f_3_6mag    ! ? Flag on 3.6mag (2)
      real*4        v4_5mag     ! (mag) [2.89/15.1]? Spitzer/IRAC 4.5um magnitude
      real*4        e_4_5mag    ! (mag) [0.02/0.47]? Error in 4.5mag
      character*5   f_4_5mag    ! ? Flag on 4.5mag (2)
      real*4        v5_8mag     ! (mag) [4.88/13.7]? Spitzer/IRAC 5.8um magnitude
      real*4        e_5_8mag    ! (mag) [0.02/0.1]? Error in 5.8mag
      character*5   f_5_8mag    ! ? Flag on 5.8mag (2)
      real*4        v8_0mag     ! (mag) [3.02/13.3]? Spitzer/IRAC 8.0um magnitude
      real*4        e_8_0mag    ! (mag) [0.02/0.2]? Error in 8.0mag
      character*5   f_8_0mag    ! ? Flag on 8.0mag (2)
      real*4        v24mag      ! (mag) [0.47/11.5]? Spitzer/IRAC 24um magnitude
      real*4        e_24mag     ! (mag) [0.02/0.4]? Error in 24mag
      character*11  f_24mag     ! ? Flag on 24mag (2)
      real*4        W1mag       ! (mag) [-1.23/16.6]? WISE W1 magnitude
      real*4        e_W1mag     ! (mag) [0.01/0.43]? Error in W1mag
      character*5   f_W1mag     ! ? Flag on W1mag (2)
      real*4        W2mag       ! (mag) [6.02/15.8]? WISE W2 magnitude
      real*4        e_W2mag     ! (mag) [0.01/0.22]? Error in W2mag
      character*5   f_W2mag     ! ? Flag on W2mag (2)
      real*4        W3mag       ! (mag) [1.11/12.3]? WISE W3 magnitude
      real*4        e_W3mag     ! (mag) [0.01/0.48]? Error in W3mag
      character*11  f_W3mag     ! ? Flag on W3mag (2)
      real*4        W4mag       ! (mag) [-4.28/8.9]? WISE W4 magnitude
      real*4        e_W4mag     ! (mag) [0/0.52]? Error in W4mag
      character*11  f_W4mag     ! ? Flag on W4mag (2)
      character*2   Exc4_5      ! ? Excess present in 4.5?
      character*5   Exc8_0      ! ? Excess present in 8.0?
      character*2   Exc24       ! ? Excess present in 24?
      character*2   ExcW2       ! ? Excess present in W2?
      character*2   ExcW3       ! ? Excess present in W3?
      character*2   ExcW4       ! ? Excess present in W4?
      character*16  DiskType    ! ? Disk Type
*Note (1): References as follows :
*          2MASS = 2MASS Point Source Catalog; II/246 (6 occurrences)
*          Gaia = Gaia = Gaia DR2; I/345 (1682 occurrences)
*          UKIDSS = UKIDSS Galactic Clusters Survey DR10 (70 occurrences)
*          bry = Bryan+, 2016ApJ...827..100B (1 occurrence)
*          ire = Ireland+, 2011ApJ...726..113I (1 occurrence)
*          lac = Lachapelle+, 2015ApJ...802...61L (1 occurrence)
*Note (2): Flags as follows:
*     nodet = non-detection;
*       sat = saturated;
*       out = outside of the camera's field of view;
*        bl = photometry may be affected by blending with a nearby star;
*       ext = photometry is known or suspected to be contaminated by
*             extended emission (no data given when extended emission dominates);
*       dif = photometery may be affected by a diffraction spike;
*       bin = includes an unresolved binary companion;
*     unres = too close to a brighter star to be detected;
*     false = detection from WISE catalog appears false or unreliable based
*             on visual inspection;
*       off = W3 and/or W4 detection appears offset from the W1/W2 detection
*             (no data given when offset is due to a known source);
*       err = W2 magnitudes brighter than ~6 are erroneous.

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

C  Declarations for 'tablec1.dat'	! Intrinsic colors of young stars and brown dwarfs

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

      character*2   SpType_3    ! Spectral Type
      real*4        BP_RP       ! (mag) ? GBP-GRP for Gaia DR2 bands
      real*4        G_RP        ! (mag) ? G-GRP for Gaia DR2 bands
      real*4        RP_K        ! (mag) ? GRP-Ks  for Gaia DR2 bands
      real*4        J_H         ! (mag) J-H for 2MASS bands
      real*4        H_K         ! (mag) H-Ks for 2MASS bands
      real*4        K_3_6       ! (mag) ? Ks-3.6 for 2MASS and Spitzer bands
      real*4        K_4_5       ! (mag) Ks-4.5 for 2MASS and Spitzer bands (1)
      real*4        K_5_8       ! (mag) Ks-5.8 for 2MASS and Spitzer bands
      real*4        K_8_0       ! (mag) Ks-8.0 for 2MASS and Spitzer bands
      real*4        K_24        ! (mag) ? Ks-24 for 2MASS and Spitzer bands (2)
      real*4        K_W1        ! (mag) ? Ks-W1 for 2MASS and WISE bands
      real*4        K_W2        ! (mag) Ks-W2 for 2MASS and WISE bands (1)
      real*4        K_W3        ! (mag) ? Ks-W3 for 2MASS and WISE bands
      real*4        K_W4        ! (mag) ? Ks-W4 for 2MASS and WISE bands (2)
*Note (1): Same values listed for Ks-4.5 and Ks-W2.
*Note (2): Same values listed for Ks-[24] and Ks-W4.

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

C  Declarations for 'ref.dat'	! References used for spectral types

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

      integer*4     Ref         ! [1/93] Reference number
      character*21  Autor       ! Main autor of the reference
      character*19  BIB         ! BIBcode
      character*63  Cat         ! Catalogue or comment

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

C  Loading file 'table1.dat'	! Candidate members of Upper Sco

C  Format for file interpretation

    1 format(
     +  A17,1X,A19,1X,A19,1X,I19,1X,F10.6,1X,F10.6,1X,A6,1X,A5,1X,I2,
     +  F7.3,1X,F5.3,1X,F7.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,
     +  1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F5.2,1X,F5.2,1X,F4.2,1X,
     +  F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,A3,1X,A49,1X,F4.2,1X,I19)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,567
        read(1,'(A312)')ar__
        read(ar__,1)
     +  v2MASS,WISEA,UGCS,Gaia,RAdeg,DEdeg,r_RAdeg,SpType,r_SpType,
     +  pmRA,e_pmRA,pmDec,e_pmDec,plx,e_plx,Gmag,e_Gmag,GBPmag,
     +  e_GBPmag,GRPmag,e_GRPmag,RUWE,Jmag,e_Jmag,Hmag,e_Hmag,Ksmag,
     +  e_Ksmag,r_JHKmag,Sel,Sep,CompGaia
        if(ar__(59:77) .EQ. '') Gaia = iNULL__
        if(ar__(114:115) .EQ. '') r_SpType = iNULL__
        if(ar__(116:122) .EQ. '') pmRA = rNULL__
        if(ar__(124:128) .EQ. '') e_pmRA = rNULL__
        if(ar__(130:136) .EQ. '') pmDec = rNULL__
        if(ar__(138:142) .EQ. '') e_pmDec = rNULL__
        if(ar__(144:149) .EQ. '') plx = rNULL__
        if(ar__(151:155) .EQ. '') e_plx = rNULL__
        if(ar__(157:162) .EQ. '') Gmag = rNULL__
        if(ar__(164:168) .EQ. '') e_Gmag = rNULL__
        if(ar__(170:175) .EQ. '') GBPmag = rNULL__
        if(ar__(177:181) .EQ. '') e_GBPmag = rNULL__
        if(ar__(183:188) .EQ. '') GRPmag = rNULL__
        if(ar__(190:194) .EQ. '') e_GRPmag = rNULL__
        if(ar__(196:200) .EQ. '') RUWE = rNULL__
        if(ar__(202:206) .EQ. '') Jmag = rNULL__
        if(ar__(208:211) .EQ. '') e_Jmag = rNULL__
        if(ar__(213:217) .EQ. '') Hmag = rNULL__
        if(ar__(219:222) .EQ. '') e_Hmag = rNULL__
        if(ar__(224:228) .EQ. '') Ksmag = rNULL__
        if(ar__(230:233) .EQ. '') e_Ksmag = rNULL__
        if(ar__(289:292) .EQ. '') Sep = rNULL__
        if(ar__(294:312) .EQ. '') CompGaia = iNULL__
c    ..............Just test output...........
        write(6,1)
     +  v2MASS,WISEA,UGCS,Gaia,RAdeg,DEdeg,r_RAdeg,SpType,r_SpType,
     +  pmRA,e_pmRA,pmDec,e_pmDec,plx,e_plx,Gmag,e_Gmag,GBPmag,
     +  e_GBPmag,GRPmag,e_GRPmag,RUWE,Jmag,e_Jmag,Hmag,e_Hmag,Ksmag,
     +  e_Ksmag,r_JHKmag,Sel,Sep,CompGaia
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table2.dat'	! Observing log

C  Format for file interpretation

    2 format(A14,1X,A11,1X,A9,1X,A10,1X,A9,1X,A6,1X,I3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,6
        read(1,'(A68)')ar__1
        read(ar__1,2)Tel,Inst,Disp,Aper,lambda,R,Ntarget
c    ..............Just test output...........
        write(6,2)Tel,Inst,Disp,Aper,lambda,R,Ntarget
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table3.dat'	! Spectroscopic data for candidate members of
*                              Upper Sco

C  Format for file interpretation

    3 format(
     +  A24,1X,F10.6,1X,F10.6,1X,A11,1X,A1,F4.2,A1,1X,A11,1X,I4,1X,A3,
     +  1X,I2,1X,I4,1X,A3,1X,I2,1X,A2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,431
        read(1,'(A104)')ar__2
        read(ar__2,3)
     +  Name,RAdeg_1,DEdeg_1,SpType_1,l_EWLi,EWLi,f_EWLi,Inst_1,Obs_Y,
     +  Obs_M,Obs_D,Obs_Y_2,Obs_M_2,Obs_D_2,Young
        if(ar__2(61:64) .EQ. '') EWLi = rNULL__
        if(ar__2(91:94) .EQ. '') Obs_Y_2 = iNULL__
        if(ar__2(100:101) .EQ. '') Obs_D_2 = iNULL__
c    ..............Just test output...........
        write(6,3)
     +  Name,RAdeg_1,DEdeg_1,SpType_1,l_EWLi,EWLi,f_EWLi,Inst_1,Obs_Y,
     +  Obs_M,Obs_D,Obs_Y_2,Obs_M_2,Obs_D_2,Young
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table4.dat'	! Young stars toward Upper Sco

C  Format for file interpretation

    4 format(
     +  A18,1X,A20,1X,A19,1X,I19,1X,A54,1X,F10.6,1X,F10.6,1X,A6,1X,
     +  A23,1X,A13,1X,A6,1X,F5.3,1X,A5,1X,F7.3,1X,F5.3,1X,F7.3,1X,
     +  F5.3,1X,F6.3,1X,F5.3,1X,F7.3,1X,F5.3,A1,1X,I2,1X,F6.2,1X,F4.2,
     +  1X,F6.2,1X,F4.2,1X,F6.2,1X,F4.2,1X,F6.3,1X,F5.3,1X,F6.3,1X,
     +  F5.3,1X,F6.3,1X,F5.3,1X,F5.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,
     +  1X,F5.2,1X,F4.2,1X,A3,1X,A5,1X,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table4.dat')
      write(6,*) '....Loading file: table4.dat'
      do i__=1,2020
        read(1,'(A405)')ar__3
        read(ar__3,4)
     +  v2MASS_1,WISEA_1,UGCS_1,Gaia_1,Name_1,RAdeg_2,DEdeg_2,
     +  r_RAdeg_1,SpType_2,r_SpType_1,SpT,Ak,f_Ak,pmRA_1,e_pmRA_1,
     +  pmDE,e_pmDE,plx_1,e_plx_1,RVel,e_RVel,f_RVel,r_RVel,U,e_U,V,
     +  e_V,W,e_W,Gmag_1,e_Gmag_1,GBPmag_1,e_GBPmag_1,GRPmag_1,
     +  e_GRPmag_1,RUWE_1,Jmag_1,eJmag,Hmag_1,eHmag,Ksmag_1,eKsmag,
     +  r_JHKmag_1,Pops,memb_
        if(ar__3(61:79) .EQ. '') Gaia_1 = iNULL__
        if(ar__3(210:214) .EQ. '') Ak = rNULL__
        if(ar__3(222:228) .EQ. '') pmRA_1 = rNULL__
        if(ar__3(230:234) .EQ. '') e_pmRA_1 = rNULL__
        if(ar__3(236:242) .EQ. '') pmDE = rNULL__
        if(ar__3(244:248) .EQ. '') e_pmDE = rNULL__
        if(ar__3(250:255) .EQ. '') plx_1 = rNULL__
        if(ar__3(257:261) .EQ. '') e_plx_1 = rNULL__
        if(ar__3(263:269) .EQ. '') RVel = rNULL__
        if(ar__3(271:275) .EQ. '') e_RVel = rNULL__
        if(ar__3(278:279) .EQ. '') r_RVel = iNULL__
        if(ar__3(281:286) .EQ. '') U = rNULL__
        if(ar__3(288:291) .EQ. '') e_U = rNULL__
        if(ar__3(293:298) .EQ. '') V = rNULL__
        if(ar__3(300:303) .EQ. '') e_V = rNULL__
        if(ar__3(305:310) .EQ. '') W = rNULL__
        if(ar__3(312:315) .EQ. '') e_W = rNULL__
        if(ar__3(317:322) .EQ. '') Gmag_1 = rNULL__
        if(ar__3(324:328) .EQ. '') e_Gmag_1 = rNULL__
        if(ar__3(330:335) .EQ. '') GBPmag_1 = rNULL__
        if(ar__3(337:341) .EQ. '') e_GBPmag_1 = rNULL__
        if(ar__3(343:348) .EQ. '') GRPmag_1 = rNULL__
        if(ar__3(350:354) .EQ. '') e_GRPmag_1 = rNULL__
        if(ar__3(356:360) .EQ. '') RUWE_1 = rNULL__
        if(ar__3(362:366) .EQ. '') Jmag_1 = rNULL__
        if(ar__3(368:371) .EQ. '') eJmag = rNULL__
        if(ar__3(373:377) .EQ. '') Hmag_1 = rNULL__
        if(ar__3(379:382) .EQ. '') eHmag = rNULL__
        if(ar__3(384:388) .EQ. '') Ksmag_1 = rNULL__
        if(ar__3(390:393) .EQ. '') eKsmag = rNULL__
c    ..............Just test output...........
        write(6,4)
     +  v2MASS_1,WISEA_1,UGCS_1,Gaia_1,Name_1,RAdeg_2,DEdeg_2,
     +  r_RAdeg_1,SpType_2,r_SpType_1,SpT,Ak,f_Ak,pmRA_1,e_pmRA_1,
     +  pmDE,e_pmDE,plx_1,e_plx_1,RVel,e_RVel,f_RVel,r_RVel,U,e_U,V,
     +  e_V,W,e_W,Gmag_1,e_Gmag_1,GBPmag_1,e_GBPmag_1,GRPmag_1,
     +  e_GRPmag_1,RUWE_1,Jmag_1,eJmag,Hmag_1,eHmag,Ksmag_1,eKsmag,
     +  r_JHKmag_1,Pops,memb_
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table5.dat'	! Mid-IR photometry for adopted Members of Upper Sco

C  Format for file interpretation

    5 format(
     +  A18,1X,A20,1X,A19,1X,I19,1X,A54,1X,F10.6,1X,F10.6,1X,A6,1X,
     +  F5.2,1X,F4.2,1X,A5,1X,F5.2,1X,F4.2,1X,A5,1X,F5.2,1X,F4.2,1X,
     +  A5,1X,F5.2,1X,F4.2,1X,A5,1X,F5.2,1X,F4.2,1X,A11,1X,F5.2,1X,
     +  F4.2,1X,A5,1X,F5.2,1X,F4.2,1X,A5,1X,F5.2,1X,F4.2,1X,A11,1X,
     +  F5.2,1X,F4.2,1X,A11,1X,A2,1X,A5,1X,A2,1X,A2,1X,A2,1X,A2,1X,A16)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table5.dat')
      write(6,*) '....Loading file: table5.dat'
      do i__=1,1761
        read(1,'(A372)')ar__4
        read(ar__4,5)
     +  v2MASS_2,WISEA_2,UGCS_2,Gaia_2,Name_2,RAdeg_3,DEdeg_3,
     +  r_RAdeg_2,v3_6mag,e_3_6mag,f_3_6mag,v4_5mag,e_4_5mag,f_4_5mag,
     +  v5_8mag,e_5_8mag,f_5_8mag,v8_0mag,e_8_0mag,f_8_0mag,v24mag,
     +  e_24mag,f_24mag,W1mag,e_W1mag,f_W1mag,W2mag,e_W2mag,f_W2mag,
     +  W3mag,e_W3mag,f_W3mag,W4mag,e_W4mag,f_W4mag,Exc4_5,Exc8_0,
     +  Exc24,ExcW2,ExcW3,ExcW4,DiskType
        if(ar__4(61:79) .EQ. '') Gaia_2 = iNULL__
        if(ar__4(165:169) .EQ. '') v3_6mag = rNULL__
        if(ar__4(171:174) .EQ. '') e_3_6mag = rNULL__
        if(ar__4(182:186) .EQ. '') v4_5mag = rNULL__
        if(ar__4(188:191) .EQ. '') e_4_5mag = rNULL__
        if(ar__4(199:203) .EQ. '') v5_8mag = rNULL__
        if(ar__4(205:208) .EQ. '') e_5_8mag = rNULL__
        if(ar__4(216:220) .EQ. '') v8_0mag = rNULL__
        if(ar__4(222:225) .EQ. '') e_8_0mag = rNULL__
        if(ar__4(233:237) .EQ. '') v24mag = rNULL__
        if(ar__4(239:242) .EQ. '') e_24mag = rNULL__
        if(ar__4(256:260) .EQ. '') W1mag = rNULL__
        if(ar__4(262:265) .EQ. '') e_W1mag = rNULL__
        if(ar__4(273:277) .EQ. '') W2mag = rNULL__
        if(ar__4(279:282) .EQ. '') e_W2mag = rNULL__
        if(ar__4(290:294) .EQ. '') W3mag = rNULL__
        if(ar__4(296:299) .EQ. '') e_W3mag = rNULL__
        if(ar__4(313:317) .EQ. '') W4mag = rNULL__
        if(ar__4(319:322) .EQ. '') e_W4mag = rNULL__
c    ..............Just test output...........
        write(6,5)
     +  v2MASS_2,WISEA_2,UGCS_2,Gaia_2,Name_2,RAdeg_3,DEdeg_3,
     +  r_RAdeg_2,v3_6mag,e_3_6mag,f_3_6mag,v4_5mag,e_4_5mag,f_4_5mag,
     +  v5_8mag,e_5_8mag,f_5_8mag,v8_0mag,e_8_0mag,f_8_0mag,v24mag,
     +  e_24mag,f_24mag,W1mag,e_W1mag,f_W1mag,W2mag,e_W2mag,f_W2mag,
     +  W3mag,e_W3mag,f_W3mag,W4mag,e_W4mag,f_W4mag,Exc4_5,Exc8_0,
     +  Exc24,ExcW2,ExcW3,ExcW4,DiskType
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'tablec1.dat'	! Intrinsic colors of young stars and brown dwarfs

C  Format for file interpretation

    6 format(
     +  A2,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,
     +  F5.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3,1X,F6.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'tablec1.dat')
      write(6,*) '....Loading file: tablec1.dat'
      do i__=1,59
        read(1,'(A99)')ar__5
        read(ar__5,6)
     +  SpType_3,BP_RP,G_RP,RP_K,J_H,H_K,K_3_6,K_4_5,K_5_8,K_8_0,K_24,
     +  K_W1,K_W2,K_W3,K_W4
        if(ar__5(4:9) .EQ. '') BP_RP = rNULL__
        if(ar__5(11:16) .EQ. '') G_RP = rNULL__
        if(ar__5(18:23) .EQ. '') RP_K = rNULL__
        if(ar__5(39:44) .EQ. '') K_3_6 = rNULL__
        if(ar__5(66:71) .EQ. '') K_24 = rNULL__
        if(ar__5(73:78) .EQ. '') K_W1 = rNULL__
        if(ar__5(87:92) .EQ. '') K_W3 = rNULL__
        if(ar__5(94:99) .EQ. '') K_W4 = rNULL__
c    ..............Just test output...........
        write(6,6)
     +  SpType_3,BP_RP,G_RP,RP_K,J_H,H_K,K_3_6,K_4_5,K_5_8,K_8_0,K_24,
     +  K_W1,K_W2,K_W3,K_W4
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'ref.dat'	! References used for spectral types

C  Format for file interpretation

    7 format(I2,1X,A21,1X,A19,1X,A63)

C  Effective file loading

      open(unit=1,status='old',file=
     +'ref.dat')
      write(6,*) '....Loading file: ref.dat'
      do i__=1,93
        read(1,'(A108)')ar__6
        read(ar__6,7)Ref,Autor,BIB,Cat
c    ..............Just test output...........
        write(6,7)Ref,Autor,BIB,Cat
c    .......End.of.Just test output...........
      end do
      close(1)

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