FORTRAN Generation
(/./ftp/cats/J/ApJS/270/12)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJS/270/12 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-Oct-12
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/270/12  The UNCOVER SPS cat. for robust sources up to z~15  (Wang+, 2024)
*================================================================================
*The UNCOVER survey: a first-look HST+JWST catalog of galaxy redshifts and
*stellar population properties spanning 0.2<~z<~15.
*    Wang B., Leja J., Labbe I., Bezanson R., Whitaker K.E., Brammer G.,
*    Furtak L.J., Weaver J.R., Price S.H., Zitrin A., Atek H., Coe D.,
*    Cutler S.E., Dayal P., van Dokkum P., Feldmann R., Marchesini D.,
*    Franx M., Forster Schreiber N., Fujimoto S., Geha M., Glazebrook K.,
*    de Graaff A., Greene J.E., Juneau S., Kassin S., Kriek M., Khullar G.,
*    Maseda M., Mowla L.A., Muzzin A., Nanayakkara T., Nelson E.J., Oesch P.A.,
*    Pacifici C., Pan R., Papovich C., Setton D.J., Shapley A.E., Smit R.,
*    Stefanon M., Suess K.A., Taylor E.N., Williams C.C.
*   <Astrophys. J. Suppl. Ser., 270, 12 (2024)>
*   =2024ApJS..270...12W
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'spscat.dat'	! Stellar population catalog based on the DR2 SUPER
                               photometric catalog

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

C  J2000 position composed of: RAdeg DEdeg
      integer*4     Seq        (nr__) ! [1/61648] Running sequence number as in
*                                      UNCOVER DR2 phot. cat.
*                                      (Weaver+ 2024, J/ApJS/270/7) (id)
      real*8        RAdeg      (nr__) ! (deg) [3.4/3.7] Right ascension (J2000) (ra)
      real*8        DEdeg      (nr__) ! (deg) [-30.5/-30.2] Declination (J2000) (dec)
      integer*4     UsePhot    (nr__) ! 1 if photometry is reliable (use_phot)
      real*8        zsp        (nr__) ! [2e-6/6.4]? Spectroscopic redshift (z_spec)
      real*8        z16        (nr__) ! [0.0017/18.8] Redshift posterior, 16%
*                                      (z_16)
      real*8        z50        (nr__) ! [0.0018/19.6] Redshift posterior, 50%
*                                      (z_50)
      real*8        z84        (nr__) ! [0.002/20] Redshift posterior, 84% (z_84)
      real*8        logMs16    (nr__) ! ([Msun]) [5.7/11.4] Log of stellar mass, 16%
*                                      (mstar_16)
      real*8        logMs50    (nr__) ! ([Msun]) [5.7/11.4] Log of stellar mass, 50%
*                                      (mstar_50)
      real*8        logMs84    (nr__) ! ([Msun]) Log of stellar mass, 84% (mstar_84)
      real*8        logZ16     (nr__) ! ([Sun]) [-2/0.2] Stellar metallicity, 16% (met_16)
      real*8        logZ50     (nr__) ! ([Sun]) [-2/0.2] Stellar metallicity, 50% (met_50)
      real*8        logZ84     (nr__) ! ([Sun]) Stellar metallicity, 84% (met_84)
      real*8        wAge16     (nr__) ! (Gyr) [0.01/12.6] Mass-weighted age, 16% (mwa_16)
      real*8        wAge50     (nr__) ! (Gyr) [0.017/13] Mass-weighted age, 50% (mwa_50)
      real*8        wAge84     (nr__) ! (Gyr) Mass-weighted age, 84% (mwa_84)
      real*8        dustd16    (nr__) ! [7.5e-5/4] Optical depth of diffuse dust,
*                                      16% (dust2_16)
      real*8        dustd50    (nr__) ! [0.00029/4] Optical depth of diffuse dust,
*                                      50% (dust2_50)
      real*8        dustd84    (nr__) ! Optical depth of diffuse dust, 84%
*                                      (dust2_84)
      real*8        dusti16    (nr__) ! [-1/0.4] Power law index for a
*                                      Calzetti+2000 attenuation curve, 16%
*                                      (dust_index_16)
      real*8        dusti50    (nr__) ! [-1/0.4] Power law index for a
