FORTRAN Generation
(/./ftp/cats/J/ApJ/779/87)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/779/87 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-Aug-15
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/779/87           MAGIICAT. III. Virial masses           (Churchill+, 2013)
*================================================================================
*MAGIICAT.
*III. Interpreting self-similarity of the circumgalactic medium with virial mass
*using MgII absorption.
*    Churchill C.W., Trujillo-Gomez S., Nielsen N.M., Kacprzak G.G.
*   <Astrophys. J., 779, 87 (2013)>
*   =2013ApJ...779...87C
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! Properties of the 182 "isolated" galaxies in
                             the "MgII Absorber-Galaxy Catalog"

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

      character*8   Field      (nr__) ! Quasar Field identifier (HHMM+DDd; B1950 or
*                                   "SDSS")
      character*20  JName      (nr__) ! Quasar J2000 identifier (JHHMMSS.ss+DDMMSS.s)
      real*4        z          (nr__) ! [0.07/1.2] Galaxy spectroscopic redshift from
*                                    Paper I (Nielsen+, 2013, J/ApJ/776/114)
      real*4        D          (nr__) ! (kpc) [5/194] Galaxy-quasar impact parameter
      real*4        rMag       (nr__) ! (mag) [-24.2/-16.6] Absolute r band AB magnitude
      real*4        logM       (nr__) ! ([Msun]) [10.7/13.9] Log of the virial mass
      real*4        e_logM     (nr__) ! ([Msun]) [0.1/0.2] Lower uncertainty in logM
      real*4        E_logM_1   (nr__) ! ([Msun]) [0.1/0.7] Upper uncertainty in logM
      integer*4     Vmax       (nr__) ! (km/s) [81/855] Maximum circular velocity (1)
      integer*4     e_Vmax     (nr__) ! (km/s) [11/59] Lower uncertainty in Vmax
      integer*4     E_Vmax_1   (nr__) ! (km/s) [30/64] Upper uncertainty in Vmax
      integer*4     Rvir       (nr__) ! (kpc) [69/845] Virial radius (1)
      integer*4     e_Rvir     (nr__) ! (kpc) [10/64] Lower uncertainty in Rvir
      integer*4     E_Rvir_1   (nr__) ! (kpc) [27/65] Upper uncertainty in Rvir
      real*4        etav       (nr__) ! [0.02/1.6] {eta}_V_ = D/Rvir (1)
      real*4        e_etav     (nr__) ! [0/0.6] Lower uncertainty in etav (1)
      real*4        E_etav_1   (nr__) ! [0/0.3] Upper uncertainty in etav (1)
      integer*4     Rc         (nr__) ! (kpc) [3/175] Theoretical cooling radius (2)
      integer*4     e_Rc       (nr__) ! (kpc) [0/18] Lower uncertainty in Rc (2)
      integer*4     E_Rc_1     (nr__) ! (kpc) [2/7] Upper uncertainty in Rc (2)
      real*4        etac       (nr__) ! [0.03/19] {eta}_C_ = D/Rc (1)
      real*4        e_etac     (nr__) ! [0/6.8] Lower uncertainty in etac (1)
      real*4        E_etac_1   (nr__) ! [0/10] Upper uncertainty in etac (1)
      real*4        Rc_Rv      (nr__) ! [0/2.6] Ratio Rc/Rvir (1)
      real*4        e_Rc_Rv    (nr__) ! [0/1.7] Lower uncertainty in Rc/Rv (1)
      real*4        E_Rc_Rv_1  (nr__) ! [0/0.5] Upper uncertainty in Rc/Rv (1)
      real*4        Wr         (nr__) ! (0.1nm) [0.003/4.5] MgII (2796{AA}) rest-frame
*                                    equivalent width in {AA}
      real*4        e_Wr       (nr__) ! (0.1nm) [0.001/0.3]?=-1 Uncertainty in Wr in {AA}
*Note (1): Uncertainties are based upon uncertainties in the virial masses. For
*          some quantities a larger (smaller) virial mass results in smaller
*          (larger) value such that the uncertainties anti-correlate.
*Note (2): Because the slope of Rc changes sign as function of virial mass,
*          where the slope is positive the uncertainties correlate and where the
*          slope is negative they anti-correlate (see Figure B1). In the narrow
*          virial mass ranges where the slope of Rc changes sign, it is possible
*          that both the upward and downward uncertainties in virial mass can
*          result in an upward (or downward) uncertainty in Rc.

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

C  Loading file 'table1.dat'	! Properties of the 182 "isolated" galaxies in
*                             the "MgII Absorber-Galaxy Catalog"

C  Format for file interpretation

    1 format(
     +  A8,1X,A20,1X,F6.4,1X,F5.1,1X,F5.1,1X,F4.1,1X,F3.1,1X,F3.1,1X,
     +  I3,1X,I2,1X,I2,1X,I3,1X,I2,1X,I2,1X,F4.2,1X,F4.2,1X,F4.2,1X,
     +  I3,1X,I2,1X,I1,1X,F5.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,F4.2,1X,
     +  F4.2,1X,F5.3,1X,F6.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,182
        read(1,'(A149)')ar__
        read(ar__,1)
     +  Field(i__),JName(i__),z(i__),D(i__),rMag(i__),logM(i__),
     +  e_logM(i__),E_logM_1(i__),Vmax(i__),e_Vmax(i__),E_Vmax_1(i__),
     +  Rvir(i__),e_Rvir(i__),E_Rvir_1(i__),etav(i__),e_etav(i__),
     +  E_etav_1(i__),Rc(i__),e_Rc(i__),E_Rc_1(i__),etac(i__),
     +  e_etac(i__),E_etac_1(i__),Rc_Rv(i__),e_Rc_Rv(i__),
     +  E_Rc_Rv_1(i__),Wr(i__),e_Wr(i__)
c    ..............Just test output...........
        write(6,1)
     +  Field(i__),JName(i__),z(i__),D(i__),rMag(i__),logM(i__),
     +  e_logM(i__),E_logM_1(i__),Vmax(i__),e_Vmax(i__),E_Vmax_1(i__),
     +  Rvir(i__),e_Rvir(i__),E_Rvir_1(i__),etav(i__),e_etav(i__),
     +  E_etav_1(i__),Rc(i__),e_Rc(i__),E_Rc_1(i__),etac(i__),
     +  e_etac(i__),E_etac_1(i__),Rc_Rv(i__),e_Rc_Rv(i__),
     +  E_Rc_Rv_1(i__),Wr(i__),e_Wr(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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