FORTRAN Generation
(/./ftp/cats/J/AJ/123/2206)

Conversion of standardized ReadMe file for file /./ftp/cats/J/AJ/123/2206 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-29
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/AJ/123/2206       BRI photometry of 3 QSO fields           (Prochaska+, 2002)
*================================================================================
*Galaxies associated with z~4 damped Ly {alpha} systems.
*I. Imaging and photometric selection.
*    Prochaska J.X., Gawiser E., Wolfe A.M., Quirrenbach A., Lanzetta K.M.,
*    Chen H.-W., Cooke J., Yahata N.
*   <Astron. J. 123, 2206 (2002)>
*   =2002AJ....123.2206P
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table4.dat'	! Photometry for PSS J0132+1341

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

C  J2000 position composed of: RAdeg DEdeg
      integer*4     v_PGW2002_ (nr__) ! Identification (1)
      real*4        Xpix       (nr__) ! (pix) The X pixel coordinate
      real*4        Ypix       (nr__) ! (pix) The Y pixel coordinate
      real*8        RAdeg      (nr__) ! (deg) Decimal degree of Right Ascension (J2000)
      real*8        DEdeg      (nr__) ! (deg) Decimal degree of Declination (J2000)
      character*1   l_Bmag     (nr__) ! Limit flag on Bmag (2)
      real*4        Bmag       (nr__) ! (mag) The B_AB_ band magnitude
      real*4        e_Bmag     (nr__) ! (mag) The 1{sigma} error in Bmag (2)
      real*4        Rmag       (nr__) ! (mag) The R_AB_ band magnitude
      real*4        e_Rmag     (nr__) ! (mag) The 1{sigma} error in Rmag (2)
      character*1   l_Imag     (nr__) ! Limit flag on Imag (2)
      real*4        Imag       (nr__) ! (mag) The I_AB_ band magnitude
      real*4        e_Imag     (nr__) ! (mag) The 1{sigma} error in Imag (2)
      real*4        z          (nr__) ! The most likely photometric redshift
      integer*4     Tmpl       (nr__) ! Best fit photometry data template number (3)
      real*4        Pz         (nr__) ! Confidence level, as the probability that the 
*                                  actual redshift is in interval [z-0.1,z+0.1]
      real*4        Pz4        (nr__) ! Probability of an actual redshift in interval
*                                  [3.5,4.5]
*Note (1): Identification:
*     Table 4: [PGW2002] PSS 0132+1341 NNNN (Nos 1-3035).
*     Table 5: [PGW2002] BR 0951-0450  NNNN (Nos 1-3804).
*     Table 6: [PGW2002] PSS 1443+2724 NNNN (Nos 1-2724).
*Note (2): Objects with fluxes <3{sigma}_sky_ have reported upper limits
*    equal to the maximum of the flux + 1{sigma}_sky_ for positive fluxes,
*    or setting its value to 1{sigma}_sky_ for negative fluxes
*    (see section 4 of the paper). 
*    These objects have the error in the magnitude set to 9.99.
*Note (3): Best fit photometry data template number:
*    1 = Elliptical
*    2 = Sab
*    3 = Scd
*    4 = Irr
*    5 = SB1
*    6 = SB2
*    7 = QSO
*   >7 = Stars

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

C  Declarations for 'table5.dat'	! Photometry for BR 0951-0450

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

C  J2000 position composed of: RAdeg DEdeg
      integer*4     v_PGW2002__1(nr__1) ! Identification (1)
      real*4        Xpix_1     (nr__1) ! (pix) The X pixel coordinate
      real*4        Ypix_1     (nr__1) ! (pix) The Y pixel coordinate
      real*8        RAdeg_1    (nr__1) ! (deg) Decimal degree of Right Ascension (J2000)
      real*8        DEdeg_1    (nr__1) ! (deg) Decimal degree of Declination (J2000)
      character*1   l_Bmag_1   (nr__1) ! Limit flag on Bmag (2)
      real*4        Bmag_1     (nr__1) ! (mag) The B_AB_ band magnitude
      real*4        e_Bmag_1   (nr__1) ! (mag) The 1{sigma} error in Bmag (2)
      real*4        Rmag_1     (nr__1) ! (mag) The R_AB_ band magnitude
      real*4        e_Rmag_1   (nr__1) ! (mag) The 1{sigma} error in Rmag (2)
      character*1   l_Imag_1   (nr__1) ! Limit flag on Imag (2)
      real*4        Imag_1     (nr__1) ! (mag) The I_AB_ band magnitude
      real*4        e_Imag_1   (nr__1) ! (mag) The 1{sigma} error in Imag (2)
      real*4        z_1        (nr__1) ! The most likely photometric redshift
      integer*4     Tmpl_1     (nr__1) ! Best fit photometry data template number (3)
      real*4        Pz_1       (nr__1) ! Confidence level, as the probability that the 