*                                      Calzetti+ 2000ApJ...533..682C
*                                      attenuation curve, 50% (dust_index_50)
      real*8        dusti84    (nr__) ! Power law index for a Calzetti+2000
*                                      attenuation curve, 84% (dust_index_84)
      real*4        dustf16    (nr__) ! [0.001/1.8] Ratio between the optical depth
*                                      of birth cloud dust and diffuse dust, 16%
*                                      (dust1_fraction_16)
      real*4        dustf50    (nr__) ! [0.006/2] Ratio between the optical depth
*                                      of birth cloud dust and diffuse dust, 50%
*                                      (dust1_fraction_50)
      real*4        dustf84    (nr__) ! Ratio between the optical depth of birth
*                                      cloud dust and diffuse dust, 84%
*                                      (dust1_fraction_84)
      real*8        logfagn16  (nr__) ! [-5/0.5] Logarithm of the ratio between
*                                      the object's AGN luminosity and its
*                                      bolometric luminosity, 16% (logfagn_16)
      real*8        logfagn50  (nr__) ! [-4.8/0.5] Log of the ratio between the
*                                      object's AGN luminosity and its
*                                      bolometric luminosity, 50% (logfagn_50)
      real*8        logfagn84  (nr__) ! Logarithm of the ratio between the
*                                      object's AGN luminosity and its
*                                      bolometric luminosity, 84% (logfagn_84)
      real*4        sfr10p16   (nr__) ! (Msun/yr) [1.4e-9/1972] Star formation rate
*                                      averaged over 10Myr, 16% (sfr10_16)
      real*4        sfr10p50   (nr__) ! (Msun/yr) [4.9e-8/2776] Star formation rate averaged
*                                      over 10Myr, 50% (sfr10_50)
      real*4        sfr10p84   (nr__) ! (Msun/yr) [6.3e-7/4020] Star formation rate
*                                      averaged over 10Myr, 84% (sfr10_84)
      real*4        sfr30p16   (nr__) ! (Msun/yr) Star formation rate averaged over 30Myr,
*                                      16% (sfr30_16)
      real*4        sfr30p50   (nr__) ! (Msun/yr) [7e-8/2047] Star formation rate averaged
*                                      over 30Myr, 50% (sfr30_50)
      real*4        sfr30p84   (nr__) ! (Msun/yr) Star formation rate averaged over 30Myr,
*                                      84% (sfr30_84)
      real*4        sfr100p16  (nr__) ! (Msun/yr) [6.5e-9/694] Star formation rate averaged
*                                      over 100Myr, 16% (sfr100_16)
      real*4        sfr100p50  (nr__) ! (Msun/yr) [1e-7/793] Star formation rate averaged
*                                      over 100Myr, 50% (sfr100_50)
      real*4        sfr100p84  (nr__) ! (Msun/yr) Star formation rate averaged over 100Myr,
*                                      84% (sfr100_84)
      real*4        ssfr10p16  (nr__) ! (yr-1) [8e-17/3e-8] Specific star formation rate
*                                      averaged over 10Myr, 16% (ssfr10_16)
      real*4        ssfr10p50  (nr__) ! (yr-1) [3.6e-15/3.5e-8] Specific star formation
*                                      rate averaged over 10Myr, 50% (ssfr10_50)
      real*4        ssfr10p84  (nr__) ! (yr-1) Specific star formation rate averaged
*                                      over 10Myr, 84% (ssfr10_84)
      real*4        ssfr30p16  (nr__) ! (yr-1) Specific star formation rate averaged
*                                      over 30Myr, 16% (ssfr30_16)
      real*4        ssfr30p50  (nr__) ! (yr-1) [4e-15/3e-8] Specific star formation rate
*                                      averaged over 30Myr, 50% (ssfr30_50)
      real*4        ssfr30p84  (nr__) ! (yr-1) Specific star formation rate averaged
*                                      over 30Myr, 84% (ssfr30_84)
      real*4        ssfr100p16 (nr__) ! (yr-1) Specific star formation rate averaged
*                                      over 100Myr, 16% (ssfr100_16)
      real*4        ssfr100p50 (nr__) ! (yr-1) [6.5e-15/9.3e-9] Specific star formation
