FORTRAN Generation
(/./ftp/cats/J/ApJS/230/6)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJS/230/6 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/ApJS/230/6 HST -COS & -STIS absorption-line spectroscopy. II. (Keeney+, 2017)
*================================================================================
*Characterizing the circumgalactic medium of nearby galaxies with HST/COS and
*HST/STIS absorption-line spectroscopy.
*II. Methods and models.
*    Keeney B.A., Stocke J.T., Danforth C.W., Shull J.M., Pratt C.T.,
*    Froning C.S., Green J.C., Penton S.V., Savage B.D.
*   <Astrophys. J. Suppl. Ser., 230, 6-6 (2017)>
*   =2017ApJS..230....6K    (SIMBAD/NED BibCode)
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! Targeted sample summary

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

      character*17  Line       (nr__) ! QSO sight line identifier
      character*24  Gal        (nr__) ! Associated galaxy identifier
      integer*4     cz1        (nr__) ! (km/s) [654/52914] First value of absorber velocity
*                              (cz_abs_)
      integer*4     cz2        (nr__) ! (km/s) ? Second value of cz_abs_
      integer*4     cz3        (nr__) ! (km/s) ? Third value of cz_abs_
      integer*4     cz4        (nr__) ! (km/s) ? Fourth value of cz_abs_
      integer*4     cz5        (nr__) ! (km/s) ? Fifth value of cz_abs_
      integer*4     cz6        (nr__) ! (km/s) ? Sixth value of cz_abs_
      character*10  Inst       (nr__) ! Instrument(s) used for the absorber detections
      integer*4     czgal      (nr__) ! (km/s) [649/52873] Galaxy's recession velocity (cz_gal_)
      real*4        Lum        (nr__) ! (Lsun) [0.003/3.5] Rest-frame g-band galaxy luminosity (1)
      character*1   u_Lum      (nr__) ! Uncertainty flag on Lum
      integer*4     rho        (nr__) ! (kpc) [7/505] Impact parameter ({rho})
      integer*4     Rvir       (nr__) ! (kpc) [47/275] Virial radius (R_vir_)
      character*1   u_Rvir     (nr__) ! Uncertainty flag on Rvir
      real*4        r_R        (nr__) ! [0.1/2.1] Ratio of impact parameter to R_vir_
      character*1   u_r_R      (nr__) ! Uncertainty flag on r/R
      character*1   l_eta      (nr__) ! Limit flag on eta
      real*4        eta        (nr__) ! [0.02/4.3] {eta}_{rho}_ ratio (2)
      real*4        eVir       (nr__) ! [0.02/4]? {eta}_vir_ ratio (3)
      character*1   u_eVir     (nr__) ! Uncertainty flag on eVir
      real*4        eDelv      (nr__) ! [0.01/6.8]? {eta}_{Delta}v_ ratio (4)
      character*1   u_eDelv    (nr__) ! Uncertainty flag on eDelv
*Note (1): The rest-frame g-band galaxy luminosity in L^*^ units (we adopt
*          M_g_^*^=-20.3, the value of Montero-Dorta & Prada 2009MNRAS.399.1106M
*          at h=0.7).
*Note (2): The ratio of the impact parameter (in kpc) of the tabulated galaxy to
*          that of the next-nearest galaxy:
*          {eta}_{rho}_={rho}_ng_/{rho}_nng_.
*Note (3): The ratio of the normalized impact parameter of the tabulated galaxy
*          to that of the next-nearest galaxy:
*          {eta}_vir_=({rho}/R_vir_)_ng_/({rho}/R_vir_)_nng_.
*Note (4): The ratio of the normalized absorber-galaxy velocity difference of
*          the tabulated galaxy to that of the next-nearest galaxy:
*          {eta}_{Delta}v_=(|{Delta}v|/v_esc_)_ng_/(|{Delta}v|/v_esc_)_nng_.

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

C  Declarations for 'table2.dat'	! Serendipitous sample summary

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

      character*17  Line_1     (nr__1) ! QSO sight line identifier
      character*24  Gal_1      (nr__1) ! Associated galaxy identifier
      integer*4     cz1_1      (nr__1) ! (km/s) [654/52914] First value of absorber velocity
