FORTRAN Generation
(/./ftp/cats/J/ApJS/237/14)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJS/237/14 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-Sep-09
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/ApJS/237/14  KASI-Yonsei Deep Imaging Survey of Clusters (KYDISC)  (Oh+, 2018)
*================================================================================
*KYDISC: galaxy morphology, quenching, and mergers in the cluster environment.
*    Oh S., Kim K., Lee J.H., Sheen Y.-K., Kim M., Ree C.H., Ho L.C., Kyeong J.,
*    Sung E.-C., Park B.-G., Yi S.K.
*   <Astrophys. J. Suppl. Ser., 237, 14 (2018)>
*   =2018ApJS..237...14O
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'table3.dat'	! Cluster characteristics

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

      character*10  Cl         (nr__) ! Cluster name
      integer*4     Vcen       (nr__) ! (km/s) [4709/43255] Central radial velocity (1)
      integer*4     sigma      (nr__) ! (km/s) [433/992] Radial velocity dispersion
*                                   ({sigma}_obs_) (1)
      integer*4     e_sigma    (nr__) ! (km/s) [30/110] sigma uncertainty
      real*4        R200       (nr__) ! (Mpc) [1/2.3] Virial radius R_200_ (see Equation 1)
      real*4        e_R200     (nr__) ! (Mpc) [0.08/0.3] R200 uncertainty
      real*4        M200       (nr__) ! (10+14Msun) [1/12.4] Cluster mass M_200_ (see Equation 2)
      real*4        e_M200     (nr__) ! (10+14Msun) [0.2/4.3] Negative uncertainty on M200
      real*4        E_M200_1   (nr__) ! (10+14Msun) [0.2/3.4] Positive uncertainty on M200
      integer*4     Nm         (nr__) ! [14/250] Number of spectroscopically
*                                   confirmed member galaxies (2)
      character*5   ACl        (nr__) ! Abbreviated cluster name as used in Table 10;
*                                   column added by CDS
*Note (1): The radial velocity dispersion ({sigma}_obs_) and the central radial
*          velocity (Vcen) of clusters are estimated from a Gaussian fit of the
*          distribution of radial velocities (Figure 3). We select galaxies as
*          possible members of each cluster when they have a radial velocity
*          (v_rv_) within +/-3{sigma}_obs_ from the central radial velocity
*          (Vcen).
*Note (2): In summary, we selected galaxies with |v_rv_-Vcen|<3{sigma}_obs_,
*          R<2R_200_, and M_r'_{<}-19.8 as the members of each cluster.
*          See section 3 for further explanations.

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

C  Declarations for 'table10.dat'	! KYDISC catalog

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

C  J2000 position composed of: RAdeg DEdeg
      character*5   ACl_1      (nr__1) ! Abbreviated cluster name
      integer*4     ID         (nr__1) ! [1/2894] Galaxy ID within the cluster
      real*8        RAdeg      (nr__1) ! (deg) Right accession in decimal degrees (J2000)
      real*8        DEdeg      (nr__1) ! (deg) Declination in decimal degrees (J2000)
      real*8        z          (nr__1) ! [0.011/0.152] Spectroscopic redshift
      integer*4     r_z        (nr__1) ! [0/6]? Origin of z (1)
      real*8        e_z        (nr__1) ! [0.00001/0.002]?=-99 Error in z
      real*8        umag       (nr__1) ! (mag) [15/25.2]?=-99 Extinction-corrected
*                                 u' band magnitude
      real*8        gmag       (nr__1) ! (mag) [13/20.7] Extinction-corrected g' band magnitude
      real*8        rmag       (nr__1) ! (mag) [12/19.7] Extinction-corrected r' band magnitude
      real*8        e_umag     (nr__1) ! (mag) [0.001/1.1]?=-99 Error in umag
      real*4        e_gmag     (nr__1) ! (mag) [0.0003/0.04] Error in gmag
      real*4        e_rmag     (nr__1) ! (mag) [0.0002/0.02] Error in rmag
      real*8        ukcor      (nr__1) ! (mag) [-4/0.6]?=-99 k-correction for u' band magnitude
      real*8        gkcor      (nr__1) ! (mag) [-0.02/0.5] k-correction for g' band magnitude
      real*8        rkcor      (nr__1) ! (mag) [-0.05/0.2] k-correction for r' band magnitude
      real*4        re         (nr__1) ! (arcsec) [0.3/29.6] Circularized effective radius
*                                 in r' band
      real*4        ee         (nr__1) ! [0.004/0.8] Ellipticity at re
      real*8        PAe        (nr__1) ! (deg) Position angle at re; N=90, E=0
      real*8        r90        (nr__1) ! (arcsec) [0.7/125.2] Circularized radius (2)
      real*4        e90        (nr__1) ! [0.004/0.8] Ellipticity at r90
      real*8        PA90       (nr__1) ! (deg) Position angle at r90; N=90, E=0
      real*4        R_R200     (nr__1) ! [0/2] Cluster-centric distance
      real*4        logD10     (nr__1) ! ([-]) [-1/3.1] log local density
      real*4        B_T        (nr__1) ! [0.01/1] Bulge-to-total ratio in r' band
