FORTRAN Generation
(/./ftp/cats/J/ApJ/700/1299)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/700/1299 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-28
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/700/1299       Gas-phase element depletions in the ISM     (Jenkins, 2009)
*================================================================================
*A unified representation of gas-phase element depletions in the interstellar
*medium.
*    Jenkins E.B.
*   <Astrophys. J., 700, 1299-1348 (2009)>
*   =2009ApJ...700.1299J
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table2.dat'	! Stellar data and line of sight information

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

C  Position composed of: GLON GLAT
      integer*4     HD         (nr__) ! ? HD number
      integer*4     Vcomp      (nr__) ! (km/s) ? Specific Heliocentric velocity component
      character*16  Name       (nr__) ! Name (1)
      real*4        GLON       (nr__) ! (deg) Galactic longitude
      real*4        GLAT       (nr__) ! (deg) Galactic latitude
      real*4        Vmag       (nr__) ! (mag) The Simbad V band magnitude (not used for
*                                  calculating the reddenings or Dist)
      character*16  SpType     (nr__) ! MK spectral type
      character*6   r_SpType   (nr__) ! Reference for SpType; see refs.dat file
      real*4        E_B_V      (nr__) ! (mag) The (B-V) extinction
      real*4        Dist       (nr__) ! (kpc) Heliocentric distance (1)
      real*4        Z          (nr__) ! (kpc) Distance from Galactic plane
      real*4        b_logNHI   (nr__) ! ([cm-2]) ? Lower limit on logNHI
      real*4        logNHI     (nr__) ! ([cm-2]) ? Log best observed HI column density
      character*1   n_logNHI   (nr__) ! Flag on logNHI (see notes.dat file)
      real*4        B_logNHI_1 (nr__) ! ([cm-2]) ? Upper limit on logNHI
      character*6   r_logNHI   (nr__) ! Reference for logNHI; see refs.dat file
      real*4        b_logNH2   (nr__) ! ([cm-2]) ? Lower limit on logNH2
      real*4        logNH2     (nr__) ! ([cm-2]) ? Log best observed H_2_ column density
      real*4        B_logNH2_1 (nr__) ! ([cm-2]) ? Upper limit on logNH2
      character*6   r_logNH2   (nr__) ! Reference for logNH2; see refs.dat file
*Note (1): Stars for which neither the Hipparcos parallaxes (Cat. I/311)
*     nor the determinations by Bowen et al. 2008, Cat. J/ApJS/176/59 were
*     applicable had their spectroscopic parallaxes determined with the
*     standard formula for the distance d (in kpc):
*     logd=(m_V_-A_V_-M_V_-10)/5,  where m_V_ is the apparent V magnitude of
*     the star, A_V_ is the extinction by dust in the visible (assumed to be
*     3.1E(B-V)), and M_V_ is the absolute V magnitude of the star. 
*     See section 4.5.1 for further details.

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

C  Declarations for 'table4.dat'	! Element depletion parameters

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

      character*2   El         (nr__1) ! Element identification "X" (1)
      real*4        v_X_H_     (nr__1) ! ([-]) Adopted reference abundance relative to H (2)
      real*4        e__X_H_    (nr__1) ! ([-]) Uncertainty in X/H (2)
      real*4        Ax         (nr__1) ! The A_X_ parameter; defined in Eq. 10 (3)
      real*4        e_Ax       (nr__1) ! Uncertainty in Ax
      real*4        Bx         (nr__1) ! The B_X_ parameter; defined in Eq. 10 (3)
      real*4        e_Bx       (nr__1) ! Uncertainty in Bx (4)
      real*4        zx         (nr__1) ! The z_X_ parameter; defined in Eq. 11 (3)
      real*4        Xgas_H0    (nr__1) ! ([-]) The [X_gas_/H]_0_ parameter; defined in Eq. 13
      real*4        e_Xgas_H0  (nr__1) ! ([-]) Uncertainty in Xgas/H0; defined in Eq. 14 (3)(4)
      real*4        Xgas_H1    (nr__1) ! ([-]) The [X_gas_/H]_1_ parameter; defined in Eq. 15
      real*4        e_Xgas_H1  (nr__1) ! ([-]) Uncertainty in Xgas/H1; defined in Eq. 16 (4)
      real*4        chi2       (nr__1) ! The {chi}^2^ value
      integer*4     dof        (nr__1) ! Degrees of freedom
      real*4        Prob       (nr__1) ! Probability of worse fit
