FORTRAN Generation
(/./ftp/cats/J/ApJ/689/1295)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/689/1295 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-Apr-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/689/1295   Lithium test implications for brown dwarfs (Kirkpatrick+, 2008)
*================================================================================
*A sample of very young field L dwarfs and implications for the brown dwarf
*"Lithium test" at early ages.
*    Kirkpatrick J.D., Cruz K.L., Barman T.S., Burgasser A.J., Looper D.L.,
*    Tinney C.G., Gelino C.R., Lowrance P.J., Liebert J., Carpenter J.M.,
*    Hillenbrand L.A., Stauffer J.R.
*   <Astrophys. J., 689, 1295-1326 (2008)>
*   =2008ApJ...689.1295K
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! L dwarf candidates sample from 2MASS

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

      character*2   S          (nr__) ! [Ld/OT] Sample (1)
      character*17  v2MASS     (nr__) ! 2MASS name (JHHMMSSss+DDMMSSs)
      real*4        Jmag       (nr__) ! (mag) 2MASS J-band magnitude
      real*4        e_Jmag     (nr__) ! (mag) Jmag uncertainty
      real*4        Hmag       (nr__) ! (mag) 2MASS H-band magnitude
      real*4        e_Hmag     (nr__) ! (mag) Hmag uncertainty
      real*4        Ksmag      (nr__) ! (mag) 2MASS Ks-band magnitude
      real*4        e_Ksmag    (nr__) ! (mag) Ksmag uncertainty
      real*4        J_Ks       (nr__) ! (mag) J-Ks color index
      real*4        e_J_Ks     (nr__) ! (mag) J-Ks uncertainty
      character*8   SpType     (nr__) ! MK optical spectral type (for Ld sample)
      character*1   f_SpType   (nr__) ! [c] Flag for 2MASS J15394189-0520428 (2)
      character*11  Type       (nr__) ! Object type (for OT sample)
      integer*4     r_Type     (nr__) ! ? Type reference; see refs.dat file
      character*6   Ref        (nr__) ! Discovery reference(s); see refs.dat file
      character*29  OName      (nr__) ! Other designation(s) or comment
      character*1   f_OName    (nr__) ! [b] Flag for 2MASS J14413716-0945590 (3)
*Note (1): Objects verified spectroscopically to be late-M or L dwarfs are
*     "Ld"; those objects known not to be late-M or L dwarfs are "OT".
*Note (2): DENIS J1539-0520 is typed in the optical as L3.5 by Reid et al.
*     (2008, Cat. J/AJ/136/1290).
*Note (3): Common proper motion for DENIS-P J1441-0945 and G 124-62A
*     confirmed by Seifahrt et al. (2005A&A...440..967S).

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

C  Declarations for 'table3.dat'	! New spectroscopic observations of late-M or L
                             dwarfs from table 1

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

      character*27  Name       (nr__1) ! Object designation
      character*1   f_Name     (nr__1) ! [e] Flag for Gl 337CD (1)
      integer*4     Ref_1      (nr__1) ! ? Discovery reference (only for table 4);
*                                    see refs.dat file
      character*10  ObsDate    (nr__1) ! ("YYYY/MM/DD") Observation date (UT)
      character*1   f_ObsDate  (nr__1) ! [m] for multidates (2)
      integer*4     tint       (nr__1) ! (s) Integration time
      real*4        CrH_a      (nr__1) ! CrH-a spectral index (3)
      character*1   l_CrT      (nr__1) ! Limit flag on CrT
      character*3   CrT        (nr__1) ! Spectral type from CrHa line
      real*4        Rb_TiO     (nr__1) ! ? Rb-b/TiO-b spectral index (3)
      character*1   l_RbT      (nr__1) ! Limit flag on RbT
      character*3   RbT        (nr__1) ! Spectral type from Rb/TiO
      real*4        Cs_VO      (nr__1) ! ? Cs-a/VO-b spectral index (3)
      character*1   l_CsT      (nr__1) ! Limit flag on CsT
      character*3   CsT        (nr__1) ! Spectral type from Cs/VO
      real*4        Col        (nr__1) ! ? Color-d spectral index (3)
      character*1   l_ColT     (nr__1) ! Limit flag on ColT
      character*3   ColT       (nr__1) ! Spectral type from col
      character*1   l_KI       (nr__1) ! [<~ ] Limit flag on KI
      character*3   KI         (nr__1) ! KI fit spectral type
      character*3   Ox         (nr__1) ! Oxyde fit spectral type
      character*6   SpType_1   (nr__1) ! Final optical spectral type
      character*1   f_SpType_1 (nr__1) ! [b] Flag for 2MASS J00145575-4844171 (4)
      integer*4     Dist       (nr__1) ! (pc) ? Distance (only for table 3)
      character*1   u_Dist     (nr__1) ! Uncertainty flag on Dist
      character*1   f_Dist     (nr__1) ! [p] Dist from trigonometric parallax (5)
      character*2   l_W_Ha     (nr__1) ! [<= ] Limit flag on W(Ha)
      real*4        W_Ha       (nr__1) ! (0.1nm) H{alpha} equivalent width
      character*1   f_W_Ha     (nr__1) ! [cd] Flag on W(Ha) (6)
      character*2   l_W_Li     (nr__1) ! [<= ] Limit flag on W(Li)
      real*4        W_Li       (nr__1) ! (0.1nm) LiI equivalent width
      character*1   u_W_Li     (nr__1) ! Uncertainty flag on W(Li)
