FORTRAN Generation
(/./ftp/cats/J/ApJ/702/1472)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/702/1472 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-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/702/1472   Column densities for HI, AlIII, SiIV, CIV, OVI  (Savage+, 2009)
*================================================================================
*The extension of the transition temperature plasma into the lower Galactic halo.
*    Savage B.D., Wakker B.P.
*   <Astrophys. J., 702, 1472-1489 (2009)>
*   =2009ApJ...702.1472S
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table2.dat'	! Column densities

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

C  Position composed of: GLON GLAT
      character*16  Name       (nr__) ! Star HD number or alternate name
      character*11  Class      (nr__) ! Object type for AGN or MK spectral type (1)
      real*4        GLON       (nr__) ! (deg) Galactic longitude
      real*4        GLAT       (nr__) ! (deg) Galactic latitude
      real*4        Dist       (nr__) ! (kpc) ? Stellar distance, d (1)
      real*4        Sig        (nr__) ! ? The estimated logarithmic error (dex) in
*                                   the distance (1)
      character*4   HII        (nr__) ! Prominent HII regions associated with the
*                                   star as identified via the footnotes
      character*1   f_HII      (nr__) ! [a-g] Note related with line of sight (2)
      integer*4     v_R_       (nr__) ! [0,3]? The 0.25-1.0keV brightness of the
*                                   X-ray sky in the direction of each object (3)
      real*4        lgN_H1     (nr__) ! ([cm-2]) ? Logarithmic HI column density (4)
      real*4        e_lgN_H1   (nr__) ! ([cm-2]) ? lgN(H1) uncertainty
      character*1   l_lgN_Al3  (nr__) ! Limit flag on lgN(Al3)
      real*4        lgN_Al3    (nr__) ! ([cm-2]) ? Logarithmic AlIII column density
      character*1   u_lgN_Al3  (nr__) ! [:] Uncertainty flag on lgN(Al3)
      real*4        e_lgN_Al3  (nr__) ! ([cm-2]) ? lgN(Al3) uncertainty
      character*1   l_lgN_Si4  (nr__) ! Limit flag on lgN(Si4)
      real*4        lgN_Si4    (nr__) ! ([cm-2]) ? Logarithmic SiIV column density
      real*4        e_lgN_Si4  (nr__) ! ([cm-2]) ? lgN(Si4) uncertainty
      character*1   l_lgN_C4   (nr__) ! Limit flag on lgN(C4)
      real*4        lgN_C4     (nr__) ! ([cm-2]) ? Logarithmic CIV column density
      character*1   u_lgN_C4   (nr__) ! [:] Uncertainty flag on lgN(C4)
      real*4        e_lgN_C4   (nr__) ! ([cm-2]) ? lgN(C4) uncertainty
      character*1   l_lgN_O6   (nr__) ! Limit flag on lgN(O6)
      real*4        lgN_O6     (nr__) ! ([cm-2]) ? Logarithmic OVI column density
      real*4        e_lgN_O6   (nr__) ! ([cm-2]) ? lgN(O6) uncertainty
      character*7   Ref        (nr__) ! AlIII, SiIV & CIV and OVI sources (5)