*                                     rate averaged over 100Myr, 50% (ssfr100_50)
      real*4        ssfr100p84 (nr__) ! (yr-1) [9e-14/9.5e-9] Specific star formation
*                                      rate averaged over 100Myr, 84%
*                                      (ssfr100_84)
      real*8        muz16      (nr__) ! [0.1/3546] Magnification based z16,
*                                      calculated from the public maps
*                                      (mu_num_16)
      real*8        muz50      (nr__) ! [0.06/3689] Magnification based z50,
*                                     calculated from the public maps (mu_num_50)
      real*8        muz84      (nr__) ! Magnification based z84, calculated from
*                                      the public maps (mu_num_84)
      real*8        restU16    (nr__) ! (mag) [17.6/38] Rest-frame U-band AB mag, 16%
*                                      (rest_U_16)
      real*8        restU50    (nr__) ! (mag) [17.6/39] Rest-frame U-band AB mag, 50%
*                                      (rest_U_50)
      real*8        restU84    (nr__) ! (mag) [17.7/40.1] Rest-frame U-band AB mag, 84%
*                                      (rest_U_84)
      real*8        restV16    (nr__) ! (mag) [17.4/36.4] Rest-frame V-band AB mag, 16%
*                                      (rest_V_16)
      real*8        restV50    (nr__) ! (mag) [17.4/37.1] Rest-frame V-band AB mag, 50%
*                                      (rest_V_50)
      real*8        restV84    (nr__) ! (mag) Rest-frame V-band AB mag, 84% (rest_V_84)
      real*8        restJ16    (nr__) ! (mag) [17/35.1] Rest-frame J-band AB mag, 16%
*                                      (rest_J_16)
      real*8        restJ50    (nr__) ! (mag) [17.2/35.7] Rest-frame J-band AB mag, 50%
*                                      (rest_J_50)
      real*8        restJ84    (nr__) ! (mag) Rest-frame J-band AB mag, 84% (rest_J_84)
      real*8        restu16_1  (nr__) ! (mag) [17.6/41] Rest-frame u-band AB mag, 16%
*                                      (rest_up_16)
      real*8        restu50_1  (nr__) ! (mag) [17.7/42] Rest-frame u-band AB mag, 50%
*                                      (rest_up_50)
      real*8        restu84_1  (nr__) ! (mag) [17.7/43] Rest-frame u-band AB mag, 84%
*                                      (rest_up_84)
      real*8        restg16    (nr__) ! (mag) [17.2/37] Rest-frame g-band AB mag, 16%
*                                      (rest_g_16)
      real*8        restg50    (nr__) ! (mag) [17.2/37.8] Rest-frame g-band AB mag, 50%
*                                      (rest_g_50)
      real*8        restg84    (nr__) ! (mag) Rest-frame g-band AB mag, 84% (rest_g_84)
      real*8        resti16    (nr__) ! (mag) [17.6/35.8] Rest-frame i-band AB mag, 16%
*                                      (rest_i_16)
      real*8        resti50    (nr__) ! (mag) [17.6/36.4] Rest-frame i-band AB mag, 50%
*                                      (rest_i_50)
      real*8        resti84    (nr__) ! (mag) Rest-frame i-band AB mag, 84% (rest_i_84)
      real*8        U_V16      (nr__) ! (mag) [-0.8/6.3] Rest-frame U-V color, 16%
*                                      (UV_16)
      real*8        U_V50      (nr__) ! (mag) [-0.7/6.4] Rest-frame U-V color, 50%
*                                      (UV_50)
      real*8        U_V84      (nr__) ! (mag) Rest-frame U-V color, 84% (UV_84)
      real*8        V_J16      (nr__) ! (mag) [-1.4/4] Rest-frame V-J color, 16% (VJ_16)
      real*8        V_J50      (nr__) ! (mag) [-1.4/4.1] Rest-frame V-J color, 50%
*                                      (VJ_50)
      real*8        V_J84      (nr__) ! (mag) Rest-frame V-J color, 84% (VJ_84)
      real*8        g_i16      (nr__) ! (mag) [-0.7/5.4] Rest-frame synthetic g-i, 16%
*                                      (gi_16)
      real*8        g_i50      (nr__) ! (mag) [-0.7/5.5] Rest-frame synthetic g-i, 50%