*                              (cz_abs_)
      integer*4     cz2_1      (nr__1) ! (km/s) ? Second value of cz_abs_
      integer*4     cz3_1      (nr__1) ! (km/s) ? Third value of cz_abs_
      integer*4     cz4_1      (nr__1) ! (km/s) ? Fourth value of cz_abs_
      integer*4     cz5_1      (nr__1) ! (km/s) ? Fifth value of cz_abs_
      integer*4     cz6_1      (nr__1) ! (km/s) ? Sixth value of cz_abs_
      character*10  Inst_1     (nr__1) ! Instrument(s) used for the absorber detections
      integer*4     czgal_1    (nr__1) ! (km/s) [649/52873] Galaxy's recession velocity (cz_gal_)
      real*4        Lum_1      (nr__1) ! (Lsun) [0.003/3.5] Rest-frame g-band galaxy luminosity (1)
      character*1   u_Lum_1    (nr__1) ! Uncertainty flag on Lum
      integer*4     rho_1      (nr__1) ! (kpc) [7/505] Impact parameter ({rho})
      integer*4     Rvir_1     (nr__1) ! (kpc) [47/275] Virial radius (R_vir_)
      character*1   u_Rvir_1   (nr__1) ! Uncertainty flag on Rvir
      real*4        r_R_1      (nr__1) ! [0.1/2.1] Ratio of impact parameter to R_vir_
      character*1   u_r_R_1    (nr__1) ! Uncertainty flag on r/R
      character*1   l_eta_1    (nr__1) ! Limit flag on eta
      real*4        eta_1      (nr__1) ! [0.02/4.3] {eta}_{rho}_ ratio (2)
      real*4        eVir_1     (nr__1) ! [0.02/4]? {eta}_vir_ ratio (3)
      character*1   u_eVir_1   (nr__1) ! Uncertainty flag on eVir
      real*4        eDelv_1    (nr__1) ! [0.01/6.8]? {eta}_{Delta}v_ ratio (4)
      character*1   u_eDelv_1  (nr__1) ! Uncertainty flag on eDelv
*Note (1): The rest-frame g-band galaxy luminosity in L^*^ units (we adopt
*          M_g_^*^=-20.3, the value of Montero-Dorta & Prada 2009MNRAS.399.1106M
*          at h=0.7).
*Note (2): The ratio of the impact parameter (in kpc) of the tabulated galaxy to
*          that of the next-nearest galaxy:
*          {eta}_{rho}_={rho}_ng_/{rho}_nng_.
*Note (3): The ratio of the normalized impact parameter of the tabulated galaxy
*          to that of the next-nearest galaxy:
*          {eta}_vir_=({rho}/R_vir_)_ng_/({rho}/R_vir_)_nng_.
*Note (4): The ratio of the normalized absorber-galaxy velocity difference of
*          the tabulated galaxy to that of the next-nearest galaxy:
*          {eta}_{Delta}v_=(|{Delta}v|/v_esc_)_ng_/(|{Delta}v|/v_esc_)_nng_.

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

C  Declarations for 'table3.dat'	! Multi-component absorption line fits

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

      integer*4     ID         (nr__2) ! [1/187] System identifier
      character*17  Line_2     (nr__2) ! Sight line identifier
      character*6   Ion        (nr__2) ! Ion identifier
      integer*4     cz         (nr__2) ! (km/s) [643/52945]? Redshift from absorption line
      real*4        b          (nr__2) ! (km/s) [4/128]? Doppler b-value
      character*1   u_b        (nr__2) ! Uncertainty flag on b
      real*4        e_b        (nr__2) ! (km/s) [0/57]? Uncertainty in b
      character*1   l_logN     (nr__2) ! Limit flag on logN
      real*4        logN       (nr__2) ! ([cm-2]) [11.6/18.4]? Log number abundance
      character*1   u_logN     (nr__2) ! Uncertainty flag in logN
      real*4        e_logN     (nr__2) ! ([cm-2]) [0.01/3]? Uncertainty in logN
      character*61  Comm       (nr__2) ! Additional comments

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

C  Loading file 'table1.dat'	! Targeted sample summary

