FORTRAN Generation
(/./ftp/cats/J/ApJS/221/24)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJS/221/24 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 2026-Sep-17
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/221/24    SDSS-III APOGEE H-band spectral line lists    (Shetrone+, 2015)
*================================================================================
*The SDSS-III APOGEE spectral line list for H-band spectroscopy.
*    Shetrone M., Bizyaev D., Lawler J.E., Allende Prieto C., Johnson J.A.,
*    Smith V.V., Cunha K., Holtzman J., Garcia Perez A.E., Meszaros SZ.,
*    Sobeck J., Zamora O., Garcia-Hernandez D.A., Souto D., Chojnowski D.,
*    Koesterke L., Majewski S., Zasowski G.
*   <Astrophys. J. Suppl. Ser., 221, 24 (2015)>
*   =2015ApJS..221...24S    (SIMBAD/NED BibCode)
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table6.dat'	! Missing lines

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

      real*8        lambda     (nr__) ! (0.1nm) [15174.1/16945.8] Vacuum wavelength in Angstroms
      integer*4     EW         (nr__) ! (10-13m) [10/219]? Equivalent width in milli-Angstroms
      character*1   f_EW       (nr__) ! Flag on EqWid (e=emission or N=N/A)
      character*10  Notes      (nr__) ! [ ,a-k] Additional notes (1)
*Note (1): Note as follows:
*    a = Seen in the cool APOGEE giant spectra residual;
*    b = Seen in the hot APOGEE spectra residual;
*    c = Between APOGEE detectors;
*    d = all FTS stars;
*    e = only M FTS stars;
*    f = found only in HD199799;
*    g = see Hasselquist et al. 2016 in prep.;
*    h = Seen in cool APOGEE dwarf spectra residuals
*    i = Seen in APOGEE Be spectra, see Chojnowski et al. (2015, J/AJ/149/7)
*    j = This may be a DIB feature, see Zasowski et al. (2015ApJ...798...35Z)
*    k = This may be a DIB feature, see Geballe et al. (2011Natur.479..200G).

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

C  Declarations for 'table7.dat'	! *ASPCAP DR12 (20131216) line list

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

      real*8        lambda_1   (nr__1) ! (nm) [1500.4/1700.5] Vacuum wavelength
      real*8        orggf      (nr__1) ! ([-]) [-16.8/1.2]? Original log(gf) value
      real*8        newgf      (nr__1) ! ([-]) [-9.5/1.1]? Improved literature or
*                              laboratory log(gf)
      real*4        e_newgf    (nr__1) ! ([-]) [0.01/1]? Error in newgf, when available
      character*3   r_newgf    (nr__1) ! Source for newgf
      real*8        astgf      (nr__1) ! ([-]) [-14.7/1.3]? Astrophysical log(gf)
      character*3   r_astgf    (nr__1) ! Source for astrogf
      real*8        specid     (nr__1) ! [1/608.18] Species identifier
      real*8        EP1        (nr__1) ! (cm-1) [0/771723.1] Lower Energy Level
      real*4        J1         (nr__1) ! [0/159] J value for EP1
      character*11  EP1id      (nr__1) ! EP1 level identification
      real*8        EP2        (nr__1) ! (cm-2) [6200.7/778194.1] Upper Energy Level
      real*4        J2         (nr__1) ! [0/160] J value for EP2
      character*11  EP2id      (nr__1) ! EP2 level identification
      real*4        Rad        (nr__1) ! [0/11.1]? Damping Rad
      real*4        Sta        (nr__1) ! [-8/0]? Damping Stark
      real*4        vdW        (nr__1) ! [-8.4/0]? Damping vdW
      integer*4     unlte      (nr__1) ! [0/99]? NLTE level number upper
      integer*4     lnlte      (nr__1) ! [0/99]? NLTE level number lower
      integer*4     iso1       (nr__1) ! [0/65]? First isotope number
      real*4        hyp        (nr__1) ! ([-]) [-9.999/0]? Hyperfine component log fractional
*                              strength
      integer*4     iso2       (nr__1) ! [0/65]? Second isotope number
      real*4        isof       (nr__1) ! ([-]) [-4/0]? Log isotopic abundance fraction
      integer*4     hE1        (nr__1) ! (mK) [-520/425]? Hyperfine shift for 1st level to be