*                                      (gi_50)
      real*8        g_i84      (nr__) ! (mag) Rest-frame synthetic g-i, 84% (gi_84)
      real*8        u_g16      (nr__) ! (mag) [-0.3/8] Rest-frame synthetic u-g, 16%
*                                      (ug_16)
      real*8        u_g50      (nr__) ! (mag) [-0.3/8.1] Rest-frame synthetic u-g, 50%
*                                      (ug_50)
      real*8        u_g84      (nr__) ! (mag) Rest-frame synthetic u-g, 84% (ug_84)
      real*4        chi2       (nr__) ! [6e-7/4502] Best-fit chi^2^, assuming a
*                                      minimum error of 5% (chi2)
      integer*4     Nb         (nr__) ! [2/15] Number of bands used in the fit
*                                      (nbands)
      real*8        zEazy      (nr__) ! [0.01/20]? Peak of p(z) from fitting with
*                                      eazy-sfhz-corr (z_eazy)
      real*8        zEazy16    (nr__) ! ? P(z), 16% (z_eazy_16)
      real*8        zEazy50    (nr__) ! [0.01/19.8]? P(z), 50% (z_eazy_50)
      real*8        zEazy84    (nr__) ! ? P(z), 84% (z_eazy_84)
      integer*4     IDDR1      (nr__) ! [1/50365]? ID of the source in the DR1
*                                      cat. (id_DR1)
      integer*4     MSA        (nr__) ! [2008/60157]? ID of the source in the
*                                      team internal MSA catalog (Price+
*                                      2023arXiv231002500P); all within <0.24"
*                                      radius (id_msa)
      integer*4     ALMA       (nr__) ! [1/69]? ID of ALMA source from
*                                      Fujimoto+ 2023, in prep. (id_alma)
      real*4        UseAper    (nr__) ! (arcsec) [0.3/1.4] Diameter adopted for color
*                                      aperture (use_aper)
      integer*4     f_kron     (nr__) ! [0/1]? 1 for systematically
*                                      underestimated photometry (i.e., the
*                                       kron radius is not robust,thus the total
*                                       photometry assumes point-lke light
*                                       distribution) (flag_kron)

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

C  Loading file 'spscat.dat'	! Stellar population catalog based on the DR2 SUPER
*                               photometric catalog