*Note (1): Coefficients for S do not appear in this table because a
*     nonstandard approach was required. The coefficients are given 
*     in the text of Section 9.
*Note (2): On a logarithmic scale with H=12. Values and their errors
*     taken from the recommended solar abundances of Lodders K.
*     (2003ApJ...591.1220L).
*Note (3): The slope of a best-fit line for the observed depletions
*     [X_gas_/H]_obs_=A_X_(F*-z_X_)+B_X_ for any given element X against
*     F_*_ (the line-of-sight depletion strength factor). The parameter A_X_
*     represents the propensity of that element to increase (the absolute
*     value of) its particular depletion level as F_*_ becomes larger. For
*     even the smallest observed values of F_*_, most elements still show
*     some depletion (see Figure 1). We denote this level of depletion as
*     [X_gas_/H]_0_. See section 3 for further details.
*Note (4): Unlike the convention for listing errors in the fit outcomes in
*     Table 7, the uncertainties with the terms listed here include both the
*     formal errors of the fit coefficients and the error in the adopted
*     value of (X/H)_{sun}_, added together in quadrature.

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

C  Declarations for 'table5.dat'	! Observed and synthetic log N(H) and F_*_

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

      integer*4     HD_1       (nr__2) ! ? HD number
      integer*4     Vcomp_1    (nr__2) ! (km/s) ? Specific Heliocentric velocity component (1)
      character*16  Name_1     (nr__2) ! Name
      real*4        b_logNHo   (nr__2) ! ([cm-2]) ? Lower limit on logNHOb
      real*4        logNHo     (nr__2) ! ([cm-2]) ? Log observed H column density
      real*4        B_logNHo_1 (nr__2) ! ([cm-2]) ? Upper limit on logNHOb
      real*4        F_o        (nr__2) ! ? Observed line-of-sight depletion strength
*                                  factor
      real*4        e_F_o      (nr__2) ! ? Uncertainty in F*o
      real*4        logNHs     (nr__2) ! ([cm-2]) ? Log synthetic H column density
      real*4        e_logNHs   (nr__2) ! ([cm-2]) ? Uncertainty in logNHs
      real*4        F_s        (nr__2) ! ? Synthetic line-of-sight depletion strength
*                                  factor
      real*4        e_F_s      (nr__2) ! ? Uncertainty in F*s
      real*4        Prob_1     (nr__2) ! ? Probability of a worse fit
      character*41  El_1       (nr__2) ! Elements considered in fit
*Note (1): Separate velocity components, explicitly identified and not
*     grouped together; see section 4.3.

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

C  Declarations for 'table6.dat'	! Synthetic log N(H) and F_*_ for stars in the
                             Local Bubble

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

C  Position composed of: GLON GLAT
      character*8   WD         (nr__3) ! White Dwarf identification (HHMM+DDd; B1950)
      character*15  AName      (nr__3) ! Alternate name
      real*4        GLON_1     (nr__3) ! (deg) Galactic longitude
      real*4        GLAT_1     (nr__3) ! (deg) Galactic latitude
      real*4        Vmag_1     (nr__3) ! (mag) The Simbad V band magnitude
      integer*4     Dist_1     (nr__3) ! (pc) Heliocentric distance (1)
      real*4        logNHs_1   (nr__3) ! ([-]) Log synthetic H column density
      real*4        e_logNHs_1 (nr__3) ! ([-]) Uncertainty in logNHs
      real*4        F_s_1      (nr__3) ! Synthetic line-of-sight depletion strength factor
      real*4        e_F_s_1    (nr__3) ! Uncertainty in F*s
      real*4        chi2_1     (nr__3) ! ? The {chi}^2^ value
      character*9   El_2       (nr__3) ! Element(s) considered (2)
*Note (1): Taken from the papers that described the FUSE observations.
*Note (2): All column density data from Lehner et al. (2003ApJ...595..858L),
*          except for HZ 21 (Oliveira et al. 2003ApJ...587..235O).

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