C  Format for file interpretation

    1 format(
     +  A17,1X,A24,1X,I5,1X,I5,1X,I5,1X,I5,1X,I5,1X,I5,1X,A10,1X,I5,
     +  1X,F5.3,A1,1X,I3,1X,I3,A1,1X,F4.2,A1,1X,A1,F4.2,1X,F4.2,A1,1X,
     +  F4.2,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,13
        read(1,'(A135)')ar__
        read(ar__,1)
     +  Line(i__),Gal(i__),cz1(i__),cz2(i__),cz3(i__),cz4(i__),
     +  cz5(i__),cz6(i__),Inst(i__),czgal(i__),Lum(i__),u_Lum(i__),
     +  rho(i__),Rvir(i__),u_Rvir(i__),r_R(i__),u_r_R(i__),l_eta(i__),
     +  eta(i__),eVir(i__),u_eVir(i__),eDelv(i__),u_eDelv(i__)
        if(ar__(50:54) .EQ. '') cz2(i__) = iNULL__
        if(ar__(56:60) .EQ. '') cz3(i__) = iNULL__
        if(ar__(62:66) .EQ. '') cz4(i__) = iNULL__
        if(ar__(68:72) .EQ. '') cz5(i__) = iNULL__
        if(ar__(74:78) .EQ. '') cz6(i__) = iNULL__
        if(ar__(125:128) .EQ. '') eVir(i__) = rNULL__
        if(ar__(131:134) .EQ. '') eDelv(i__) = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  Line(i__),Gal(i__),cz1(i__),cz2(i__),cz3(i__),cz4(i__),
     +  cz5(i__),cz6(i__),Inst(i__),czgal(i__),Lum(i__),u_Lum(i__),
     +  rho(i__),Rvir(i__),u_Rvir(i__),r_R(i__),u_r_R(i__),l_eta(i__),
     +  eta(i__),eVir(i__),u_eVir(i__),eDelv(i__),u_eDelv(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table2.dat'	! Serendipitous sample summary

C  Format for file interpretation

    2 format(
     +  A17,1X,A24,1X,I5,1X,I5,1X,I5,1X,I5,1X,I5,1X,I5,1X,A10,1X,I5,
     +  1X,F5.3,A1,1X,I3,1X,I3,A1,1X,F4.2,A1,1X,A1,F4.2,1X,F4.2,A1,1X,
     +  F4.2,A1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,35
        read(1,'(A135)')ar__1
        read(ar__1,2)
     +  Line_1(i__),Gal_1(i__),cz1_1(i__),cz2_1(i__),cz3_1(i__),
     +  cz4_1(i__),cz5_1(i__),cz6_1(i__),Inst_1(i__),czgal_1(i__),
     +  Lum_1(i__),u_Lum_1(i__),rho_1(i__),Rvir_1(i__),u_Rvir_1(i__),
     +  r_R_1(i__),u_r_R_1(i__),l_eta_1(i__),eta_1(i__),eVir_1(i__),
     +  u_eVir_1(i__),eDelv_1(i__),u_eDelv_1(i__)
        if(ar__1(50:54) .EQ. '') cz2_1(i__) = iNULL__
        if(ar__1(56:60) .EQ. '') cz3_1(i__) = iNULL__
        if(ar__1(62:66) .EQ. '') cz4_1(i__) = iNULL__
        if(ar__1(68:72) .EQ. '') cz5_1(i__) = iNULL__
        if(ar__1(74:78) .EQ. '') cz6_1(i__) = iNULL__
        if(ar__1(125:128) .EQ. '') eVir_1(i__) = rNULL__
        if(ar__1(131:134) .EQ. '') eDelv_1(i__) = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  Line_1(i__),Gal_1(i__),cz1_1(i__),cz2_1(i__),cz3_1(i__),
     +  cz4_1(i__),cz5_1(i__),cz6_1(i__),Inst_1(i__),czgal_1(i__),
     +  Lum_1(i__),u_Lum_1(i__),rho_1(i__),Rvir_1(i__),u_Rvir_1(i__),
     +  r_R_1(i__),u_r_R_1(i__),l_eta_1(i__),eta_1(i__),eVir_1(i__),
     +  u_eVir_1(i__),eDelv_1(i__),u_eDelv_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table3.dat'	! Multi-component absorption line fits

C  Format for file interpretation

    3 format(
     +  I3,1X,A17,1X,A6,1X,I5,1X,F5.1,A1,1X,F4.1,1X,A1,1X,F5.2,A1,1X,
     +  F4.2,1X,A61)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,467
        read(1,'(A122)')ar__2
        read(ar__2,3)
     +  ID(i__),Line_2(i__),Ion(i__),cz(i__),b(i__),u_b(i__),e_b(i__),
     +  l_logN(i__),logN(i__),u_logN(i__),e_logN(i__),Comm(i__)
        if(ar__2(30:34) .EQ. '') cz(i__) = iNULL__
        if(ar__2(36:40) .EQ. '') b(i__) = rNULL__
        if(ar__2(43:46) .EQ. '') e_b(i__) = rNULL__
        if(ar__2(50:54) .EQ. '') logN(i__) = rNULL__
        if(ar__2(57:60) .EQ. '') e_logN(i__) = rNULL__
c    ..............Just test output...........
        write(6,3)
     +  ID(i__),Line_2(i__),Ion(i__),cz(i__),b(i__),u_b(i__),e_b(i__),
     +  l_logN(i__),logN(i__),u_logN(i__),e_logN(i__),Comm(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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