FORTRAN Generation
(/./ftp/cats/J/ApJ/595/685)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/595/685 into FORTRAN code for reading data files line by line.

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-Apr-20
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/595/685       HST imaging in Chandra Deep Field-South    (Grogin+, 2003)
*================================================================================
*Hubble Space Telescope imaging in the Chandra Deep Field-South.
*III. Quantitative morphology of the 1 million second Chandra counterparts and
*comparison with the field population.
*    Grogin N.A., Koekemoer A.M., Schreier E.J., Bergeron J., Giacconi R.,
*    Hasinger G., Kewley L., Norman C., Rosati P., Tozzi P., Zirm A.
*   <Astrophys. J. 595, 685 (2003)>
*   =2003ApJ...595..685G
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table1.dat'	! WFPC2 Source Catalog: Coordinates, Magnitudes,
                           and Shape Parameters

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

C  J2000 position composed of: RAdeg DEdeg
      real*8        RAdeg       ! (deg) Right Ascension in decimal degrees (J2000)
      real*8        DEdeg       ! (deg) Declination in decimal degrees (J2000)
      character*1   l_F606W     ! Limit flag on Vmag
      real*4        F606W       ! (mag) The WFPC2 (F606W) V band magnitude
      character*1   l_F814W     ! Limit flag on Imag
      real*4        F814W       ! (mag) The WFPC2 (F814W) I band magnitude
      real*4        R0_5        ! (arcsec) The I band half-light radius
      real*4        etaV        ! The V band stellarity measure (G1)
      real*4        etaI        ! The I band stellarity measure (G1)
      real*4        AI          ! The I band asymmetry index
      character*4   CV          ! The V band concentration index
      character*4   CI          ! The I band concentration index

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

C  Declarations for 'table2.dat'	! WFPC2 CDFS Counterparts: X-ray and optical properties

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

      integer*4     XID         ! Giacconi et al. (2002,
*                                    Cat. <J/ApJS/139/369>) designation,
*                                    [GZW2002] XID NNN in Simbad
      character*16  IAU         ! IAU designation (JHHMMSS.s-DDMMSS)
      character*1   l_logF2_10keV ! Limit flag on logF2-10keV (1)
      real*4        logF2_10keV ! ([aW/m2]) X-ray flux in the hard (2-10keV) band
*                                    (exprssed in 10-15erg/cm2/s)
      real*4        HR          ! Hardness ratio (2)
      real*4        F606W_1     ! (mag) F606W (V) magnitude
      real*4        F814W_1     ! (mag) F814W (I) magnitude
      real*4        R0_5_1      ! (arcsec) The I band half-light radius
      real*4        etaV_1      ! The V band stellarity measure (G1)
      real*4        etaI_1      ! The I band stellarity measure (G1)
      character*1   l_Fn_FtV    ! [~] Limit flag on Fn/FtV (3)
      real*8        Fn_FtV      ! ? Nuclear-to-total flux ratio in V band
      character*1   l_E_Fn_FtV  ! Limit flag on E_Fn/FtV (4)
      real*8        E_Fn_FtV    ! ? Error in Fn/FtV (upper limit)
      real*4        e_Fn_FtV_1  ! ? Error in Fn/FtV (lower limit)
      character*1   l_Fn_FtI    ! [~] Limit flag on Fn/FtI
      real*4        Fn_FtI      ! ? Nuclear-to-total flux ratio in I band
      character*1   l_E_Fn_FtI  ! Limit flag on E_Fn/FtI (4)
      real*4        E_Fn_FtI    ! ? Error in Fn/FtI (upper limit)
      real*4        e_Fn_FtI_1  ! ? Error in Fn/FtI (lower limit)
      real*4        AI_1        ! The I band asymmetry index
      character*1   u_AI        ! Uncertainty flag on AI
      real*4        e_AI        ! rms uncertainty on AI
      real*4        CI_1        ! The I band concentration index
      real*4        e_CI        ! rms uncertainty on CI
      character*30  Comm        ! Comments