*                                  actual redshift is in interval [z-0.1,z+0.1]
      real*4        Pz4_1      (nr__1) ! Probability of an actual redshift in interval
*                                  [3.5,4.5]
*Note (1): Identification:
*     Table 4: [PGW2002] PSS 0132+1341 NNNN (Nos 1-3035).
*     Table 5: [PGW2002] BR 0951-0450  NNNN (Nos 1-3804).
*     Table 6: [PGW2002] PSS 1443+2724 NNNN (Nos 1-2724).
*Note (2): Objects with fluxes <3{sigma}_sky_ have reported upper limits
*    equal to the maximum of the flux + 1{sigma}_sky_ for positive fluxes,
*    or setting its value to 1{sigma}_sky_ for negative fluxes
*    (see section 4 of the paper). 
*    These objects have the error in the magnitude set to 9.99.
*Note (3): Best fit photometry data template number:
*    1 = Elliptical
*    2 = Sab
*    3 = Scd
*    4 = Irr
*    5 = SB1
*    6 = SB2
*    7 = QSO
*   >7 = Stars

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

C  Declarations for 'table6.dat'	! Photometry for PSS J1443+2724

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

C  J2000 position composed of: RAdeg DEdeg
      integer*4     v_PGW2002__2(nr__2) ! Identification (1)
      real*4        Xpix_2     (nr__2) ! (pix) The X pixel coordinate
      real*4        Ypix_2     (nr__2) ! (pix) The Y pixel coordinate
      real*8        RAdeg_2    (nr__2) ! (deg) Decimal degree of Right Ascension (J2000)
      real*8        DEdeg_2    (nr__2) ! (deg) Decimal degree of Declination (J2000)
      character*1   l_Bmag_2   (nr__2) ! Limit flag on Bmag (2)
      real*4        Bmag_2     (nr__2) ! (mag) The B_AB_ band magnitude
      real*4        e_Bmag_2   (nr__2) ! (mag) The 1{sigma} error in Bmag (2)
      real*4        Rmag_2     (nr__2) ! (mag) The R_AB_ band magnitude
      real*4        e_Rmag_2   (nr__2) ! (mag) The 1{sigma} error in Rmag (2)
      character*1   l_Imag_2   (nr__2) ! Limit flag on Imag (2)
      real*4        Imag_2     (nr__2) ! (mag) The I_AB_ band magnitude
      real*4        e_Imag_2   (nr__2) ! (mag) The 1{sigma} error in Imag (2)
      real*4        z_2        (nr__2) ! The most likely photometric redshift
      integer*4     Tmpl_2     (nr__2) ! Best fit photometry data template number (3)
      real*4        Pz_2       (nr__2) ! Confidence level, as the probability that the 
*                                  actual redshift is in interval [z-0.1,z+0.1]
      real*4        Pz4_2      (nr__2) ! Probability of an actual redshift in interval
*                                  [3.5,4.5]
*Note (1): Identification:
*     Table 4: [PGW2002] PSS 0132+1341 NNNN (Nos 1-3035).
*     Table 5: [PGW2002] BR 0951-0450  NNNN (Nos 1-3804).
*     Table 6: [PGW2002] PSS 1443+2724 NNNN (Nos 1-2724).
*Note (2): Objects with fluxes <3{sigma}_sky_ have reported upper limits
*    equal to the maximum of the flux + 1{sigma}_sky_ for positive fluxes,
*    or setting its value to 1{sigma}_sky_ for negative fluxes
*    (see section 4 of the paper). 
*    These objects have the error in the magnitude set to 9.99.
*Note (3): Best fit photometry data template number:
*    1 = Elliptical
*    2 = Sab
*    3 = Scd
*    4 = Irr
*    5 = SB1
*    6 = SB2
*    7 = QSO
*   >7 = Stars

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

C  Loading file 'table4.dat'	! Photometry for PSS J0132+1341

