Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/702/A187 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-Jul-13
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/A+A/702/A187 GLACE survey: The mass--SFR relation (De Daniloff+, 2025)
*================================================================================
*GLACE survey: OSIRIS/GTC multi-object spectroscopy of the rich galaxy cluster
*ZwCl 0024.0+1652 at z~0.4.
*III. The mass-SFR relation and the quenching of cluster galaxies.
* De Daniloff S.B., Bongiovanni A., Sanchez-Portal M., Cedres B.,
* Padilla-Torres C.P., Perez-Garcia A.M., Perez-Martinez R., Espada D.,
* de la Casa C.C., Torres-Rios G., Gonzalez-Otero M., de Diego J.A.,
* Rodriguez M.I., Cervino M., Lara-Lopez M.A., Cepa J., Valtchanov I.,
* Gonzalez-Serrano J.I., Cruz-Gonzalez I., Negrete C.A., Beyoro-Amado Z.,
* Castillo-Fraile M., Mancillas B., Povic M.
* <Astron. Astrophys. 702, A187 (2025)>
* =2025A&A...702A.187D (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'catalog.dat' ! Catalog of the properties of the galaxies
integer*4 nr__
parameter (nr__=102) ! Number of records
character*3288 ar__ ! Full-size record
C J2000 position composed of: RAdeg DEdeg
character*7 ID_unique (nr__) ! GLACE ID with pointing label,
* NNNNN_a (id_unique) (1)(2)
integer*4 ID (nr__) ! GLACE ID without pointing label
* (ID) (1)
character*5 ID_MOS (nr__) ! ID of the MOS spectrum after
* reduction (id_MOS)
real*8 RAdeg (nr__) ! (deg) Right ascension (J2000) from TF
* (tunable filters) calibrated data
* (RA_J2000) (1)
real*8 DEdeg (nr__) ! (deg) Declination (J2000) from TF
* calibrated data (DEC_J2000) (1)
real*8 RAGdeg (nr__) ! (deg) Right ascension from GAIA-DR1
* (ICRS) at Ep=2015.0 (RA_GAIA) (1)
real*8 DEGdeg (nr__) ! (deg) Declination from GAIA- DR1 (ICRS)
* at Ep=2015.0 (DEC_GAIA) (1)
real*8 zSP15 (nr__) ! Redshift from TF calibrated data
* (Sanchez-Portal et al.,
* 2015A&A...578A..30S,) (z_SP15)
real*8 zHaIC (nr__) ! Redshift of Halpha emission line
* obtained from inverse convolution
* method (z_Halpha_InvConv) (3)
real*8 b_zHaIC (nr__) ! Lower value of redshift of Halpha
* emission line
* (z_Halpha_InvConv_d) (3)
real*8 B_zHaIC_1 (nr__) ! Upper value of redshift of Halpha
* emission line
* (z_Halpha_InvConv_u) (3)
real*8 zOII (nr__) ! ? Redshift of [OII] emission
* line (z_OII)
real*8 e_zOII (nr__) ! ? Error in redshift of [OII]
* emission line (e_z_OII)
real*8 zHb (nr__) ! ? Redshift of Hbeta emission line
* (z_Hbeta)
real*8 e_zHb (nr__) ! ? Error in redshift of Hbeta
* emission line (e_z_Hbeta)
real*8 zOIII (nr__) ! ? Redshift of [OIII]5007 emission
* line (z_OIII)
real*8 e_zOIII (nr__) ! ? Error in redshift of [OIII]5007
* emission line (e_z_OIII)
real*8 zMOSavg (nr__) ! ? Average redshift of [OII],
* Hbeta, and [OIII]5007 emission
* lines (z_MOS_avg)
real*8 e_zMOSavg (nr__) ! ? Error in average redshift of
* [OII], Hbeta, and [OIII]5007
* emission lines (e_z_MOS_avg)
character*10 Structure (nr__) ! Kinematical structure in which
* the galaxy is located
* (Structure) (4)
real*8 Bmag (nr__) ! (mag) ? B-band (CFHT/CFH12k) AUTO
* magnitude (Vega system);
* automatic aperture from
* SExtractor; proxy for total
* magnitude (B_AUTO) (1)
real*8 e_Bmag (nr__) ! (mag) ? Error in B-band AUTO magnitude
* (e_B_AUTO) (1)
real*8 Vmag (nr__) ! (mag) ? V-band (CFHT/CFH12k) AUTO
* magnitude V_AUTO) (1)
real*8 e_Vmag (nr__) ! (mag) ? Error in V-band AUTO magnitude
* (e_V_AUTO) (1)
real*8 Rmag (nr__) ! (mag) ? R-band (CFHT/CFH12k) AUTO
* magnitude (R_AUTO) (1)
real*8 e_Rmag (nr__) ! (mag) ? Error in R-band AUTO magnitude
* (e_R_AUTO) (1)
real*8 Imag (nr__) ! (mag) ? I-band (CFHT/CFH12k) AUTO
* magnitude (I_AUTO) (1)
real*8 e_Imag (nr__) ! (mag) ? Error in I-band AUTO magnitude
* (e_I_AUTO) (1)
real*8 Jmag (nr__) ! (mag) ? J-band (Palomar200/WIRC) AUTO
* magnitude (J_AUTO) (1)
real*8 e_Jmag (nr__) ! (mag) ? Error in J-band AUTO magnitude
* (e_J_AUTO) (1)
real*8 Ksmag (nr__) ! (mag) ? Ks-band (Palomar200/WIRC) AUTO
* magnitude (Ks_AUTO) (1)
real*8 e_Ksmag (nr__) ! (mag) ? Error in Ks-band AUTO magnitude
* (e_Ks_AUTO) (1)
real*8 F814Wmag (nr__) ! (mag) ? HST/WFPC2 814nm wide band AUTO
* magnitude (F814W_AUTO) (1)
real*8 e_F814Wmag (nr__) ! (mag) ? Error in HST/WFPC2 814nm wide
* band AUTO magnitude
* (e_F814W_AUTO) (1)
real*4 FHaSP15 (nr__) ! (mW/m2) Flux of Halpha emission line
* (F_Halpha_SP15) (1)
real*4 e_FHaSP15 (nr__) ! (mW/m2) Error in flux of Halpha emission
* line (e_F_Halpha_SP15) (1)
real*4 FNIISP15 (nr__) ! (mW/m2) ? Flux of [NII]6548 emission line
* (F_NII_SP15) (1)
real*4 e_FNIISP15 (nr__) ! (mW/m2) ? Error in flux of [NII]6548
* emission line (e_FNII_SP15) (1)
real*4 FHacSP15 (nr__) ! (mW/m2) Flux of Halpha adjusted and
* corrected for absorptions and
* extinction
* (FHalpha_corrected_SP15) (1)
real*4 e_FHacSP15 (nr__) ! (mW/m2) Error in flux of Halpha adjusted
* and corrected for absorptions and
* extinction
* (eFHalpha_corrected_SP15) (1)
real*4 FNIIcSP15 (nr__) ! (mW/m2) ? Flux of [NII]6548 adjusted and
* corrected for absorptions and
* extinction
* (FNII_corrected_SP15) (1)
real*4 e_FNIIcSP15(nr__) ! (mW/m2) ? Error in flux of [NII]6548
* adjusted and corrected for
* absorptions and extinction
* (eFNII_corrected_SP15) (1)
real*8 EWHaSP15 (nr__) ! (0.1nm) Equivalent width of Halpha
* emission line obtained from
* tunable filters
* (EW_Halpha_SP15) (1)
real*8 e_EWHaSP15 (nr__) ! (0.1nm) Error in equivalent width of
* Halpha emission line obtained
* from tunable filters
* (e_EW_Halpha_SP15) (1)
real*4 ContHaIC (nr__) ! (10mW/m2/nm) Continuum level of Halpha line
* obtained from inverse convolution
* (Cont_Halpha_InvConv) (3)
real*4 b_ContHaIC (nr__) ! (10mW/m2/nm) Lower value of continuum level of
* Halpha line obtained from
* inverse convolution
* (Cont_Halpha_InvConv_d) (3)
real*4 B_ContHaIC_1(nr__) ! (10mW/m2/nm) Upper value of continuum level of
* Halpha line obtained from
* inverse convolution
* (Cont_Halpha_InvConv_u) (3)
real*8 sigmaHaIC (nr__) ! (0.1nm) Sigma of the Gaussian fit for
* Halpha line obtained from
* inverse convolution
* (sigma_Halpha_InvConv) (3)
real*8 b_sigmaHaIC(nr__) ! (0.1nm) Lower value of sigma of the
* Gaussian fit for Halpha line
* obtained from inverse convolution
* (sigma_Halpha_InvConv_d) (3)
real*8 B_sigmaHaIC_1(nr__) ! (0.1nm) Upper value of sigma of the
* Gaussian fit for Halpha line
* obtained from inverse convolution
* (sigma_Halpha_InvConv_u) (3)
real*4 FNII6548IC (nr__) ! (mW/m2) Flux of [NII]6548 emission line
* obtained from inverse convolution
* (F_NII6548_InvConv) (3)
real*4 b_FNII6548IC(nr__) ! (mW/m2) Lower value of flux of [NII]6548
* emission line obtained from
* inverse convolution
* (F_NII6548_InvConv_d) (3)
real*4 B_FNII6548IC_1(nr__) ! (mW/m2) Upper value of flux of [NII]6548
* emission line obtained from
* inverse convolution
* (F_NII6548_InvConv_u) (3)
real*4 FHaIC (nr__) ! (mW/m2) Flux of Halpha emission line
* obtained from inverse convolution
* (F_Halpha_InvConv) (3)
real*4 b_FHaIC (nr__) ! (mW/m2) Lower value of flux of Halpha
* obtained from inverse convolution
* (F_Halpha_InvConv_d) (3)
real*4 B_FHaIC_1 (nr__) ! (mW/m2) Upper value of flux of Halpha
* obtained from inverse convolution
* (F_Halpha_InvConv_u) (3)
real*4 FNII6583IC (nr__) ! (mW/m2) Flux of [NII]6583 emission line
* obtained from inverse convolution
* (F_NII6583_InvConv) (3)
real*4 b_FNII6583IC(nr__) ! (mW/m2) Lower value of flux of [NII]6583
* obtained from inverse convolution
* (F_NII6583_InvConv_d) (3)
real*4 B_FNII6583IC_1(nr__) ! (mW/m2) Upper value of flux of [NII]6583
* obtained from inverse convolution
* (F_NII6583_InvConv_u) (3)
real*4 EWNII6548IC(nr__) ! (0.1nm) Equivalent width of [NII]6548
* emission line obtained from
* inverse convolution
* (EW_NII6548_InvConv) (3)
real*4 b_EWNII6548IC(nr__) ! (0.1nm) Lower value of equivalent width of
* [NII]6548 emission line obtained
* from inverse convolution
* (EW_NII6548_InvConv_d) (3)
real*8 B_EWNII6548IC_1(nr__) ! (0.1nm) Upper value of equivalent width of
* [NII]6548 emission line
* obtained from inverse convolution
* (EW_NII6548_InvConv_u) (3)
real*8 EWHaIC (nr__) ! (0.1nm) Equivalent width of Halpha
* emission line obtained from
* inverse convolution
* (EW_Halpha_InvConv) (3)
real*4 b_EWHaIC (nr__) ! (0.1nm) Lower value of equivalent width of
* Halpha emission line obtained
* from inverse convolution
* (EW_Halpha_InvConv_d) (3)
real*8 B_EWHaIC_1 (nr__) ! (0.1nm) Upper value of equivalent width of
* Halpha emission line obtained
* from inverse convolution
* (EW_Halpha_InvConv_u) (3)
real*8 EWNII6583IC(nr__) ! (0.1nm) Equivalent width of [NII]6583
* emission line obtained from
* inverse convolution
* (EW_NII6583_InvConv) (3)
real*4 b_EWNII6583IC(nr__) ! (0.1nm) Lower value of equivalent width of
* [NII]6583 emission line obtained
* from inverse convolution
* (EW_NII6583_InvConv_d) (3)
real*8 B_EWNII6583IC_1(nr__) ! (0.1nm) Upper value of equivalent width of
* [NII]6583 emission line obtained
* from inverse convolution
* (EW_NII6583_InvConv_u) (3)
real*8 EWOII (nr__) ! (0.1nm) ? Equivalent width of [OII]
* emission line (EW_OII)
real*8 e_EWOII (nr__) ! (0.1nm) ? Error in equivalent width of
* [OII] emission line (e_EW_OII)
real*8 EWHb (nr__) ! (0.1nm) ? Equivalent width of Hbeta
* emission line (EW_Hbeta)
real*8 e_EWHb (nr__) ! (0.1nm) ? Error in equivalent width of
* Hbeta emission line (e_EW_Hbeta)
real*8 EWOIII (nr__) ! (0.1nm) ? Equivalent width of [OIII]5007
* emission line (EW_OIII)
real*8 e_EWOIII (nr__) ! (0.1nm) ? Error in equivalent width of
* [OIII]5007 emission line
* (e_EW_OIII)
real*8 contOII (nr__) ! (10mW/m2/nm) ? Continuum level of [OII]
* emission line (cont_OII)
real*8 e_contOII (nr__) ! (10mW/m2/nm) ? Error in continuum level of
* [OII] emission line (e_cont_OII)
real*8 contHb (nr__) ! (10mW/m2/nm) ? Continuum level of Hbeta
* emission line (cont_Hbeta)
real*8 e_contHb (nr__) ! (10mW/m2/nm) ? Error in continuum level of
* Hbeta emission line
* (e_cont_Hbeta)
real*8 contOIII (nr__) ! (10mW/m2/nm) ? Continuum level of [OIII]5007
* emission line (cont_OIII)
real*8 e_contOIII (nr__) ! (10mW/m2/nm) ? Error in continuum level of
* [OIII]5007 emission line
* (e_cont_OIII)
real*8 coreOII (nr__) ! (10mW/m2/nm) ? Core (amplitude) of [OII]
* emission line (core_OII)
real*8 e_coreOII (nr__) ! (10mW/m2/nm) ? Error in core (amplitude) of
* [OII] emission line (e_core_OII)
real*8 coreHb (nr__) ! (10mW/m2/nm) ? Core (amplitude) of Hbeta
* emission line (core_Hbeta)
real*8 e_coreHb (nr__) ! (10mW/m2/nm) ? Error in core (amplitude) of
* Hbeta emission line
* (e_core_Hbeta)
real*8 coreOIII (nr__) ! (10mW/m2/nm) ? Core (amplitude) of [OIII]5007
* emission line (core_OIII)
real*8 e_coreOIII (nr__) ! (10mW/m2/nm) ? Error in core (amplitude) of
* [OIII]5007 emission line
* (e_core_OIII)
real*8 sigmaOII (nr__) ! (0.1nm) ? sigma of the Gaussian fit for
* [OII] emission line (sigma_OII)
real*8 e_sigmaOII (nr__) ! (0.1nm) ? Error in sigma of the Gaussian
* fit for [OII] emission line
* (e_sigma_OII)
real*8 sigmaHb (nr__) ! (0.1nm) ? sigma of the Gaussian fit for
* Hbeta emission line (sigma_Hbeta)
real*8 e_sigmaHb (nr__) ! (0.1nm) ? Error in sigma of the Gaussian
* fit for Hbeta emission line
* (e_sigma_Hbeta)
real*8 sigmaOIII (nr__) ! (0.1nm) ? sigma of the Gaussian fit for
* [OIII]5007 emission line
* (sigma_OIII)
real*8 e_sigmaOIII(nr__) ! (0.1nm) ? Error in sigma of the Gaussian
* fit for [OIII]5007 emission line
* (e_sigma_OIII)
real*8 FWHMOII (nr__) ! (0.1nm) ? Full-Width at Half Maximum
* (FWHM) of the Gaussian fit for
* [OII] emission line (FWHM_OII)
real*8 e_FWHMOII (nr__) ! (0.1nm) ? Error in FWHM of the Gaussian
* fit for [OII] emission line
* (e_FWHM_OII)
real*8 FWHMHb (nr__) ! (0.1nm) ? Full-Width at Half Maximum
* (FWHM) of the Gaussian fit for
* Hbeta emission line (FWHM_Hbeta)
real*8 e_FWHMHb (nr__) ! (0.1nm) ? Error in FWHM of the Gaussian
* fit for Hbeta emission line
* (e_FWHM_Hbeta)
real*8 FWHMOIII (nr__) ! (0.1nm) ? Full-Width at Half Maximum
* (FWHM) of the Gaussian fit for
* [OIII]5007 emission line
* (FWHM_OIII)
real*8 e_FWHMOIII (nr__) ! (0.1nm) ? Error in FWHM of the Gaussian
* fit for [OIII]5007 emission line
* (e_FWHM_OIII)
real*8 SNRlineOII (nr__) ! ? Signal-to-Noise Ratio (S/N) of
* the [OII] emission line
* (SNR_line_OII)
real*8 SNRlineHb (nr__) ! ? Signal-to-Noise Ratio (S/N) of
* the Hbeta emission line
* (SNR_line_Hbeta)
real*8 SNRlineOIII(nr__) ! ? Signal-to-Noise Ratio (S/N) of
* the [OIII]5007 emission line
* (SNR_line_OIII)
real*4 FOII (nr__) ! (mW/m2) ? Flux of [OII] emission line
* (F_OII)
real*4 e_FOII (nr__) ! (mW/m2) ? Error in flux of [OII] emission
* line (e_F_OII)
real*4 FHb (nr__) ! (mW/m2) ? Flux of Hbeta emission line
* (F_Hbeta)
real*4 e_FHb (nr__) ! (mW/m2) ? Error in flux of Hbeta emission
* line (e_F_Hbeta)
real*4 FHblim (nr__) ! (mW/m2) ? Upper-limit flux of Hbeta
* (if no fitting was performed)
* (F_Hbeta_lim)
real*4 FOIII (nr__) ! (mW/m2) ? Flux of [OIII]5007 emission line
* (F_OIII)
real*4 e_FOIII (nr__) ! (mW/m2) ? Error in flux of [OIII]5007
* emission line (e_F_OIII)
real*4 FOIIIlim (nr__) ! (mW/m2) ? Upper-limit flux of [OIII]5007
* (if no fitting was performed)
* (F_OIII_lim)
real*4 FOIIadj (nr__) ! (mW/m2) ? Flux of [OII] emission line
* adjusted for aperture effect
* (F_OII_adj) (5)
real*4 e_FOIIadj (nr__) ! (mW/m2) ? Error in flux of [OII] emission
* line adjusted for aperture effect
* (e_F_OII_adj) (5)
real*4 FHbadj (nr__) ! (mW/m2) ? Flux of Hbeta emission line
* adjusted for aperture effect
* (F_Hbeta_adj) (5)
real*4 e_FHbadj (nr__) ! (mW/m2) ? Error in flux of Hbeta emission
* line adjusted for aperture effect
* (e_F_Hbeta_adj) (5)
real*4 FHbadjlim (nr__) ! (mW/m2) ? Upper-limit flux of Hbeta
* adjusted for aperture effect
* (if no fitting was performed)
* (F_Hbeta_lim_adj) (4)
real*4 FOIIIadj (nr__) ! (mW/m2) ? Flux of [OIII]5007 emission line
* adjusted for aperture effect
* (F_OIII_adj) (5)
real*4 e_FOIIIadj (nr__) ! (mW/m2) ? Error in flux of [OIII]5007
* emission line adjusted for
* aperture effect
* (e_F_OI )5)II_adj) (5)
real*4 FOIIIadjlim(nr__) ! (mW/m2) ? Upper-limit flux of [OIII]5007
* adjusted for aperture effect
* (if no fitting was performed)
* (F_OIII_lim_adj) (4)
real*4 FHaabs (nr__) ! (mW/m2) ? Flux of Halpha emission line
* adjusted and corrected for
* absorption (F_Halpha_abs) (3)
real*4 e_FHaabs (nr__) ! (mW/m2) ? Error in flux of Halpha emission
* line adjusted and corrected for
* absorption (e_F_Halpha_abs) (3)
real*4 FHbabs (nr__) ! (mW/m2) ? Flux of Hbeta emission line
* adjusted and corrected for
* absorption (F_Hbeta_abs)
real*4 e_FHbabs (nr__) ! (mW/m2) ? Error in flux of Hbeta emission
* line adjusted and corrected for
* absorption (e_F_Hbeta_abs)
real*4 FHbabslim (nr__) ! (mW/m2) ? Upper-limit flux of Hbeta
* adjusted and corrected for
* absorption (if no fitting was
* performed) (F_Hbeta_lim_abs)
real*4 FNII6548ec (nr__) ! (mW/m2) ? Flux of [NII]6548 emission line
* adjusted and corrected for
* absorptions and extinction
* (F_NII6548_extinc) (3)
real*4 e_FNII6548ec(nr__) ! (mW/m2) ? Error in flux of [NII]6548
* emission line adjusted and
* corrected for absorptions and
* extinction
* (e_F_NII6548_extinc) (3)
real*4 FHaec (nr__) ! (mW/m2) ? Flux of Halpha emission line
* adjusted and corrected for
* absorptions and extinction
* (F_Halpha_extinc)
real*4 e_FHaec (nr__) ! (mW/m2) ? Error in flux of Halpha emission
* line adjusted and corrected for
* absorptions and extinction
* (e_F_Halpha_extinc) (3)
real*4 FNII6583ec (nr__) ! (mW/m2) Flux of [NII]6583 emission line
* adjusted and corrected for
* absorptions and extinction
* (F_NII6583_extinc) (3)
real*4 e_FNII6583ec(nr__) ! (mW/m2) Error in flux of [NII]6583
* emission line adjusted and
* corrected for absorptions and
* extinction
* (e_F_NII6583_extinc) (3)
real*4 FOIIec (nr__) ! (mW/m2) ? Flux of [OII] emission line
* adjusted and corrected for
* absorptions and extinction
* (F_OII_extinc)
real*4 e_FOIIec (nr__) ! (mW/m2) ? Error in flux of [OII] emission
* line adjusted and corrected for
* absorptions and extinction
* (e_F_OII_extinc)
real*4 FHbec (nr__) ! (mW/m2) ? Flux of Hbeta emission line
* adjusted and corrected for
* absorptions and extinction
* (F_Hbeta_extinc)
real*4 e_FHbec (nr__) ! (mW/m2) ? Error in flux of Hbeta emission
* line adjusted and corrected for
* absorptions and extinction
* (e_F_Hbeta_extinc)
real*4 FHbeclim (nr__) ! (mW/m2) ? Upper-limit flux of Hbeta
* adjusted and corrected for
* absorptions and extinction
* (if no fitting was performed)
* (F_Hbeta_lim_extinc)
real*4 FOIIIec (nr__) ! (mW/m2) ? Flux of [OIII]5007 emission line
* adjusted and corrected for
* absorptions and extinction
* (F_OIII_extinc)
real*4 e_FOIIIec (nr__) ! (mW/m2) ? Error in flux of [OIII]5007
* emission line adjusted and
* corrected for absorptions and
* extinction (e_F_OIII_extinc)
real*4 FOIIIeclim (nr__) ! (mW/m2) ? Upper-limit flux of [OIII]5007
* adjusted and corrected for
* absorptions and extinction
* (if no fitting was performed)
* (F_OIII_lim_extinc)
real*8 FluxadjfacOII(nr__) ! ? Flux adjustment factor for
* [OII] emission line
* (flux_adj_factor_OII) (5)
real*8 FluxadjfacHb(nr__) ! ? Flux adjustment factor for Hbeta
* emission line
* (flux_adj_factor_Hbeta) (5)
real*8 FluxadjfacOIII(nr__) ! ? Flux adjustment factor for
* [OIII]5007 emission line
* (flux_adj_factor_OIII) (5)
real*4 EWabsHaCSP (nr__) ! (0.1nm) ? Equivalent width of Halpha line
* obtained from CSP
* (EW_abs_Halpha_CSP)
real*4 EWabsHbCSP (nr__) ! (0.1nm) ? Equivalent width of Hbeta line
* obtained from CSP
* (EW_abs_Hbeta_CSP)
real*8 E_B_V_Balmer(nr__) ! ? Balmer decrement- derived color
* excess of the stellar continuum
* (E_B_V_Balmer)
real*8 E_B_V_CSP (nr__) ! ? Color excess of the stellar
* continuum obtained from CSP
* (E_B_V_CSP)
real*8 E_B_V_used (nr__) ! ? Color excess of the stellar
* continuum used to apply
* correction of extinction
* (E_B_V_used)
character*6 FlagE_B_V (nr__) ! Flag indicating the source of
* E(B_V_used: "Balmer" for Balmer
* decrement, "CSP" for CSP
* (flag_EBV)
real*8 logM_ (nr__) ! ([Msun]) ? Logarithm of the stellar mass
* obtained from CSP
* (logStellar_mass)
real*8 e_logM_ (nr__) ! ([Msun]) ? Error in logarithm of the
* stellar mass obtained from
* CSP (e_logStellar_mass)
real*8 D4000 (nr__) ! ? D4000 break index (D4000)
real*8 e_D4000 (nr__) ! ? Error in D4000 break index
* (e_D4000)
real*8 MetN2IC (nr__) ! ? N2 metallicity (12+log(O/H))
* from inverse convolution
* (Met_N2_InvConv) (3)
real*8 b_MetN2IC (nr__) ! Lower value of N2 metallicity
* (12+log(O/H)) from inverse
* convolution
* (Met_N2_InvConv_d) (3)
real*8 B_MetN2IC_1(nr__) ! Upper value of N2 metallicity
* (12+log(O/H)) from inverse
* convolution
* (Met_N2_InvConv_u) (3)
real*8 MetP05 (nr__) ! ? P05 metallicity (12+log(O/H))
* (Met_P05)
real*8 e_MetP05 (nr__) ! ? Error in P05 metallicity
* (12+log(O/H)) (e_Met_P05)
real*8 MetZ94 (nr__) ! ? Z94 metallicity (12+log(O/H))
* (Met_Z94)
real*8 e_MetZ94 (nr__) ! ? Error in Z94 metallicity
* (12+log(O/H)) (e_Met_Z94)
real*8 MetM91 (nr__) ! ? M91 metallicity (12+log(O/H))
* (Met_M91)
real*8 e_MetM91 (nr__) ! ? Error in M91 metallicity
* (12+log(O/H)) (e_Met_M91)
real*8 MetPP04 (nr__) ! ? PP04 metallicity (12+log(O/H))
* (Met_PP04)
real*8 e_MetPP04 (nr__) ! ? Error in PP04 metallicity
* (12+log(O/H)) (e_Met_PP04)
real*8 MetT04 (nr__) ! ? T04 metallicity (12+log(O/H))
* (Met_T04)
real*8 e_MetT04 (nr__) ! ? Error in T04 metallicity
* (12+log(O/H)) (e_Met_T04)
real*8 SFRHa (nr__) ! (Msun/yr) ? Halpha star formation rate
* (Kennicutt & Evans,
* 2012ARA&A..50..531K) (SFR_Ha)
real*8 e_SFRHa (nr__) ! (Msun/yr) ? Error in Halpha star formation
* rate based on Halpha (Kennicutt &
* Evans, 2012ARA&A..50..531K)
* (e_SFR_Ha)
real*8 SFRHb (nr__) ! (Msun/yr) ? Hbeta star formation rate based
* (Kennicutt & Evans,
* 2012ARA&A..50..531K) (SFR_Hbeta)
real*8 e_SFRHb (nr__) ! (Msun/yr) ? Error in Hbeta star formation
* rate based (Kennicutt & Evans,
* 2012ARA&A..50..531K)
* (e_SFR_Hbeta)
real*8 SFROIIZ94 (nr__) ! (Msun/yr) ? [OII] star formation rate
* corrected by Z94 (Kewley+2004)
* (SFR_OII_Z94)
real*8 e_SFROIIZ94(nr__) ! (Msun/yr) ? Error in [OII] star formation
* rate corrected by Z94
* (Kewley+2004) (e_SFR_OII_Z94)
real*8 SFROIIM91 (nr__) ! (Msun/yr) ? [OII] star formation rate
* corrected by M91 (Kewley et al.
* 2004AJ....127.2002K)
* (SFR_OII_M91)
real*8 e_SFROIIM91(nr__) ! (Msun/yr) ? Error in [OII] star formation
* rate corrected by M91
* (Kewley, 2004AJ....127.2002K)
* (e_SFR_OII_M91)
real*8 Dclusteras (nr__) ! (arcsec) Cluster-centric projected distance
* in arcsec (D_cluster_arcsec) (6)
real*8 DclusterMpc(nr__) ! (Mpc) Cluster-centric projected distance
* in Mpc (D_cluster_Mpc) (6)
real*8 sigma5 (nr__) ! (1/Mpc2) Local density estimated to the
* 5th neighbor (sigma5) (6)
real*8 sigma10 (nr__) ! (1/Mpc2) Local density estimated to the
* 10th neighbor (sigma10) (6)
real*4 isAGNSP15 (nr__) ! True/False flag for AGN activity
* (isAGN_SP15) (1)
real*4 isBLAGNSP15(nr__) ! True/False flag for broad-line AGN
* (isBLAGN_SP15) (1)
real*4 isAGNBPT (nr__) ! [0/1]? True/False flag for AGN
* activity based on BPT
* (Kauffmann+2003); NaN if not
* conclusive (isAGN_BPT)
real*4 isAGNWHAN (nr__) ! [0/1]? True/False flag for AGN
* activity based on EWalphaN2
* (Ho+1997); NaN if not conclusive
* (isAGN_WHAN)
real*4 isSFGMOS (nr__) ! [0/1]True/False flag for
* star-forming galaxies based on
* BPT and EWalphaN2; NaN if not
* conclusive (isSFG_MOS)
real*4 isAGNMOS (nr__) ! [0/1]True/False flag for AGN
* activity based on BPT and
* EWalphaN2; NaN if not conclusive
* (isAGN_MOS)
real*4 isCompMOS (nr__) ! [0/1]? True/False flag for
* composite galaxies based on BPT
* and EWalphaN2; NaN if not
* conclusive (isComposite_MOS)
real*4 isOutlier (nr__) ! [0/1] True/False flag for outlier
* galaxies (TF redshift not
* confident) (is_Outlier)
integer*4 QualIC (nr__) ! [1/6] Quality criterion of
* inverse convolution
* (Quality_InvConv) (3)
*Note (1): Values from Sanchez-Portal et al., 2015, Cat. J/A+A/578/A30
*Note (2): 'a' or 'b'; see Fig. 4 in Sanchez-Portal et al . 2015A&A...578A..30S,
* 'C' and 'O', respectively, in Table 5 of Sanchez-Portal et al.,
* 2015A&A...578A..30S, Cat. J/A+A/578/A30
*Note (3): Following the methodology and data described in Cedres et al.,
* 2024A&A...686A..60C
*Note (4): Based on z_avg or, if unavailable, z_Halpha_InvConv
* (if Quality_InvConv<4) or z_SP15 (if Quality_InvConv>=4).
* If 0.366<z<0.3876, it is "B", else if z>0.3876, it is "A";
* if z<0.366, it is "Interloper"
*Note (5): Fluxes were adjusted due to the use of MOS slits, as described in
* Appendix B.1.
*Note (6): Information provided by Perez-Martinez et al. (in prep.).
C=============================================================================
C Loading file 'catalog.dat' ! Catalog of the properties of the galaxies
C Format for file interpretation
1 format(
+ A7,1X,I5,1X,A5,1X,F9.7,1X,F10.7,1X,F8.6,1X,F8.5,1X,F8.6,1X,
+ F8.6,1X,F8.6,1X,F8.6,1X,F19.17,1X,F22.20,1X,F19.17,1X,F22.20,
+ 1X,F19.17,1X,F22.20,1X,F19.17,1X,F22.20,1X,A10,1X,F18.15,1X,
+ F21.19,1X,F17.14,1X,F21.19,1X,F18.15,1X,F20.18,1X,F18.15,1X,
+ F18.16,1X,F18.15,1X,F18.16,1X,F18.15,1X,F18.16,1X,F18.15,1X,
+ F18.16,1X,E22.17,1X,E22.17,1X,E22.17,1X,E21.16,1X,E22.17,1X,
+ E22.17,1X,E22.17,1X,E22.17,1X,F10.6,1X,F9.6,1X,E22.17,1X,
+ E22.17,1X,E22.17,1X,F9.6,1X,F9.6,1X,F9.6,1X,E22.17,1X,E22.17,
+ 1X,E22.17,1X,E22.17,1X,E11.6,1X,E22.17,1X,E22.17,1X,E22.17,1X,
+ E22.17,1X,E8.6,1X,E11.6,1X,F9.6,1X,F10.6,1X,E8.6,1X,F9.5,1X,
+ F9.6,1X,E11.9,1X,F9.6,1X,F20.15,1X,F20.16,1X,F20.16,1X,F19.16,
+ 1X,F20.16,1X,F19.16,1X,F19.16,1X,F18.16,1X,F19.16,1X,F18.16,
+ 1X,F19.16,1X,F18.16,1X,F19.16,1X,F18.16,1X,F19.16,1X,F18.16,
+ 1X,F20.16,1X,F18.16,1X,F18.16,1X,F18.16,1X,F18.16,1X,F19.16,
+ 1X,F18.16,1X,F18.16,1X,F18.15,1X,F18.16,1X,F18.15,1X,F19.16,
+ 1X,F18.15,1X,F18.16,1X,F19.15,1X,F20.16,1X,F19.16,1X,E22.17,
+ 1X,E22.17,1X,E22.17,1X,E22.17,1X,E22.17,1X,E22.17,1X,E22.17,
+ 1X,E22.17,1X,E22.17,1X,E22.17,1X,E22.17,1X,E22.17,1X,E22.17,
+ 1X,E22.17,1X,E22.17,1X,E22.17,1X,E22.17,1X,E22.17,1X,E22.17,
+ 1X,E22.17,1X,E21.16,1X,E22.17,1X,E22.17,1X,E22.17,1X,E22.17,
+ 1X,E22.17,1X,E22.17,1X,E22.17,1X,E22.17,1X,E22.17,1X,E22.17,
+ 1X,E22.17,1X,E22.17,1X,E22.17,1X,E22.17,1X,F18.16,1X,F18.16,
+ 1X,F18.16,1X,F4.2,1X,F4.2,1X,F19.16,1X,F18.16,1X,F18.16,1X,A6,
+ 1X,F18.15,1X,F20.18,1X,F18.16,1X,F18.16,1X,F7.5,1X,F7.5,1X,
+ F7.5,1X,F18.16,1X,F22.10,1X,F17.15,1X,F18.16,1X,F17.15,1X,
+ F18.16,1X,F17.15,1X,F18.16,1X,F17.15,1X,F18.16,1X,F19.16,1X,
+ F18.16,1X,F19.16,1X,F18.16,1X,F19.16,1X,F18.16,1X,F18.16,1X,
+ F18.16,1X,F18.14,1X,F18.16,1X,F20.16,1X,F19.16,1X,F3.1,1X,
+ F3.1,1X,F3.1,1X,F3.1,1X,F3.1,1X,F3.1,1X,F3.1,1X,F3.1,1X,I1)
C Effective file loading
open(unit=1,status='old',file=
+'catalog.dat')
write(6,*) '....Loading file: catalog.dat'
do i__=1,102
read(1,'(A3288)')ar__
read(ar__,1)
+ ID_unique(i__),ID(i__),ID_MOS(i__),RAdeg(i__),DEdeg(i__),
+ RAGdeg(i__),DEGdeg(i__),zSP15(i__),zHaIC(i__),b_zHaIC(i__),
+ B_zHaIC_1(i__),zOII(i__),e_zOII(i__),zHb(i__),e_zHb(i__),
+ zOIII(i__),e_zOIII(i__),zMOSavg(i__),e_zMOSavg(i__),
+ Structure(i__),Bmag(i__),e_Bmag(i__),Vmag(i__),e_Vmag(i__),
+ Rmag(i__),e_Rmag(i__),Imag(i__),e_Imag(i__),Jmag(i__),
+ e_Jmag(i__),Ksmag(i__),e_Ksmag(i__),F814Wmag(i__),
+ e_F814Wmag(i__),FHaSP15(i__),e_FHaSP15(i__),FNIISP15(i__),
+ e_FNIISP15(i__),FHacSP15(i__),e_FHacSP15(i__),FNIIcSP15(i__),
+ e_FNIIcSP15(i__),EWHaSP15(i__),e_EWHaSP15(i__),ContHaIC(i__),
+ b_ContHaIC(i__),B_ContHaIC_1(i__),sigmaHaIC(i__),
+ b_sigmaHaIC(i__),B_sigmaHaIC_1(i__),FNII6548IC(i__),
+ b_FNII6548IC(i__),B_FNII6548IC_1(i__),FHaIC(i__),b_FHaIC(i__),
+ B_FHaIC_1(i__),FNII6583IC(i__),b_FNII6583IC(i__),
+ B_FNII6583IC_1(i__),EWNII6548IC(i__),b_EWNII6548IC(i__),
+ B_EWNII6548IC_1(i__),EWHaIC(i__),b_EWHaIC(i__),
+ B_EWHaIC_1(i__),EWNII6583IC(i__),b_EWNII6583IC(i__),
+ B_EWNII6583IC_1(i__),EWOII(i__),e_EWOII(i__),EWHb(i__),
+ e_EWHb(i__),EWOIII(i__),e_EWOIII(i__),contOII(i__),
+ e_contOII(i__),contHb(i__),e_contHb(i__),contOIII(i__),
+ e_contOIII(i__),coreOII(i__),e_coreOII(i__),coreHb(i__),
+ e_coreHb(i__),coreOIII(i__),e_coreOIII(i__),sigmaOII(i__),
+ e_sigmaOII(i__),sigmaHb(i__),e_sigmaHb(i__),sigmaOIII(i__),
+ e_sigmaOIII(i__),FWHMOII(i__),e_FWHMOII(i__),FWHMHb(i__),
+ e_FWHMHb(i__),FWHMOIII(i__),e_FWHMOIII(i__),SNRlineOII(i__),
+ SNRlineHb(i__),SNRlineOIII(i__),FOII(i__),e_FOII(i__),
+ FHb(i__),e_FHb(i__),FHblim(i__),FOIII(i__),e_FOIII(i__),
+ FOIIIlim(i__),FOIIadj(i__),e_FOIIadj(i__),FHbadj(i__),
+ e_FHbadj(i__),FHbadjlim(i__),FOIIIadj(i__),e_FOIIIadj(i__),
+ FOIIIadjlim(i__),FHaabs(i__),e_FHaabs(i__),FHbabs(i__),
+ e_FHbabs(i__),FHbabslim(i__),FNII6548ec(i__),
+ e_FNII6548ec(i__),FHaec(i__),e_FHaec(i__),FNII6583ec(i__),
+ e_FNII6583ec(i__),FOIIec(i__),e_FOIIec(i__),FHbec(i__),
+ e_FHbec(i__),FHbeclim(i__),FOIIIec(i__),e_FOIIIec(i__),
+ FOIIIeclim(i__),FluxadjfacOII(i__),FluxadjfacHb(i__),
+ FluxadjfacOIII(i__),EWabsHaCSP(i__),EWabsHbCSP(i__),
+ E_B_V_Balmer(i__),E_B_V_CSP(i__),E_B_V_used(i__),
+ FlagE_B_V(i__),logM_(i__),e_logM_(i__),D4000(i__),
+ e_D4000(i__),MetN2IC(i__),b_MetN2IC(i__),B_MetN2IC_1(i__),
+ MetP05(i__),e_MetP05(i__),MetZ94(i__),e_MetZ94(i__),
+ MetM91(i__),e_MetM91(i__),MetPP04(i__),e_MetPP04(i__),
+ MetT04(i__),e_MetT04(i__),SFRHa(i__),e_SFRHa(i__),SFRHb(i__),
+ e_SFRHb(i__),SFROIIZ94(i__),e_SFROIIZ94(i__),SFROIIM91(i__),
+ e_SFROIIM91(i__),Dclusteras(i__),DclusterMpc(i__),sigma5(i__),
+ sigma10(i__),isAGNSP15(i__),isBLAGNSP15(i__),isAGNBPT(i__),
+ isAGNWHAN(i__),isSFGMOS(i__),isAGNMOS(i__),isCompMOS(i__),
+ isOutlier(i__),QualIC(i__)
if(ar__(96:114) .EQ. '') zOII(i__) = rNULL__
if(ar__(116:137) .EQ. '') e_zOII(i__) = rNULL__
if(ar__(139:157) .EQ. '') zHb(i__) = rNULL__
if(ar__(159:180) .EQ. '') e_zHb(i__) = rNULL__
if(ar__(182:200) .EQ. '') zOIII(i__) = rNULL__
if(ar__(202:223) .EQ. '') e_zOIII(i__) = rNULL__
if(ar__(225:243) .EQ. '') zMOSavg(i__) = rNULL__
if(ar__(245:266) .EQ. '') e_zMOSavg(i__) = rNULL__
if(ar__(279:296) .EQ. '') Bmag(i__) = rNULL__
if(ar__(298:318) .EQ. '') e_Bmag(i__) = rNULL__
if(ar__(320:336) .EQ. '') Vmag(i__) = rNULL__
if(ar__(338:358) .EQ. '') e_Vmag(i__) = rNULL__
if(ar__(360:377) .EQ. '') Rmag(i__) = rNULL__
if(ar__(379:398) .EQ. '') e_Rmag(i__) = rNULL__
if(ar__(400:417) .EQ. '') Imag(i__) = rNULL__
if(ar__(419:436) .EQ. '') e_Imag(i__) = rNULL__
if(ar__(438:455) .EQ. '') Jmag(i__) = rNULL__
if(ar__(457:474) .EQ. '') e_Jmag(i__) = rNULL__
if(ar__(476:493) .EQ. '') Ksmag(i__) = rNULL__
if(ar__(495:512) .EQ. '') e_Ksmag(i__) = rNULL__
if(ar__(514:531) .EQ. '') F814Wmag(i__) = rNULL__
if(ar__(533:550) .EQ. '') e_F814Wmag(i__) = rNULL__
if(ar__(598:619) .EQ. '') FNIISP15(i__) = rNULL__
if(ar__(621:641) .EQ. '') e_FNIISP15(i__) = rNULL__
if(ar__(689:710) .EQ. '') FNIIcSP15(i__) = rNULL__
if(ar__(712:733) .EQ. '') e_FNIIcSP15(i__) = rNULL__
if(ar__(1144:1163) .EQ. '') EWOII(i__) = rNULL__
if(ar__(1165:1184) .EQ. '') e_EWOII(i__) = rNULL__
if(ar__(1186:1205) .EQ. '') EWHb(i__) = rNULL__
if(ar__(1207:1225) .EQ. '') e_EWHb(i__) = rNULL__
if(ar__(1227:1246) .EQ. '') EWOIII(i__) = rNULL__
if(ar__(1248:1266) .EQ. '') e_EWOIII(i__) = rNULL__
if(ar__(1268:1286) .EQ. '') contOII(i__) = rNULL__
if(ar__(1288:1305) .EQ. '') e_contOII(i__) = rNULL__
if(ar__(1307:1325) .EQ. '') contHb(i__) = rNULL__
if(ar__(1327:1344) .EQ. '') e_contHb(i__) = rNULL__
if(ar__(1346:1364) .EQ. '') contOIII(i__) = rNULL__
if(ar__(1366:1383) .EQ. '') e_contOIII(i__) = rNULL__
if(ar__(1385:1403) .EQ. '') coreOII(i__) = rNULL__
if(ar__(1405:1422) .EQ. '') e_coreOII(i__) = rNULL__
if(ar__(1424:1442) .EQ. '') coreHb(i__) = rNULL__
if(ar__(1444:1461) .EQ. '') e_coreHb(i__) = rNULL__
if(ar__(1463:1482) .EQ. '') coreOIII(i__) = rNULL__
if(ar__(1484:1501) .EQ. '') e_coreOIII(i__) = rNULL__
if(ar__(1503:1520) .EQ. '') sigmaOII(i__) = rNULL__
if(ar__(1522:1539) .EQ. '') e_sigmaOII(i__) = rNULL__
if(ar__(1541:1558) .EQ. '') sigmaHb(i__) = rNULL__
if(ar__(1560:1578) .EQ. '') e_sigmaHb(i__) = rNULL__
if(ar__(1580:1597) .EQ. '') sigmaOIII(i__) = rNULL__
if(ar__(1599:1616) .EQ. '') e_sigmaOIII(i__) = rNULL__
if(ar__(1618:1635) .EQ. '') FWHMOII(i__) = rNULL__
if(ar__(1637:1654) .EQ. '') e_FWHMOII(i__) = rNULL__
if(ar__(1656:1673) .EQ. '') FWHMHb(i__) = rNULL__
if(ar__(1675:1693) .EQ. '') e_FWHMHb(i__) = rNULL__
if(ar__(1695:1712) .EQ. '') FWHMOIII(i__) = rNULL__
if(ar__(1714:1731) .EQ. '') e_FWHMOIII(i__) = rNULL__
if(ar__(1733:1751) .EQ. '') SNRlineOII(i__) = rNULL__
if(ar__(1753:1772) .EQ. '') SNRlineHb(i__) = rNULL__
if(ar__(1774:1792) .EQ. '') SNRlineOIII(i__) = rNULL__
if(ar__(1794:1815) .EQ. '') FOII(i__) = rNULL__
if(ar__(1817:1838) .EQ. '') e_FOII(i__) = rNULL__
if(ar__(1840:1861) .EQ. '') FHb(i__) = rNULL__
if(ar__(1863:1884) .EQ. '') e_FHb(i__) = rNULL__
if(ar__(1886:1907) .EQ. '') FHblim(i__) = rNULL__
if(ar__(1909:1930) .EQ. '') FOIII(i__) = rNULL__
if(ar__(1932:1953) .EQ. '') e_FOIII(i__) = rNULL__
if(ar__(1955:1976) .EQ. '') FOIIIlim(i__) = rNULL__
if(ar__(1978:1999) .EQ. '') FOIIadj(i__) = rNULL__
if(ar__(2001:2022) .EQ. '') e_FOIIadj(i__) = rNULL__
if(ar__(2024:2045) .EQ. '') FHbadj(i__) = rNULL__
if(ar__(2047:2068) .EQ. '') e_FHbadj(i__) = rNULL__
if(ar__(2070:2091) .EQ. '') FHbadjlim(i__) = rNULL__
if(ar__(2093:2114) .EQ. '') FOIIIadj(i__) = rNULL__
if(ar__(2116:2137) .EQ. '') e_FOIIIadj(i__) = rNULL__
if(ar__(2139:2160) .EQ. '') FOIIIadjlim(i__) = rNULL__
if(ar__(2162:2183) .EQ. '') FHaabs(i__) = rNULL__
if(ar__(2185:2206) .EQ. '') e_FHaabs(i__) = rNULL__
if(ar__(2208:2229) .EQ. '') FHbabs(i__) = rNULL__
if(ar__(2231:2252) .EQ. '') e_FHbabs(i__) = rNULL__
if(ar__(2254:2274) .EQ. '') FHbabslim(i__) = rNULL__
if(ar__(2276:2297) .EQ. '') FNII6548ec(i__) = rNULL__
if(ar__(2299:2320) .EQ. '') e_FNII6548ec(i__) = rNULL__
if(ar__(2322:2343) .EQ. '') FHaec(i__) = rNULL__
if(ar__(2345:2366) .EQ. '') e_FHaec(i__) = rNULL__
if(ar__(2414:2435) .EQ. '') FOIIec(i__) = rNULL__
if(ar__(2437:2458) .EQ. '') e_FOIIec(i__) = rNULL__
if(ar__(2460:2481) .EQ. '') FHbec(i__) = rNULL__
if(ar__(2483:2504) .EQ. '') e_FHbec(i__) = rNULL__
if(ar__(2506:2527) .EQ. '') FHbeclim(i__) = rNULL__
if(ar__(2529:2550) .EQ. '') FOIIIec(i__) = rNULL__
if(ar__(2552:2573) .EQ. '') e_FOIIIec(i__) = rNULL__
if(ar__(2575:2596) .EQ. '') FOIIIeclim(i__) = rNULL__
if(ar__(2598:2615) .EQ. '') FluxadjfacOII(i__) = rNULL__
if(ar__(2617:2634) .EQ. '') FluxadjfacHb(i__) = rNULL__
if(ar__(2636:2653) .EQ. '') FluxadjfacOIII(i__) = rNULL__
if(ar__(2655:2658) .EQ. '') EWabsHaCSP(i__) = rNULL__
if(ar__(2660:2663) .EQ. '') EWabsHbCSP(i__) = rNULL__
if(ar__(2665:2683) .EQ. '') E_B_V_Balmer(i__) = rNULL__
if(ar__(2685:2702) .EQ. '') E_B_V_CSP(i__) = rNULL__
if(ar__(2704:2721) .EQ. '') E_B_V_used(i__) = rNULL__
if(ar__(2730:2747) .EQ. '') logM_(i__) = rNULL__
if(ar__(2749:2768) .EQ. '') e_logM_(i__) = rNULL__
if(ar__(2770:2787) .EQ. '') D4000(i__) = rNULL__
if(ar__(2789:2806) .EQ. '') e_D4000(i__) = rNULL__
if(ar__(2808:2814) .EQ. '') MetN2IC(i__) = rNULL__
if(ar__(2832:2849) .EQ. '') MetP05(i__) = rNULL__
if(ar__(2851:2872) .EQ. '') e_MetP05(i__) = rNULL__
if(ar__(2874:2890) .EQ. '') MetZ94(i__) = rNULL__
if(ar__(2892:2909) .EQ. '') e_MetZ94(i__) = rNULL__
if(ar__(2911:2927) .EQ. '') MetM91(i__) = rNULL__
if(ar__(2929:2946) .EQ. '') e_MetM91(i__) = rNULL__
if(ar__(2948:2964) .EQ. '') MetPP04(i__) = rNULL__
if(ar__(2966:2983) .EQ. '') e_MetPP04(i__) = rNULL__
if(ar__(2985:3001) .EQ. '') MetT04(i__) = rNULL__
if(ar__(3003:3020) .EQ. '') e_MetT04(i__) = rNULL__
if(ar__(3022:3040) .EQ. '') SFRHa(i__) = rNULL__
if(ar__(3042:3059) .EQ. '') e_SFRHa(i__) = rNULL__
if(ar__(3061:3079) .EQ. '') SFRHb(i__) = rNULL__
if(ar__(3081:3098) .EQ. '') e_SFRHb(i__) = rNULL__
if(ar__(3100:3118) .EQ. '') SFROIIZ94(i__) = rNULL__
if(ar__(3120:3137) .EQ. '') e_SFROIIZ94(i__) = rNULL__
if(ar__(3139:3156) .EQ. '') SFROIIM91(i__) = rNULL__
if(ar__(3158:3175) .EQ. '') e_SFROIIM91(i__) = rNULL__
if(ar__(3264:3266) .EQ. '') isAGNBPT(i__) = rNULL__
if(ar__(3268:3270) .EQ. '') isAGNWHAN(i__) = rNULL__
if(ar__(3280:3282) .EQ. '') isCompMOS(i__) = rNULL__
c ..............Just test output...........
write(6,1)
+ ID_unique(i__),ID(i__),ID_MOS(i__),RAdeg(i__),DEdeg(i__),
+ RAGdeg(i__),DEGdeg(i__),zSP15(i__),zHaIC(i__),b_zHaIC(i__),
+ B_zHaIC_1(i__),zOII(i__),e_zOII(i__),zHb(i__),e_zHb(i__),
+ zOIII(i__),e_zOIII(i__),zMOSavg(i__),e_zMOSavg(i__),
+ Structure(i__),Bmag(i__),e_Bmag(i__),Vmag(i__),e_Vmag(i__),
+ Rmag(i__),e_Rmag(i__),Imag(i__),e_Imag(i__),Jmag(i__),
+ e_Jmag(i__),Ksmag(i__),e_Ksmag(i__),F814Wmag(i__),
+ e_F814Wmag(i__),FHaSP15(i__),e_FHaSP15(i__),FNIISP15(i__),
+ e_FNIISP15(i__),FHacSP15(i__),e_FHacSP15(i__),FNIIcSP15(i__),
+ e_FNIIcSP15(i__),EWHaSP15(i__),e_EWHaSP15(i__),ContHaIC(i__),
+ b_ContHaIC(i__),B_ContHaIC_1(i__),sigmaHaIC(i__),
+ b_sigmaHaIC(i__),B_sigmaHaIC_1(i__),FNII6548IC(i__),
+ b_FNII6548IC(i__),B_FNII6548IC_1(i__),FHaIC(i__),b_FHaIC(i__),
+ B_FHaIC_1(i__),FNII6583IC(i__),b_FNII6583IC(i__),
+ B_FNII6583IC_1(i__),EWNII6548IC(i__),b_EWNII6548IC(i__),
+ B_EWNII6548IC_1(i__),EWHaIC(i__),b_EWHaIC(i__),
+ B_EWHaIC_1(i__),EWNII6583IC(i__),b_EWNII6583IC(i__),
+ B_EWNII6583IC_1(i__),EWOII(i__),e_EWOII(i__),EWHb(i__),
+ e_EWHb(i__),EWOIII(i__),e_EWOIII(i__),contOII(i__),
+ e_contOII(i__),contHb(i__),e_contHb(i__),contOIII(i__),
+ e_contOIII(i__),coreOII(i__),e_coreOII(i__),coreHb(i__),
+ e_coreHb(i__),coreOIII(i__),e_coreOIII(i__),sigmaOII(i__),
+ e_sigmaOII(i__),sigmaHb(i__),e_sigmaHb(i__),sigmaOIII(i__),
+ e_sigmaOIII(i__),FWHMOII(i__),e_FWHMOII(i__),FWHMHb(i__),
+ e_FWHMHb(i__),FWHMOIII(i__),e_FWHMOIII(i__),SNRlineOII(i__),
+ SNRlineHb(i__),SNRlineOIII(i__),FOII(i__),e_FOII(i__),
+ FHb(i__),e_FHb(i__),FHblim(i__),FOIII(i__),e_FOIII(i__),
+ FOIIIlim(i__),FOIIadj(i__),e_FOIIadj(i__),FHbadj(i__),
+ e_FHbadj(i__),FHbadjlim(i__),FOIIIadj(i__),e_FOIIIadj(i__),
+ FOIIIadjlim(i__),FHaabs(i__),e_FHaabs(i__),FHbabs(i__),
+ e_FHbabs(i__),FHbabslim(i__),FNII6548ec(i__),
+ e_FNII6548ec(i__),FHaec(i__),e_FHaec(i__),FNII6583ec(i__),
+ e_FNII6583ec(i__),FOIIec(i__),e_FOIIec(i__),FHbec(i__),
+ e_FHbec(i__),FHbeclim(i__),FOIIIec(i__),e_FOIIIec(i__),
+ FOIIIeclim(i__),FluxadjfacOII(i__),FluxadjfacHb(i__),
+ FluxadjfacOIII(i__),EWabsHaCSP(i__),EWabsHbCSP(i__),
+ E_B_V_Balmer(i__),E_B_V_CSP(i__),E_B_V_used(i__),
+ FlagE_B_V(i__),logM_(i__),e_logM_(i__),D4000(i__),
+ e_D4000(i__),MetN2IC(i__),b_MetN2IC(i__),B_MetN2IC_1(i__),
+ MetP05(i__),e_MetP05(i__),MetZ94(i__),e_MetZ94(i__),
+ MetM91(i__),e_MetM91(i__),MetPP04(i__),e_MetPP04(i__),
+ MetT04(i__),e_MetT04(i__),SFRHa(i__),e_SFRHa(i__),SFRHb(i__),
+ e_SFRHb(i__),SFROIIZ94(i__),e_SFROIIZ94(i__),SFROIIM91(i__),
+ e_SFROIIM91(i__),Dclusteras(i__),DclusterMpc(i__),sigma5(i__),
+ sigma10(i__),isAGNSP15(i__),isBLAGNSP15(i__),isAGNBPT(i__),
+ isAGNWHAN(i__),isSFGMOS(i__),isAGNMOS(i__),isCompMOS(i__),
+ isOutlier(i__),QualIC(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end