FORTRAN Generation
(/./ftp/cats/J/ApJ/701/587)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/701/587 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/701/587 A deep HST H-band imaging survey. II. QUEST QSOs (Veilleux+, 2009)
*================================================================================
*A deep Hubble space telescope H-band imaging survey of massive gas-rich mergers.
*II. The QUEST QSOs.
*    Veilleux S., Kim D.-C., Rupke D.S.N., Peng C.Y., Tacconi L.J., Genzel R.,
*    Lutz D., Sturm E., Contursi A., Schweitzer M., Dasyra K.M., Ho L.C.,
*    Sanders D.B., Burkert A.
*   <Astrophys. J., 701, 587-606 (2009)>
*   =2009ApJ...701..587V
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! PG QSOs sample

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

      character*11  Name       (nr__) ! Object name
      character*10  OName      (nr__) ! Other name
      real*4        z          (nr__) ! Redshift
      real*4        logLB      (nr__) ! ([Lsun]) Blue luminosity
      character*1   l_logLIR   (nr__) ! Limit flag on logLIR
      real*4        logLIR     (nr__) ! ([Lsun]) ? Infrared luminosity (1)
      real*4        logLbol    (nr__) ! ([Lsun]) Bolometric luminosity (2)
      character*1   l_LIR_LB   (nr__) ! Limit flag on LIR/LB
      real*4        LIR_LB     (nr__) ! ? Ratio of infrared to blue luminosities
      character*1   l_LIR_Lbol (nr__) ! Limit flag on LIR/Lbol
      real*4        LIR_Lbol   (nr__) ! ? Ratio of infrared to bolometric luminosities
      character*1   l_F25_F60  (nr__) ! Limit flag on F25/F60
      real*4        F25_F60    (nr__) ! 25/60{mu}m flux ratio
      character*1   u_F25_F60  (nr__) ! [:] Uncertainty flag on F25/F60
      character*1   RL         (nr__) ! [LQ] Radio loudness (Q=quiet vs L=loud) (3)
      character*10  Class      (nr__) ! FIR strength: Strong, Weak or Undetected (4)
*Note (1): Infrared luminosity calculated from prescription of
*     Sanders & Mirabel (1996ARA&A..34..749S).
*Note (2): Bolometric luminosity calculated from 7xL(5100{AA})+L(IR)
*     (Netzer et al. 2007, Cat. J/ApJ/666/806)
*Note (3): From Kellermann et al. (1994, Cat. J/AJ/108/1163) except for two
*     objects (B2 2201+31A and PG2349-014 from Becker et al. 1991, Cat. VIII/13).
*Note (4): far-IR strength according to Netzer et al. (2007, Cat. J/ApJ/666/806)

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

C  Declarations for 'table5.dat'	! *Summary results from the GALFIT analysis

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

      character*11  Name_1     (nr__1) ! Object name
      character*1   u_Name     (nr__1) ! [b] uncertain values (1)
      real*4        MAGt       (nr__1) ! (mag) Total absolute magnitude (host+PSF-comp.)
      real*4        MAGp       (nr__1) ! (mag) Absolute magnitude of PSF component
      real*4        MAGh       (nr__1) ! (mag) Absolute magnitude of host galaxy (including
*                                    tidal features)
      real*4        MAGm       (nr__1) ! (mag) Absolute magnitude of best-fitting galaxy
*                                    host model
      real*4        Ip_Ih      (nr__1) ! PSF-to-host intensity ratio
      real*4        Im_Ih      (nr__1) ! Model-to-host intensity ratio
      real*4        r1_2       (nr__1) ! (kpc) ? Half-light radius of n=4 Sersic component
      real*4        v_mu1_2_   (nr__1) ! (mag/arcsec2) ? Mean surface brightness of n=4 Sersic
*                                    component within half-light radius
*                                    in Hmag/arcsec^2^
      character*4   Bar        (nr__1) ! [Yes?/No ] Presence of a stellar bar?
      character*4   Arms       (nr__1) ! [Yes?/No ] Presence of spiral arms?
      character*4   Dis        (nr__1) ! [Yes?/No/Str ] Presence of merger-induced
*                                    disturbance (Str=strongly disturbed)
      character*3   MType      (nr__1) ! Morphological class (2)
      character*3   IC         (nr__1) ! Interaction class (IVb or V) (3)
*Note (1): Entries for this object are uncertain because the host galaxy
*     fills most of the field of view, so the sky background is not well
*     determined.
*Note (2): Morphological classes are:
*   E   = elliptical,
*   B+D = bulge+disk,
*   D   = disk,
*   A   = ambiguous.
*   Question marks ("?") indicate uncertain classification 
*   due to PSF mismatch (see sections 3 and 4.1).
*Note (3): IVb = compact merger, V = old merger; see Veilleux et al. 2002,
*   Cat. J/ApJS/143/315, section 3.4, for details on the definitions.

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