*Note (1): From Bowen et al. (2008, Cat. J/ApJS/176/59) and references therein.
*Note (2): Flag as follows:
*  a = Line of sight samples gas in the Scutum Supershell (see Savage et al.
*      2001ApJ...547..907S).
*  b = Line of sight samples gas in the Gum nebula.
*  c = Line of sight samples gas in the Vela SNR.
*  d = Line of sight samples gas in the Carina HII region.
*  e = Line of sight samples gas in Gum 42 which is also R62.
*  f = Line of sight samples gas in Gum 56 which is within R113.
*  g = Line of sight samples gas in Gum 57a which is within S2.
*Note (3): [R], a measure of the 0.25-1.0keV brightness of the X-ray sky in
*     the direction of each object from Bowen et al. with R=0, 1, 2, and 3
*     implying an increasing brightness of the X-ray sky in the immediate
*     vicinity of the star or AGN.
*Note (4): The HI column densities for the stars are from the Ly{alpha}
*     observations of Diplas & Savage (1994ApJ...427..274D). For the AGNs,
*     the HI column densities are from the 21cm results compiled by Wakker
*     et al. (2003, Cat. J/ApJS/146/1).
*Note (5): Up to 3 references are given, the last number referencing the
*  source of OVI column density, and PF=Profile Fit method. 
*  The sources are as follows:
*  1  = Bowen et al. 2008, Cat. J/ApJS/176/59;
*  2  = Savage et al. 2001, Cat. J/ApJS/136/631;
*  3  = Zsargo et al. 2003ApJ...586.1019Z;
*  4  = Savage & Sembach 1994ApJ...434..145S;
*  5  = Fox et al. 2003ApJ...582..793F;
*  6  = Savage et al. 2001ApJ...547..907S;
*  7  = Jenkins 1978ApJ...219..845J;
*  8  = Edgar & Savage 1992ApJ...396..124E;
*  9  = Brandt et al. 1999AJ....117..400B;
*  10 = Fitzpatrick & Spitzer 1997ApJ...475..623F;
*  11 = Huang et al. 1995ApJ...450..163H;
*  12 = Howk & Savage 1999ApJ...517..746H;
*  13 = AGN results from B. P. Wakker et al. (2010, in preparation);
*  14 = LMC results from Lehner & Howk, 2007MNRAS.377..687L;
*  15 = Sembach & Savage, 1992ApJS...83..147S;
*  16 = Hoopes et al. 2002ApJ...569..233H.

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

C  Loading file 'table2.dat'	! Column densities

C  Format for file interpretation

    1 format(
     +  A16,1X,A11,1X,F6.2,1X,F6.2,1X,F4.1,1X,F4.2,1X,A4,A1,1X,I1,1X,
     +  F5.2,1X,F4.2,1X,A1,F5.2,A1,1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,A1,
     +  F5.2,A1,1X,F4.2,1X,A1,F5.2,1X,F4.2,1X,A7)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,140
        read(1,'(A129)')ar__
        read(ar__,1)
     +  Name(i__),Class(i__),GLON(i__),GLAT(i__),Dist(i__),Sig(i__),
     +  HII(i__),f_HII(i__),v_R_(i__),lgN_H1(i__),e_lgN_H1(i__),
     +  l_lgN_Al3(i__),lgN_Al3(i__),u_lgN_Al3(i__),e_lgN_Al3(i__),
     +  l_lgN_Si4(i__),lgN_Si4(i__),e_lgN_Si4(i__),l_lgN_C4(i__),
     +  lgN_C4(i__),u_lgN_C4(i__),e_lgN_C4(i__),l_lgN_O6(i__),
     +  lgN_O6(i__),e_lgN_O6(i__),Ref(i__)
        if(ar__(44:47) .EQ. '') Dist(i__) = rNULL__
        if(ar__(49:52) .EQ. '') Sig(i__) = rNULL__
        if(ar__(60:60) .EQ. '') v_R_(i__) = iNULL__
        if(ar__(62:66) .EQ. '') lgN_H1(i__) = rNULL__
        if(ar__(68:71) .EQ. '') e_lgN_H1(i__) = rNULL__
        if(ar__(74:78) .EQ. '') lgN_Al3(i__) = rNULL__
        if(ar__(81:84) .EQ. '') e_lgN_Al3(i__) = rNULL__
        if(ar__(87:91) .EQ. '') lgN_Si4(i__) = rNULL__
        if(ar__(93:96) .EQ. '') e_lgN_Si4(i__) = rNULL__
        if(ar__(99:103) .EQ. '') lgN_C4(i__) = rNULL__
        if(ar__(106:109) .EQ. '') e_lgN_C4(i__) = rNULL__
        if(ar__(112:116) .EQ. '') lgN_O6(i__) = rNULL__
        if(ar__(118:121) .EQ. '') e_lgN_O6(i__) = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  Name(i__),Class(i__),GLON(i__),GLAT(i__),Dist(i__),Sig(i__),
     +  HII(i__),f_HII(i__),v_R_(i__),lgN_H1(i__),e_lgN_H1(i__),
     +  l_lgN_Al3(i__),lgN_Al3(i__),u_lgN_Al3(i__),e_lgN_Al3(i__),
     +  l_lgN_Si4(i__),lgN_Si4(i__),e_lgN_Si4(i__),l_lgN_C4(i__),
     +  lgN_C4(i__),u_lgN_C4(i__),e_lgN_C4(i__),l_lgN_O6(i__),
     +  lgN_O6(i__),e_lgN_O6(i__),Ref(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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