C  Declarations for 'table7.dat'	! Observations and fits for Carbon - Krypton

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

      character*2   El_3       (nr__4) ! Element "X"
      integer*4     HD_2       (nr__4) ! ? HD number
      integer*4     Vcomp_2    (nr__4) ! (km/s) ? Specific Heliocentric velocity component
      character*16  Name_2     (nr__4) ! Name (1)
      real*4        b_logNx    (nr__4) ! ([cm-2]) ? Lower limit on logNx
      real*4        logNx      (nr__4) ! ([cm-2]) ? Log best observed column density (2)
      character*1   f_logNx    (nr__4) ! [if] infinite or fixed value of logNx  (3)
      real*4        B_logNx_1  (nr__4) ! ([cm-2]) ? Upper limit on logNx
      character*1   n_logNx    (nr__4) ! [if] infinite or fixed value of B_logNx  (3)
      character*6   r_logNx    (nr__4) ! Reference codes for logNx  (see refs.dat file)
      real*4        F_         (nr__4) ! ? Line-of-sight depletion strength factor
      real*4        e_F_       (nr__4) ! ? Uncertainty in F*
      real*4        v_X_H__1   (nr__4) ! ([-]) ? Log expected depletion [X_gas_/H] (4)
      real*4        e__X_H__1  (nr__4) ! ([-]) ? Uncertainty in Elgas/H
      real*4        Res        (nr__4) ! ? Observed [X/H] minus computed from Eq. 10
      real*4        sigma      (nr__4) ! ? Standard deviation (5)
*Note (1): In table 7, {omicron}^2^ Sco is a misprint for {omicron}^1^ Sco;
*          corrected at CDS.
*Note (2): Column densities for CII, NI, OI, MgII, SiII, PII, SII, ClII, TiII,
*          CrII, MnII, FeII, NiII, CuII, ZnII, GeII, and KrI.
*Note (3): Flag as follows:
*    i = indicates value is infinite;
*    f = The entry 13.44 for the adopted value of N (NiII) stated by CLMS06
*        (Cartledge et al. 2006ApJ...641..327C) is incorrect; this number should
*        be 13.94, as confirmed by Cartledge (private communication). The value
*        listed here includes a correction for a change in the transition
*        f -value by -0.106dex - see Table 1.
*Note (4): Computed using Eq. 10. The listed errors do not include an overall
*          systematic uncertainty of
*   C/H = 0.04,
*   N/H = 0.11,
*   O/H = 0.05,
*  Mg/H = 0.02,
*  Si/H = 0.02,
*   P/H = 0.04,
*   S/H = 0.04,
*  Cl/H = 0.06,
*  Ti/H = 0.03,
*  Cr/H = 0.05,
*  Mn/H = 0.03,
*  Fe/H = 0.03,
*  Ni/H = 0.03,
*  Cu/H = 0.06,
*  Zn/H = 0.04,
*  Ge/H = 0.05, and
*  Kr/H = 0.08
*     in the solar abundance in order to show just the formal error that
*     arises from the uncertainties in location of line of best fit and 
*     the value of F*.
*Note (5): The difference shown in the previous column divided by the
*     expected difference due to the uncertainties in both the measured
*     column density and the coefficients that appear in Eq. 10.

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

C  Declarations for 'refs.dat'	! References

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

      character*6   Ref        (nr__5) ! Code of reference
      character*19  BibCode    (nr__5) ! Bibcode
      character*21  Auth       (nr__5) ! Author's name(s)
      character*19  Comm       (nr__5) ! Comment

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

C  Declarations for 'notes.dat'	! Notes of table 2

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

      character*1   F          (nr__6) ! [d-r] Flag code
      character*74  Note       (nr__6) ! Text of note

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

C  Loading file 'table2.dat'	! Stellar data and line of sight information

C  Format for file interpretation

    1 format(
     +  I6,4X,I3,2X,A16,1X,F6.2,1X,F6.2,1X,F5.2,1X,A16,1X,A6,1X,F5.2,
     +  1X,F5.3,1X,F5.2,1X,F5.2,3X,F5.2,1X,A1,1X,F5.2,3X,A6,2X,F5.2,
     +  1X,F5.2,1X,F5.2,1X,A6)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,276
        read(1,'(A150)')ar__
        read(ar__,1)
     +  HD(i__),Vcomp(i__),Name(i__),GLON(i__),GLAT(i__),Vmag(i__),
     +  SpType(i__),r_SpType(i__),E_B_V(i__),Dist(i__),Z(i__),
     +  b_logNHI(i__),logNHI(i__),n_logNHI(i__),B_logNHI_1(i__),
     +  r_logNHI(i__),b_logNH2(i__),logNH2(i__),B_logNH2_1(i__),
     +  r_logNH2(i__)
        if(ar__(1:6) .EQ. '') HD(i__) = iNULL__
        if(ar__(11:13) .EQ. '') Vcomp(i__) = iNULL__
        if(ar__(95:99) .EQ. '') b_logNHI(i__) = rNULL__
        if(ar__(103:107) .EQ. '') logNHI(i__) = rNULL__
        if(ar__(111:115) .EQ. '') B_logNHI_1(i__) = rNULL__
        if(ar__(127:131) .EQ. '') b_logNH2(i__) = rNULL__
        if(ar__(133:137) .EQ. '') logNH2(i__) = rNULL__
        if(ar__(139:143) .EQ. '') B_logNH2_1(i__) = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  HD(i__),Vcomp(i__),Name(i__),GLON(i__),GLAT(i__),Vmag(i__),
     +  SpType(i__),r_SpType(i__),E_B_V(i__),Dist(i__),Z(i__),
     +  b_logNHI(i__),logNHI(i__),n_logNHI(i__),B_logNHI_1(i__),
     +  r_logNHI(i__),b_logNH2(i__),logNH2(i__),B_logNH2_1(i__),
     +  r_logNH2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table4.dat'	! Element depletion parameters