*Note (1): Sources detected only in CXO soft band (0.5-2keV)
*     indicated with upper limits.
*Note (2): HR (H-S)/(H+S), where H=2-10keV counts, S=0.5-2keV counts
*Note (3): When Fn/FtV=~1, HST-unresolved sources omitted from
*    two-dimensional fitting, plotted with Fnuc~Ftot in fig.4 of the paper.
*Note (4): If 0 is not excluded at >3{sigma}, only 3{sigma} upper limit is shown.

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

C  Loading file 'table1.dat'	! WFPC2 Source Catalog: Coordinates, Magnitudes,
*                           and Shape Parameters

C  Format for file interpretation

    1 format(
     +  F9.6,1X,F10.6,1X,A1,1X,F5.2,1X,A1,1X,F5.2,1X,F5.2,1X,F4.2,1X,
     +  F4.2,1X,F5.2,1X,A4,1X,A4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table1.dat')
      write(6,*) '....Loading file: table1.dat'
      do i__=1,3030
        read(1,'(A69)')ar__
        read(ar__,1)
     +  RAdeg,DEdeg,l_F606W,F606W,l_F814W,F814W,R0_5,etaV,etaI,AI,CV,
     +  CI
c    ..............Just test output...........
        write(6,1)
     +  RAdeg,DEdeg,l_F606W,F606W,l_F814W,F814W,R0_5,etaV,etaI,AI,CV,
     +  CI
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table2.dat'	! WFPC2 CDFS Counterparts: X-ray and optical properties

C  Format for file interpretation

    2 format(
     +  I3,1X,A16,1X,A1,F6.2,1X,F5.2,1X,F5.2,1X,F5.2,1X,F4.2,1X,F4.2,
     +  1X,F4.2,1X,A1,F7.5,1X,A1,F7.5,F6.4,1X,A1,F6.4,1X,A1,F6.4,1X,
     +  F6.4,3X,F6.3,A1,1X,F5.3,1X,F4.2,1X,F4.2,1X,A30)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,37
        read(1,'(A165)')ar__1
        read(ar__1,2)
     +  XID,IAU,l_logF2_10keV,logF2_10keV,HR,F606W_1,F814W_1,R0_5_1,
     +  etaV_1,etaI_1,l_Fn_FtV,Fn_FtV,l_E_Fn_FtV,E_Fn_FtV,e_Fn_FtV_1,
     +  l_Fn_FtI,Fn_FtI,l_E_Fn_FtI,E_Fn_FtI,e_Fn_FtI_1,AI_1,u_AI,e_AI,
     +  CI_1,e_CI,Comm
        if(ar__1(64:70) .EQ. '') Fn_FtV = rNULL__
        if(ar__1(73:79) .EQ. '') E_Fn_FtV = rNULL__
        if(ar__1(80:85) .EQ. '') e_Fn_FtV_1 = rNULL__
        if(ar__1(88:93) .EQ. '') Fn_FtI = rNULL__
        if(ar__1(96:101) .EQ. '') E_Fn_FtI = rNULL__
        if(ar__1(103:108) .EQ. '') e_Fn_FtI_1 = rNULL__
c    ..............Just test output...........
        write(6,2)
     +  XID,IAU,l_logF2_10keV,logF2_10keV,HR,F606W_1,F814W_1,R0_5_1,
     +  etaV_1,etaI_1,l_Fn_FtV,Fn_FtV,l_E_Fn_FtV,E_Fn_FtV,e_Fn_FtV_1,
     +  l_Fn_FtI,Fn_FtI,l_E_Fn_FtI,E_Fn_FtI,e_Fn_FtI_1,AI_1,u_AI,e_AI,
     +  CI_1,e_CI,Comm
c    .......End.of.Just test output...........
      end do
      close(1)

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