*                              added to E1
      integer*4     hE2        (nr__1) ! (mK) [-603/469]? Hyperfine shift for 1st level to be
*                              added to E2
      character*1   F0         (nr__1) ! [F] Hyperfine F symbol
      integer*4     F1         (nr__1) ! [0/9]? Hyperfine F for the first level
      character*1   Note1      (nr__1) ! [?z] Note on character of hyperfine data
*                              for first level (z=none; ?=guessed)
      character*1   S          (nr__1) ! [-] The symbol "-" for legibility
      integer*4     F2         (nr__1) ! [0/9]? Hyperfine F' for the second level
      character*1   note2      (nr__1) ! [?z] Note on character of hyperfine data for
*                              second level (z=none; ?=guessed)
      integer*4     g1         (nr__1) ! [-720/3934]? Lande g for first level times 1000
      integer*4     g2         (nr__1) ! [-689/4006]? Lande g for second level times 1000
      character*3   vdWorg     (nr__1) ! Source for the original vdW damping
      real*4        vdWast     (nr__1) ! [-8.1/-5.2]? Astrophysical vdW damping

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

C  Declarations for 'table8.dat'	! *ASPCAP INT (20120216) line list

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

      real*8        lambda_2   (nr__2) ! (nm) [1500.4/1700.5] Vacuum wavelength
      real*8        orggf_1    (nr__2) ! ([-]) [-16.8/1.2]? Original log(gf) value
      real*8        newgf_1    (nr__2) ! ([-]) [-9.5/1.1]? Improved literature or
*                              laboratory log(gf)
      real*4        e_newgf_1  (nr__2) ! ([-]) [0.01/1]? Error in newgf, when available
      character*3   r_newgf_1  (nr__2) ! Source for newgf
      real*8        astgf_1    (nr__2) ! ([-]) [-14.7/1.3]? Astrophysical log(gf)
      character*3   r_astgf_1  (nr__2) ! Source for astrogf
      real*8        specid_1   (nr__2) ! [1/608.18] Species identifier
      real*8        EP1_1      (nr__2) ! (cm-1) [0/771723.1] Lower Energy Level
      real*4        J1_1       (nr__2) ! [0/159] J value for EP1
      character*11  EP1id_1    (nr__2) ! EP1 level identification
      real*8        EP2_1      (nr__2) ! (cm-2) [6200.7/778194.1] Upper Energy Level
      real*4        J2_1       (nr__2) ! [0/160] J value for EP2
      character*11  EP2id_1    (nr__2) ! EP2 level identification
      real*4        Rad_1      (nr__2) ! [0/11.1]? Damping Rad
      real*4        Sta_1      (nr__2) ! [-8/0]? Damping Stark
      real*4        vdW_1      (nr__2) ! [-8.4/0]? Damping vdW
      integer*4     unlte_1    (nr__2) ! [0/99]? NLTE level number upper
      integer*4     lnlte_1    (nr__2) ! [0/99]? NLTE level number lower
      integer*4     iso1_1     (nr__2) ! [0/65]? First isotope number
      real*4        hyp_1      (nr__2) ! ([-]) [-9.999/0]? Hyperfine component log fractional
*                              strength
      integer*4     iso2_1     (nr__2) ! [0/65]? Second isotope number
      real*4        isof_1     (nr__2) ! ([-]) [-4/0]? Log isotopic abundance fraction
      integer*4     hE1_1      (nr__2) ! (mK) [-520/425]? Hyperfine shift for 1st level to be
*                              added to E1
      integer*4     hE2_1      (nr__2) ! (mK) [-603/469]? Hyperfine shift for 1st level to be
*                              added to E2
      character*1   F0_1       (nr__2) ! [F] Hyperfine F symbol
      integer*4     F1_1       (nr__2) ! [0/9]? Hyperfine F for the first level
      character*1   Note1_1    (nr__2) ! [?z] Note on character of hyperfine data
*                              for first level (z=none; ?=guessed)
      character*1   S_1        (nr__2) ! [-] The symbol "-" for legibility
      integer*4     F2_1       (nr__2) ! [0/9]? Hyperfine F' for the second level
      character*1   note2_1    (nr__2) ! [?z] Note on character of hyperfine data for