C  Format for file interpretation

    1 format(
     +  I5,1X,F13.11,1X,F15.11,1X,I1,1X,F8.6,1X,F9.6,1X,F9.6,1X,F9.6,
     +  1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,
     +  F7.4,1X,F7.4,1X,F8.6,1X,F8.6,1X,F8.6,1X,F7.4,1X,F7.4,1X,F7.4,
     +  1X,F6.4,1X,F6.4,1X,F6.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,E11.6,1X,
     +  E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,
     +  1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,
     +  E11.6,1X,E11.6,1X,E11.6,1X,E11.6,1X,F9.4,1X,F9.4,1X,F9.4,1X,
     +  F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,
     +  1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,
     +  F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,
     +  1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,
     +  E10.5,1X,I2,1X,F9.6,1X,F9.6,1X,F9.6,1X,F9.6,1X,I5,1X,I5,1X,I2,
     +  1X,F4.2,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'spscat.dat')
      write(6,*) '....Loading file: spscat.dat'
      do i__=1,61648
        read(1,'(A806)')ar__
        read(ar__,1)
     +  Seq(i__),RAdeg(i__),DEdeg(i__),UsePhot(i__),zsp(i__),z16(i__),
     +  z50(i__),z84(i__),logMs16(i__),logMs50(i__),logMs84(i__),
     +  logZ16(i__),logZ50(i__),logZ84(i__),wAge16(i__),wAge50(i__),
     +  wAge84(i__),dustd16(i__),dustd50(i__),dustd84(i__),
     +  dusti16(i__),dusti50(i__),dusti84(i__),dustf16(i__),
     +  dustf50(i__),dustf84(i__),logfagn16(i__),logfagn50(i__),
     +  logfagn84(i__),sfr10p16(i__),sfr10p50(i__),sfr10p84(i__),
     +  sfr30p16(i__),sfr30p50(i__),sfr30p84(i__),sfr100p16(i__),
     +  sfr100p50(i__),sfr100p84(i__),ssfr10p16(i__),ssfr10p50(i__),
     +  ssfr10p84(i__),ssfr30p16(i__),ssfr30p50(i__),ssfr30p84(i__),
     +  ssfr100p16(i__),ssfr100p50(i__),ssfr100p84(i__),muz16(i__),
     +  muz50(i__),muz84(i__),restU16(i__),restU50(i__),restU84(i__),
     +  restV16(i__),restV50(i__),restV84(i__),restJ16(i__),
     +  restJ50(i__),restJ84(i__),restu16_1(i__),restu50_1(i__),
     +  restu84_1(i__),restg16(i__),restg50(i__),restg84(i__),
     +  resti16(i__),resti50(i__),resti84(i__),U_V16(i__),U_V50(i__),
     +  U_V84(i__),V_J16(i__),V_J50(i__),V_J84(i__),g_i16(i__),
     +  g_i50(i__),g_i84(i__),u_g16(i__),u_g50(i__),u_g84(i__),
     +  chi2(i__),Nb(i__),zEazy(i__),zEazy16(i__),zEazy50(i__),
     +  zEazy84(i__),IDDR1(i__),MSA(i__),ALMA(i__),UseAper(i__),
     +  f_kron(i__)
        if(ar__(39:46) .EQ. '') zsp(i__) = rNULL__
        if(ar__(746:754) .EQ. '') zEazy(i__) = rNULL__
        if(ar__(756:764) .EQ. '') zEazy16(i__) = rNULL__
        if(ar__(766:774) .EQ. '') zEazy50(i__) = rNULL__
        if(ar__(776:784) .EQ. '') zEazy84(i__) = rNULL__
        if(ar__(786:790) .EQ. '') IDDR1(i__) = iNULL__
        if(ar__(792:796) .EQ. '') MSA(i__) = iNULL__
        if(ar__(798:799) .EQ. '') ALMA(i__) = iNULL__
        if(ar__(806:806) .EQ. '') f_kron(i__) = iNULL__
c    ..............Just test output...........
        write(6,1)
     +  Seq(i__),RAdeg(i__),DEdeg(i__),UsePhot(i__),zsp(i__),z16(i__),
     +  z50(i__),z84(i__),logMs16(i__),logMs50(i__),logMs84(i__),
     +  logZ16(i__),logZ50(i__),logZ84(i__),wAge16(i__),wAge50(i__),
     +  wAge84(i__),dustd16(i__),dustd50(i__),dustd84(i__),
     +  dusti16(i__),dusti50(i__),dusti84(i__),dustf16(i__),
     +  dustf50(i__),dustf84(i__),logfagn16(i__),logfagn50(i__),
     +  logfagn84(i__),sfr10p16(i__),sfr10p50(i__),sfr10p84(i__),
     +  sfr30p16(i__),sfr30p50(i__),sfr30p84(i__),sfr100p16(i__),
     +  sfr100p50(i__),sfr100p84(i__),ssfr10p16(i__),ssfr10p50(i__),
     +  ssfr10p84(i__),ssfr30p16(i__),ssfr30p50(i__),ssfr30p84(i__),
     +  ssfr100p16(i__),ssfr100p50(i__),ssfr100p84(i__),muz16(i__),
     +  muz50(i__),muz84(i__),restU16(i__),restU50(i__),restU84(i__),
     +  restV16(i__),restV50(i__),restV84(i__),restJ16(i__),
     +  restJ50(i__),restJ84(i__),restu16_1(i__),restu50_1(i__),
     +  restu84_1(i__),restg16(i__),restg50(i__),restg84(i__),
     +  resti16(i__),resti50(i__),resti84(i__),U_V16(i__),U_V50(i__),
     +  U_V84(i__),V_J16(i__),V_J50(i__),V_J84(i__),g_i16(i__),
     +  g_i50(i__),g_i84(i__),u_g16(i__),u_g50(i__),u_g84(i__),
     +  chi2(i__),Nb(i__),zEazy(i__),zEazy16(i__),zEazy50(i__),
     +  zEazy84(i__),IDDR1(i__),MSA(i__),ALMA(i__),UseAper(i__),
     +  f_kron(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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