C  Declarations for 'table6.dat'	! Black Hole Mass

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

      character*7   Sample     (nr__2) ! [ULIRGs/PG QSOs]
      character*16  Name_2     (nr__2) ! Galaxy name
      real*8        Mph        (nr__2) ! (10+7Msun) ? Black hole mass from spheroid luminosity
      real*4        Mdyn       (nr__2) ! (10+7Msun) ? Black hole mass from spheroid velocity
*                                   dispersion
      real*4        Mrev       (nr__2) ! (10+7Msun) ? Black hole mass from reverberation mapping
      real*4        Mvir       (nr__2) ! (10+7Msun) ? Black hole mass from virial relation
      character*7   Ref        (nr__2) ! M_BH_ reference(s) (1)
*Note (1): References as follows:
*   1 = Dasyra et al. 2006ApJ...651..835D
*   2 = Veilleux et al. 2002, Cat. J/ApJS/143/315
*   3 = Veilleux et al. 2006ApJ...643..707V
*   4 = Dasyra et al. 2006ApJ...638..745D
*   5 = This paper;
*   6 = Peterson et al. 2004, Cat. J/ApJ/613/682
*   7 = Dasyra et al. 2007ApJ...657..102D
*   8 = Vestergaard & Peterson 2006ApJ...641..689V

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

C  Loading file 'table1.dat'	! PG QSOs sample

C  Format for file interpretation

    1 format(
     +  A11,1X,A10,1X,F5.3,1X,F6.2,1X,A1,F5.2,1X,F5.2,1X,A1,F5.2,1X,
     +  A1,F4.2,1X,A1,F5.3,A1,1X,A1,1X,A10)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,28
        read(1,'(A82)')ar__
        read(ar__,1)
     +  Name(i__),OName(i__),z(i__),logLB(i__),l_logLIR(i__),
     +  logLIR(i__),logLbol(i__),l_LIR_LB(i__),LIR_LB(i__),
     +  l_LIR_Lbol(i__),LIR_Lbol(i__),l_F25_F60(i__),F25_F60(i__),
     +  u_F25_F60(i__),RL(i__),Class(i__)
        if(ar__(38:42) .EQ. '') logLIR(i__) = rNULL__
        if(ar__(51:55) .EQ. '') LIR_LB(i__) = rNULL__
        if(ar__(58:61) .EQ. '') LIR_Lbol(i__) = rNULL__
c    ..............Just test output...........
        write(6,1)
     +  Name(i__),OName(i__),z(i__),logLB(i__),l_logLIR(i__),
     +  logLIR(i__),logLbol(i__),l_LIR_LB(i__),LIR_LB(i__),
     +  l_LIR_Lbol(i__),LIR_Lbol(i__),l_F25_F60(i__),F25_F60(i__),
     +  u_F25_F60(i__),RL(i__),Class(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table5.dat'	! *Summary results from the GALFIT analysis

C  Format for file interpretation

    2 format(
     +  A11,A1,1X,F6.2,1X,F6.2,1X,F6.2,1X,F6.2,1X,F5.2,1X,F4.2,1X,
     +  F4.2,1X,F5.2,1X,A4,1X,A4,1X,A4,1X,A3,1X,A3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table5.dat')
      write(6,*) '....Loading file: table5.dat'
      do i__=1,28
        read(1,'(A85)')ar__1
        read(ar__1,2)
     +  Name_1(i__),u_Name(i__),MAGt(i__),MAGp(i__),MAGh(i__),
     +  MAGm(i__),Ip_Ih(i__),Im_Ih(i__),r1_2(i__),v_mu1_2_(i__),
     +  Bar(i__),Arms(i__),Dis(i__),MType(i__),IC(i__)
        if(ar__1(53:56) .EQ. '') r1_2(i__) = rNULL__
        if(ar__1(58:62) .EQ. '') v_mu1_2_(i__) = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  Name_1(i__),u_Name(i__),MAGt(i__),MAGp(i__),MAGh(i__),
     +  MAGm(i__),Ip_Ih(i__),Im_Ih(i__),r1_2(i__),v_mu1_2_(i__),
     +  Bar(i__),Arms(i__),Dis(i__),MType(i__),IC(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table6.dat'	! Black Hole Mass

C  Format for file interpretation

    3 format(A7,1X,A16,1X,F7.2,1X,F6.3,1X,F6.2,1X,F6.2,1X,A7)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table6.dat')
      write(6,*) '....Loading file: table6.dat'
      do i__=1,217
        read(1,'(A61)')ar__2
        read(ar__2,3)
     +  Sample(i__),Name_2(i__),Mph(i__),Mdyn(i__),Mrev(i__),
     +  Mvir(i__),Ref(i__)
        if(ar__2(26:32) .EQ. '') Mph(i__) = rNULL__
        if(ar__2(34:39) .EQ. '') Mdyn(i__) = rNULL__
        if(ar__2(41:46) .EQ. '') Mrev(i__) = rNULL__
        if(ar__2(48:53) .EQ. '') Mvir(i__) = rNULL__
c    ..............Just test output...........
        write(6,3)
     +  Sample(i__),Name_2(i__),Mph(i__),Mdyn(i__),Mrev(i__),
     +  Mvir(i__),Ref(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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