*Note (1): Gl 337CD is also known as 2MASS J09121469+1459396 and was identified
*          as a close double by Burgasser et al. (2005AJ....129.2849B).
*Note (2): SDSSp J0423-0414 (table 3) was observed on four separate dates: 2001
*          November 13 (900s integration), 2002 January 1 (1200s integration),
*          2002 January 2 (1200s integration), and 2002 January 3 (1800s
*          integration). All analysis was performed on a co-added spectrum
*          combining all 5100s of integration.
*          SDSS J0830+4828 (table 4) was observed on two separate dates: 2002
*          January 1 (1200s integration) and 2002 January 2 (1200s integration).
*          All analysis was performed on a co-added spectrum combining all 2400s
*          of integration.
*Note (3): Primary standards and spectral indices are defined in Kirkpatrick et
*          al. (1999ApJ...519..802K). The indices are defined as flux ratios in
*          spectral ranges (wavelength in {AA}):
* CrH-a (CrH {lambda}8611)   = F(8580-8600)/F(8621-8641)
*  Rb-b (RbI {lambda}7947.6) =
*                        ((F(7922.6-7932.6)+F(7962.6-7972.6))/2/(F(7942.6-7952.6)
* TiO-b (TiO {lambda}8432)   = F(8400-8415)/F(8435-8470)
*  Cs-a (CsI {lambda}8521.1) =
*                         ((F(8496.1-8506.1)+F(8536.1-8546.1))/2/F(8516.1-8526.1)
*  VO-b (VO {lambda}7912)    = ((F(7860-7880)+F(8080-8100))/F(7960-8000)
*  Color-d = F(9675-9875)/F(7350-7550)
*Note (4): The optical spectrum of 2MASS J0014-4844 is slightly peculiar, perhaps
*          indicative of lower metallicity.
*Note (5):
*     p = distances as measured through trignonometric parallax. All other
*     entries are distances derived from spectrophotometric parallaxes. For
*     SDSSp J0423-0414, the trigonometric parallax is from Vrba et al.
*     (2004AJ....127.2948V); for Gl 618.1B the distance is from its
*     association with the primary, Gl 618.1A, whose trigonometric parallax
*     was measured by Hipparcos (Perryman et al. 1997, Cat. I/239).
*Note (6): Flag as follows:
*    c = 2MASS J0144-0716 shows variable H{alpha} emission when spectra from
*        different epochs are compared. See Liebert et al. (2003AJ....125..343L)
*        for further discussion.
*    d = 2MASS J0318-3421 shows a possible absorption trough of ~10{AA} EW
*        at the location of H{alpha}, but this is likely a data artifact
*       since it is aphysically broad.

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

C  Declarations for 'table4.dat'	! Observations of L dwarfs from litterature sample

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

      character*27  Name_1     (nr__2) ! Object designation
      character*1   f_Name_1   (nr__2) ! [e] Flag for Gl 337CD (1)
      integer*4     Ref_2      (nr__2) ! ? Discovery reference (only for table 4);
*                                    see refs.dat file
      character*10  ObsDate_1  (nr__2) ! ("YYYY/MM/DD") Observation date (UT)
      character*1   f_ObsDate_1(nr__2) ! [m] for multidates (2)
      integer*4     tint_1     (nr__2) ! (s) Integration time
      real*4        CrH_a_1    (nr__2) ! CrH-a spectral index (3)
      character*1   l_CrT_1    (nr__2) ! Limit flag on CrT
      character*3   CrT_1      (nr__2) ! Spectral type from CrHa line
      real*4        Rb_TiO_1   (nr__2) ! ? Rb-b/TiO-b spectral index (3)
      character*1   l_RbT_1    (nr__2) ! Limit flag on RbT
      character*3   RbT_1      (nr__2) ! Spectral type from Rb/TiO
      real*4        Cs_VO_1    (nr__2) ! ? Cs-a/VO-b spectral index (3)
      character*1   l_CsT_1    (nr__2) ! Limit flag on CsT
      character*3   CsT_1      (nr__2) ! Spectral type from Cs/VO
      real*4        Col_1      (nr__2) ! ? Color-d spectral index (3)
      character*1   l_ColT_1   (nr__2) ! Limit flag on ColT
      character*3   ColT_1     (nr__2) ! Spectral type from col
      character*1   l_KI_1     (nr__2) ! [<~ ] Limit flag on KI
      character*3   KI_1       (nr__2) ! KI fit spectral type
      character*3   Ox_1       (nr__2) ! Oxyde fit spectral type
      character*6   SpType_2   (nr__2) ! Final optical spectral type
      character*1   f_SpType_2 (nr__2) ! [b] Flag for 2MASS J00145575-4844171 (4)
      integer*4     Dist_1     (nr__2) ! (pc) ? Distance (only for table 3)
      character*1   u_Dist_1   (nr__2) ! Uncertainty flag on Dist
      character*1   f_Dist_1   (nr__2) ! [p] Dist from trigonometric parallax (5)
      character*2   l_W_Ha_1   (nr__2) ! [<= ] Limit flag on W(Ha)
      real*4        W_Ha_1     (nr__2) ! (0.1nm) H{alpha} equivalent width
      character*1   f_W_Ha_1   (nr__2) ! [cd] Flag on W(Ha) (6)
      character*2   l_W_Li_1   (nr__2) ! [<= ] Limit flag on W(Li)
      real*4        W_Li_1     (nr__2) ! (0.1nm) LiI equivalent width
      character*1   u_W_Li_1   (nr__2) ! Uncertainty flag on W(Li)
*Note (1): Gl 337CD is also known as 2MASS J09121469+1459396 and was identified
*          as a close double by Burgasser et al. (2005AJ....129.2849B).
*Note (2): SDSSp J0423-0414 (table 3) was observed on four separate dates: 2001
*          November 13 (900s integration), 2002 January 1 (1200s integration),
*          2002 January 2 (1200s integration), and 2002 January 3 (1800s
*          integration). All analysis was performed on a co-added spectrum
*          combining all 5100s of integration.
*          SDSS J0830+4828 (table 4) was observed on two separate dates: 2002
*          January 1 (1200s integration) and 2002 January 2 (1200s integration).
*          All analysis was performed on a co-added spectrum combining all 2400s
*          of integration.
*Note (3): Primary standards and spectral indices are defined in Kirkpatrick et
*          al. (1999ApJ...519..802K). The indices are defined as flux ratios in
*          spectral ranges (wavelength in {AA}):
* CrH-a (CrH {lambda}8611)   = F(8580-8600)/F(8621-8641)
*  Rb-b (RbI {lambda}7947.6) =
*                        ((F(7922.6-7932.6)+F(7962.6-7972.6))/2/(F(7942.6-7952.6)
* TiO-b (TiO {lambda}8432)   = F(8400-8415)/F(8435-8470)
*  Cs-a (CsI {lambda}8521.1) =
*                         ((F(8496.1-8506.1)+F(8536.1-8546.1))/2/F(8516.1-8526.1)
*  VO-b (VO {lambda}7912)    = ((F(7860-7880)+F(8080-8100))/F(7960-8000)
*  Color-d = F(9675-9875)/F(7350-7550)
*Note (4): The optical spectrum of 2MASS J0014-4844 is slightly peculiar, perhaps
*          indicative of lower metallicity.
*Note (5):
*     p = distances as measured through trignonometric parallax. All other
*     entries are distances derived from spectrophotometric parallaxes. For
*     SDSSp J0423-0414, the trigonometric parallax is from Vrba et al.
*     (2004AJ....127.2948V); for Gl 618.1B the distance is from its
*     association with the primary, Gl 618.1A, whose trigonometric parallax
*     was measured by Hipparcos (Perryman et al. 1997, Cat. I/239).
*Note (6): Flag as follows:
*    c = 2MASS J0144-0716 shows variable H{alpha} emission when spectra from
*        different epochs are compared. See Liebert et al. (2003AJ....125..343L)
*        for further discussion.
*    d = 2MASS J0318-3421 shows a possible absorption trough of ~10{AA} EW
*        at the location of H{alpha}, but this is likely a data artifact
*       since it is aphysically broad.

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

C  Declarations for 'table5.dat'	! Observations of T dwarfs from litterature sample

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

      character*25  Name_2     (nr__3) ! Object designation
      integer*4     Ref_3      (nr__3) ! Discovery reference; see refs.dat file
      character*10  ObsDate_2  (nr__3) ! ("YYYY/MM/DD") Observation date (UT)
      character*1   f_ObsDate_2(nr__3) ! [m] for multidates (1)
      integer*4     tint_2     (nr__3) ! (s) Integration time
      real*4        CsI        (nr__3) ! CsI(A) spectral index (2)
      character*1   l_CsT_2    (nr__3) ! Limit flag on CsT
      character*3   CsT_2      (nr__3) ! Spectral type from CsI
      real*4        CrH_H2O    (nr__3) ! CrH(A)/H_2_O spectral index (2)
      character*1   l_CrT_2    (nr__3) ! Limit flag on CrT
      character*3   CrT_2      (nr__3) ! Spectral type from CrH/H20
      real*4        FeH        (nr__3) ! FeH(B) spectral type (2)
      character*1   u_FeT      (nr__3) ! [~ ] Approximative flag on FeT
      integer*4     FeT        (nr__3) ! ? Spectral type from FeH
      real*4        Col_2      (nr__3) ! Color-e spectral index (2)
      character*1   l_ColT_2   (nr__3) ! Limit flag on ColT
      integer*4     ColT_2     (nr__3) ! ? Spectral type from Col
      character*5   Type_1     (nr__3) ! By-eye type
      character*6   SpType_3   (nr__3) ! Final optical spectral type
      character*1   l_W_Ha_1   (nr__3) ! Limit flag on W(Ha)
      integer*4     W_Ha_1     (nr__3) ! (0.1nm) ? H{alpha} equivalent width
      character*1   f_W_Ha_1   (nr__3) ! [c] Non-significant equivalent widths (3)
      character*1   l_W_Li_1   (nr__3) ! Limit flag on W(Li)
      integer*4     W_Li_1     (nr__3) ! (0.1nm) ? LiI equivalent width
*Note (1): SDSS J0837-0000 was observed on three separate dates:
*     2000 December 26 (7200s integration), 2000 December 27 (3600s
*     integration), and 2000 December 28 (3600s integration). All analysis
*     was performed on a co-added spectrum combining all 14400s of
*     integration. 2MASS J1209-1004 was observed on three separate dates:
*     2003 December 22 (2400s integration), 2003 December 23 (3600s
*     integration), and 2003 December 24 (3600s integration). All analysis
*     was performed on a co-added spectrum combining all 9600s of integration.
*Note (2): Indices of T dwarf spectral ratios in the red optical are defined in
*          Burgasser et al. 2003ApJ...594..510B:
*   CsI(A) (8521{AA}) = <F8496.1-8506.1>+<F8536.1-8546.1>
*   CrH(A) (8611{AA}) = {int}F8560-8600
*   H_2_O  (9250{AA}) = {int}F9220-9240
*   FeH(B) (9896{AA}) = {int}F9855-9885
*   Color-e (spectral slope) = <F9140-9240>
*Note (3): 2MASS J0415-0935 has a flux measurement at or very close to zero
*     at these wavelengths, so meaningful equivalent width measures are not
*     possible.

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

C  Declarations for 'table6.dat'	! *Ultracool field dwarfs (>=M7) with spectroscopic
                             signatures of low gravity (Youth)

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

      character*24  Name_3     (nr__4) ! Object name
      character*1   f_Name_2   (nr__4) ! [d] SSSPM J1102-3431 = 2MASS J11020983-3430355
      character*8   SpType_4   (nr__4) ! Optical spectral type
      character*1   f_SpType_3 (nr__4) ! [e] Flag for 2MASSI J1615425+495321 (1)
      integer*4     r_SpType   (nr__4) ! ? Optical reference; see refs.dat file
      character*10  IRSpT      (nr__4) ! Near-IR spectral type
      integer*4     r_IRSpT    (nr__4) ! ? Near-IR reference; see refs.dat file
      real*4        J_Ks_1     (nr__4) ! (mag) ? 2MASS J-Ks color index
      real*4        e_J_Ks_1   (nr__4) ! (mag) ? J-Ks uncertainty
      character*1   f_J_Ks     (nr__4) ! [c] SDSS J102552.43+321234.0 not detected at
*                                  2MASS J-band
      integer*4     Ref_4      (nr__4) ! Discovery reference; see refs.dat file
      character*4   Ref2       (nr__4) ! Youth reference(s); see refs.dat file
      character*45  Diag       (nr__4) ! Diagnostics
*Note (1): Spectrum has poor signal-to-noise ratio and may prove to be
*     dusty rather than low gravity.

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

C  Declarations for 'table7.dat'	! Log of additional spectroscopic observations for
                             objects shown in Figs. 7-13

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

      character*23  Name_4     (nr__5) ! Object name
      character*24  OName_1    (nr__5) ! Other designation
      character*10  Date       (nr__5) ! ("YYYY/MM/DD") Observation date (UT)
      character*1   f_Date     (nr__5) ! [b-e] Multidate (1)
      character*8   Tel        (nr__5) ! Telescope
      character*1   f_Tel      (nr__5) ! [f] Flag for Subaru (1)
      character*7   Inst       (nr__5) ! Instrument
      integer*4     tint_3     (nr__5) ! (s) Integration time
      integer*4     tint2      (nr__5) ! (s) ? Second integration time
      integer*4     r_SpType_1 (nr__5) ! ? Spectroscopic reference; see refs.dat file
      character*13  SpType_5   (nr__5) ! Spectral type
      character*4   Cor        (nr__5) ! Telluric corrected? (3)
      character*1   f_Cor      (nr__5) ! [c] Flag for G 196-3B (1)
*Note (1): Flags as follows:
*   b = 2MASS J0850+1057 was observed on five different dates; the plotted
*       spectrum is a sum of all five. Exposure times were 2400s on 1997 November
*       9 (UT), 1200s on 1997 December 7 (UT), 1200s on 1997 December 8 (UT),
*       1200s on 1997 December 9 (UT), and 4800s on 1998 January 24 (UT).
*   c = G 196-3B was observed using a blue setup (3900-8700{AA}) on one night and
*       our standard red setup (6300-10100{AA}) on two other nights. Observations
*       for the blue setup were done in superior conditions, so the plotted
*       spectrum shortward of 8500{AA} is from this set, and data longward of
*       8500{AA} come from the poorer data taken with the red setup. Data taken
*       with the blue setup were corrected for telluric absorption; data with the
*       red setup were not telluric corrected. The exposure time for the blue set
*       up was 3600s on 2001 February 19 (UT). Exposure times for the red setup
*       were 2400s on 1999 March 04 (UT) and 3600s on 1999 March 5 (UT).
*   d = 2MASS J1022+0200 was observed for 600s on two different dates, 2003 April
*       04 (UT) and 2006 January 15 (UT). The plotted spectrum is the sum of
*       these two observations.
*   e = Gl 417B was observed on three different dates; the plotted spectrum is a
*       sum of all three. Exposure times were 1200s on 1998 December 14 (UT),
*       1200s on 1999 March 4 (UT), and 3600s on 1999 March 5 (UT).
*   f = Observation and reduction procedures for Subaru-FOCAS data are
*       discussed in Looper et al. (2008ApJ...686..528L).
*Note (3): Indicates whether or not the spectrum was telluric corrected
*     using a G dwarf spectrum acquired near in time and near on sky to the
*     program object.

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

C  Declarations for 'table8.dat'	! Lithium equivalent width measurements for L0,
                             L2, L4, and L6.5 dwarfs

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

      character*4   SpT        (nr__6) ! Spectral type
      character*27  Name_5     (nr__6) ! Name
      character*1   f_Name_3   (nr__6) ! [a-g] Flag on Name (1)
      character*2   l_W_Li_1   (nr__6) ! [<= ] Limit flag on W(Li)
      real*4        W_Li_1     (nr__6) ! (0.1nm) Lithium equivalent width
      real*4        e_W_Li     (nr__6) ! (0.1nm) ? W(Li) uncertainty
      integer*4     f_SpT      (nr__6) ! [1-4] Particularity of dwarfs (2)
*Note (1): Flag as follows:
*   a = 2MASS J0141-4633: age estimate from fitting of theoretical models to
*       spectra is 1-50Myr (Kirkpatrick et al. 2006ApJ...639.1120K).
*   b = Kelu-1: age estimate from model fits is 300-800Myr (Liu & Leggett,
*       2005ApJ...634..616L).
*   c = G 196-3B: age estimates from the presence of the lithium in G 196-3B
*       itself as well as activity diagnostics in the primary is 20-300Myr
*       (Rebolo et al. 1998Sci...282.1309R); an age estimate using indicators
*       from the early-M primary alone is 60-300Myr (Kirkpatrick et al.
*       2001AJ....121.3235K).
*   d = Gl 417B: age estimate from the G dwarf primary is 80-300Myr (Kirkpatrick
*       et al. 2001AJ....121.3235K). Gl 417B is actually typed as L4.5 but is
*       included here as an L4 so that this bin of peculiar, low-gravity L4
*       dwarfs could be populated.
*   e = GD 165B: age estimate using the white dwarf primary is 1.2-5.5Gyr
*       (Kirkpatrick et al. 1999ApJ...519..834K).
*   f = 2MASS J2244+2043: This object has a normal L6.5 optical spectrum, but its
*       near-infrared spectrum is extremely peculiar. In this case, the
*       peculiarity may be attributable to low gravity (McLean et al.
*       2003ApJ...596..561M).
*   g = 2MASS J2148+4003: This object has a normal L6.5 optical spectrum, but its
*       near-infrared spectrum is extremely peculiar. In this case, the
*       peculiarity is best explained by excessive atmospheric dust (Looper et
*       al. 2008ApJ...686..528L).
*Note (2): Dwarf type as follows:
*   1 = Normal dwarf without Lithium
*   2 = Normal dwarf with Lithium
*   3 = Peculiar dwarf (low-g)
*   4 = Very peculiar dwarf (very low-g)

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

C  Declarations for 'refs.dat'	! References

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

      integer*4     Ref_5      (nr__7) ! Reference number
      character*19  BibCode    (nr__7) ! Bibcode
      character*18  Aut        (nr__7) ! Author's name
      character*38  Com        (nr__7) ! Comment

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

C  Loading file 'table1.dat'	! L dwarf candidates sample from 2MASS

C  Format for file interpretation

    1 format(
     +  A2,7X,A17,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,
     +  F5.3,1X,F5.3,1X,A8,A1,1X,A11,1X,I2,1X,A6,1X,A29,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,68
        read(1,'(A140)')ar__
        read(ar__,1)
     +  S(i__),v2MASS(i__),Jmag(i__),e_Jmag(i__),Hmag(i__),
     +  e_Hmag(i__),Ksmag(i__),e_Ksmag(i__),J_Ks(i__),e_J_Ks(i__),
     +  SpType(i__),f_SpType(i__),Type(i__),r_Type(i__),Ref(i__),
     +  OName(i__),f_OName(i__)
        if(ar__(101:102) .EQ. '') r_Type(i__) = iNULL__
c    ..............Just test output...........
        write(6,1)
     +  S(i__),v2MASS(i__),Jmag(i__),e_Jmag(i__),Hmag(i__),
     +  e_Hmag(i__),Ksmag(i__),e_Ksmag(i__),J_Ks(i__),e_J_Ks(i__),
     +  SpType(i__),f_SpType(i__),Type(i__),r_Type(i__),Ref(i__),
     +  OName(i__),f_OName(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table3.dat'	! New spectroscopic observations of late-M or L
*                             dwarfs from table 1

C  Format for file interpretation

    2 format(
     +  A27,A1,1X,I2,1X,A10,A1,1X,I4,1X,F4.2,1X,A1,A3,1X,F4.2,1X,A1,
     +  A3,1X,F5.2,1X,A1,A3,1X,F5.2,1X,A1,A3,1X,A1,A3,1X,A3,1X,A6,A1,
     +  1X,I2,A1,A1,1X,A2,F5.1,A1,1X,A2,F4.1,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,28
        read(1,'(A129)')ar__1
        read(ar__1,2)
     +  Name(i__),f_Name(i__),Ref_1(i__),ObsDate(i__),f_ObsDate(i__),
     +  tint(i__),CrH_a(i__),l_CrT(i__),CrT(i__),Rb_TiO(i__),
     +  l_RbT(i__),RbT(i__),Cs_VO(i__),l_CsT(i__),CsT(i__),Col(i__),
     +  l_ColT(i__),ColT(i__),l_KI(i__),KI(i__),Ox(i__),SpType_1(i__),
     +  f_SpType_1(i__),Dist(i__),u_Dist(i__),f_Dist(i__),l_W_Ha(i__),
     +  W_Ha(i__),f_W_Ha(i__),l_W_Li(i__),W_Li(i__),u_W_Li(i__)
        if(ar__1(30:31) .EQ. '') Ref_1(i__) = iNULL__
        if(ar__1(60:63) .EQ. '') Rb_TiO(i__) = rNULL__
        if(ar__1(70:74) .EQ. '') Cs_VO(i__) = rNULL__
        if(ar__1(81:85) .EQ. '') Col(i__) = rNULL__
        if(ar__1(109:110) .EQ. '') Dist(i__) = iNULL__
c    ..............Just test output...........
        write(6,2)
     +  Name(i__),f_Name(i__),Ref_1(i__),ObsDate(i__),f_ObsDate(i__),
     +  tint(i__),CrH_a(i__),l_CrT(i__),CrT(i__),Rb_TiO(i__),
     +  l_RbT(i__),RbT(i__),Cs_VO(i__),l_CsT(i__),CsT(i__),Col(i__),
     +  l_ColT(i__),ColT(i__),l_KI(i__),KI(i__),Ox(i__),SpType_1(i__),
     +  f_SpType_1(i__),Dist(i__),u_Dist(i__),f_Dist(i__),l_W_Ha(i__),
     +  W_Ha(i__),f_W_Ha(i__),l_W_Li(i__),W_Li(i__),u_W_Li(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table4.dat'	! Observations of L dwarfs from litterature sample

C  Format for file interpretation

    3 format(
     +  A27,A1,1X,I2,1X,A10,A1,1X,I4,1X,F4.2,1X,A1,A3,1X,F4.2,1X,A1,
     +  A3,1X,F5.2,1X,A1,A3,1X,F5.2,1X,A1,A3,1X,A1,A3,1X,A3,1X,A6,A1,
     +  1X,I2,A1,A1,1X,A2,F5.1,A1,1X,A2,F4.1,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table4.dat')
      write(6,*) '....Loading file: table4.dat'
      do i__=1,6
        read(1,'(A129)')ar__2
        read(ar__2,3)
     +  Name_1(i__),f_Name_1(i__),Ref_2(i__),ObsDate_1(i__),
     +  f_ObsDate_1(i__),tint_1(i__),CrH_a_1(i__),l_CrT_1(i__),
     +  CrT_1(i__),Rb_TiO_1(i__),l_RbT_1(i__),RbT_1(i__),Cs_VO_1(i__),
     +  l_CsT_1(i__),CsT_1(i__),Col_1(i__),l_ColT_1(i__),ColT_1(i__),
     +  l_KI_1(i__),KI_1(i__),Ox_1(i__),SpType_2(i__),f_SpType_2(i__),
     +  Dist_1(i__),u_Dist_1(i__),f_Dist_1(i__),l_W_Ha_1(i__),
     +  W_Ha_1(i__),f_W_Ha_1(i__),l_W_Li_1(i__),W_Li_1(i__),
     +  u_W_Li_1(i__)
        if(ar__2(30:31) .EQ. '') Ref_2(i__) = iNULL__
        if(ar__2(60:63) .EQ. '') Rb_TiO_1(i__) = rNULL__
        if(ar__2(70:74) .EQ. '') Cs_VO_1(i__) = rNULL__
        if(ar__2(81:85) .EQ. '') Col_1(i__) = rNULL__
        if(ar__2(109:110) .EQ. '') Dist_1(i__) = iNULL__
c    ..............Just test output...........
        write(6,3)
     +  Name_1(i__),f_Name_1(i__),Ref_2(i__),ObsDate_1(i__),
     +  f_ObsDate_1(i__),tint_1(i__),CrH_a_1(i__),l_CrT_1(i__),
     +  CrT_1(i__),Rb_TiO_1(i__),l_RbT_1(i__),RbT_1(i__),Cs_VO_1(i__),
     +  l_CsT_1(i__),CsT_1(i__),Col_1(i__),l_ColT_1(i__),ColT_1(i__),
     +  l_KI_1(i__),KI_1(i__),Ox_1(i__),SpType_2(i__),f_SpType_2(i__),
     +  Dist_1(i__),u_Dist_1(i__),f_Dist_1(i__),l_W_Ha_1(i__),
     +  W_Ha_1(i__),f_W_Ha_1(i__),l_W_Li_1(i__),W_Li_1(i__),
     +  u_W_Li_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table5.dat'	! Observations of T dwarfs from litterature sample

C  Format for file interpretation

    4 format(
     +  A25,1X,I2,1X,A10,A1,1X,I5,1X,F4.2,1X,A1,A3,1X,F4.2,1X,A1,A3,
     +  1X,F4.2,1X,A1,I1,1X,F4.2,1X,A1,I1,1X,A5,1X,A6,1X,A1,I2,A1,1X,
     +  A1,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table5.dat')
      write(6,*) '....Loading file: table5.dat'
      do i__=1,5
        read(1,'(A103)')ar__3
        read(ar__3,4)
     +  Name_2(i__),Ref_3(i__),ObsDate_2(i__),f_ObsDate_2(i__),
     +  tint_2(i__),CsI(i__),l_CsT_2(i__),CsT_2(i__),CrH_H2O(i__),
     +  l_CrT_2(i__),CrT_2(i__),FeH(i__),u_FeT(i__),FeT(i__),
     +  Col_2(i__),l_ColT_2(i__),ColT_2(i__),Type_1(i__),
     +  SpType_3(i__),l_W_Ha_1(i__),W_Ha_1(i__),f_W_Ha_1(i__),
     +  l_W_Li_1(i__),W_Li_1(i__)
        if(ar__3(74:74) .EQ. '') FeT(i__) = iNULL__
        if(ar__3(82:82) .EQ. '') ColT_2(i__) = iNULL__
        if(ar__3(98:99) .EQ. '') W_Ha_1(i__) = iNULL__
        if(ar__3(103:103) .EQ. '') W_Li_1(i__) = iNULL__
c    ..............Just test output...........
        write(6,4)
     +  Name_2(i__),Ref_3(i__),ObsDate_2(i__),f_ObsDate_2(i__),
     +  tint_2(i__),CsI(i__),l_CsT_2(i__),CsT_2(i__),CrH_H2O(i__),
     +  l_CrT_2(i__),CrT_2(i__),FeH(i__),u_FeT(i__),FeT(i__),
     +  Col_2(i__),l_ColT_2(i__),ColT_2(i__),Type_1(i__),
     +  SpType_3(i__),l_W_Ha_1(i__),W_Ha_1(i__),f_W_Ha_1(i__),
     +  l_W_Li_1(i__),W_Li_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table6.dat'	! *Ultracool field dwarfs (>=M7) with spectroscopic
*                             signatures of low gravity (Youth)

C  Format for file interpretation

    5 format(
     +  A24,A1,1X,A8,A1,1X,I2,1X,A10,1X,I2,1X,F4.2,1X,F4.2,A1,1X,I2,
     +  1X,A4,1X,A45)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table6.dat')
      write(6,*) '....Loading file: table6.dat'
      do i__=1,20
        read(1,'(A117)')ar__4
        read(ar__4,5)
     +  Name_3(i__),f_Name_2(i__),SpType_4(i__),f_SpType_3(i__),
     +  r_SpType(i__),IRSpT(i__),r_IRSpT(i__),J_Ks_1(i__),
     +  e_J_Ks_1(i__),f_J_Ks(i__),Ref_4(i__),Ref2(i__),Diag(i__)
        if(ar__4(37:38) .EQ. '') r_SpType(i__) = iNULL__
        if(ar__4(51:52) .EQ. '') r_IRSpT(i__) = iNULL__
        if(ar__4(54:57) .EQ. '') J_Ks_1(i__) = rNULL__
        if(ar__4(59:62) .EQ. '') e_J_Ks_1(i__) = rNULL__
c    ..............Just test output...........
        write(6,5)
     +  Name_3(i__),f_Name_2(i__),SpType_4(i__),f_SpType_3(i__),
     +  r_SpType(i__),IRSpT(i__),r_IRSpT(i__),J_Ks_1(i__),
     +  e_J_Ks_1(i__),f_J_Ks(i__),Ref_4(i__),Ref2(i__),Diag(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table7.dat'	! Log of additional spectroscopic observations for
*                             objects shown in Figs. 7-13

C  Format for file interpretation

    6 format(
     +  A23,1X,A24,1X,A10,1X,A1,1X,A8,A1,1X,A7,1X,I5,1X,I4,1X,I2,1X,
     +  A13,1X,A4,1X,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table7.dat')
      write(6,*) '....Loading file: table7.dat'
      do i__=1,31
        read(1,'(A114)')ar__5
        read(ar__5,6)
     +  Name_4(i__),OName_1(i__),Date(i__),f_Date(i__),Tel(i__),
     +  f_Tel(i__),Inst(i__),tint_3(i__),tint2(i__),r_SpType_1(i__),
     +  SpType_5(i__),Cor(i__),f_Cor(i__)
        if(ar__5(87:90) .EQ. '') tint2(i__) = iNULL__
        if(ar__5(92:93) .EQ. '') r_SpType_1(i__) = iNULL__
c    ..............Just test output...........
        write(6,6)
     +  Name_4(i__),OName_1(i__),Date(i__),f_Date(i__),Tel(i__),
     +  f_Tel(i__),Inst(i__),tint_3(i__),tint2(i__),r_SpType_1(i__),
     +  SpType_5(i__),Cor(i__),f_Cor(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table8.dat'	! Lithium equivalent width measurements for L0,
*                             L2, L4, and L6.5 dwarfs

C  Format for file interpretation

    7 format(A4,1X,A27,1X,A1,1X,A2,F4.1,1X,F3.1,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table8.dat')
      write(6,*) '....Loading file: table8.dat'
      do i__=1,41
        read(1,'(A47)')ar__6
        read(ar__6,7)
     +  SpT(i__),Name_5(i__),f_Name_3(i__),l_W_Li_1(i__),W_Li_1(i__),
     +  e_W_Li(i__),f_SpT(i__)
        if(ar__6(43:45) .EQ. '') e_W_Li(i__) = rNULL__
c    ..............Just test output...........
        write(6,7)
     +  SpT(i__),Name_5(i__),f_Name_3(i__),l_W_Li_1(i__),W_Li_1(i__),
     +  e_W_Li(i__),f_SpT(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'refs.dat'	! References

C  Format for file interpretation

    8 format(I2,1X,A19,1X,A18,1X,A38)

C  Effective file loading

      open(unit=1,status='old',file=
     +'refs.dat')
      write(6,*) '....Loading file: refs.dat'
      do i__=1,41
        read(1,'(A80)')ar__7
        read(ar__7,8)Ref_5(i__),BibCode(i__),Aut(i__),Com(i__)
c    ..............Just test output...........
        write(6,8)Ref_5(i__),BibCode(i__),Aut(i__),Com(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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