FORTRAN Generation
(/./ftp/cats/J/ApJ/925/34)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJ/925/34 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/ApJ/925/34        The COSMOS-DASH morphological catalog        (Cutler+, 2022)
*================================================================================
*Diagnosing DASH: a catalog of structural properties for the COSMOS-DASH survey.
*    Cutler S.E., Whitaker K.E., Mowla L.A., Brammer G.B., van der Wel A.,
*    Marchesini D., van Dokkum P.G., Momcheva I.G., Song M., Akhshik M.,
*    Nelson E.J., Bezanson R., Franx M., Kriek M., Lange-Vagle D., Leja J.,
*    MacKenty J.W., Muzzin A., Shipley H.
*   <Astrophys. J., 925, 34 (2022)>
*   =2022ApJ...925...34C
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'morphcat.dat'	! Morphological catalog of galaxies, derived from
                                the version 1.3 mosaic --updated by the author

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

C  J2000 position composed of: RAdeg DEdeg
      integer*4     v_MMS2013_ (nr__) ! [1/262615] Object identifier from UVISTA
*                                 catalog of Muzzin+ 2013, J/ApJS/206/8 (ID)
      real*8        RAdeg      (nr__) ! (deg) [149.3/151] Right ascension (J2000)
      real*8        DEdeg      (nr__) ! (deg) [1.6/3] Declination (J2000)
      integer*4     Flag       (nr__) ! [0/4] GALFIT flag (0=Good) (1)
      integer*4     Use        (nr__) ! [0/1] General use flag (2)
      real*8        mag        (nr__) ! (mag) [12.5/29.2]?=-99 GALFIT best-fit magnitude
      real*8        e_mag      (nr__) ! (mag) [0.0003/1.8]?=-99 Uncertainty in GALFIT
*                                 magnitude (see Section 3.5) (dmag)
      real*8        re         (nr__) ! (arcsec) [0.03/40]?=-99 GALFIT best-fit effective
*                                  (half-light) radius
      real*8        e_re       (nr__) ! (arcsec) [0.0002/7]?=-99 Uncertainty in re (dre)
      real*8        n          (nr__) ! [0.2/6]?=-99 GALFIT best-fit Sersic index
      real*8        e_n        (nr__) ! [0.001/12]?=-99 Uncertainty in n (dn)
      real*8        q          (nr__) ! [0/1]?=-99 GALFIT best-fit axis ratio
      real*8        e_q        (nr__) ! [0.0002/1.2]?=-99 Uncertainty in q (dq)
      real*8        PA         (nr__) ! (deg) [-90/90]?=-99 GALFIT best-fit
*                                  position angle
      real*8        e_PA       (nr__) ! (deg) [0.01/133]?=-99 Uncertainty in PA (dpa)
      real*8        KRad       (nr__) ! (arcsec) [0/22.6]?=-99 Kron radius from SExtractor (kron)
      real*8        F160W      (nr__) ! [0.04/96430]?=-99 F160W AUTO flux (SExtractor
*                                  measured, within Kron radius); ZP=25
*                                  (f_F160W_auto)
      real*8        e_F160W    (nr__) ! [0.03/128]?=-99 Error in F160W AUTO flux; ZP=25
*                                  (e_F160W_auto)
      real*8        F160Wt     (nr__) ! [0.2/96430]?=-99 F160W total flux (scaled from
*                                  AUTO flux); ZP=25 (f_F160W_tot)
      real*8        e_F160Wt   (nr__) ! [0.04/128]?=-99 Error in F160Wt (e_F160W_tot)
      real*8        SNR        (nr__) ! [0.39/8856]?=-99 Total scaled signal to noise
*                                  ratio (F160Wt/e_F160Wt)
      integer*4     f_re       (nr__) ! [0/1] Parameter robustness flag for re
*                                  (1=robust; 49,440 occurrences) (see Section 4)
*                                   (flag_limit_r)
      integer*4     f_n        (nr__) ! [0/1] Parameter robustness flag for n
*                                  (1=robust; 23,587 occurrences) (see Section 4)
*                                  (flag_limit_n)
      integer*4     Deb        (nr__) ! [0/4] Deblending flag (0=not blended)
*                                  (flag_deb) (3)
      real*8        Mcorr      (nr__) ! [0.3/1] Mass correction for deblended galaxies
*                                  (see Section 2.2)
      real*8        Depth      (nr__) ! (mag) [-0.23/27]?=-99 Total 5-sigma depth measured
*                                  from COSMOS-DASH mosaic ( 5_sigma_depth) (4)
      real*8        Chi2       (nr__) ! [-99/5.6e+06] GALFIT chi-squared (chi)
      real*8        Chi_nu     (nr__) ! [-99/141] GALFIT chi-squared per degree of
*                                  freedom (chi_nu)
*Note (1): GALFIT flags (see Section 3.1) as follows:
*   0 = good (14757 occurrences)
*   1 = suspicious (30008 occurrences)
*   2 = bad (6567 occurrences)
*   3 = failed (254 occurrences)
*   4 = no coverage (211029 occurrences)
*Note (2): General use flag (see Section 4);
*   1 = GALFIT flag<2, re>FWHM (38026 occurrences)
*Note (3): Deblending flag (see Section 2.2) as follows:
*   0 = not blended (49726 occurrences)
*   1 = deblended (1689 occurrences)
*   2 = blended (171 occurrences)
*   4 = no coverage (211029 occurrences)
*Note (4): Total 5-sigma depth, measured with empty aperture method assuming
*    0.7" diameter (see Section 3.6).

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

C  Loading file 'morphcat.dat'	! Morphological catalog of galaxies, derived from
*                                the version 1.3 mosaic --updated by the author

C  Format for file interpretation

    1 format(
     +  I6,1X,F9.5,1X,F9.7,1X,I1,1X,I1,1X,F10.6,1X,F10.6,1X,F10.6,1X,
     +  F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,
     +  1X,F10.6,1X,F12.6,1X,F10.6,1X,F12.6,1X,F10.6,1X,F11.6,1X,I1,
     +  1X,I1,1X,I1,1X,F8.6,1X,F10.6,1X,F14.6,1X,F10.6)

C  Effective file loading

      open(unit=1,status='old',file=
     +'morphcat.dat')
      write(6,*) '....Loading file: morphcat.dat'
      do i__=1,262615
        read(1,'(A263)')ar__
        read(ar__,1)
     +  v_MMS2013_(i__),RAdeg(i__),DEdeg(i__),Flag(i__),Use(i__),
     +  mag(i__),e_mag(i__),re(i__),e_re(i__),n(i__),e_n(i__),q(i__),
     +  e_q(i__),PA(i__),e_PA(i__),KRad(i__),F160W(i__),e_F160W(i__),
     +  F160Wt(i__),e_F160Wt(i__),SNR(i__),f_re(i__),f_n(i__),
     +  Deb(i__),Mcorr(i__),Depth(i__),Chi2(i__),Chi_nu(i__)
c    ..............Just test output...........
        write(6,1)
     +  v_MMS2013_(i__),RAdeg(i__),DEdeg(i__),Flag(i__),Use(i__),
     +  mag(i__),e_mag(i__),re(i__),e_re(i__),n(i__),e_n(i__),q(i__),
     +  e_q(i__),PA(i__),e_PA(i__),KRad(i__),F160W(i__),e_F160W(i__),
     +  F160Wt(i__),e_F160Wt(i__),SNR(i__),f_re(i__),f_n(i__),
     +  Deb(i__),Mcorr(i__),Depth(i__),Chi2(i__),Chi_nu(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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