*                              second level (z=none; ?=guessed)
      integer*4     g1_1       (nr__2) ! [-720/3934]? Lande g for first level times 1000
      integer*4     g2_1       (nr__2) ! [-689/4006]? Lande g for second level times 1000
      character*3   vdWorg_1   (nr__2) ! Source for the original vdW damping
      real*4        vdWast_1   (nr__2) ! [-8.1/-5.2]? Astrophysical vdW damping

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

C  Declarations for 'table9.dat'	! *ASPCAP DR10 (20110510) line list

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

      real*8        lambda_3   (nr__3) ! (nm) [1500.4/1700.5] Vacuum wavelength
      real*8        orggf_2    (nr__3) ! ([-]) [-16.8/1.2]? Original log(gf) value
      real*8        newgf_2    (nr__3) ! ([-]) [-9.5/1.1]? Improved literature or
*                              laboratory log(gf)
      real*4        e_newgf_2  (nr__3) ! ([-]) [0.01/1]? Error in newgf, when available
      character*3   r_newgf_2  (nr__3) ! Source for newgf
      real*8        astgf_2    (nr__3) ! ([-]) [-14.7/1.3]? Astrophysical log(gf)
      character*3   r_astgf_2  (nr__3) ! Source for astrogf
      real*8        specid_2   (nr__3) ! [1/608.18] Species identifier
      real*8        EP1_2      (nr__3) ! (cm-1) [0/771723.1] Lower Energy Level
      real*4        J1_2       (nr__3) ! [0/159] J value for EP1
      character*11  EP1id_2    (nr__3) ! EP1 level identification
      real*8        EP2_2      (nr__3) ! (cm-2) [6200.7/778194.1] Upper Energy Level
      real*4        J2_2       (nr__3) ! [0/160] J value for EP2
      character*11  EP2id_2    (nr__3) ! EP2 level identification
      real*4        Rad_2      (nr__3) ! [0/11.1]? Damping Rad
      real*4        Sta_2      (nr__3) ! [-8/0]? Damping Stark
      real*4        vdW_2      (nr__3) ! [-8.4/0]? Damping vdW
      integer*4     unlte_2    (nr__3) ! [0/99]? NLTE level number upper
      integer*4     lnlte_2    (nr__3) ! [0/99]? NLTE level number lower
      integer*4     iso1_2     (nr__3) ! [0/65]? First isotope number
      real*4        hyp_2      (nr__3) ! ([-]) [-9.999/0]? Hyperfine component log fractional
*                              strength
      integer*4     iso2_2     (nr__3) ! [0/65]? Second isotope number
      real*4        isof_2     (nr__3) ! ([-]) [-4/0]? Log isotopic abundance fraction
      integer*4     hE1_2      (nr__3) ! (mK) [-520/425]? Hyperfine shift for 1st level to be
*                              added to E1
      integer*4     hE2_2      (nr__3) ! (mK) [-603/469]? Hyperfine shift for 1st level to be
*                              added to E2
      character*1   F0_2       (nr__3) ! [F] Hyperfine F symbol
      integer*4     F1_2       (nr__3) ! [0/9]? Hyperfine F for the first level
      character*1   Note1_2    (nr__3) ! [?z] Note on character of hyperfine data
*                              for first level (z=none; ?=guessed)
      character*1   S_2        (nr__3) ! [-] The symbol "-" for legibility
      integer*4     F2_2       (nr__3) ! [0/9]? Hyperfine F' for the second level
      character*1   note2_2    (nr__3) ! [?z] Note on character of hyperfine data for
*                              second level (z=none; ?=guessed)
      integer*4     g1_2       (nr__3) ! [-720/3934]? Lande g for first level times 1000
      integer*4     g2_2       (nr__3) ! [-689/4006]? Lande g for second level times 1000
      character*3   vdWorg_2   (nr__3) ! Source for the original vdW damping
      real*4        vdWast_2   (nr__3) ! [-8.1/-5.2]? Astrophysical vdW damping

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

C  Loading file 'table6.dat'	! Missing lines