*                                 using galfit
      real*4        A          (nr__1) ! [0.02/0.8] CAS asymmetry parameter (detailed in
*                                 Conselice+, 2013, J/ApJS/147/1) in r' band
      integer*4     HType      (nr__1) ! [0/3] Hubble type; E=0, S0=1, SE=2, and SL=3
      integer*4     SigH       (nr__1) ! [2/5] Significance of the Hubble classification
      integer*4     MType      (nr__1) ! [0/2]? Merger type; Normal=0, OM=1, and PM=2 (3)
      integer*4     SigM       (nr__1) ! [2/5] Significance of the Merger classification
*Note (1): Origin of redshift as follows:
*    0 = the 2.5m du Pont Telescopes located at Las Campanas Observatory, Chile
*         (see section 2.4.2);
*    1 = the Inamori Magellan Areal Camera and Spectrograph (IMACS) on the 6.5m
*         Magellan Baade telescope (see section 2.4.1);
*    2 = the Hydra on the WIYN 3.5m telescope (see section 2.4.3);
*    3 = HeCS (Hectospec Cluster Survey; Rines+ 2013, J/ApJ/767/15);
*    4 = the NASA/IPAC Extragalactic Database (NED);
*    5 = SDSS (see V/147);
*    6 = 6dF (see VII/259).
*          See also the "Description" section above.
*Note (2): Containing 90% of total light in r' band.
*Note (3): Type as follows:
*    1 = ongoing mergers (OM) with close pairs within four times <r_90_>
*         (53 occurrences)
*    2 = galaxies with features seemingly related to mergers in the past without
*         relevant pairs for the features ("PM"; 279 occurrences).
*    0 = Normal (1077 galaxies).
*          See section 4.2.2.

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

C  Loading file 'table3.dat'	! Cluster characteristics

C  Format for file interpretation

    1 format(
     +  A10,1X,I5,1X,I3,1X,I3,1X,F4.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F4.2,
     +  1X,I3,1X,A5)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,14
        read(1,'(A60)')ar__
        read(ar__,1)
     +  Cl(i__),Vcen(i__),sigma(i__),e_sigma(i__),R200(i__),
     +  e_R200(i__),M200(i__),e_M200(i__),E_M200_1(i__),Nm(i__),
     +  ACl(i__)
c    ..............Just test output...........
        write(6,1)
     +  Cl(i__),Vcen(i__),sigma(i__),e_sigma(i__),R200(i__),
     +  e_R200(i__),M200(i__),e_M200(i__),E_M200_1(i__),Nm(i__),
     +  ACl(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'table10.dat'	! KYDISC catalog

C  Format for file interpretation

    2 format(
     +  A5,1X,I5,1X,F11.7,1X,F11.7,1X,F7.5,1X,I1,1X,F9.5,1X,F8.4,1X,
     +  F7.4,1X,F7.4,1X,F9.4,1X,F6.4,1X,F6.4,1X,F9.4,1X,F7.4,1X,F7.4,
     +  1X,F6.3,1X,F5.3,1X,F7.3,1X,F7.3,1X,F5.3,1X,F7.3,1X,F5.3,1X,
     +  F5.2,1X,F5.3,1X,F5.2,1X,I1,1X,I1,1X,I1,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table10.dat')
      write(6,*) '....Loading file: table10.dat'
      do i__=1,1409
        read(1,'(A205)')ar__1
        read(ar__1,2)
     +  ACl_1(i__),ID(i__),RAdeg(i__),DEdeg(i__),z(i__),r_z(i__),
     +  e_z(i__),umag(i__),gmag(i__),rmag(i__),e_umag(i__),
     +  e_gmag(i__),e_rmag(i__),ukcor(i__),gkcor(i__),rkcor(i__),
     +  re(i__),ee(i__),PAe(i__),r90(i__),e90(i__),PA90(i__),
     +  R_R200(i__),logD10(i__),B_T(i__),A(i__),HType(i__),SigH(i__),
     +  MType(i__),SigM(i__)
        if(ar__1(45:45) .EQ. '') r_z(i__) = iNULL__
        if(ar__1(203:203) .EQ. '') MType(i__) = iNULL__
c    ..............Just test output...........
        write(6,2)
     +  ACl_1(i__),ID(i__),RAdeg(i__),DEdeg(i__),z(i__),r_z(i__),
     +  e_z(i__),umag(i__),gmag(i__),rmag(i__),e_umag(i__),
     +  e_gmag(i__),e_rmag(i__),ukcor(i__),gkcor(i__),rkcor(i__),
     +  re(i__),ee(i__),PAe(i__),r90(i__),e90(i__),PA90(i__),
     +  R_R200(i__),logD10(i__),B_T(i__),A(i__),HType(i__),SigH(i__),
     +  MType(i__),SigM(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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