FORTRAN Generation
(/./ftp/cats/J/ApJ/723/251)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/723/251 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-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/723/251    Keck/LRIS confirmation of Coma membership    (Chiboucas+, 2010)
*================================================================================
*Keck/LRIS spectroscopic confirmation of Coma cluster dwarf galaxy membership
*assignments.
*    Chiboucas K., Tully R.B., Marzke R.O., Trentham N., Ferguson H.C.,
*    Hammer D., Carter D., Khosroshahi H.
*   <Astrophys. J., 723, 251-266 (2010)>
*   =2010ApJ...723..251C
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table2.dat'	! LRIS (membership test sample) redshift results

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

C  J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs
      real*8        RAdeg       ! (deg) Right Ascension J2000
      real*8        DEdeg       ! (deg)     Declination J2000
C  ---------------------------------- ! (position vector(s) in degrees)

      character*2   S           ! Status (Mm, Bg, Un, PM, or PB) (1)
      integer*4     ID          ! [10162/4042826] Identification number
      character*1   f_ID        ! [he] Note on the observation (2)
      integer*4     PT          ! ?=x [1/5] Membership probability defined by
*                                  N.T. on HST/ACS image (3)
      integer*4     PF          ! ?=x [1/4] Membership probability defined by
*                                  H.F. on HST/ACS image (3)
      real*4        Rmag        ! (mag) [18.68/23.41]? R-band magnitude (4)
      character*1   f_Rmag      ! [f] Blended photometry (2)
      real*4        F814Wc      ! (mag) ? HST/ACS F814W magnitude; corrected (5)
      real*4        F814W       ! (mag) ? HST/ACS F814W magnitude (5)
      real*4        B_V         ! (mag) ? (B-V) color index (3" aperture) (4)
      character*1   f_B_V       ! [f] Blended photometry (2)
      real*4        g_I         ! (mag) ? (F475W-F814W) color index (2.25" aperture)
      real*4        mu0         ! (mag/arcsec2) ? Central surface brightness ({mu}_0_(R))
      real*4        mue         ! (mag/arcsec2) ? Mean effective surface brightness
*                                  (<{mu}_e,F814W_>)
      integer*4     RAh         ! (h) [12/13] Hour of right ascension (J2000)
      integer*4     RAm         ! (min) Minute of right ascension (J2000)
      real*4        RAs         ! (s) Second of right ascension (J2000)
      character*1   DE_         ! [+] Sign of declination (J2000)
      integer*4     DEd         ! (deg) [27/28] Degree of declination (J2000)
      integer*4     DEm         ! (arcmin) Arcminute of declination (J2000)
      real*4        DEs         ! (arcsec) Arcsecond of declination (J2000)
      integer*4     czA         ! (km/s) ? Absorption redshift
      integer*4     e_czA       ! (km/s) ? czA uncertainty (6)
      real*4        Rfx         ! (aW/m2/nm) ? Relative flux (in 10^-16^erg/s/cm2/{AA})
      integer*4     czE         ! (km/s) ? Emission redshift
      character*1   f_czE       ! [g] redshift from red-side (2)
      integer*4     e_czE       ! (km/s) ? czE uncertainty
      real*4        S_N         ! [3.7/61.1]? Signal to noise ratio per
*                                  Angstrom around 5000{AA}
      character*3   MType       ! Morphological type
*Note (1): Membership status as follows:
*    Mm = Member
*    Bg = Background
*    Un = Uncertain
*    PM = Possible member (failed)
*    PB = Possible background (failed)
*Note (2): Flag as follows:
*     e = Object was observed in a separate run in 2009.
*     f = Photometry is for two blended components including the UCD target and
*         another source centered 0.9 arcsec west of it.
*     g = Redshift measurement from red-side data.
*     h = Object was observed in two masks.
*Note (3): Membership probability classes of the authors (N.Trenton and
*          H.Ferguson) as follows:
*     1 = expected member
*     2 = probable member
*     3 = possible member
*     4 = expected background
*     5 = the galaxy was an extra object that landed in the slit
*     x = no membership probability assigned
*          See section 2.1.1 for further explanations.
*Note (4): BVR photometry from catalog of Adami et al. (2006,
*          Cat. J/A+A/451/1159).
*Note (5): F814W photometry comes from SExtractor measurements from our ACS data
*          (Hammer et al. 2010ApJS..191..143H).
*Note (6): Errors include contributions from the cross-correlation measurement
*          uncertainty and from the uncertainties in the sky line wavelength
*          calibration shift.

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

C  Loading file 'table2.dat'	! LRIS (membership test sample) redshift results

C  Format for file interpretation

    1 format(
     +  A2,1X,I7,A1,1X,I1,1X,I1,1X,F5.2,A1,1X,F5.2,1X,F5.2,1X,F4.2,A1,
     +  1X,F5.2,1X,F5.2,1X,F5.2,1X,I2,1X,I2,1X,F5.2,1X,A1,I2,1X,I2,1X,
     +  F5.2,1X,I6,1X,I3,1X,F4.1,1X,I6,A1,1X,I3,1X,F4.1,1X,A3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table2.dat')
      write(6,*) '....Loading file: table2.dat'
      do i__=1,101
        read(1,'(A120)')ar__
        read(ar__,1)
     +  S,ID,f_ID,PT,PF,Rmag,f_Rmag,F814Wc,F814W,B_V,f_B_V,g_I,mu0,
     +  mue,RAh,RAm,RAs,DE_,DEd,DEm,DEs,czA,e_czA,Rfx,czE,f_czE,e_czE,
     +  S_N,MType
        if (PT .EQ. 120) PT =  iNULL__
        if (PF .EQ. 120) PF =  iNULL__
        if(ar__(17:21) .EQ. '') Rmag = rNULL__
        if(ar__(24:28) .EQ. '') F814Wc = rNULL__
        if(ar__(30:34) .EQ. '') F814W = rNULL__
        if(ar__(36:39) .EQ. '') B_V = rNULL__
        if(ar__(42:46) .EQ. '') g_I = rNULL__
        if(ar__(48:52) .EQ. '') mu0 = rNULL__
        if(ar__(54:58) .EQ. '') mue = rNULL__
        if(ar__(85:90) .EQ. '') czA = iNULL__
        if(ar__(92:94) .EQ. '') e_czA = iNULL__
        if(ar__(96:99) .EQ. '') Rfx = rNULL__
        if(ar__(101:106) .EQ. '') czE = iNULL__
        if(ar__(109:111) .EQ. '') e_czE = iNULL__
        if(ar__(113:116) .EQ. '') S_N = rNULL__
        RAdeg = rNULL__
        DEdeg = rNULL__
c  Derive coordinates RAdeg and DEdeg from input data
c  (RAdeg and DEdeg are set to rNULL__ when unknown)
        if(RAh .GT. -180) RAdeg=RAh*15.
        if(RAm .GT. -180) RAdeg=RAdeg+RAm/4.
        if(RAs .GT. -180) RAdeg=RAdeg+RAs/240.
        if(DEd .GE. 0) DEdeg=DEd
        if(DEm .GE. 0) DEdeg=DEdeg+DEm/60.
        if(DEs .GE. 0) DEdeg=DEdeg+DEs/3600.
        if(DE_.EQ.'-'.AND.DEdeg.GE.0) DEdeg=-DEdeg
c    ..............Just test output...........
        write(6,1)
     +  S,ID,f_ID,PT,PF,Rmag,f_Rmag,F814Wc,F814W,B_V,f_B_V,g_I,mu0,
     +  mue,RAh,RAm,RAs,DE_,DEd,DEm,DEs,czA,e_czA,Rfx,czE,f_czE,e_czE,
     +  S_N,MType
        write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
c    .......End.of.Just test output...........
      end do
      close(1)

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