C  Format for file interpretation

    1 format(F7.1,1X,I3,A1,1X,A10)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table6.dat')
      write(6,*) '....Loading file: table6.dat'
      do i__=1,445
        read(1,'(A23)')ar__
        read(ar__,1)lambda(i__),EW(i__),f_EW(i__),Notes(i__)
        if(ar__(9:11) .EQ. '') EW(i__) = iNULL__
c    ..............Just test output...........
        write(6,1)lambda(i__),EW(i__),f_EW(i__),Notes(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table7.dat'	! *ASPCAP DR12 (20131216) line list

C  Format for file interpretation

    2 format(
     +  F9.4,1X,F7.3,1X,F7.3,1X,F4.2,1X,A3,F7.3,1X,A3,1X,F8.2,F12.3,
     +  F5.1,A11,F12.3,F5.1,A11,F6.2,F6.2,F6.2,1X,I2,I2,I3,F6.3,I3,
     +  F6.3,I5,I5,A1,I1,A1,A1,I1,A1,I5,I5,A3,F6.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table7.dat')
      write(6,*) '....Loading file: table7.dat'
      do i__=1,134457
        read(1,'(A186)')ar__1
        read(ar__1,2)
     +  lambda_1(i__),orggf(i__),newgf(i__),e_newgf(i__),r_newgf(i__),
     +  astgf(i__),r_astgf(i__),specid(i__),EP1(i__),J1(i__),
     +  EP1id(i__),EP2(i__),J2(i__),EP2id(i__),Rad(i__),Sta(i__),
     +  vdW(i__),unlte(i__),lnlte(i__),iso1(i__),hyp(i__),iso2(i__),
     +  isof(i__),hE1(i__),hE2(i__),F0(i__),F1(i__),Note1(i__),S(i__),
     +  F2(i__),note2(i__),g1(i__),g2(i__),vdWorg(i__),vdWast(i__)
        if(ar__1(11:17) .EQ. '') orggf(i__) = rNULL__
        if(ar__1(19:25) .EQ. '') newgf(i__) = rNULL__
        if(ar__1(27:30) .EQ. '') e_newgf(i__) = rNULL__
        if(ar__1(35:41) .EQ. '') astgf(i__) = rNULL__
        if(ar__1(111:116) .EQ. '') Rad(i__) = rNULL__
        if(ar__1(117:122) .EQ. '') Sta(i__) = rNULL__
        if(ar__1(123:128) .EQ. '') vdW(i__) = rNULL__
        if(ar__1(130:131) .EQ. '') unlte(i__) = iNULL__
        if(ar__1(132:133) .EQ. '') lnlte(i__) = iNULL__
        if(ar__1(134:136) .EQ. '') iso1(i__) = iNULL__
        if(ar__1(137:142) .EQ. '') hyp(i__) = rNULL__
        if(ar__1(143:145) .EQ. '') iso2(i__) = iNULL__
        if(ar__1(146:151) .EQ. '') isof(i__) = rNULL__
        if(ar__1(152:156) .EQ. '') hE1(i__) = iNULL__
        if(ar__1(157:161) .EQ. '') hE2(i__) = iNULL__
        if(ar__1(163:163) .EQ. '') F1(i__) = iNULL__
        if(ar__1(166:166) .EQ. '') F2(i__) = iNULL__
        if(ar__1(168:172) .EQ. '') g1(i__) = iNULL__
        if(ar__1(173:177) .EQ. '') g2(i__) = iNULL__
        if(ar__1(181:186) .EQ. '') vdWast(i__) = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  lambda_1(i__),orggf(i__),newgf(i__),e_newgf(i__),r_newgf(i__),
     +  astgf(i__),r_astgf(i__),specid(i__),EP1(i__),J1(i__),
     +  EP1id(i__),EP2(i__),J2(i__),EP2id(i__),Rad(i__),Sta(i__),
     +  vdW(i__),unlte(i__),lnlte(i__),iso1(i__),hyp(i__),iso2(i__),
     +  isof(i__),hE1(i__),hE2(i__),F0(i__),F1(i__),Note1(i__),S(i__),
     +  F2(i__),note2(i__),g1(i__),g2(i__),vdWorg(i__),vdWast(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table8.dat'	! *ASPCAP INT (20120216) line list

C  Format for file interpretation

    3 format(
     +  F9.4,1X,F7.3,1X,F7.3,1X,F4.2,1X,A3,F7.3,1X,A3,1X,F8.2,F12.3,
     +  F5.1,A11,F12.3,F5.1,A11,F6.2,F6.2,F6.2,1X,I2,I2,I3,F6.3,I3,
     +  F6.3,I5,I5,A1,I1,A1,A1,I1,A1,I5,I5,A3,F6.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table8.dat')
      write(6,*) '....Loading file: table8.dat'
      do i__=1,134230
        read(1,'(A186)')ar__2
        read(ar__2,3)
     +  lambda_2(i__),orggf_1(i__),newgf_1(i__),e_newgf_1(i__),
     +  r_newgf_1(i__),astgf_1(i__),r_astgf_1(i__),specid_1(i__),
     +  EP1_1(i__),J1_1(i__),EP1id_1(i__),EP2_1(i__),J2_1(i__),
     +  EP2id_1(i__),Rad_1(i__),Sta_1(i__),vdW_1(i__),unlte_1(i__),
     +  lnlte_1(i__),iso1_1(i__),hyp_1(i__),iso2_1(i__),isof_1(i__),
     +  hE1_1(i__),hE2_1(i__),F0_1(i__),F1_1(i__),Note1_1(i__),
     +  S_1(i__),F2_1(i__),note2_1(i__),g1_1(i__),g2_1(i__),
     +  vdWorg_1(i__),vdWast_1(i__)
        if(ar__2(11:17) .EQ. '') orggf_1(i__) = rNULL__
        if(ar__2(19:25) .EQ. '') newgf_1(i__) = rNULL__
        if(ar__2(27:30) .EQ. '') e_newgf_1(i__) = rNULL__
        if(ar__2(35:41) .EQ. '') astgf_1(i__) = rNULL__
        if(ar__2(111:116) .EQ. '') Rad_1(i__) = rNULL__
        if(ar__2(117:122) .EQ. '') Sta_1(i__) = rNULL__
        if(ar__2(123:128) .EQ. '') vdW_1(i__) = rNULL__
        if(ar__2(130:131) .EQ. '') unlte_1(i__) = iNULL__
        if(ar__2(132:133) .EQ. '') lnlte_1(i__) = iNULL__
        if(ar__2(134:136) .EQ. '') iso1_1(i__) = iNULL__
        if(ar__2(137:142) .EQ. '') hyp_1(i__) = rNULL__
        if(ar__2(143:145) .EQ. '') iso2_1(i__) = iNULL__
        if(ar__2(146:151) .EQ. '') isof_1(i__) = rNULL__
        if(ar__2(152:156) .EQ. '') hE1_1(i__) = iNULL__
        if(ar__2(157:161) .EQ. '') hE2_1(i__) = iNULL__
        if(ar__2(163:163) .EQ. '') F1_1(i__) = iNULL__
        if(ar__2(166:166) .EQ. '') F2_1(i__) = iNULL__
        if(ar__2(168:172) .EQ. '') g1_1(i__) = iNULL__
        if(ar__2(173:177) .EQ. '') g2_1(i__) = iNULL__
        if(ar__2(181:186) .EQ. '') vdWast_1(i__) = rNULL__
c    ..............Just test output...........
        write(6,3)
     +  lambda_2(i__),orggf_1(i__),newgf_1(i__),e_newgf_1(i__),
     +  r_newgf_1(i__),astgf_1(i__),r_astgf_1(i__),specid_1(i__),
     +  EP1_1(i__),J1_1(i__),EP1id_1(i__),EP2_1(i__),J2_1(i__),
     +  EP2id_1(i__),Rad_1(i__),Sta_1(i__),vdW_1(i__),unlte_1(i__),
     +  lnlte_1(i__),iso1_1(i__),hyp_1(i__),iso2_1(i__),isof_1(i__),
     +  hE1_1(i__),hE2_1(i__),F0_1(i__),F1_1(i__),Note1_1(i__),
     +  S_1(i__),F2_1(i__),note2_1(i__),g1_1(i__),g2_1(i__),
     +  vdWorg_1(i__),vdWast_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table9.dat'	! *ASPCAP DR10 (20110510) line list

C  Format for file interpretation

    4 format(
     +  F9.4,1X,F7.3,1X,F7.3,1X,F4.2,1X,A3,F7.3,1X,A3,1X,F8.2,F12.3,
     +  F5.1,A11,F12.3,F5.1,A11,F6.2,F6.2,F6.2,1X,I2,I2,I3,F6.3,I3,
     +  F6.3,I5,I5,A1,I1,A1,A1,I1,A1,I5,I5,A3,F6.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table9.dat')
      write(6,*) '....Loading file: table9.dat'
      do i__=1,134052
        read(1,'(A186)')ar__3
        read(ar__3,4)
     +  lambda_3(i__),orggf_2(i__),newgf_2(i__),e_newgf_2(i__),
     +  r_newgf_2(i__),astgf_2(i__),r_astgf_2(i__),specid_2(i__),
     +  EP1_2(i__),J1_2(i__),EP1id_2(i__),EP2_2(i__),J2_2(i__),
     +  EP2id_2(i__),Rad_2(i__),Sta_2(i__),vdW_2(i__),unlte_2(i__),
     +  lnlte_2(i__),iso1_2(i__),hyp_2(i__),iso2_2(i__),isof_2(i__),
     +  hE1_2(i__),hE2_2(i__),F0_2(i__),F1_2(i__),Note1_2(i__),
     +  S_2(i__),F2_2(i__),note2_2(i__),g1_2(i__),g2_2(i__),
     +  vdWorg_2(i__),vdWast_2(i__)
        if(ar__3(11:17) .EQ. '') orggf_2(i__) = rNULL__
        if(ar__3(19:25) .EQ. '') newgf_2(i__) = rNULL__
        if(ar__3(27:30) .EQ. '') e_newgf_2(i__) = rNULL__
        if(ar__3(35:41) .EQ. '') astgf_2(i__) = rNULL__
        if(ar__3(111:116) .EQ. '') Rad_2(i__) = rNULL__
        if(ar__3(117:122) .EQ. '') Sta_2(i__) = rNULL__
        if(ar__3(123:128) .EQ. '') vdW_2(i__) = rNULL__
        if(ar__3(130:131) .EQ. '') unlte_2(i__) = iNULL__
        if(ar__3(132:133) .EQ. '') lnlte_2(i__) = iNULL__
        if(ar__3(134:136) .EQ. '') iso1_2(i__) = iNULL__
        if(ar__3(137:142) .EQ. '') hyp_2(i__) = rNULL__
        if(ar__3(143:145) .EQ. '') iso2_2(i__) = iNULL__
        if(ar__3(146:151) .EQ. '') isof_2(i__) = rNULL__
        if(ar__3(152:156) .EQ. '') hE1_2(i__) = iNULL__
        if(ar__3(157:161) .EQ. '') hE2_2(i__) = iNULL__
        if(ar__3(163:163) .EQ. '') F1_2(i__) = iNULL__
        if(ar__3(166:166) .EQ. '') F2_2(i__) = iNULL__
        if(ar__3(168:172) .EQ. '') g1_2(i__) = iNULL__
        if(ar__3(173:177) .EQ. '') g2_2(i__) = iNULL__
        if(ar__3(181:186) .EQ. '') vdWast_2(i__) = rNULL__
c    ..............Just test output...........
        write(6,4)
     +  lambda_3(i__),orggf_2(i__),newgf_2(i__),e_newgf_2(i__),
     +  r_newgf_2(i__),astgf_2(i__),r_astgf_2(i__),specid_2(i__),
     +  EP1_2(i__),J1_2(i__),EP1id_2(i__),EP2_2(i__),J2_2(i__),
     +  EP2id_2(i__),Rad_2(i__),Sta_2(i__),vdW_2(i__),unlte_2(i__),
     +  lnlte_2(i__),iso1_2(i__),hyp_2(i__),iso2_2(i__),isof_2(i__),
     +  hE1_2(i__),hE2_2(i__),F0_2(i__),F1_2(i__),Note1_2(i__),
     +  S_2(i__),F2_2(i__),note2_2(i__),g1_2(i__),g2_2(i__),
     +  vdWorg_2(i__),vdWast_2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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