C  Format for file interpretation

    1 format(
     +  I4,1X,F6.1,1X,F6.1,1X,F10.6,1X,F10.6,1X,A1,1X,F5.2,1X,F4.2,1X,
     +  F5.2,1X,F4.2,1X,A1,1X,F5.2,1X,F4.2,1X,F4.2,1X,I3,1X,F4.2,1X,
     +  F4.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table4.dat')
      write(6,*) '....Loading file: table4.dat'
      do i__=1,3035
        read(1,'(A97)')ar__
        read(ar__,1)
     +  v_PGW2002_(i__),Xpix(i__),Ypix(i__),RAdeg(i__),DEdeg(i__),
     +  l_Bmag(i__),Bmag(i__),e_Bmag(i__),Rmag(i__),e_Rmag(i__),
     +  l_Imag(i__),Imag(i__),e_Imag(i__),z(i__),Tmpl(i__),Pz(i__),
     +  Pz4(i__)
c    ..............Just test output...........
        write(6,1)
     +  v_PGW2002_(i__),Xpix(i__),Ypix(i__),RAdeg(i__),DEdeg(i__),
     +  l_Bmag(i__),Bmag(i__),e_Bmag(i__),Rmag(i__),e_Rmag(i__),
     +  l_Imag(i__),Imag(i__),e_Imag(i__),z(i__),Tmpl(i__),Pz(i__),
     +  Pz4(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table5.dat'	! Photometry for BR 0951-0450

C  Format for file interpretation

    2 format(
     +  I4,1X,F6.1,1X,F6.1,1X,F10.6,1X,F10.6,1X,A1,1X,F5.2,1X,F4.2,1X,
     +  F5.2,1X,F4.2,1X,A1,1X,F5.2,1X,F4.2,1X,F4.2,1X,I3,1X,F4.2,1X,
     +  F4.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table5.dat')
      write(6,*) '....Loading file: table5.dat'
      do i__=1,3804
        read(1,'(A97)')ar__1
        read(ar__1,2)
     +  v_PGW2002__1(i__),Xpix_1(i__),Ypix_1(i__),RAdeg_1(i__),
     +  DEdeg_1(i__),l_Bmag_1(i__),Bmag_1(i__),e_Bmag_1(i__),
     +  Rmag_1(i__),e_Rmag_1(i__),l_Imag_1(i__),Imag_1(i__),
     +  e_Imag_1(i__),z_1(i__),Tmpl_1(i__),Pz_1(i__),Pz4_1(i__)
c    ..............Just test output...........
        write(6,2)
     +  v_PGW2002__1(i__),Xpix_1(i__),Ypix_1(i__),RAdeg_1(i__),
     +  DEdeg_1(i__),l_Bmag_1(i__),Bmag_1(i__),e_Bmag_1(i__),
     +  Rmag_1(i__),e_Rmag_1(i__),l_Imag_1(i__),Imag_1(i__),
     +  e_Imag_1(i__),z_1(i__),Tmpl_1(i__),Pz_1(i__),Pz4_1(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table6.dat'	! Photometry for PSS J1443+2724

C  Format for file interpretation

    3 format(
     +  I4,1X,F6.1,1X,F6.1,1X,F10.6,1X,F10.6,1X,A1,1X,F5.2,1X,F4.2,1X,
     +  F5.2,1X,F4.2,1X,A1,1X,F5.2,1X,F4.2,1X,F4.2,1X,I3,1X,F4.2,1X,
     +  F4.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table6.dat')
      write(6,*) '....Loading file: table6.dat'
      do i__=1,2667
        read(1,'(A97)')ar__2
        read(ar__2,3)
     +  v_PGW2002__2(i__),Xpix_2(i__),Ypix_2(i__),RAdeg_2(i__),
     +  DEdeg_2(i__),l_Bmag_2(i__),Bmag_2(i__),e_Bmag_2(i__),
     +  Rmag_2(i__),e_Rmag_2(i__),l_Imag_2(i__),Imag_2(i__),
     +  e_Imag_2(i__),z_2(i__),Tmpl_2(i__),Pz_2(i__),Pz4_2(i__)
c    ..............Just test output...........
        write(6,3)
     +  v_PGW2002__2(i__),Xpix_2(i__),Ypix_2(i__),RAdeg_2(i__),
     +  DEdeg_2(i__),l_Bmag_2(i__),Bmag_2(i__),e_Bmag_2(i__),
     +  Rmag_2(i__),e_Rmag_2(i__),l_Imag_2(i__),Imag_2(i__),
     +  e_Imag_2(i__),z_2(i__),Tmpl_2(i__),Pz_2(i__),Pz4_2(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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