C  Format for file interpretation

    2 format(
     +  A2,1X,F4.2,1X,F4.2,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F5.3,1X,
     +  F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F5.1,1X,I3,1X,F5.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table4.dat')
      write(6,*) '....Loading file: table4.dat'
      do i__=1,16
        read(1,'(A86)')ar__1
        read(ar__1,2)
     +  El(i__),v_X_H_(i__),e__X_H_(i__),Ax(i__),e_Ax(i__),Bx(i__),
     +  e_Bx(i__),zx(i__),Xgas_H0(i__),e_Xgas_H0(i__),Xgas_H1(i__),
     +  e_Xgas_H1(i__),chi2(i__),dof(i__),Prob(i__)
c    ..............Just test output...........
        write(6,2)
     +  El(i__),v_X_H_(i__),e__X_H_(i__),Ax(i__),e_Ax(i__),Bx(i__),
     +  e_Bx(i__),zx(i__),Xgas_H0(i__),e_Xgas_H0(i__),Xgas_H1(i__),
     +  e_Xgas_H1(i__),chi2(i__),dof(i__),Prob(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table5.dat'	! Observed and synthetic log N(H) and F_*_

C  Format for file interpretation

    3 format(
     +  I6,4X,I3,2X,A16,1X,F5.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F4.2,1X,
     +  F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.3,1X,A41)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table5.dat')
      write(6,*) '....Loading file: table5.dat'
      do i__=1,239
        read(1,'(A130)')ar__2
        read(ar__2,3)
     +  HD_1(i__),Vcomp_1(i__),Name_1(i__),b_logNHo(i__),logNHo(i__),
     +  B_logNHo_1(i__),F_o(i__),e_F_o(i__),logNHs(i__),e_logNHs(i__),
     +  F_s(i__),e_F_s(i__),Prob_1(i__),El_1(i__)
        if(ar__2(1:6) .EQ. '') HD_1(i__) = iNULL__
        if(ar__2(11:13) .EQ. '') Vcomp_1(i__) = iNULL__
        if(ar__2(33:37) .EQ. '') b_logNHo(i__) = rNULL__
        if(ar__2(39:43) .EQ. '') logNHo(i__) = rNULL__
        if(ar__2(45:49) .EQ. '') B_logNHo_1(i__) = rNULL__
        if(ar__2(51:55) .EQ. '') F_o(i__) = rNULL__
        if(ar__2(57:60) .EQ. '') e_F_o(i__) = rNULL__
        if(ar__2(62:66) .EQ. '') logNHs(i__) = rNULL__
        if(ar__2(68:71) .EQ. '') e_logNHs(i__) = rNULL__
        if(ar__2(73:77) .EQ. '') F_s(i__) = rNULL__
        if(ar__2(79:82) .EQ. '') e_F_s(i__) = rNULL__
        if(ar__2(84:88) .EQ. '') Prob_1(i__) = rNULL__
c    ..............Just test output...........
        write(6,3)
     +  HD_1(i__),Vcomp_1(i__),Name_1(i__),b_logNHo(i__),logNHo(i__),
     +  B_logNHo_1(i__),F_o(i__),e_F_o(i__),logNHs(i__),e_logNHs(i__),
     +  F_s(i__),e_F_s(i__),Prob_1(i__),El_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table6.dat'	! Synthetic log N(H) and F_*_ for stars in the
*                             Local Bubble

C  Format for file interpretation

    4 format(
     +  A8,1X,A15,1X,F6.2,1X,F6.2,1X,F4.1,1X,I3,1X,F5.2,1X,F4.2,1X,
     +  F5.2,1X,F4.2,1X,F4.2,1X,A9)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table6.dat')
      write(6,*) '....Loading file: table6.dat'
      do i__=1,29
        read(1,'(A84)')ar__3
        read(ar__3,4)
     +  WD(i__),AName(i__),GLON_1(i__),GLAT_1(i__),Vmag_1(i__),
     +  Dist_1(i__),logNHs_1(i__),e_logNHs_1(i__),F_s_1(i__),
     +  e_F_s_1(i__),chi2_1(i__),El_2(i__)
        if(ar__3(71:74) .EQ. '') chi2_1(i__) = rNULL__
c    ..............Just test output...........
        write(6,4)
     +  WD(i__),AName(i__),GLON_1(i__),GLAT_1(i__),Vmag_1(i__),
     +  Dist_1(i__),logNHs_1(i__),e_logNHs_1(i__),F_s_1(i__),
     +  e_F_s_1(i__),chi2_1(i__),El_2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table7.dat'	! Observations and fits for Carbon - Krypton

C  Format for file interpretation

    5 format(
     +  A2,1X,I6,2X,I3,2X,A16,1X,F5.2,3X,F5.2,1X,A1,1X,F5.2,1X,A1,1X,
     +  A6,1X,F5.2,1X,F4.2,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table7.dat')
      write(6,*) '....Loading file: table7.dat'
      do i__=1,1475
        read(1,'(A108)')ar__4
        read(ar__4,5)
     +  El_3(i__),HD_2(i__),Vcomp_2(i__),Name_2(i__),b_logNx(i__),
     +  logNx(i__),f_logNx(i__),B_logNx_1(i__),n_logNx(i__),
     +  r_logNx(i__),F_(i__),e_F_(i__),v_X_H__1(i__),e__X_H__1(i__),
     +  Res(i__),sigma(i__)
        if(ar__4(4:9) .EQ. '') HD_2(i__) = iNULL__
        if(ar__4(12:14) .EQ. '') Vcomp_2(i__) = iNULL__
        if(ar__4(34:38) .EQ. '') b_logNx(i__) = rNULL__
        if(ar__4(42:46) .EQ. '') logNx(i__) = rNULL__
        if(ar__4(50:54) .EQ. '') B_logNx_1(i__) = rNULL__
        if(ar__4(65:69) .EQ. '') F_(i__) = rNULL__
        if(ar__4(71:74) .EQ. '') e_F_(i__) = rNULL__
        if(ar__4(76:81) .EQ. '') v_X_H__1(i__) = rNULL__
        if(ar__4(83:87) .EQ. '') e__X_H__1(i__) = rNULL__
        if(ar__4(89:94) .EQ. '') Res(i__) = rNULL__
        if(ar__4(96:100) .EQ. '') sigma(i__) = rNULL__
c    ..............Just test output...........
        write(6,5)
     +  El_3(i__),HD_2(i__),Vcomp_2(i__),Name_2(i__),b_logNx(i__),
     +  logNx(i__),f_logNx(i__),B_logNx_1(i__),n_logNx(i__),
     +  r_logNx(i__),F_(i__),e_F_(i__),v_X_H__1(i__),e__X_H__1(i__),
     +  Res(i__),sigma(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'refs.dat'	! References

C  Format for file interpretation

    6 format(A6,1X,A19,1X,A21,2X,A19)

C  Effective file loading

      open(unit=1,status='old',file=
     +'refs.dat')
      write(6,*) '....Loading file: refs.dat'
      do i__=1,138
        read(1,'(A69)')ar__5
        read(ar__5,6)Ref(i__),BibCode(i__),Auth(i__),Comm(i__)
c    ..............Just test output...........
        write(6,6)Ref(i__),BibCode(i__),Auth(i__),Comm(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'notes.dat'	! Notes of table 2

C  Format for file interpretation

    7 format(A1,1X,A74)

C  Effective file loading

      open(unit=1,status='old',file=
     +'notes.dat')
      write(6,*) '....Loading file: notes.dat'
      do i__=1,44
        read(1,'(A76)')ar__6
        read(ar__6,7)F(i__),Note(i__)
c    ..............Just test output...........
        write(6,7)F(i__),Note(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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