FORTRAN Generation
(/./ftp/cats/J/ApJS/262/36)

Conversion of standardized ReadMe file for file /./ftp/cats/J/ApJS/262/36 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 2026-Aug-14
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/262/36  SDSS-IV MaNGA: pyPipe3D data for 10000 galaxies  (Sanchez+, 2022)
*================================================================================
*SDSS-IV MaNGA: pyPipe3D analysis release for 10000 galaxies.
*    Sanchez S.F., Barrera-Ballesteros J.K., Lacerda E., Mejia-Narvaez A.,
*    Camps-farina A., Bruzual G., Espinosa-Ponce C., Rodriguez-Puebla A.,
*    Calette A.R., Ibarra-Medel H., Avila-Reese V., Hernandez-Toledo H.,
*    Bershady M.A., Cano-Diaz M., Munguia-Cordova A.M.
*   <Astrophys. J. Suppl. Ser., 262, 36 (2022)>
*   =2022ApJS..262...36S    (SIMBAD/NED BibCode)
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'catalog.dat'	! Integrated and characteristic parameters delivered
                               for each analyzed data cube

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

C  J2000 position composed of: RAdeg DEdeg
      character*17  Name        ! Manga-plate-ifudsgn unique name
      integer*4     Plate       ! [7443/12772]? Plate ID of the
*                                                  MaNGA cube
      character*5   IFUdsgn     ! IFU bundle ID of the MaNGA cube
      character*11  PlateIFU    ! Code formed by the plate and
*                                                  ifu names
      character*9   MaNGAId     ! MaNGA ID
      real*8        RAdeg       ! (deg) Right ascension of the object
*                                                  (J2000) (objra)
      real*8        DEdeg       ! (deg) [-9.4/65.3] Declination of
*                                                  the object (J2000) (objdec)
      real*8        logSFRHa    ! ([Msun/yr]) [-5.4/3.9]? Integrated SFR
*                                                  derived from the integrated
*                                                  H{alpha} luminosity
*                                                  (log_SFR_Ha)
      real*8        FoV         ! [-0.003/32] Ratio between the
*                                                  diagonal radius of the cube
*                                                  and effective radius
      real*8        Reff        ! (kpc) [-25635/137] Effective radius
*                                                  in kiloparsecs (Re_kpc)
      real*4        e_logMass   ! ([Msun]) [0.05/0.15] logMass
*                                                  uncertainty (e_log_Mass)
      real*4        e_logSFRHa  ! ([Msun/yr]) [0.05/4.5] logSFRHa
*                                                  uncertainty (e_log_SFR_Ha)
      real*8        logMass     ! ([Msun]) [5.9/17.8] Integrated stellar
*                                                  mass in units of the solar
*                                                  mass in log (log_Mass)
      real*8        logSFRs     ! ([Msun/yr]) [-3.9/2.4]? Integrated SFR
*                                                  derived from the SSP analysis
*                                                  t<32Myr (log_SFR_ssp)
      real*8        logNII_Hac  ! [-1.4/5]? Log of the
*                                                  [NII]6583/H{alpha} line ratio
*                                                  in the central aperture
*                                                  (log_NII_Ha_cen)
      real*4        e_logNII_Hac ! [3.5e-5/1.3]? logNII/Hac
*                                                  uncertainty (e_log_NII_Ha_cen)
      real*8        logOIII_Hbc ! [-1.3/4.6]? Log of the
*                                                  [OIII]5007/H{beta} line ratio
*                                                  in the central aperture
*                                                  (log_OIII_Hb_cen)
      real*8        e_logOIII_Hbc ! [0.0015/5]? logOIII/Hbc
*                                                  uncertainty
*                                                  (e_log_OIII_Hb_cen)
      real*8        logSII_Hac  ! [-1/6.3]? Log of the
*                                                  [SII]6717+6731/H{alpha} line
*                                                  ratio in the central aperture
*                                                  (log_SII_Ha_cen)
      real*8        e_logSII_Hac ! [0.0013/3.5]? logSII/Hac
*                                                  uncertainty (e_log_SII_Ha_cen)
      real*8        logOII_Hbc  ! [-0.94/4]? Log of the
*                                                  [OII]3727/H{beta} line ratio
*                                                  in the central aperture
*                                                  (log_OII_Hb_cen)
      real*8        e_logOII_Hbc ! [0.0027/5.3]? logOII/Hbc
*                                                  uncertainty (e_log_OII_Hb_cen)
      real*8        EWHac       ! (0.1nm) [-465.1/-0.001]? Equivalent
*                                                  width of H{alpha} in the
*                                                  central aperture in {AA}
*                                                  (EW_Ha_cen)
      real*8        e_EWHac     ! (0.1nm) [0.006/158]? EWHac
*                                                  uncertainty (e_EW_Ha_cen)
      real*8        logZReLW    ! ([Sun]) [-1.6/0.2] Log of
*                                                  luminosity-weighted (LW)
*                                                  metallicity of the stellar
*                                                  population at the effective
*                                                  radius (ZH_LW_Re_fit)
      real*4        e_logZReLW  ! ([Sun]) [0/0.3]? logZReLW uncertainty
*                                                  (e_ZH_LW_Re_fit)
      real*8        alphaZReLW  ! [-0.9/0.6]? Slope of the
*                                                  gradient of the LW
*                                                  metallicity of the stellar
*                                                  population in dex/Re
*                                                  (alpha_ZH_LW_Re_fit)
      real*4        e_alphaZReLW ! [0/0.3]? alphaZReLW
*                                                  uncertainty
*                                                  (e_alpha_ZH_LW_Re_fit)
      real*8        logZReMW    ! ([Sun]) [-1.7/0.4]? Log of
*                                                  mass-weighted (MW)
*                                                  metallicity of the stellar
*                                                  population at the effective
*                                                  radius (ZH_MW_Re_fit)
      real*4        e_logZReMW  ! ([-]) [0.0029/0.7]? logZReMW
*                                                  uncertainty (e_ZH_MW_Re_fit)
      real*8        alphaZReMW  ! [-1.8/1.6]? Slope of the
*                                                  gradient of the MW
*                                                  log-metallicity of the
*                                                  stellar population in dex/Re
*                                                  (alpha_ZH_MW_Re_fit)
      real*4        e_alphaZReMW ! [0.003/0.7]? alphaZReMW
*                                                  uncertainty
*                                                  (e_alpha_ZH_MW_Re_fit)
      real*8        logAgeReLW  ! ([yr]) [7.4/10]? Log of LW age of
*                                                  the stellar population
*                                                  (Age_LW_Re_fit)
      real*8        e_logAgeReLW ! ([yr]) [0.004/0.8]? logAgeReLW
*                                                  uncertainty (e_Age_LW_Re_fit)
      real*8        alphaAgeReLW ! [-1.9/1.6]? Slope of the
*                                                  gradient of the LW log-age of
*                                                  the stellar population in
*                                                  log(yr)/Re
*                                                  (alpha_Age_LW_Re_fit)
      real*8        e_alphaAgeReLW ! [0.004/1.3]? alphaAgeReLW
*                                                  uncertainty
*                                                  (e_alpha_Age_LW_Re_fit)
      real*8        logAgeReMW  ! ([yr]) [8.9/10.1]? Log of MW age of
*                                                  the stellar population
*                                                  (Age_MW_Re_fit)
      real*8        e_logAgeReMW ! ([yr]) [0.002/0.5]? logAgeReMW
*                                                  uncertainty (e_Age_MW_Re_fit)
      real*8        alphaAgeReMW ! [-1/0.8]? Slope of the
*                                                  gradient of the MW log-age of
*                                                  the stellar population in
*                                                  log(yr)/Re
*                                                  (alpha_Age_MW_Re_fit)
      real*8        e_alphaAgeReMW ! [0.0039/0.8]? alphaAgeReMW
*                                                  uncertainty
*                                                  (e_alpha_Age_MW_Re_fit)
      real*8        ReffA       ! (arcsec) [0.46/162]? Adopted effective
*                                                  radius in arcseconds (Re_arc)
      real*4        DL          ! (Mpc) [3e+24/2.2e+27] Adopted
*                                                  luminosity distance
      real*8        DA          ! (Mpc) [0.0047/2.6] Adopted angular
*                                                  diameter distance
      real*4        PA          ! (deg) [90/270] Adopted position angle
      real*8        ell         ! [0/0.95] Adopted ellipticity
*                                                  (ellip)
      real*8        logMassgas  ! ([Msun]) [5/195]? Log of integrated
*                                                  gas mass (log_Mass_gas)
      real*4        vel_sig     ! [0.017/1.1]?
*                                                  Velocity/dispersion ratio for
*                                                  the stellar populations
*                                                  within 1.5Re (vel_sigma_Re)
      real*4        e_vel_sig   ! [0/0.2]? vel/sig uncertainty
*                                                  (e_vel_sigma_Re)
      real*8        logSFRsf    ! ([Msun/yr]) [-5.5/4]? Integrated SFR
*                                                  using only the spaxels
*                                                  compatible with SF
*                                                  (log_SFR_SF)
      real*8        logSFRdc    ! ([Msun/yr]) [-7.1/4]? Integrated SFR
*                                                  diffuse corrected
*                                                  (log_SFR_D_C)
      real*8        OH_O3N2c    ! [8/9.5]? Oxygen abundance
*                                                  using the calibrator O3N2 at
*                                                  the central aperture
*                                                  (OH_O3N2_cen) (1)
      real*4        e_OH_O3N2c  ! [1.7e-6/0.3]? OH-O3N2c
*                                                  uncertainty (e_OH_O3N2_cen)
      real*8        OH_N2c      ! [7.9/10.2]? Oxygen abundance
*                                                  using the calibrator N2 at
*                                                  the central aperture
*                                                  (OH_N2_cen) (1)
      real*8        e_OH_N2c    ! [0.0002/0.4]? OH-N2c
*                                                  uncertainty (e_OH_N2_cen)
      real*8        OH_ONSc     ! [6.4/9.1]? Oxygen abundance
*                                                  using the calibrator ONS at
*                                                  the central aperture
*                                                  (OH_ONS_cen) (1)
      real*4        e_OH_ONSc   ! [1.8e-5/1.4]? OH-ONSc
*                                                  uncertainty (e_OH_ONS_cen)
      real*8        OH_R23c     ! [-26.4/13.1]? Oxygen
*                                                  abundance using the
*                                                  calibrator R23 at the central
*                                                  aperture (OH_R23_cen) (1)
      real*8        e_OH_R23c   ! [0.00017/57]? OH-R23c
*                                                  uncertainty (e_OH_R23_cen)
      real*8        OH_pyqzc    ! [7.7/9.3]? Oxygen abundance
*                                                  using the calibrator pyqz at
*                                                  the central aperture
*                                                  (OH_pyqz_cen) (1)
      real*8        e_OH_pyqzc  ! [0.0015/0.8]? OH-pyqzc
*                                                  uncertainty (e_OH_pyqz_cen)
      real*8        OH_t2c      ! [8.2/10]? Oxygen abundance
*                                                  using the calibrator t2 at
*                                                  the central aperture
*                                                  (OH_t2_cen) (1)
      real*8        e_OH_t2c    ! [0.00047/0.8]? OH-t2c
*                                                  uncertainty (e_OH_t2_cen)
      real*8        OH_M08c     ! [8/74]? Oxygen abundance
*                                                  using the calibrator M08 at
*                                                  the central aperture
*                                                  (OH_M08_cen) (1)
      real*8        e_OH_M08c   ! [0/100]? OH-M08c uncertainty
*                                                  (e_OH_M08_cen)
      real*8        OH_T04c     ! [-44/9.8]? Oxygen abundance
*                                                  using the calibrator T04 at
*                                                  the central aperture
*                                                  (OH_T04_cen) (1)
      real*8        e_OH_T04c   ! [0/83]? OH-T04c uncertainty
*                                                  (e_OH_T04_cen)
      real*8        OH_dopc     ! [-870/18795]? Oxygen
*                                                  abundance using the
*                                                  calibrator dop at the central
*                                                  aperture (OH_dop_cen) (1)
      real*8        e_OH_dopc   ! [0.00037/33338]? OH-dopc
*                                                  uncertainty (e_OH_dop_cen)
      real*8        OH_O3N2_EPc ! [7/11]? Oxygen abundance
*                                                  using the calibrator
*                                                  O3N2_EPM09 at the central
*                                                  aperture
*                                                  (OH_O3N2_EPM09_cen) (1)
      real*8        e_OH_O3N2_EPc ! [0.0003/2.4]? OH-O3N2-EPc
*                                                  uncertainty
*                                                  (e_OH_O3N2_EPM09_cen)
      real*8        logOI_Hac   ! [-2.8/4.8]? Log of the
*                                                  [OI]6301/H{alpha} within a
*                                                  central aperture
*                                                  (log_OI_Ha_cen)
      real*8        e_logOI_Hac ! [0.00049/3.4]? logOI/Hac
*                                                  uncertainty (e_log_OI_Ha_cen)
      real*8        Ha_Hbc      ! [-29421/334448]? Ratio
*                                                  between the flux of H{alpha}
*                                                  and H{beta} within a central
*                                                  aperture (Ha_Hb_cen)
      real*8        e_Ha_Hbc    ! [-4531/9527]? Ha/Hbc
*                                                  uncertainty (e_Ha_Hb_cen)
      real*8        logNII_HaRe ! [-1.4/4.3]? Log of the
*                                                  [NII]6583/H{alpha} line ratio
*                                                  at 1Re (log_NII_Ha_Re)
      real*8        e_logNII_HaRe ! [0.002/2]? logNII/HaRe
*                                                  uncertainty (e_log_NII_Ha_Re)
      real*8        logOIII_HbRe ! [-1.3/6.6]? Log of the
*                                                  [OIII]5007/H{beta} line ratio
*                                                  at 1Re (log_OIII_Hb_Re)
      real*8        e_logOIII_HbRe ! [0.00013/3.5]? logOIII/HbRe
*                                                  uncertainty (e_log_OIII_Hb_Re)
      real*8        logSII_HaRe ! [-0.99/4.8]? Log of the
*                                                  [SII]6717+6731/H{alpha} ratio
*                                                  at 1Re (log_SII_Ha_Re)
      real*8        e_logSII_HaRe ! [0.0002/3.9]? logSII/HaRe
*                                                  uncertainty (e_log_SII_Ha_Re)
      real*8        logOII_HbRe ! [-0.94/4.5]? Log of the
*                                                  [OII]3727/H{beta} ratio at
*                                                  1Re (log_OII_Hb_Re)
      real*8        e_logOII_HbRe ! [0.0004/3.8]? logOII/HbRe
*                                                  uncertainty (e_log_OII_Hb_Re)
      real*8        logOI_HaRe  ! [-2.6/6]? Log of the
*                                                  [OI]6301/H{alpha} ratio at
*                                                  1Re (log_OI_Ha_Re)
      real*8        e_logOI_HaRe ! [0.00013/5]? logOI/HaRe
*                                                  uncertainty (e_log_OI_Ha_Re)
      real*8        EWHaRe      ! (0.1nm) [-488/-0]? EW of H{alpha} at
*                                                  1Re in {AA} (EW_Ha_Re)
      real*8        e_EWHaRe    ! (0.1nm) [0.012/122]? EWHaRe
*                                                  uncertainty (e_EW_Ha_Re)
      real*8        Ha_HbRe     ! [2.8/123]? Ratio between the
*                                                  flux of H{alpha} and H{beta}
*                                                  within at 1Re (Ha_Hb_Re)
      real*8        e_Ha_HbRe   ! [0/309]? Ha/HbRe uncertainty
*                                                  (e_Ha_Hb_Re)
      real*8        logNII_Ha   ! [-1.3/7.3]? Log of the
*                                                  [NII]6583/H{alpha} line ratio
*                                                  within the entire FOV
*                                                  (log_NII_Ha_ALL)
      real*8        e_logNII_Ha ! [0.00037/2.8]? logNII/Ha
*                                                  uncertainty (e_log_NII_Ha_ALL)
      real*8        logOIII_Hb  ! [-1.2/6.4]? Log of the
*                                                  [OIII]5007/H{beta} line ratio
*                                                  within the FOV
*                                                  (log_OIII_Hb_ALL)
      real*8        e_logOIII_Hb ! [0.00047/6]? logOIII/Hb
*                                                  uncertainty
*                                                  (e_log_OIII_Hb_ALL)
      real*8        logSII_Ha   ! [-0.99/7.8]? Log of the
*                                                  [SII]6717+6731/H{alpha} ratio
*                                                  within the entire FOV
*                                                  (log_SII_Ha_ALL)
      real*8        e_logSII_Ha ! [0.0008/3.5]? logSII/Ha
*                                                  uncertainty (e_log_SII_Ha_ALL)
      real*8        logOII_Hb   ! [-0.99/5.2]? Log of the
*                                                  [OII]3727/H{beta} ratio
*                                                  within the entire FOV
*                                                  (log_OII_Hb_ALL)
      real*8        e_logOII_Hb ! [0.0004/4.4]? logOII/Hb
*                                                  uncertainty (e_log_OII_Hb_ALL)
      real*8        logOI_Ha    ! [-2.6/5]? Log of the
*                                                  [OI]6301/H{alpha} ratio
*                                                  within the entire FOV
*                                                  (log_OI_Ha_ALL)
      real*8        e_logOI_Ha  ! [0.0018/5]? logOI/Ha
*                                                  uncertainty (e_log_OI_Ha_ALL)
      real*8        EWHa        ! (0.1nm) [-285/0]? EW of H{alpha}
*                                                  within the entire FOV in {AA}
*                                                  (EW_Ha_ALL)
      real*8        e_EWHa      ! (0.1nm) [0/159]? EWHa uncertainty
*                                                  (e_EW_Ha_ALL)
      real*4        Ha_Hb       ! [2.8/58]? Ratio between the
*                                                  flux of H{alpha} and H{beta}
*                                                  within the entire FOV
*                                                  (Ha_Hb_ALL)
      real*8        logsigMc    ! ([Msun/pc2]) [0.6/4.5]? Stellar-mass
*                                                  surface density in the
*                                                  central aperture
*                                                  (Sigma_Mass_cen)
      real*8        e_logsigMc  ! ([Msun/pc2]) [0/0.9]? logsigMc uncertainty
*                                                  (e_Sigma_Mass_cen)
      real*8        logsigMRe   ! ([Msun/pc2]) [0.37/4.1]? Stellar-mass
*                                                  surface density at 1Re
*                                                  (Sigma_Mass_Re)
      real*8        e_logsigMRe ! ([Msun/pc2]) [0.0009/0.9]? logsigMRe
*                                                  uncertainty (e_Sigma_Mass_Re)
      real*8        logsigM     ! ([Msun/pc2]) [0.5/3.8]? Average
*                                                  stellar-mass surface density
*                                                  within the entire FOV
*                                                  (Sigma_Mass_ALL)
      real*8        e_logsigM   ! ([Msun/pc2]) [0.07/1.3]? logsigM
*                                                  uncertainty (e_Sigma_Mass_ALL)
      real*8        T30         ! (Gyr) [2.3/13] Lookback time at
*                                                  which the galaxy has formed
*                                                  30% of its stellar mass
      real*8        Z_T30       ! ([Sun]) [-2.31/0.3] Stellar
*                                                  metallicity at T30 time
*                                                  (ZH_T30)
      real*8        Z_T30Re     ! ([Sun]) [-9.4/1.1]? Stellar
*                                                  metallicity at Re at T30 time
*                                                  (ZH_Re_T30)
      real*8        aZ_T30      ! [-1.94/2.1]? Slope of the ZH
*                                                  gradient at T30 time in
*                                                  dex/Re (a_ZH_T30)
      real*8        T40         ! (Gyr) [1.9/12.8] Lookback time at
*                                                  which the galaxy has formed
*                                                  40% of its stellar mass
      real*8        Z_T40       ! ([Sun]) [-2.31/0.3] Stellar
*                                                  metallicity at T40 time
*                                                  (ZH_T40)
      real*8        Z_T40Re     ! ([Sun]) [-9/2.3]? Stellar metallicity
*                                                  at Re at T40 time (ZH_Re_T40)
      real*8        aZ_T40      ! [-2.5/2.1]? Slope of the ZH
*                                                  gradient at T40 time in
*                                                  dex/Re (a_ZH_T40)
      real*8        T50         ! (Gyr) [1.49/12.7] Lookback time at
*                                                  which the galaxy has formed
*                                                  50% of its stellar mass
      real*8        Z_T50       ! ([Sun]) [-2.31/0.3] Stellar
*                                                  metallicity at T50 time
*                                                  (ZH_T50)
      real*8        Z_T50Re     ! ([Sun]) [-5.8/0.8]? Stellar
*                                                  metallicity at Re at T50 time
*                                                  (ZH_Re_T50)
      real*8        aZ_T50      ! [-2.5/2.1]? Slope of the ZH
*                                                  gradient at T50 time in
*                                                  dex/Re (a_ZH_T50)
      real*8        T60         ! (Gyr) [1.26/12.5] Lookback time at
*                                                  which the galaxy has formed
*                                                  60% of its stellar mass
      real*8        Z_T60       ! ([Sun]) [-2.31/0.3] Stellar
*                                                  metallicity at T60 time
*                                                  (ZH_T60)
      real*8        Z_T60Re     ! ([Sun]) [-4.3/2]? Stellar metallicity
*                                                  at Re at T60 time (ZH_Re_T60)
      real*8        aZ_T60      ! [-2.4/2.1]? Slope of the ZH
*                                                  gradient at T60 time in
*                                                  dex/Re (a_ZH_T60)
      real*8        T70         ! (Gyr) [1/12.3] Lookback time at
*                                                  which the galaxy has formed
*                                                  70% of its stellar mass
      real*8        Z_T70       ! ([Sun]) [-2.31/0.3] Stellar
*                                                  metallicity at T70 time
*                                                  (ZH_T70)
      real*8        Z_T70Re     ! ([Sun]) [-3.4/0.7]? Stellar
*                                                  metallicity at Re at T70 time
*                                                  (ZH_Re_T70)
      real*8        aZH_T70     ! [-2.2/2.1]? Slope of the ZH
*                                                  gradient at T70 time in
*                                                  dex/Re (a_ZH_T70)
      real*8        T80         ! (Gyr) [0.88/12.1] Lookback time at
*                                                  which the galaxy has formed
*                                                  80% of its stellar mass
      real*8        Z_T80       ! ([Sun]) [-2.31/0.3] Stellar
*                                                  metallicity at T80 time
*                                                  (ZH_T80)
      real*8        Z_T80Re     ! ([Sun]) [-3.4/0.5]? Stellar
*                                                  metallicity at Re at T80 time
*                                                  (ZH_Re_T80)
      real*8        aZ_T80      ! [-1.95/1.5]? Slope of the ZH
*                                                  gradient at T80 time in
*                                                  dex/Re (a_ZH_T80)
      real*8        T90         ! (Gyr) [0.66/10.6] Lookback time at
*                                                  which the galaxy has formed
*                                                  90% of its stellar mass
      real*8        Z_T90       ! ([Sun]) [-2.1/0.3] Stellar
*                                                  metallicity at T90 time
*                                                  (ZH_T90)
      real*8        Z_T90Re     ! ([Sun]) [-2.7/1.71]? Stellar
*                                                  metallicity at Re at T90 time
*                                                  (ZH_Re_T90)
      real*8        aZ_T90      ! [-1.5/1.2]? Slope of the ZH
*                                                  gradient at T90 time in
*                                                  dex/Re (a_ZH_T90)
      real*8        T95         ! (Gyr) [0.25/7.7] Lookback time at
*                                                  which the galaxy has formed
*                                                  95% of its stellar mass
      real*8        Z_T95       ! ([Sun]) [-2.2/0.3] Stellar
*                                                  metallicity at T95 time
*                                                  (ZH_T95)
      real*8        Z_T95Re     ! ([Sun]) [-1.81/2.8]? Stellar
*                                                  metallicity at Re at T95 time
*                                                  (ZH_Re_T95)
      real*8        aZ_T95      ! [-1.5/1.3]? Slope of the ZH
*                                                  gradient at T95 time in
*                                                  dex/Re (a_ZH_T95)
      real*4        T99         ! (Gyr) [0.0074/4.2] Lookback time at
*                                                  which the galaxy has formed
*                                                  99% of its stellar mass
      real*8        Z_T99       ! ([Sun]) [-1.8/0.3] Stellar
*                                                  metallicity at T99 time
*                                                  (ZH_T99)
      real*8        Z_T99Re     ! ([Sun]) [-2.2/0.4]? Stellar
*                                                  metallicity at Re at T99 time
*                                                  (ZH_Re_T99)
      real*8        aZ_T99      ! [-1.4/1.3]? Slope of the ZH
*                                                  gradient at T99 time in
*                                                  dex/Re (a_ZH_T99)
      real*4        logMgasAvg  ! ([Msun]) [5.7/10.7]? Log of integrated
*                                                  gas mass using the calibrator
*                                                  log_Mass_gas_Av_gas_OH
*                                                  (log_Mass_gas_Av_gas_OH)
      real*4        logMgasAvs  ! ([Msun]) [6.4/10.8]? Log of integrated
*                                                  gas mass using the calibrator
*                                                  log_Mass_gas_Av_ssp_OH
*                                                  (log_Mass_gas_Av_ssp_OH)
      real*8        Vels2Re     ! (km/s) [0.09/1042]? Stellar velocity
*                                                  at 2Re (vel_ssp_2)
      real*8        e_Vels2Re   ! (km/s) [0.029/696]? Vels2Re
*                                                  uncertainty (e_vel_ssp_2)
      real*8        VelHa2Re    ! (km/s) [0/1542]? H{alpha} velocity
*                                                  at 2Re (vel_Ha_2)
      real*8        e_VelHa2Re  ! (km/s) [0/919]? VelHa2Re uncertainty
*                                                  (e_vel_Ha_2)
      real*8        Vels1Re     ! (km/s) [0.028/1580]? Stellar
*                                                  velocity at 1Re (vel_ssp_1)
      real*8        e_Vels1Re   ! (km/s) [0/875]? Vels1Re uncertainty
*                                                  (e_vel_ssp_1)
      real*8        VelHa1Re    ! (km/s) [0.01/2244]? H{alpha}
*                                                  velocity at 1Re (vel_Ha_1)
      real*8        e_VelHa1Re  ! (km/s) [0/865]? VelHa1Re uncertainty
*                                                  (e_vel_Ha_1)
      real*8        logSFRs100  ! ([Msun/yr]) [-4/1.9]? Integrated SFR
*                                                  derived from the SSP analysis
*                                                  for t<100Myr
*                                                  (log_SFR_ssp_100Myr)
      real*8        logSFRs10   ! ([Msun/yr]) [-3.5/2.7]? Integrated SFR
*                                                  derived from the SSP analysis
*                                                  for t<10Myr
*                                                  (log_SFR_ssp_10Myr)
      real*8        VdispHac    ! (km/s) [12/823]? H{alpha} velocity
*                                                  dispersion at the central
*                                                  aperture (vel_disp_Ha_cen)
      real*8        Vdispsc     ! (km/s) [0.28/466]? Stellar velocity
*                                                  dispersion at the central
*                                                  aperture (vel_disp_ssp_cen)
      real*8        VdispHa1re  ! (km/s) [0.09/703]? H{alpha} velocity
*                                                  dispersion at 1Re
*                                                  (vel_disp_Ha_1Re)
      real*8        Vdisps1Re   ! (km/s) [6/490]? Stellar velocity
*                                                  dispersion at 1Re
*                                                  (vel_disp_ssp_1Re)
      real*4        logM1Re     ! ([Msun]) [5.9/16]? Integrated stellar
*                                                  mass within 1 optical Re
*                                                  (log_Mass_in_Re)
      real*4        MLint       ! (Msun/Lsun) [0.2/4.4e+7] V-band
*                                                  mass-to-light ratio from
*                                                  integrated quantities (ML_int)
      real*8        MLavg       ! (Msun/Lsun) [0.00029/547] V-band
*                                                  mass-to-light ratio averaged
*                                                  across the FOV (ML_avg)
      real*8        FHac        ! (10-19W/m2) [0/4454] Flux intensity of
*                                                  H{alpha} in the central
*                                                  aperture in
*                                                  10^-16^erg/s/cm^2^ (F_Ha_cen)
      real*8        e_FHac      ! (10-19W/m2) [0/7021] FHac uncertainty
*                                                  (e_F_Ha_cen)
      real*8        R50V        ! (kpc) [0.1/45]? Radius at which
*                                                  half of the light in the
*                                                  V-band within the FOV is
*                                                  integrated (R50_kpc_V)
      real*8        e_R50V      ! (kpc) [0.02/14]? R50V uncertainty
*                                                  (e_R50_kpc_V)
      real*8        R50M        ! (kpc) [0.1/44.1]? Radius at which
*                                                  half of the stellar mass
*                                                  within the FOV is integrated
*                                                  (R50_kpc_Mass)
      real*8        e_R50M      ! (kpc) [0.02/14]? R50M uncertainty
*                                                  (e_R50_kpc_Mass)
      real*8        logMR50V    ! ([Msun]) [5.8/15.8]? Integrated
*                                                  stellar mass within R50V
*                                                  (log_Mass_corr_in_R50_V)
      real*4        e_logMR50V  ! ([Msun]) [0/1] logMR50V uncertainty
*                                                  (e_log_Mass_corr_in_R50_V)
      real*8        logMgasAvgl ! ([Msun]) [2.59/11.7]? Molecular gas
*                                                  mass derived from the dust
*                                                  extinction using the log-log
*                                                  calibration
*                                                  (log_Mass_gas_Av_gas_log_log)
      real*4        Avgas1Re    ! (mag) [4.4e-5/3.4]? Dust extinction
*                                                  in the V band derived from
*                                                  the H{alpha}/H{beta} line
*                                                  ratio at 1Re (Av_gas_Re)
      real*8        e_Avgas1Re  ! (mag) [0/1.3]? Avgas1Re uncertainty
*                                                  (e_Av_gas_Re)
      real*4        Avs         ! (mag) [5.9e-5/1.8]? Dust extinction
*                                                  in the V band derived from
*                                                  the analysis of the stellar
*                                                  population (Av_ssp_Re)
      real*8        e_Avs       ! (mag) [0/0.6]? Avs uncertainty
*                                                  (e_Av_ssp_Re)
      real*8        Lambda1Re   ! [0.014/1]? Specific angular
*                                                  momentum ({lambda}-parameter)
*                                                  for the stellar populations
*                                                  within 1Re (Lambda_Re)
      real*4        e_Lambda1Re ! [7e-6/0.5]? Lambda1Re
*                                                  uncertainty (e_Lambda_Re)
      real*8        zNSA        ! [0.0002/0.15]? NSA redshift
*                                                  (nsa_redshift)
      real*8        MNSA        ! ([Msun]) [6.2/13.2]? Stellar mass
*                                                  derived by the NSA
*                                                  collaboration (nsa_mstar)
      real*8        IncNSA      ! (deg) [4/86]? Inclination derived
*                                                  by the NSA collaboration
*                                                  (nsa_inclination)
      real*8        F_OII_26Re  ! (10-19W/m2) [-2.4/16]? Flux intensity of
*                                                  [OII]3726.03 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_OII3726_Re_fit)
      real*8        e_F_OII_26Re ! (10-19W/m2) [0.00037/4]? F[OII]26Re
*                                                  uncertainty
*                                                  (e_flux_OII3726_Re_fit)
      real*8        F_OII_26a   ! (10-19W/m2) [-11/1.3]? Flux intensity of
*                                                  [OII]3726.03 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_OII3726_alpha_fit)
      real*8        e_F_OII_26a ! (10-19W/m2) [0.0003/3.2]? F[OII]26a
*                                                  uncertainty
*                                                  (e_flux[OII]3726.03_alpha_fit)
      real*8        F_OII_28Re  ! (10-19W/m2) [-1.3/16]? Flux intensity of
*                                                  [OII]3728.82 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_OII3728_Re_fit)
      real*8        e_F_OII_28Re ! (10-19W/m2) [0.0005/3]? F[OII]28Re
*                                                  uncertainty
*                                                  (e_flux_OII3728_Re_fit)
      real*8        F_OII_28a   ! (10-19W/m2) [-11/2.3]? Flux intensity of
*                                                  [OII]3728.82 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_OII3728_alpha_fit)
      real*8        e_F_OII_28a ! (10-19W/m2) [0.0004/3]? F[OII]28a
*                                                  uncertainty
*                                                  (e_flux[OII]3728.82_alpha_fit)
      real*8        FHI34Re     ! (10-19W/m2) [-2.5/16]? Flux intensity of
*                                                  HI3734.37 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_HI3734_Re_fit)
      real*8        e_FHI34Re   ! (10-19W/m2) [0.0003/2.3]? FHI34Re
*                                                  uncertainty
*                                                  (e_flux_HI3734_Re_fit)
      real*8        FHI34a      ! (10-19W/m2) [-11/1.4]? Flux intensity of
*                                                  HI3734.37 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_HI3734_alpha_fit)
      real*8        e_FHI34a    ! (10-19W/m2) [0.00027/2]? FHI34a
*                                                  uncertainty
*                                                  (e_flux_HI3734_alpha_fit)
      real*8        FHI97Re     ! (10-19W/m2) [-2.1/4.8]? Flux intensity of
*                                                  HI3797.9 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_HI3797.9_Re_fit)
      real*8        e_FHI97Re   ! (10-19W/m2) [0.00024/5]? FHI97Re
*                                                  uncertainty
*                                                  (e_flux_HI3797_Re_fit)
      real*8        FHI97a      ! (10-19W/m2) [-5.82/1.1]? Flux intensity
*                                                  of HI3797.9 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_HI3797_alpha_fit)
      real*8        e_FHI97a    ! (10-19W/m2) [0.0002/4.5]? FHI97a
*                                                  uncertainty
*                                                  (e_flux_HI3797_alpha_fit)
      real*8        FHeI38Re    ! (10-19W/m2) [-0.3/8]? Flux intensity of
*                                                  HeI3888.65 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_HeI3888_Re_fit)
      real*8        e_FHeI38Re  ! (10-19W/m2) [0.00035/2.2]? FHeI38Re
*                                                  uncertainty
*                                                  (e_flux_HeI3888_Re_fit)
      real*8        FHeI38a     ! (10-19W/m2) [-4.8/0.7]? Flux intensity of
*                                                  HeI3888.65 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_HeI3888_alpha_fit)
      real*8        e_FHeI38a   ! (10-19W/m2) [0.0003/1.7]? FHeI38a
*                                                  uncertainty
*                                                  (e_flux_HeI3888_alpha_fit)
      real*8        FHI89Re     ! (10-19W/m2) [-0.31/8]? Flux intensity of
*                                                  HI3889.05 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_HI3889_Re_fit)
      real*8        e_FHI89Re   ! (10-19W/m2) [0.00027/5]? FHI89Re
*                                                  uncertainty
*                                                  (e_flux_HI3889_Re_fit)
      real*8        FHI89a      ! (10-19W/m2) [-4.8/2]? Flux intensity of
*                                                  HI3889.05 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_HI3889_alpha_fit)
      real*8        e_FHI89a    ! (10-19W/m2) [0.0002/4]? FHI89a
*                                                  uncertainty
*                                                  (e_flux_HI3889_alpha_fit)
      real*8        FHeI39Re    ! (10-19W/m2) [-0.83/11]? Flux intensity of
*                                                  HeI3964.73 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_HeI3964_Re_fit)
      real*8        e_FHeI39Re  ! (10-19W/m2) [0.00028/1.2]? FHeI39Re
*                                                  uncertainty
*                                                  (e_flux_HeI3964_Re_fit)
      real*8        FHeI39a     ! (10-19W/m2) [-7/1.4]? Flux intensity of
*                                                  HeI3964.73 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_HeI3964_alpha_fit)
      real*8        e_FHeI39a   ! (10-19W/m2) [0.00029/1.1]? FHeI39a
*                                                  uncertainty
*                                                  (e_flux_HeI3964_alpha_fit)
      real*8        F_NeIII_Re  ! (10-19W/m2) [-0.75/10.3]? Flux intensity
*                                                  of [NeIII]3967.46 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_[NeIII]3967_Re_fit)
      real*8        e_F_NeIII_Re ! (10-19W/m2) [0.00028/2.8]? F[NeIII]Re
*                                                  uncertainty
*                                                  (e_flux_[NeIII]3967_Re_fit)
      real*8        F_NeIII_a   ! (10-19W/m2) [-7.3/1.2]? Flux intensity of
*                                                  [NeIII]3967.46 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux[NeIII]3967.46_alpha_fit)
      real*8        e_F_NeIII_a ! (10-19W/m2) [0.00024/2.4]? F[NeIII]a
*                                                  uncertainty
*                                                  (e_fluxNeIII3967_alpha_fit)
      real*8        FCaIIRe     ! (10-19W/m2) [-0.72/10.4]? Flux intensity
*                                                  of CaII3968.47 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_CaII3968_Re_fit)
      real*8        e_FCaIIRe   ! (10-19W/m2) [0.00026/2.6]? FCaIIRe
*                                                  uncertainty
*                                                  (e_flux_CaII3968_Re_fit)
      real*8        FCaIIa      ! (10-19W/m2) [-7.4/1.2]? Flux intensity of
*                                                  CaII3968.47 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_CaII3968_alpha_fit)
      real*8        e_FCaIIa    ! (10-19W/m2) [0.0002/2.2]? FCaIIa
*                                                  uncertainty
*                                                  (e_flux_CaII3968_alpha_fit)
      real*8        FHepsRe     ! (10-19W/m2) [-0.7/10.4]? Flux intensity
*                                                  of H{epsilon}3970.07 at 1Re
*                                                  in 10^-16^erg/s/cm^2^
*                                                  (flux_Hepsilon3970_Re_fit)
      real*8        e_FHepsRe   ! (10-19W/m2) [0.00025/2.4]? FHepsRe
*                                                  uncertainty
*                                                  (e_fluxHepsilon3970_Re_fit)
      real*8        FHepsa      ! (10-19W/m2) [-8.81/1.6]? Flux intensity
*                                                  of H{epsilon}3970.07 at 1Re
*                                                  in 10^-16^erg/s/cm^2^
*                                                  (flux_Heps3970_alpha_fit)
      real*8        e_FHepsa    ! (10-19W/m2) [0.0002/2.1]? FHepsa
*                                                  uncertainty
*                                                  (e_flux_Heps3970.07_alpha_fit)
      real*8        FHdRe       ! (10-19W/m2) [-0.099/9.7]? Flux intensity
*                                                  of H{delta}4101.77 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_Hdelta4101_Re_fit)
      real*8        e_FHdRe     ! (10-19W/m2) [0.0003/4.3]? FHdRe
*                                                  uncertainty
*                                                  (e_flux_Hdelta4101_Re_fit)
      real*8        FHda        ! (10-19W/m2) [-5.4/2]? Flux intensity of
*                                                  H{delta}4101.77 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_Hdelta4101_alpha_fit)
      real*8        e_FHda      ! (10-19W/m2) [0.00028/3.6]? FHda
*                                                  uncertainty
*                                                  (e_flux_Hdel4101_alpha_fit)
      real*8        FHgRe       ! (10-19W/m2) [-1.25/11.9]? Flux intensity
*                                                  of H{gamma}4340.49 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_Hgamma4340_Re_fit)
      real*8        e_FHgRe     ! (10-19W/m2) [0.0002/5.3]? FHgRe
*                                                  uncertainty
*                                                  (e_flux_Hgamma4340_Re_fit)
      real*8        FHga        ! (10-19W/m2) [-7/6.6]? Flux intensity of
*                                                  H{gamma}4340.49 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_Hgamma4340_alpha_fit)
      real*8        e_FHga      ! (10-19W/m2) [0.00018/4]? FHga uncertainty
*                                                  (e_flux_Hgam4340_alpha_fit)
      real*8        FHbRe       ! (10-19W/m2) [-0.24/14.7]? Flux intensity
*                                                  of H{beta}4861.36 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_Hbeta4861_Re_fit)
      real*8        e_FHbRe     ! (10-19W/m2) [0.00024/8.1]? FHbRe
*                                                  uncertainty
*                                                  (e_flux_Hbeta4861_Re_fit)
      real*8        FHba        ! (10-19W/m2) [-8.1/0.7]? Flux intensity of
*                                                  H{beta}4861.36 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_Hbeta4861_alpha_fit)
      real*8        e_FHba      ! (10-19W/m2) [0.0002/5.6]? FHba
*                                                  uncertainty
*                                                  (e_flux_Hb4861_alpha_fit)
      real*8        F_OIII_49Re ! (10-19W/m2) [-0.9/16]? Flux intensity of
*                                                  [OIII]4958.91 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_[OIII]4958_Re_fit)
      real*8        e_F_OIII_49Re ! (10-19W/m2) [0.0002/17]? F[OIII]49Re
*                                                  uncertainty
*                                                  (e_flux_OIII4958_Re_fit)
      real*8        F_OIII_49a  ! (10-19W/m2) [-9.3/1.3]? Flux intensity of
*                                                  [OIII]4958.91 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_OIII4958_alpha_fit)
      real*8        e_F_OIII_49a ! (10-19W/m2) [0.00017/11]? F[OIII]49a
*                                                  uncertainty
*                                                  (e_fluxOIII4958_alpha_fit)
      real*8        F_OIII_50Re ! (10-19W/m2) [-1.8/12]? Flux intensity of
*                                                  [OIII]5006.84 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_OIII5006_Re_fit)
      real*8        e_F_OIII_50Re ! (10-19W/m2) [0.0002/3.5]? F[OIII]50Re
*                                                  uncertainty
*                                                  (e_flux_OIII5006_Re_fit)
      real*8        F_OIII_50a  ! (10-19W/m2) [-12.6/1.1]? Flux intensity
*                                                  of [OIII]5006.84 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_OIII5006_alpha_fit)
      real*8        e_F_OIII_50a ! (10-19W/m2) [0.00018/3]? F[OIII]50a
*                                                  uncertainty
*                                                  (e_fluxOIII5006_alpha_fit)
      real*8        FHeI50Re    ! (10-19W/m2) [-2.2/8]? Flux intensity of
*                                                  HeI5015.68 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_HeI5015_Re_fit)
      real*8        e_FHeI50Re  ! (10-19W/m2) [0.0002/3.5]? FHeI50Re
*                                                  uncertainty
*                                                  (e_flux_HeI5015_Re_fit)
      real*8        FHeI50a     ! (10-19W/m2) [-11/2]? Flux intensity of
*                                                  HeI5015.68 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_HeI5015_alpha_fit)
      real*8        e_FHeI50a   ! (10-19W/m2) [0.00018/2.5]? FHeI50a
*                                                  uncertainty
*                                                  (e_flux_HeI5015_alpha_fit)
      real*8        F_NI_51Re   ! (10-19W/m2) [-0.22/4]? Flux intensity of
*                                                  [NI]5197.9 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_NI5197_Re_fit)
      real*8        e_F_NI_51Re ! (10-19W/m2) [0.00016/4.2]? F[NI]51Re
*                                                  uncertainty
*                                                  (e_flux_NI5197_Re_fit)
      real*8        F_NI_51a    ! (10-19W/m2) [-4.1/0.7]? Flux intensity of
*                                                  [NI]5197.9 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_NI5197_alpha_fit)
      real*8        e_F_NI_51a  ! (10-19W/m2) [0.00014/3.4]? F[NI]51a
*                                                  uncertainty
*                                                  (e_flux_NI5197_alpha_fit)
      real*8        F_NI_52Re   ! (10-19W/m2) [-0.34/1.7]? Flux intensity
*                                                  of [NI]5200.26 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_NI5200_Re_fit)
      real*8        e_F_NI_52Re ! (10-19W/m2) [0.00018/3.1]? F[NI]52Re
*                                                  uncertainty
*                                                  (e_flux_NI5200_Re_fit)
      real*8        F_NI_52a    ! (10-19W/m2) [-2.6/1.1]? Flux intensity of
*                                                  [NI]5200.26 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_NI5200_alpha_fit)
      real*8        e_F_NI_52a  ! (10-19W/m2) [0.00015/3]? F[NI]52a
*                                                  uncertainty
*                                                  (e_flux_NI5200_alpha_fit)
      real*8        FHeI58Re    ! (10-19W/m2) [-0.6/2.2]? Flux intensity of
*                                                  HeI5876.0 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_HeI5876_Re_fit)
      real*8        e_FHeI58Re  ! (10-19W/m2) [0.00017/1.7]? FHeI58Re
*                                                  uncertainty
*                                                  (e_flux_HeI5876_Re_fit)
      real*8        FHeI58a     ! (10-19W/m2) [-2.8/0.95]? Flux intensity
*                                                  of HeI5876.0 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_HeI5876_alpha_fit)
      real*8        e_FHeI58a   ! (10-19W/m2) [0.00014/1.5]? FHeI58a
*                                                  uncertainty
*                                                  (e_flux_HeI5876_alpha_fit)
      real*8        FNaI89Re    ! (10-19W/m2) [-1/4.6]? Flux intensity of
*                                                  NaI5889.95 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_NaI5889_Re_fit)
      real*8        e_FNaI89Re  ! (10-19W/m2) [0.0002/4.6]? FNaI89Re
*                                                  uncertainty
*                                                  (e_flux_NaI5889_Re_fit)
      real*8        FNaI89a     ! (10-19W/m2) [-5.81/1.4]? Flux intensity
*                                                  of NaI5889.95 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_NaI5889_alpha_fit)
      real*8        e_FNaI89a   ! (10-19W/m2) [0.0002/3.4]? FNaI89a
*                                                  uncertainty
*                                                  (e_flux_NaI5889_alpha_fit)
      real*8        FNaI95Re    ! (10-19W/m2) [-1.01/4.4]? Flux intensity
*                                                  of NaI5895.92 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_NaI5895_Re_fit)
      real*8        e_FNaI95Re  ! (10-19W/m2) [0.00019/4]? FNaI95Re
*                                                  uncertainty
*                                                  (e_flux_NaI5895_Re_fit)
      real*8        FNaI95a     ! (10-19W/m2) [-6.22/2.7]? Flux intensity
*                                                  of NaI5895.92 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_NaI5895_alpha_fit)
      real*8        e_FNaI95a   ! (10-19W/m2) [0.00016/3.2]? FNaI95a
*                                                  uncertainty
*                                                  (e_flux_NaI5895_alpha_fit)
      real*8        F_OI_Re     ! (10-19W/m2) [-1.13/2.6]? Flux intensity
*                                                  of [OI]6300.3 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_OI6300_Re_fit)
      real*8        e_F_OI_Re   ! (10-19W/m2) [0.00014/3.3]? F[OI]Re
*                                                  uncertainty
*                                                  (e_flux_OI6300_Re_fit)
      real*8        F_OI_a      ! (10-19W/m2) [-5.4/0.7]? Flux intensity of
*                                                  [OI]6300.3 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_OI6300_alpha_fit)
      real*8        e_F_OI_a    ! (10-19W/m2) [0.00012/2.7]? F[OI]a
*                                                  uncertainty
*                                                  (e_flux_OI6300_alpha_fit)
      real*8        F_NII_48Re  ! (10-19W/m2) [-0.3/12.4]? Flux intensity
*                                                  of [NII]6548.05 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_NII6548_Re_fit)
      real*8        e_F_NII_48Re ! (10-19W/m2) [0.00016/3.4]? F[NII]48Re
*                                                  uncertainty
*                                                  (e_flux_NII6548_Re_fit)
      real*8        F_NII_48a   ! (10-19W/m2) [-10.1/0.6]? Flux intensity
*                                                  of [NII]6548.05 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_NII6548_alpha_fit)
      real*8        e_F_NII_48a ! (10-19W/m2) [0.00014/2.5]? F[NII]48a
*                                                  uncertainty
*                                                  (e_flux_NII6548_alpha_fit)
      real*8        FHaRe       ! (10-19W/m2) [-0.22/20.2]? Flux intensity
*                                                  of H{alpha}6562.85 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_Halpha6562_Re_fit)
      real*8        e_FHaRe     ! (10-19W/m2) [0.00017/17.6]? FHaRe
*                                                  uncertainty
*                                                  (e_flux_Halpha6562_Re_fit)
      real*8        FHaa        ! (10-19W/m2) [-13/5.1]? Flux intensity of
*                                                  H{alpha}6562.85 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_Halpha6562_alpha_fit)
      real*8        e_FHaa      ! (10-19W/m2) [0.00015/10.6]? FHaa
*                                                  uncertainty
*                                                  (e_flux_Ha6562_alpha_fit)
      real*8        F_NII_83Re  ! (10-19W/m2) [-0.3/16]? Flux intensity of
*                                                  [NII]6583.45 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_NII6583_Re_fit)
      real*8        e_F_NII_83Re ! (10-19W/m2) [0.0002/4.5]? F[NII]83Re
*                                                  uncertainty
*                                                  (e_flux_NII6583_Re_fit)
      real*8        F_NII_83a   ! (10-19W/m2) [-13/1]? Flux intensity of
*                                                  [NII]6583.45 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_NII6583_alpha_fit)
      real*8        e_F_NII_83a ! (10-19W/m2) [0.00017/3.2]? F[NII]83a
*                                                  uncertainty
*                                                  (e_flux_NII6583_alpha_fit)
      real*8        F_SII_16Re  ! (10-19W/m2) [-0.6/5.5]? Flux intensity of
*                                                  [SII]6716.44 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_SII6716_Re_fit)
      real*8        e_F_SII_16Re ! (10-19W/m2) [0.00014/2]? F[SII]16Re
*                                                  uncertainty
*                                                  (e_flux_SII6716_Re_fit)
      real*8        F_SII_16a   ! (10-19W/m2) [-4.6/1]? Flux intensity of
*                                                  [SII]6716.44 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_SII6716_alpha_fit)
      real*8        e_F_SII_16a ! (10-19W/m2) [0.00012/1.5]? F[SII]16a
*                                                  uncertainty
*                                                  (e_flux_SII6716_alpha_fit)
      real*8        F_SII_30Re  ! (10-19W/m2) [-0.6/5.3]? Flux intensity of
*                                                  [SII]6730.82 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_[SII]6730.82_Re_fit)
      real*8        e_F_SII_30Re ! (10-19W/m2) [0.00017/2.1]? F[SII]30Re
*                                                  uncertainty
*                                                  (e_flux_SII6730_Re_fit)
      real*8        F_SII_30a   ! (10-19W/m2) [-5/1]? Flux intensity of
*                                                  [SII]6730.82 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_SII6730_alpha_fit)
      real*8        e_F_SII_30a ! (10-19W/m2) [0.00014/1.8]? F[SII]30a
*                                                  uncertainty
*                                                  (e_flux_SII6730_alpha_fit)
      real*8        F_ArIII_Re  ! (10-19W/m2) [-0.95/1.7]? Flux intensity
*                                                  of [ArIII]7135.8 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_ArIII7135_Re_fit)
      real*8        e_F_ArIII_Re ! (10-19W/m2) [0.00016/2.5]? F[ArIII]Re
*                                                  uncertainty
*                                                  (e_flux_ArIII7135_Re_fit)
      real*8        F_ArIII_a   ! (10-19W/m2) [-3/0.7]? Flux intensity of
*                                                  [ArIII]7135.8 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_ArIII7135_alpha_fit)
      real*8        e_F_ArIII_a ! (10-19W/m2) [0.00015/2]? F[ArIII]a
*                                                  uncertainty
*                                                  (e_flux_ArIII7135_alpha_fit)
      real*8        FHI14Re     ! (10-19W/m2) [-1.6/2.2]? Flux intensity of
*                                                  HI9014.91 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_HI9014_Re_fit)
      real*8        e_FHI14Re   ! (10-19W/m2) [0.00012/3]? FHI14Re
*                                                  uncertainty
*                                                  (e_flux_HI9014_Re_fit)
      real*8        FHI14a      ! (10-19W/m2) [-3.1/1.3]? Flux intensity of
*                                                  HI9014.91 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_HI9014_alpha_fit)
      real*8        e_FHI14a    ! (10-19W/m2) [0.0001/2.5]? FHI14a
*                                                  uncertainty
*                                                  (e_flux_HI9014_alpha_fit)
      real*8        F_SIII_90Re ! (10-19W/m2) [-0.5/3.4]? Flux intensity of
*                                                  [SIII]9069.0 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_SIII9069_Re_fit)
      real*8        e_F_SIII_90Re ! (10-19W/m2) [0.0001/3]? F[SIII]90Re
*                                                  uncertainty
*                                                  (e_flux_SIII9069_Re_fit)
      real*8        F_SIII_90a  ! (10-19W/m2) [-2.8/3]? Flux intensity of
*                                                  [SIII]9069.0 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_SIII9069_alpha_fit)
      real*4        e_F_SIII_90a ! (10-19W/m2) [9.6e-5/2.6]? F[SIII]90a
*                                                  uncertainty
*                                                  (e_flux_SIII9069_alpha_fit)
      real*8        F_FeII_Re   ! (10-19W/m2) [-0.91/0.71]? Flux intensity
*                                                  of [FeII]9470.93 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_FeII9470_Re_fit)
      real*8        e_F_FeII_Re ! (10-19W/m2) [0.00018/1.4]? F[FeII]Re
*                                                  uncertainty
*                                                  (e_flux_FeII9470_Re_fit)
      real*8        F_FeII_a    ! (10-19W/m2) [-1.2/0.5]? Flux intensity of
*                                                  [FeII]9470.93 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_FeII9470_alpha_fit)
      real*8        e_F_FeII_a  ! (10-19W/m2) [0.00016/1.1]? F[FeII]a
*                                                  uncertainty
*                                                  (e_flux_FeII9470_alpha_fit)
      real*8        F_SIII_95Re ! (10-19W/m2) [-0.92/8]? Flux intensity of
*                                                  [SIII]9531.1 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_SIII9531_Re_fit)
      real*8        e_F_SIII_95Re ! (10-19W/m2) [0.0003/1.6]? F[SIII]95Re
*                                                  uncertainty
*                                                  (e_flux_SIII9531_Re_fit)
      real*8        F_SIII_95a  ! (10-19W/m2) [-5.3/2]? Flux intensity of
*                                                  [SIII]9531.1 at 1Re in
*                                                  10^-16^erg/s/cm^2^
*                                                  (flux_SIII9531_alpha_fit)
      real*8        e_F_SIII_95a ! (10-19W/m2) [0.00027/1.4]? F[SIII]95a
*                                                  uncertainty
*                                                  (e_flux_SIII9531_alpha_fit)
      real*4        OH_Mar13_N2Re ! [8/8.8]? Oxygen abundance
*                                                 using the calibrator Mar13_N2
*                                                 at 1Re (OH_Mar13_N2_Re_fit) (1)
      real*4        e_OH_Mar13_N2Re ! [0.001/0.4]? OH-Mar13-N2Re
*                                                  uncertainty
*                                                  (e_OH_Mar13_N2_Re_fit)
      real*8        OH_Mar13_N2a ! [-0.5/0.5]? Slope of the O/H
*                                                  gradient using the calibrator
*                                                  Mar13_N2 in dex/Re
*                                                  (OH_Mar13_N2_alpha_fit)
      real*8        e_OH_Mar13_N2a ! [0.0008/0.21]? OH-Mar13-N2a
*                                                  uncertainty
*                                                  (e_OH_Mar13_N2_alpha_fit)
      real*8        OH_Mar13_O3N2Re ! [8.14/8.7]? Oxygen abundance
*                                                  using the calibrator
*                                                  Mar13_O3N2 at 1Re
*                                                  (OH_Mar13_O3N2_Re_fit) (1)
      real*8        e_OH_Mar13_O3N2Re ! [0.00068/0.4]?
*                                                  OH-Mar13-O3N2Re uncertainty
*                                                  (e_OH_Mar13_O3N2_Re_fit)
      real*8        OH_Mar13_O3N2a ! [-0.42/0.3]? Slope of the O/H
*                                                  gradient using the calibrator
*                                                  Mar13_O3N2 in dex/Re
*                                                  (OH_Mar13_O3N2_alpha_fit)
      real*8        e_OH_Mar13_O3N2a ! [0.00058/0.22]?
*                                                  OH-Mar13-O3N2a uncertainty
*                                                  (e_OH_Mar13_O3N2_alpha_fit)
      real*8        OH_T04Re    ! [7/11]? Oxygen abundance
*                                                  using the calibrator T04 at
*                                                  1Re (OH_T04_Re_fit) (1)
      real*4        e_OH_T04Re  ! [0.003/5.1]? OH-T04Re
*                                                  uncertainty (e_OH_T04_Re_fit)
      real*8        OH_T04a     ! [-3.9/2.8]? Slope of the O/H
*                                                  gradient using the calibrator
*                                                  T04 in dex/Re
*                                                  (OH_T04_alpha_fit)
      real*4        e_OH_T04a   ! [0.0029/3.1]? OH-T04a
*                                                  uncertainty
*                                                  (e_OH_T04_alpha_fit)
      real*4        OH_Pet04_N2linRe ! [7.9/9]? Oxygen abundance
*                                                  using the calibrator
*                                                  Pet04_N2_lin at 1Re
*                                                  (OH_Pet04_N2_lin_Re_fit) (1)
      real*4        e_OH_Pet04_N2linRe ! [0.001/0.5]? OH-Pet04-N2linRe
*                                                  uncertainty
*                                                  (e_OH_Pet04_N2_lin_Re_fit)
      real*8        OH_Pet04_N2lina ! [-0.9/0.6]? Slope of the O/H
*                                                  gradient using the calibrator
*                                                  Pet04_N2_lin in dex/Re
*                                                  (OH_Pet04_N2_lin_alpha_fit)
      real*8        e_OH_Pet04_N2lina ! [0.0008/0.3]? OH-Pet04-N2lina
*                                                  uncertainty
*                                                  (e_OH_Pet04_N2_lin_alpha_fit)
      real*4        OH_Pet04_N2polRe ! [7.9/9.1]? Oxygen abundance
*                                                  using the calibrator
*                                                  Pet04_N2_poly at 1Re
*                                                  (OH_Pet04_N2_poly_Re_fit) (1)
      real*4        e_OH_Pet04_N2polRe ! [0.002/0.6]? OH-Pet04-N2polRe
*                                                  uncertainty
*                                                  (e_OH_Pet04_N2_poly_Re_fit)
      real*8        OH_Pet04_N2pola ! [-0.9/0.6]? Slope of the O/H
*                                                  gradient using the calibrator
*                                                  Pet04_N2_poly in dex/Re
*                                                  (OH_Pet04_N2_poly_alpha_fit)
      real*4        e_OH_Pet04_N2pola ! [0.0017/0.4]? OH-Pet04-N2pola
*                                                  uncertainty
*                                                  (e_OH_Pet04_N2_poly_alpha_fit)
      real*4        OH_Pet04_O3N2Re ! [8/9]? Oxygen abundance using
*                                                  the calibrator Pet04_O3N2 at
*                                                  1Re (OH_Pet04_O3N2_Re_fit) (1)
      real*4        e_OH_Pet04_O3N2Re ! [0.001/0.6]? OH-Pet04-O3N2Re
*                                                  uncertainty
*                                                  (e_OH_Pet04_O3N2_Re_fit)
      real*8        OH_Pet04_O3N2a ! [-0.7/0.5]? Slope of the O/H
*                                                  gradient using the calibrator
*                                                  Pet04_O3N2 in dex/Re
*                                                  (OH_Pet04_O3N2_alpha_fit)
      real*8        e_OH_Pet04_O3N2a ! [0.0008/0.4]? OH-Pet04-O3N2a
*                                                  uncertainty
*                                                  (e_OH_Pet04_O3N2_alpha_fit)
      real*4        OH_Kew02_N2O2Re ! [8/9.5]? Oxygen abundance
*                                                  using the calibrator
*                                                  Kew02_N2O2 at 1Re
*                                                  (OH_Kew02_N2O2_Re_fit) (1)
      real*4        e_OH_Kew02_N2O2Re ! [0.0018/1.3]? OH-Kew02-N2O2Re
*                                                  uncertainty
*                                                  (e_OH_Kew02_N2O2_Re_fit)
      real*8        OH_Kew02_N2O2a ! [-1.1/1.1]? Slope of the O/H
*                                                  gradient using the calibrator
*                                                  Kew02_N2O2 in dex/Re
*                                                  (OH_Kew02_N2O2_alpha_fit)
      real*4        e_OH_Kew02_N2O2a ! [0.0015/1]? OH-Kew02-N2O2a
*                                                  uncertainty
*                                                  (e_OH_Kew02_N2O2_alpha_fit)
      real*8        OH_Pil10_ONSRe ! [7.5/9]? Oxygen abundance
*                                                  using the calibrator
*                                                  Pil10_ONS at 1Re
*                                                  (OH_Pil10_ONS_Re_fit) (1)
      real*8        e_OH_Pil10_ONSRe ! [0.0009/2]? OH-Pil10-ONSRe
*                                                  uncertainty
*                                                  (e_OH_Pil10_ONS_Re_fit)
      real*8        OH_Pil10_ONSa ! [-1.6/1.3]? Slope of the O/H
*                                                  gradient using the calibrator
*                                                  Pil10_ONS
*                                                  (OH_Pil10_ONS_alpha_fit)
      real*8        e_OH_Pil10_ONSa ! [0.0008/2.3]? OH-Pil10-ONSa
*                                                  uncertainty
*                                                  (e_OH_Pil10_ONS_alpha_fit)
      real*4        OH_Pil10_ONRe ! [7.5/9]? Oxygen abundance
*                                                 using the calibrator Pil10_ON
*                                                 at 1Re (OH_Pil10_ON_Re_fit) (1)
      real*4        e_OH_Pil10_ONRe ! [0.0014/2]? OH-Pil10-ONRe
*                                                  uncertainty
*                                                  (e_OH_Pil10_ON_Re_fit)
      real*8        OH_Pil10_ONa ! [-1.7/1.3]? Slope of the O/H
*                                                  gradient using the calibrator
*                                                  Pil10_ON in dex/Re
*                                                  (OH_Pil10_ON_alpha_fit)
      real*4        e_OH_Pil10_ONa ! [0.0012/2.2]? OH-Pil10-ONa
*                                                  uncertainty
*                                                  (e_OH_Pil10_ON_alpha_fit)
      real*4        OH_Pil11_NSRe ! [7/9.3]? Oxygen abundance
*                                                 using the calibrator Pil11_NS
*                                                 at 1Re (OH_Pil11_NS_Re_fit) (1)
      real*4        e_OH_Pil11_NSRe ! [0.001/2.6]? OH-Pil11-NSRe
*                                                  uncertainty
*                                                  (e_OH_Pil11_NS_Re_fit)
      real*8        OH_Pil11_NSa ! [-1.7/3.4]? Slope of the O/H
*                                                  gradient using the calibrator
*                                                  Pil11_NS in dex/Re
*                                                  (OH_Pil11_NS_alpha_fit)
      real*8        e_OH_Pil11_NSa ! [0.0009/2.2]? OH-Pil11-NSa
*                                                  uncertainty
*                                                  (e_OH_Pil11_NS_alpha_fit)
      real*4        OH_Cur20_RS32Re ! [7.9/9]? Oxygen abundance
*                                                  using the calibrator
*                                                  Cur20_RS32 at 1Re
*                                                  (OH_Cur20_RS32_Re_fit) (1)
      real*4        e_OH_Cur20_RS32Re ! [0.001/0.6]? OH-Cur20-RS32Re
*                                                  uncertainty
*                                                  (e_OH_Cur20_RS32_Re_fit)
      real*8        OH_Cur20_RS32a ! [-0.6/0.8]? Slope of the O/H
*                                                  gradient using the calibrator
*                                                  Cur20_RS32 in dex/Re
*                                                  (OH_Cur20_RS32_alpha_fit)
      real*8        e_OH_Cur20_RS32a ! [0.0008/0.5]? OH-Cur20-RS32a
*                                                  uncertainty
*                                                  (e_OH_Cur20_RS32_alpha_fit)
      real*4        OH_Cur20_R3Re ! [7.9/9]? Oxygen abundance
*                                                 using the calibrator Cur20_R3
*                                                 at 1Re (OH_Cur20_R3_Re_fit) (1)
      real*4        e_OH_Cur20_R3Re ! [0.001/0.6]? OH-Cur20-R3Re
*                                                  uncertainty
*                                                  (e_OH_Cur20_R3_Re_fit)
      real*8        OH_Cur20_R3a ! [-0.6/0.8]? Slope of the O/H
*                                                  gradient using the calibrator
*                                                  Cur20_R3 in dex/Re
*                                                  (OH_Cur20_R3_alpha_fit)
      real*8        e_OH_Cur20_R3a ! [0.0009/0.4]? OH-Cur20-R3a
*                                                  uncertainty
*                                                  (e_OH_Cur20_R3_alpha_fit)
      real*4        OH_Cur20_O3O2Re ! [7.9/9]? Oxygen abundance
*                                                  using the calibrator
*                                                  Cur20_O3O2 at 1Re
*                                                  (OH_Cur20_O3O2_Re_fit) (1)
      real*4        e_OH_Cur20_O3O2Re ! [0.0013/0.7]? OH-Cur20-O3O2Re
*                                                  uncertainty
*                                                  (e_OH_Cur20_O3O2_Re_fit)
      real*8        OH_Cur20_O3O2a ! [-0.6/0.5]? Slope of the O/H
*                                                  gradient using the calibrator
*                                                  Cur20_O3O2 in dex/Re
*                                                  (OH_Cur20_O3O2_alpha_fit)
      real*4        e_OH_Cur20_O3O2a ! [0.001/0.6]? OH-Cur20-O3O2a
*                                                  uncertainty
*                                                  (e_OH_Cur20_O3O2_alpha_fit)
      real*4        OH_Cur20_S2Re ! [7.8/9]? Oxygen abundance
*                                                 using the calibrator Cur20_S2
*                                                 at 1Re (OH_Cur20_S2_Re_fit) (1)
      real*4        e_OH_Cur20_S2Re ! [0.0017/1.2]? OH-Cur20-S2Re
*                                                  uncertainty
*                                                  (e_OH_Cur20_S2_Re_fit)
      real*8        OH_Cur20_S2a ! [-0.8/0.9]? Slope of the O/H
*                                                  gradient using the calibrator
*                                                  Cur20_S2 in dex/Re
*                                                  (OH_Cur20_S2_alpha_fit)
      real*4        e_OH_Cur20_S2a ! [0.0016/0.8]? OH-Cur20-S2a
*                                                  uncertainty
*                                                  (e_OH_Cur20_S2_alpha_fit)
      real*4        OH_Cur20_R2Re ! [7.9/9]? Oxygen abundance
*                                                 using the calibrator Cur20_R2
*                                                 at 1Re (OH_Cur20_R2_Re_fit) (1)
      real*4        e_OH_Cur20_R2Re ! [0.0019/0.5]? OH-Cur20-R2Re
*                                                  uncertainty
*                                                  (e_OH_Cur20_R2_Re_fit)
      real*8        OH_Cur20_R2a ! [-0.5/1.1]? Slope of the O/H
*                                                  gradient using the calibrator
*                                                  Cur20_R2 in dex/Re
*                                                  (OH_Cur20_R2_alpha_fit)
      real*4        e_OH_Cur20_R2a ! [0.0018/0.4]? OH-Cur20-R2a
*                                                  uncertainty
*                                                  (e_OH_Cur20_R2_alpha_fit)
      real*4        OH_Cur20_N2Re ! [7.9/9]? Oxygen abundance
*                                                 using the calibrator Cur20_N2
*                                                 at 1Re (OH_Cur20_N2_Re_fit) (1)
      real*4        e_OH_Cur20_N2Re ! [0.0014/0.6]? OH-Cur20-N2Re
*                                                  uncertainty
*                                                  (e_OH_Cur20_N2_Re_fit)
      real*8        OH_Cur20_N2a ! [-0.8/0.6]? Slope of the O/H
*                                                  gradient using the calibrator
*                                                  Cur20_N2 in dex/Re
*                                                  (OH_Cur20_N2_alpha_fit)
      real*4        e_OH_Cur20_N2a ! [0.0012/0.4]? OH-Cur20-N2a
*                                                  uncertainty
*                                                  (e_OH_Cur20_N2_alpha_fit)
      real*4        OH_Cur20_R23Re ! [7.8/9]? Oxygen abundance
*                                                  using the calibrator
*                                                  Cur20_R23 at 1Re
*                                                  (OH_Cur20_R23_Re_fit) (1)
      real*4        e_OH_Cur20_R23Re ! [0.0025/0.8]? OH-Cur20-R23Re
*                                                  uncertainty
*                                                  (e_OH_Cur20_R23_Re_fit)
      real*8        OH_Cur20_R23a ! [-1.2/0.7]? Slope of the O/H
*                                                  gradient using the calibrator
*                                                  Cur20_R23 in dex/Re
*                                                  (OH_Cur20_R23_alpha_fit)
      real*4        e_OH_Cur20_R23a ! [0.002/0.7]? OH-Cur20-R23a
*                                                  uncertainty
*                                                  (e_OH_Cur20_R23_alpha_fit)
      real*8        OH_Cur20_O3N2Re ! [8/9]? Oxygen abundance using
*                                                  the calibrator Cur20_O3N2 at
*                                                  1Re (OH_Cur20_O3N2_Re_fit) (1)
      real*8        e_OH_Cur20_O3N2Re ! [0.0007/0.5]? OH-Cur20-O3N2Re
*                                                  uncertainty
*                                                  (e_OH_Cur20_O3N2_Re_fit)
      real*8        OH_Cur20_O3N2a ! [-0.6/0.4]? Slope of the O/H
*                                                  gradient using the calibrator
*                                                  Cur20_O3N2 in dex/Re
*                                                  (OH_Cur20_O3N2_alpha_fit)
      real*8        e_OH_Cur20_O3N2a ! [0.0006/0.3]? OH-Cur20-O3N2a
*                                                  uncertainty
*                                                  (e_OH_Cur20_O3N2_alpha_fit)
      real*4        OH_Cur20_O3S2Re ! [8/9]? Oxygen abundance using
*                                                  the calibrator Cur20_O3S2 at
*                                                  1Re (OH_Cur20_O3S2_Re_fit) (1)
      real*4        e_OH_Cur20_O3S2Re ! [0.002/0.9]? OH-Cur20-O3S2Re
*                                                  uncertainty
*                                                  (e_OH_Cur20_O3S2_Re_fit)
      real*8        OH_Cur20_O3S2a ! [-0.8/0.6]? Slope of the O/H
*                                                  gradient using the calibrator
*                                                  Cur20_O3S2 in dex/Re
*                                                  (OH_Cur20_O3S2_alpha_fit)
      real*4        e_OH_Cur20_O3S2a ! [0.0019/0.7]? OH-Cur20-O3S2a
*                                                  uncertainty
*                                                  (e_OH_Cur20_O3S2_alpha_fit)
      real*8        OH_KK04Re   ! [7.8/9.8]? Oxygen abundance
*                                                  using the calibrator KK04 at
*                                                  1Re (OH_KK04_Re_fit) (1)
      real*8        e_OH_KK04Re ! [0.0008/2.1]? OH-KK04Re
*                                                  uncertainty (e_OH_KK04_Re_fit)
      real*8        OH_KK04a    ! [-1.8/3.4]? Slope of the O/H
*                                                  gradient using the calibrator
*                                                  KK04 in dex/Re
*                                                  (OH_KK04_alpha_fit)
      real*8        e_OH_KK04a  ! [0.0007/2.1]? OH-KK04a
*                                                  uncertainty
*                                                  (e_OH_KK04_alpha_fit)
      real*4        OH_Pil16_RRe ! [7/9.1]? Oxygen abundance
*                                                  using the calibrator Pil16_R
*                                                  at 1Re (OH_Pil16_R_Re_fit) (1)
      real*4        e_OH_Pil16_RRe ! [0.0019/2.4]? OH-Pil16-RRe
*                                                  uncertainty
*                                                  (e_OH_Pil16_R_Re_fit)
      real*8        OH_Pil16_Ra ! [-1.6/1.5]? Slope of the O/H
*                                                  gradient using the calibrator
*                                                  Pil16_R in dex/Re
*                                                  (OH_Pil16_R_alpha_fit)
      real*4        e_OH_Pil16_Ra ! [0.0016/1.7]? OH-Pil16-Ra
*                                                  uncertainty
*                                                  (e_OH_Pil16_R_alpha_fit)
      real*4        OH_Pil16_SRe ! [7.2/9]? Oxygen abundance
*                                                  using the calibrator Pil16_S
*                                                  at 1Re (OH_Pil16_S_Re_fit) (1)
      real*4        e_OH_Pil16_SRe ! [0.0014/1.4]? OH-Pil16-SRe
*                                                  uncertainty
*                                                  (e_OH_Pil16_S_Re_fit)
      real*8        OH_Pil16_Sa ! [-1.8/1.5]? Slope of the O/H
*                                                  gradient using the calibrator
*                                                  Pil16_S in dex/Re
*                                                  (OH_Pil16_S_alpha_fit)
      real*4        e_OH_Pil16_Sa ! [0.0012/1.4]? OH-Pil16-Sa
*                                                  uncertainty
*                                                  (e_OH_Pil16_S_alpha_fit)
      real*4        OH_HoRe     ! [7.7/9.6]? Oxygen abundance
*                                                  using the calibrator Ho at
*                                                  1Re (OH_Ho_Re_fit) (1)
      real*4        e_OH_HoRe   ! [0.0015/1.1]? OH-HoRe
*                                                  uncertainty (e_OH_Ho_Re_fit)
      real*8        OH_Hoa      ! [-0.9/1.2]? Slope of the O/H
*                                                  gradient using the calibrator
*                                                  Ho in dex/Re (OH_Ho_alpha_fit)
      real*4        e_OH_Hoa    ! [0.0013/1.2]? OH-Hoa
*                                                  uncertainty
*                                                  (e_OH_Ho_alpha_fit)
      real*8        U_Dors_O32Re ! [-4.5/-1]? Log(U) ionization
*                                                  parameter using the
*                                                  calibrator Dors_O32 at 1Re
*                                                  (U_Dors_O32_Re_fit) (1)
      real*4        e_U_Dors_O32Re ! [0.0048/3]? U-Dors-O32Re
*                                                  uncertainty
*                                                  (e_U_Dors_O32_Re_fit)
      real*8        U_Dors_O32a ! [-1.93/2.4]? Slope of the
*                                                  log(U) gradient using the
*                                                  calibrator Dors_O32 in dex/Re
*                                                  (U_Dors_O32_alpha_fit)
      real*4        e_U_Dors_O32a ! [0.004/2]? U-Dors-O32a
*                                                  uncertainty
*                                                  (e_U_Dors_O32_alpha_fit)
      real*8        U_Dors_SRe  ! [-5.21/1.8]? Log(U)
*                                                  ionization parameter using
*                                                  the calibrator Dors_S at 1Re
*                                                  (U_Dors_S_Re_fit) (1)
      real*4        e_U_Dors_SRe ! [0.004/4]? U-Dors-SRe
*                                                  uncertainty
*                                                  (e_U_Dors_S_Re_fit)
      real*8        U_Dors_Sa   ! [-3.1/2.4]? Slope of the
*                                                  log(U) gradient using the
*                                                  calibrator Dors_S in dex/Re
*                                                  (U_Dors_S_alpha_fit)
      real*4        e_U_Dors_Sa ! [0.003/2.5]? U-Dors-Sa
*                                                  uncertainty
*                                                  (e_U_Dors_S_alpha_fit)
      real*8        U_Mor16_O23fsRe ! [-7/-0.3]? Log(U) ionization
*                                                 parameter using the
*                                                 calibrator Mor16_O23_fs at
*                                                 1Re (U_Mor16_O23_fs_Re_fit) (1)
      real*4        e_U_Mor16_O23fsRe ! [0.009/6]? U-Mor16-O23fsRe
*                                                  uncertainty
*                                                  (e_U_Mor16_O23_fs_Re_fit)
      real*8        U_Mor16_O23fsa ! [-3.8/4.5]? Slope of the
*                                                  log(U) gradient using the
*                                                  calibrator Mor16_O23_fs in
*                                                  dex/Re
*                                                  (U_Mor16_O23_fs_alpha_fit)
      real*4        e_U_Mor16_O23fsa ! [0.007/4]? U-Mor16-O23fsa
*                                                  uncertainty
*                                                  (e_U_Mor16_O23_fs_alpha_fit)
      real*8        U_Mor16_O23tsRe ! [-4.7/-1.8]? Log(U)
*                                                  ionization parameter using
*                                                  the calibrator Mor16_O23_ts
*                                                  at 1Re
*                                                  (U_Mor16_O23_ts_Re_fit) (1)
      real*4        e_U_Mor16_O23tsRe ! [0.0039/2.4]? U-Mor16-O23tsRe
*                                                  uncertainty
*                                                  (e_U_Mor16_O23_ts_Re_fit)
      real*8        U_Mor16_O23tsa ! [-1.6/2]? Slope of the log(U)
*                                                  gradient using the calibrator
*                                                  Mor16_O23_ts in dex/Re
*                                                  (U_Mor16_O23_ts_alpha_fit)
      real*4        e_U_Mor16_O23tsa ! [0.003/1.6]? U-Mor16-O23tsa
*                                                  uncertainty
*                                                  (e_U_Mor16_O23_ts_alpha_fit)
      real*4        NH_Pil16RRe ! [5.3/9.1]? Nitrogen abundance
*                                                  using the calibrator Pil16_R
*                                                  at 1Re (NH_Pil16_R_Re_fit) (1)
      real*4        e_NH_Pil16RRe ! [0.003/2]? NH-Pil16RRe
*                                                  uncertainty
*                                                  (e_NH_Pil16_R_Re_fit)
      real*8        NH_Pil16Ra  ! [-2.7/2.6]? Slope of the N/H
*                                                  gradient using the calibrator
*                                                  Pil16_R in dex/Re
*                                                  (NH_Pil16_R_alpha_fit)
      real*4        e_NH_Pil16Ra ! [0.003/2]? NH-Pil16Ra
*                                                  uncertainty
*                                                  (e_NH_Pil16_R_alpha_fit)
      real*8        NO_Pil16RRe ! [-2.7/0.5]? N/O abundance
*                                                  using the calibrator Pil16_R
*                                                  at 1Re (NO_Pil16_R_Re_fit) (1)
      real*4        e_NO_Pil16RRe ! [0.003/3.1]? NO-Pil16RRe
*                                                  uncertainty
*                                                  (e_NO_Pil16_R_Re_fit)
      real*8        NO_Pil16Ra  ! [-2.4/1.9]? Slope of the N/O
*                                                  gradient using the calibrator
*                                                  Pil16_R in dex/Re
*                                                  (NO_Pil16_R_alpha_fit)
      real*4        e_NO_Pil16Ra ! [0.0029/2]? NO-Pil16Ra
*                                                  uncertainty
*                                                  (e_NO_Pil16_R_alpha_fit)
      real*8        NO_Pil16_HoRRe ! [-3.5/0.5]? N/O abundance
*                                                  using the calibrator
*                                                  Pil16_Ho_R at 1Re
*                                                  (NO_Pil16_Ho_R_Re_fit) (1)
      real*4        e_NO_Pil16_HoRRe ! [0.0028/2.3]? NO-Pil16-HoRRe
*                                                  uncertainty
*                                                  (e_NO_Pil16_Ho_R_Re_fit)
      real*8        NO_Pil16_HoRa ! [-2.9/2.7]? Slope of the N/O
*                                                  gradient using the calibrator
*                                                  Pil16_Ho_R in dex/Re
*                                                  (NO_Pil16_Ho_R_alpha_fit)
      real*4        e_NO_Pil16_HoRa ! [0.002/2.2]? NO-Pil16-HoRa
*                                                  uncertainty
*                                                  (e_NO_Pil16_Ho_R_alpha_fit)
      real*8        NO_Pil16_N2R2Re ! [-1.61/-0.27]? N/O abundance
*                                                  using the calibrator
*                                                  Pil16_N2_R2 at 1Re
*                                                  (NO_Pil16_N2_R2_Re_fit) (1)
      real*4        e_NO_Pil16_N2R2Re ! [0.0029/1.2]? NO-Pil16-N2R2Re
*                                                  uncertainty
*                                                  (e_NO_Pil16_N2_R2_Re_fit)
      real*8        NO_Pil16_N2R2a ! [-1.2/1.1]? Slope of the N/O
*                                                  gradient using the calibrator
*                                                  Pil16_N2_R2 in dex/Re
*                                                  (NO_Pil16_N2_R2_alpha_fit) (1)
      real*4        e_NO_Pil16_N2R2a ! [0.002/0.8]? NO-Pil16-N2R2a
*                                                  uncertainty
*                                                  (e_NO_Pil16_N2_R2_alpha_fit)
      real*4        Ne_Oster_SRe ! [1.39/3.6]? Electron density
*                                                  using the Oster_S estimator
*                                                  at 1Re (Ne_Oster_S_Re_fit)
      real*4        e_Ne_Oster_SRe ! [0.015/1]? Ne-Oster-SRe
*                                                  uncertainty
*                                                  (e_Ne_Oster_S_Re_fit)
      real*4        Ne_Oster_Sa ! [-0.98/1.6]? Slope of the n_e
*                                                  gradient using the Oster_S
*                                                  estimator in dex/Re
*                                                  (Ne_Oster_S_alpha_fit)
      real*4        e_Ne_Oster_Sa ! [0.013/0.7]? Ne-Oster-Sa
*                                                  uncertainty
*                                                  (e_Ne_Oster_S_alpha_fit)
      real*4        HdRe        ! (0.1nm) [-3.61/11]? Value of the Hd
*                                                  stellar index at 1Re in {AA}
*                                                  (Hd_Re_fit)
      real*4        e_HdRe      ! (0.1nm) [0.02/4.4]? HdRe uncertainty
*                                                  (e_Hd_Re_fit)
      real*8        Hda         ! [-32.3/17]? Slope of the
*                                                  gradient of the Hd index in
*                                                  {AA}/Re (Hd_alpha_fit)
      real*4        e_Hda       ! [0.03/6]? Hda uncertainty
*                                                  (e_Hd_alpha_fit)
      real*4        HbRe        ! (0.1nm) [-2.02/16.2]? Value of the
*                                                  H{beta} stellar index at 1Re
*                                                  in {AA} (Hb_Re_fit)
      real*4        e_HbRe      ! (0.1nm) [0.018/6]? HbRe uncertainty
*                                                  (e_Hb_Re_fit)
      real*8        Hba         ! [-7.3/14]? Slope of the
*                                                  gradient of the H{beta} index
*                                                  in {AA}/Re (Hb_alpha_fit)
      real*4        e_Hba       ! [0.015/7.5]? Hba uncertainty
*                                                  (e_Hb_alpha_fit)
      real*4        MgbRe       ! (0.1nm) [-2.8/16.2]? Value of the Mgb
*                                                  stellar index at 1Re in {AA}
*                                                  (Mgb_Re_fit)
      real*4        e_MgbRe     ! (0.1nm) [0.019/6]? MgbRe uncertainty
*                                                  (e_Mgb_Re_fit)
      real*4        Mgba        ! [-8.8/17]? Slope of the
*                                                  gradient of the Mgb index in
*                                                  {AA}/Re (Mgb_alpha_fit)
      real*4        e_Mgba      ! [0.016/10]? Mgba uncertainty
*                                                  (e_Mgb_alpha_fit)
      real*4        Fe5270Re    ! (0.1nm) [-3.7/22.2]? Value of the
*                                                  Fe5270 stellar index at 1Re
*                                                  in {AA} (Fe5270_Re_fit)
      real*4        e_Fe5270Re  ! (0.1nm) [0.017/5]? Fe5270Re
*                                                  uncertainty (e_Fe5270_Re_fit)
      real*4        Fe5270a     ! [-9.3/32]? Slope of the
*                                                  gradient of the Fe5270 index
*                                                  in {AA}/Re (Fe5270_alpha_fit)
      real*4        e_Fe5270a   ! [0.015/9]? Fe5270a
*                                                  uncertainty
*                                                  (e_Fe5270_alpha_fit)
      real*4        Fe5335Re    ! (0.1nm) [-3.8/21]? Value of the
*                                                  Fe5335 stellar index at 1Re
*                                                  in {AA} (Fe5335_Re_fit)
      real*4        e_Fe5335Re  ! (0.1nm) [0.016/18]? Fe5335Re
*                                                  uncertainty (e_Fe5335_Re_fit)
      real*8        Fe5335a     ! [-12/24]? Slope of the
*                                                  gradient of the Fe5335 index
*                                                  in {AA}/Re (Fe5335_alpha_fit)
      real*4        e_Fe5335a   ! [0.013/12]? Fe5335a
*                                                  uncertainty
*                                                  (e_Fe5335_alpha_fit)
      real*4        D4000Re     ! [0.57/4.6]? Value of the
*                                                  D4000 stellar index at 1Re
*                                                  (D4000_Re_fit1)
      real*4        e_D4000Re   ! [0.002/1.3]? D4000Re
*                                                  uncertainty (e_D4000_Re_fit)
      real*8        D4000a      ! [-0.8/6]? Slope of the
*                                                  gradient of the D4000 index
*                                                  (D4000_alpha_fit)
      real*4        e_D4000a    ! [0.002/1.5]? D4000a
*                                                  uncertainty
*                                                  (e_D4000_alpha_fit)
      real*4        HdmodRe     ! (0.1nm) [-3.2/14]? Value of the Hdmod
*                                                  stellar index at 1Re in {AA}
*                                                  (Hdmod_Re_fit)
      real*4        e_HdmodRe   ! (0.1nm) [0.05/6]? HdmodRe uncertainty
*                                                  (e_Hdmod_Re_fit)
      real*8        Hdmoda      ! [-27/24]? Slope of the
*                                                  gradient of the Hdmod index
*                                                  in {AA}/Re (Hdmod_alpha_fit)
      real*4        e_Hdmoda    ! [0.04/10]? Hdmoda uncertainty
*                                                  (e_Hdmod_alpha_fit)
      real*4        HgRe        ! (0.1nm) [-5/11]? Value of the Hg
*                                                  stellar index at 1Re in {AA}
*                                                  (Hg_Re_fit)
      real*4        e_HgRe      ! (0.1nm) [0.03/6]? HgRe uncertainty
*                                                  (e_Hg_Re_fit)
      real*4        Hga         ! [-5.4/10]? Slope of the
*                                                  gradient of the Hg index in
*                                                  {AA}/Re (Hg_alpha_fit)
      real*4        e_Hga       ! [0.027/6]? Hga uncertainty
*                                                  (e_Hg_alpha_fit)
      real*4        umag        ! (mag) [14.48/23.15] u-band
*                                                  magnitude derived from the
*                                                  original cube (u_band_mag)
      real*4        e_umag      ! (mag) [0.05/6] Uncertainty on the
*                                                  umag (u_band_mag_error)
      real*8        uMag_1      ! (mag) [-24.64/-7.9] u-band absolute
*                                                  magnitude derived from the
*                                                  original cube (u_band_abs_mag)
      real*4        e_uMag_1    ! (mag) [0.05/6] uMag uncertainty
*                                                  (u_band_abs_mag_error)
      real*4        gmag        ! (mag) [13.41/20.83] g-band
*                                                  magnitude derived from the
*                                                  original cube (g_band_mag)
      real*4        e_gmag      ! (mag) [0.05/1.1] Uncertainty on the
*                                                  gmag (g_band_mag_error)
      real*8        gMag_1      ! (mag) [-25.15/-8.61] g-band absolute
*                                                  magnitude derived from the
*                                                  original cube (g_band_abs_mag)
      real*4        e_gMag_1    ! (mag) [0.05/1.1] gMag uncertainty
*                                                  (g_band_abs_mag_error)
      real*4        rmag        ! (mag) [12.75/20.4] r-band magnitude
*                                                  derived from the original
*                                                  cube (r_band_mag)
      real*4        e_rmag      ! (mag) [0.05/1.2] Uncertainty on the
*                                                  rmag (r_band_mag_error)
      real*8        rMag_1      ! (mag) [-25.3/-8.7] r-band absolute
*                                                  magnitude derived from the
*                                                  original cube (r_band_abs_mag)
      real*4        e_rMag_1    ! (mag) [0.05/1.2] rMag uncertainty
*                                                  (r_band_abs_mag_error)
      real*4        imag        ! (mag) [12.32/21.76] i-band
*                                                  magnitude derived from the
*                                                  original cube (i_band_mag)
      real*4        e_imag      ! (mag) [0.05/17.1] Uncertainty on
*                                                  the imag (i_band_mag_error)
      real*8        iMag_1      ! (mag) [-25.4/-8.7] i-band absolute
*                                                  magnitude derived from the
*                                                  original cube (iMag)
      real*4        e_iMag_1    ! (mag) [0.05/17.1] iMag uncertainty
*                                                  (iMag_error)
      real*4        Bmag        ! (mag) [13.6/20.6] B-band magnitude
*                                                  derived from the original
*                                                  cube (B_band_mag)
      real*4        e_Bmag      ! (mag) [0.05/0.8] Uncertainty on the
*                                                  Bmag (B_band_mag_error)
      real*8        BMag_1      ! (mag) [-25.04/-8.5] B-band absolute
*                                                  magnitude derived from the
*                                                  original cube (B_band_abs_mag)
      real*4        e_BMag_1    ! (mag) [0.05/0.8] BMag uncertainty
*                                                  (B_band_abs_mag_error)
      real*4        Vmag        ! (mag) [13.12/21.3] V-band magnitude
*                                                  derived from the original
*                                                  cube (V_band_mag)
      real*4        e_Vmag      ! (mag) [0.05/2] Uncertainty on the
*                                                  Vmag (V_band_mag_error)
      real*8        VMag_1      ! (mag) [-25.25/-8.6] V-band absolute
*                                                  magnitude derived from the
*                                                  original cube (V_band_abs_mag)
      real*4        e_VMag_1    ! (mag) [0.05/2] VMag uncertainty
*                                                  (V_band_abs_mag_error)
      real*4        RJmag       ! (mag) [12.66/20.35] R-band
*                                                  magnitude derived from the
*                                                  original cube (RJ_band_mag)
      real*4        e_RJmag     ! (mag) [0.05/1.3] Uncertainty on the
*                                                  RJmag (RJ_band_mag_error)
      real*8        RJMag_1     ! (mag) [-25.3/-8.7] R-band absolute
*                                                  magnitude derived from the
*                                                  original cube
*                                                  (RJ_band_abs_mag)
      real*4        e_RJMag_1   ! (mag) [0.05/1.3] RJMag uncertainty
*                                                  (RJ_band_abs_mag_error)
      real*8        R50         ! (arcsec) [1.9/23] Radius at which half
*                                                  of the light within the FOV
*                                                  in the g-band is integrated
      real*4        e_R50       ! (arcsec) [2.7e-5/9.6] R50 uncertainty
*                                                  (error_R50)
      real*8        R90         ! (arcsec) [3.8/30.2] Radius at which
*                                                  90% of the light within the
*                                                  FOV in the g-band is
*                                                  integrated
      real*4        e_R90       ! (arcsec) [6e-7/2] R90 uncertainty
*                                                  (error_R90)
      real*8        C           ! [1/6] R90/R50 concentration
*                                                  index
      real*4        e_C         ! [5.9e-5/3.2] C uncertainty
      real*4        B_V         ! (mag) [-1/1.5] B-V color index
      real*4        e_B_V       ! (mag) [0.07/2.1] B-V uncertainty
*                                                  (error_B-V)
      real*4        B_R         ! (mag) [-0.28/2.2] B-R color index
      real*4        e_B_R       ! (mag) [0.07/1.6] B-R uncertainty
*                                                  (error_B-R)
      real*4        logMph      ! ([Msun]) [4.7/12] Stellar masses
*                                                  derived from the photometric
*                                                  data within the FOV
*                                                  (log_Mass_phot)
      real*4        e_logMph    ! ([Msun]) [0.13/4] logMph uncertainty
*                                                  (e_log_Mass_phot)
      real*4        SBV1Re      ! (mag/arcsec2) [17.8/38.9] V-band surface
*                                                  brightness at 1Re
*                                                  (V-band_SB_at_Re)
      real*4        e_SBV1Re    ! (mag/arcsec2) [0/9] SBV1Re uncertainty
*                                                  (error_V-band_SB_at_Re)
      real*4        SBVR50      ! (mag/arcsec2) [18/39] V-band surface
*                                                  brightness at R50
*                                                  (V-band_SB_at_R_50)
      real*4        e_SBVR50    ! (mag/arcsec2) [0/10] SBVR50 uncertainty
*                                                  (error_V-band_SB_at_R_50)
      real*8        nNSA        ! [0.5/6] NSA Sersic index
*                                                  (nsa_sersic_n_morph)
      real*8        u_g         ! (mag) [-3.2/9.2] u-g NSA color index
      real*8        g_r         ! (mag) [-0.94/2.5] g-r NSA color index
      real*8        r_i         ! (mag) [-2.2/2] r-i NSA color index
      real*8        i_z         ! (mag) [-7/2.8] i-z NSA color index
      real*8        P_CD        ! [0/1] Probability of being a
*                                                  cD galaxy
      real*8        P_E         ! [0/0.97] Probability of being
*                                                  an E galaxy
      real*8        P_S0        ! [0.0001/0.94] Probability of
*                                                  being an S0 galaxy
      real*8        P_Sa        ! [0.0002/0.92] Probability of
*                                                  being an Sa galaxy
      real*8        P_Sab       ! [0.0001/0.8] Probability of
*                                                  being an Sab galaxy
      real*8        P_Sb        ! [0.0005/0.9] Probability of
*                                                  being an Sb galaxy
      real*8        P_Sbc       ! [0.0003/0.85] Probability of
*                                                  being an Sbc galaxy
      real*8        P_Sc        ! [0.0004/0.92] Probability of
*                                                  being an Sc galaxy
      real*8        P_Scd       ! [0/0.93] Probability of being
*                                                  an Scd galaxy
      real*8        P_Sd        ! [0/0.94] Probability of being
*                                                  an Sd galaxy
      real*8        P_Sdm       ! [0/0.99] Probability of being
*                                                  an Sdm galaxy
      real*8        P_Sm        ! [0/0.94] Probability of being
*                                                  an Sm galaxy
      real*8        P_Irr       ! [0/0.96] Probability of being
*                                                  an Irr galaxy
      integer*4     Type        ! [-2/10]? Best morphological
*                                                  type index based on the NN
*                                                  analysis (best_type_n)
      character*3   Morph       ! [nan] Morphological type
*                                                  derived by the NN analysis;
*                                                  always "nan" (best_type)
      integer*4     NSAid       ! [58/698275]? NSA identifier
*                                                  (nsa_nsaid)
      real*8        w_Vmax      ! [0/0.0012]? Weight for the
*                                                  volume correction in volume
*                                                  in Mpc^-3^/dex (Vmax_w)
      real*8        w_Num       ! [0/13314]? Weight for the
*                                                  volume correction in number
*                                                  (Num_w)
      integer*4     Q           ! [0/7]? QC flag: 0=good,
*                                                  2=bad, >2 warning (QCFLAG)
*Note (1): List of oxygen, nitrogen, and ionization parameter calibrators
*          adopted to estimate the values at the effective radius (see Table 15):
* ------------------------------------------------------------------------------
* ID at Re                 Emission lines                Reference
*    (central aperture)     and ratios
* ------------------------------------------------------------------------------
* OH_Mar13_N2 (OH_N2)     [NII]6548/H{alpha}             Marino+ (2013,
*                                                         J/A+A/559/A114)
* OH_Mar13_O3N2 (OH_O3N2) O3N2                           Marino+ (2013,
*                                                         J/A+A/559/A114)
* OH_T04                  [NII]6548/H{alpha}, R_23_      Tremonti+
*                                                         2004ApJ...613..898T
* OH_Pet04_N2_lin         [NII]6548/H{alpha}             Pettini & Pagel
*                                                         2004MNRAS.348L..59P
* OH_Pet04_N2_poly        [NII]6548/H{alpha}             Pettini & Pagel
*                                                         2004MNRAS.348L..59P
* OH_Pet04_O3N2           O3N2                           Pettini & Pagel
*                                                         2004MNRAS.348L..59P
* OH_Kew02_N2O2           [NII]6583/[OII]3727            Kewley & Dopita
*                                                         2002ApJS..142...35K
* OH_Pil10_ONS (OH_ONS)   [NII]6583/H{beta}, R_2_, R_3_, Pilyugin et al.
*                         P, [SII]/H{alpha}               2010ApJ...720.1738P
* OH_Pil10_ON             [NII]6583/H{beta}, R_2_, R_3_, Pilyugin et al.
*                         [SII]/H{alpha}                  2010ApJ...720.1738P
* OH_Pil11_NS             [NII]6583/H{beta}, R_3_,       Pilyugin & Mattsson
*                         [SII]/H{alpha}                  2011MNRAS.412.1145P
* OH_Cur20_RS32           R_3_+[SII]/H{alpha}/H{alpha}   Curti et al.
*                                                         2020MNRAS.491..944C
* OH_Cur20_R3             R_3_                           Curti et al.
*                                                         2020MNRAS.491..944C
* OH_Cur20_O3O2           [OIII]5007, [OII]3727 + 29     Curti et al.
*                                                         2020MNRAS.491..944C
* OH_Cur20_S2             [SII]/H{alpha}                 Curti et al.
*                                                         2020MNRAS.491..944C
* OH_Cur20_R2             R_2_                           Curti et al.
*                                                         2020MNRAS.491..944C
* OH_Cur20_N2             [NII]6583/H{alpha}             Curti et al.
*                                                         2020MNRAS.491..944C
* OH_Cur20_R23            R_23_                          Curti et al.
*                                                         2020MNRAS.491..944C
* OH_Cur20_O3N2           R_3_, [NII]6583/H{alpha}       Curti et al.
*                                                         2020MNRAS.491..944C
* OH_Cur20_O3S2           R_3_/[SII]/H{alpha}/H{alpha}   Curti et al.
*                                                         2020MNRAS.491..944C
* OH_KK04 (OH_R23)        R_23_, [OIII]/[OII]            Kobulnicky & Kewley
*                                                         (2004, J/ApJ/617/240)
* OH_Pil16_R              [NII]6583/H{beta}, R_2_, R_3_  Pilyugin & Grebel
*                                                         2016MNRAS.457.3678P
* OH_Pil16_S              [NII]6583/H{beta}, R_3_,       Pilyugin & Grebel
*                         [SII]/H{alpha}                  2016MNRAS.457.3678P
* OH_Ho                   R_2_, R_3_, [NII]6583/H{beta}, Ho I.-T.
*                         [SII]/H{alpha}/H{alpha}         2019MNRAS.485.3569H
*
* U_Dors_O32              [OIII]5007, [OII]3727 + 29     Dors et al.
*                                                         2011MNRAS.415.3616D
* U_Dors_S2               [NII]6583/H{alpha},            Dors et al.
*                         [SII]/H{alpha}                  2011MNRAS.415.3616D
* U_Mor16_O32_fs          [OIII]5007, [OII]3727 + 29     Morisset et al.
*                                                         2016A&A...594A..37M
* U_Mor16_O32_ts          [OIII]5007, [OII]3727 + 29     Morisset et al.
*                                                         2016A&A...594A..37M
*
* NH_Pil16_R              [NII]6583/H{beta}, R_2_, R_3_  Pilyugin & Grebel
*                                                         2016MNRAS.457.3678P
* NO_Pil16_R              [NII]6583/H{beta}, R_2_, R_3_  Pilyugin & Grebel
*                                                         2016MNRAS.457.3678P
* NO_Pil16_Ho_R           [NII]6583/H{beta}, R_2_, R_3_, Pilyugin & Grebel
*                         [SII]/H{alpha}/H{alpha}        2016MNRAS.457.3678P and
*                                                        Ho (2019MNRAS.485.3569H)
* NO_Pil16_N2_R2          [NII]6583/H{beta}, R_2_, R_3_  Pilyugin & Grebel
*                                                         2016MNRAS.457.3678P
* ------------------------------------------------------------------------------

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

C  Loading file 'catalog.dat'	! Integrated and characteristic parameters delivered
*                               for each analyzed data cube

C  Format for file interpretation

    1 format(
     +  A17,1X,I5,1X,A5,1X,A11,1X,A9,1X,F11.7,1X,F10.7,1X,F7.4,1X,
     +  F8.5,1X,F10.3,1X,F6.4,1X,F6.4,1X,F7.4,1X,F7.4,1X,F9.6,1X,
     +  E10.5,1X,F8.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,F8.3,
     +  1X,F7.3,1X,F7.4,1X,F6.4,1X,F7.4,1X,F6.4,1X,F7.4,1X,F6.4,1X,
     +  F7.4,1X,F6.4,1X,F7.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,F8.5,1X,F7.5,
     +  1X,F8.5,1X,F7.5,1X,F9.3,1X,E9.4,1X,F7.5,1X,F6.2,1X,F7.5,1X,
     +  F8.4,1X,F6.3,1X,F6.3,1X,F7.4,1X,F7.4,1X,F8.6,1X,E11.6,1X,F8.5,
     +  1X,F7.5,1X,F7.5,1X,E10.5,1X,F9.5,1X,F8.5,1X,F7.5,1X,F7.5,1X,
     +  F7.5,1X,F7.5,1X,F8.5,1X,F8.5,1X,F9.5,1X,F8.5,1X,F11.5,1X,
     +  F11.5,1X,F8.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,F10.3,1X,F9.3,1X,
     +  F8.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,F8.5,1X,F7.5,
     +  1X,F8.5,1X,F7.5,1X,F8.3,1X,F7.3,1X,F7.3,1X,F7.3,1X,F8.5,1X,
     +  F7.5,1X,F8.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,F8.5,
     +  1X,F7.5,1X,F8.3,1X,F7.3,1X,F6.3,1X,F7.5,1X,F7.5,1X,F7.5,1X,
     +  F7.5,1X,F7.5,1X,F7.5,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,F7.4,1X,F7.4,1X,F6.4,1X,F7.4,1X,F7.4,
     +  1X,F7.4,1X,F6.3,1X,F6.3,1X,F8.3,1X,F7.3,1X,F8.3,1X,F7.3,1X,
     +  F8.3,1X,F7.3,1X,F8.3,1X,F7.3,1X,F7.4,1X,F7.4,1X,F7.3,1X,F7.3,
     +  1X,F7.3,1X,F7.3,1X,F6.3,1X,E9.4,1X,F9.5,1X,F10.5,1X,F10.5,1X,
     +  F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F7.4,1X,F6.4,1X,F7.4,1X,E10.5,
     +  1X,F7.5,1X,E10.5,1X,F7.5,1X,F7.5,1X,E10.5,1X,F8.6,1X,F7.4,1X,
     +  F7.4,1X,F8.5,1X,F7.5,1X,F9.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,F9.5,
     +  1X,F7.5,1X,F8.5,1X,F7.5,1X,F9.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,
     +  F8.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,F8.5,1X,F7.5,
     +  1X,F8.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,F8.5,1X,
     +  F7.5,1X,F8.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,F8.5,
     +  1X,F7.5,1X,F8.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,
     +  F8.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,F8.5,1X,F7.5,
     +  1X,F8.5,1X,F8.5,1X,F8.5,1X,F8.5,1X,F8.5,1X,F7.5,1X,F9.5,1X,
     +  F7.5,1X,F8.5,1X,F7.5,1X,F9.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,F8.5,
     +  1X,F7.5,1X,F8.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,
     +  F8.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,F8.5,1X,F7.5,
     +  1X,F8.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,F8.5,1X,
     +  F7.5,1X,F9.5,1X,F7.5,1X,F8.5,1X,F8.5,1X,F9.5,1X,F8.5,1X,F8.5,
     +  1X,F7.5,1X,F9.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,
     +  F8.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,F8.5,1X,F7.5,
     +  1X,F8.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,F8.5,1X,
     +  E10.5,1X,F8.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,F8.5,
     +  1X,F7.5,1X,F6.4,1X,F6.4,1X,F8.5,1X,F7.5,1X,F7.5,1X,F7.5,1X,
     +  F8.5,1X,F7.5,1X,F7.4,1X,F6.4,1X,F7.4,1X,F6.4,1X,F6.4,1X,F6.4,
     +  1X,F8.5,1X,F7.5,1X,F6.4,1X,F6.4,1X,F7.4,1X,F6.4,1X,F6.4,1X,
     +  F6.4,1X,F8.5,1X,F7.5,1X,F6.4,1X,F6.4,1X,F7.4,1X,F6.4,1X,F7.5,
     +  1X,F7.5,1X,F8.5,1X,F7.5,1X,F6.4,1X,F6.4,1X,F7.4,1X,F6.4,1X,
     +  F6.4,1X,F6.4,1X,F8.5,1X,F7.5,1X,F6.4,1X,F6.4,1X,F8.5,1X,F7.5,
     +  1X,F6.4,1X,F6.4,1X,F8.5,1X,F7.5,1X,F6.4,1X,F6.4,1X,F7.4,1X,
     +  F6.4,1X,F6.4,1X,F6.4,1X,F7.4,1X,F6.4,1X,F6.4,1X,F6.4,1X,F7.4,
     +  1X,F6.4,1X,F6.4,1X,F6.4,1X,F7.4,1X,F6.4,1X,F6.4,1X,F6.4,1X,
     +  F7.4,1X,F6.4,1X,F7.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,F6.4,1X,F6.4,
     +  1X,F7.4,1X,F6.4,1X,F7.5,1X,F7.5,1X,F8.5,1X,F7.5,1X,F6.4,1X,
     +  F6.4,1X,F7.4,1X,F6.4,1X,F6.4,1X,F6.4,1X,F7.4,1X,F6.4,1X,F6.4,
     +  1X,F6.4,1X,F7.4,1X,F6.4,1X,F7.4,1X,F6.4,1X,F7.4,1X,F6.4,1X,
     +  F7.4,1X,F6.4,1X,F7.4,1X,F6.4,1X,F7.4,1X,F6.4,1X,F7.4,1X,F6.4,
     +  1X,F7.4,1X,F6.4,1X,F7.4,1X,F6.4,1X,F6.4,1X,F6.4,1X,F7.4,1X,
     +  F6.4,1X,F7.4,1X,F6.4,1X,F7.4,1X,F6.4,1X,F7.4,1X,F6.4,1X,F7.4,
     +  1X,F6.4,1X,F7.4,1X,F6.4,1X,F7.4,1X,F6.4,1X,F5.3,1X,F5.3,1X,
     +  F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F7.3,1X,F5.3,1X,F6.3,1X,F5.3,
     +  1X,F7.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,
     +  F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F6.3,1X,F7.3,1X,F6.3,1X,F6.4,
     +  1X,F6.4,1X,F7.4,1X,F6.4,1X,F6.3,1X,F5.3,1X,F7.3,1X,F5.3,1X,
     +  F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F7.3,1X,F5.3,
     +  1X,F6.3,1X,F5.3,1X,F7.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F7.3,1X,
     +  F5.3,1X,F6.3,1X,F6.3,1X,F7.3,1X,F6.3,1X,F6.3,1X,F5.3,1X,F7.3,
     +  1X,F5.3,1X,F6.3,1X,F5.3,1X,F7.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,
     +  F7.3,1X,F5.3,1X,F8.5,1X,E10.5,1X,F8.5,1X,E10.5,1X,F7.5,1X,
     +  E10.5,1X,F6.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F5.2,1X,F4.2,1X,F6.3,
     +  1X,F5.3,1X,F6.3,1X,F5.3,1X,F8.6,1X,F9.6,1X,F9.6,1X,F9.6,1X,
     +  F9.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,
     +  1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,F8.6,1X,I2,1X,A3,
     +  1X,I6,1X,F8.6,1X,F8.2,1X,I1)

C  Effective file loading

      open(unit=1,status='old',file=
     +'catalog.dat')
      write(6,*) '....Loading file: catalog.dat'
      do i__=1,10220
        read(1,'(A4325)')ar__
        read(ar__,1)
     +  Name,Plate,IFUdsgn,PlateIFU,MaNGAId,RAdeg,DEdeg,logSFRHa,FoV,
     +  Reff,e_logMass,e_logSFRHa,logMass,logSFRs,logNII_Hac,
     +  e_logNII_Hac,logOIII_Hbc,e_logOIII_Hbc,logSII_Hac,
     +  e_logSII_Hac,logOII_Hbc,e_logOII_Hbc,EWHac,e_EWHac,logZReLW,
     +  e_logZReLW,alphaZReLW,e_alphaZReLW,logZReMW,e_logZReMW,
     +  alphaZReMW,e_alphaZReMW,logAgeReLW,e_logAgeReLW,alphaAgeReLW,
     +  e_alphaAgeReLW,logAgeReMW,e_logAgeReMW,alphaAgeReMW,
     +  e_alphaAgeReMW,ReffA,DL,DA,PA,ell,logMassgas,vel_sig,
     +  e_vel_sig,logSFRsf,logSFRdc,OH_O3N2c,e_OH_O3N2c,OH_N2c,
     +  e_OH_N2c,OH_ONSc,e_OH_ONSc,OH_R23c,e_OH_R23c,OH_pyqzc,
     +  e_OH_pyqzc,OH_t2c,e_OH_t2c,OH_M08c,e_OH_M08c,OH_T04c,
     +  e_OH_T04c,OH_dopc,e_OH_dopc,OH_O3N2_EPc,e_OH_O3N2_EPc,
     +  logOI_Hac,e_logOI_Hac,Ha_Hbc,e_Ha_Hbc,logNII_HaRe,
     +  e_logNII_HaRe,logOIII_HbRe,e_logOIII_HbRe,logSII_HaRe,
     +  e_logSII_HaRe,logOII_HbRe,e_logOII_HbRe,logOI_HaRe,
     +  e_logOI_HaRe,EWHaRe,e_EWHaRe,Ha_HbRe,e_Ha_HbRe,logNII_Ha,
     +  e_logNII_Ha,logOIII_Hb,e_logOIII_Hb,logSII_Ha,e_logSII_Ha,
     +  logOII_Hb,e_logOII_Hb,logOI_Ha,e_logOI_Ha,EWHa,e_EWHa,Ha_Hb,
     +  logsigMc,e_logsigMc,logsigMRe,e_logsigMRe,logsigM,e_logsigM,
     +  T30,Z_T30,Z_T30Re,aZ_T30,T40,Z_T40,Z_T40Re,aZ_T40,T50,Z_T50,
     +  Z_T50Re,aZ_T50,T60,Z_T60,Z_T60Re,aZ_T60,T70,Z_T70,Z_T70Re,
     +  aZH_T70,T80,Z_T80,Z_T80Re,aZ_T80,T90,Z_T90,Z_T90Re,aZ_T90,T95,
     +  Z_T95,Z_T95Re,aZ_T95,T99,Z_T99,Z_T99Re,aZ_T99,logMgasAvg,
     +  logMgasAvs,Vels2Re,e_Vels2Re,VelHa2Re,e_VelHa2Re,Vels1Re,
     +  e_Vels1Re,VelHa1Re,e_VelHa1Re,logSFRs100,logSFRs10,VdispHac,
     +  Vdispsc,VdispHa1re,Vdisps1Re,logM1Re,MLint,MLavg,FHac,e_FHac,
     +  R50V,e_R50V,R50M,e_R50M,logMR50V,e_logMR50V,logMgasAvgl,
     +  Avgas1Re,e_Avgas1Re,Avs,e_Avs,Lambda1Re,e_Lambda1Re,zNSA,MNSA,
     +  IncNSA,F_OII_26Re,e_F_OII_26Re,F_OII_26a,e_F_OII_26a,
     +  F_OII_28Re,e_F_OII_28Re,F_OII_28a,e_F_OII_28a,FHI34Re,
     +  e_FHI34Re,FHI34a,e_FHI34a,FHI97Re,e_FHI97Re,FHI97a,e_FHI97a,
     +  FHeI38Re,e_FHeI38Re,FHeI38a,e_FHeI38a,FHI89Re,e_FHI89Re,
     +  FHI89a,e_FHI89a,FHeI39Re,e_FHeI39Re,FHeI39a,e_FHeI39a,
     +  F_NeIII_Re,e_F_NeIII_Re,F_NeIII_a,e_F_NeIII_a,FCaIIRe,
     +  e_FCaIIRe,FCaIIa,e_FCaIIa,FHepsRe,e_FHepsRe,FHepsa,e_FHepsa,
     +  FHdRe,e_FHdRe,FHda,e_FHda,FHgRe,e_FHgRe,FHga,e_FHga,FHbRe,
     +  e_FHbRe,FHba,e_FHba,F_OIII_49Re,e_F_OIII_49Re,F_OIII_49a,
     +  e_F_OIII_49a,F_OIII_50Re,e_F_OIII_50Re,F_OIII_50a,
     +  e_F_OIII_50a,FHeI50Re,e_FHeI50Re,FHeI50a,e_FHeI50a,F_NI_51Re,
     +  e_F_NI_51Re,F_NI_51a,e_F_NI_51a,F_NI_52Re,e_F_NI_52Re,
     +  F_NI_52a,e_F_NI_52a,FHeI58Re,e_FHeI58Re,FHeI58a,e_FHeI58a,
     +  FNaI89Re,e_FNaI89Re,FNaI89a,e_FNaI89a,FNaI95Re,e_FNaI95Re,
     +  FNaI95a,e_FNaI95a,F_OI_Re,e_F_OI_Re,F_OI_a,e_F_OI_a,
     +  F_NII_48Re,e_F_NII_48Re,F_NII_48a,e_F_NII_48a,FHaRe,e_FHaRe,
     +  FHaa,e_FHaa,F_NII_83Re,e_F_NII_83Re,F_NII_83a,e_F_NII_83a,
     +  F_SII_16Re,e_F_SII_16Re,F_SII_16a,e_F_SII_16a,F_SII_30Re,
     +  e_F_SII_30Re,F_SII_30a,e_F_SII_30a,F_ArIII_Re,e_F_ArIII_Re,
     +  F_ArIII_a,e_F_ArIII_a,FHI14Re,e_FHI14Re,FHI14a,e_FHI14a,
     +  F_SIII_90Re,e_F_SIII_90Re,F_SIII_90a,e_F_SIII_90a,F_FeII_Re,
     +  e_F_FeII_Re,F_FeII_a,e_F_FeII_a,F_SIII_95Re,e_F_SIII_95Re,
     +  F_SIII_95a,e_F_SIII_95a,OH_Mar13_N2Re,e_OH_Mar13_N2Re,
     +  OH_Mar13_N2a,e_OH_Mar13_N2a,OH_Mar13_O3N2Re,e_OH_Mar13_O3N2Re,
     +  OH_Mar13_O3N2a,e_OH_Mar13_O3N2a,OH_T04Re,e_OH_T04Re,OH_T04a,
     +  e_OH_T04a,OH_Pet04_N2linRe,e_OH_Pet04_N2linRe,OH_Pet04_N2lina,
     +  e_OH_Pet04_N2lina,OH_Pet04_N2polRe,e_OH_Pet04_N2polRe,
     +  OH_Pet04_N2pola,e_OH_Pet04_N2pola,OH_Pet04_O3N2Re,
     +  e_OH_Pet04_O3N2Re,OH_Pet04_O3N2a,e_OH_Pet04_O3N2a,
     +  OH_Kew02_N2O2Re,e_OH_Kew02_N2O2Re,OH_Kew02_N2O2a,
     +  e_OH_Kew02_N2O2a,OH_Pil10_ONSRe,e_OH_Pil10_ONSRe,
     +  OH_Pil10_ONSa,e_OH_Pil10_ONSa,OH_Pil10_ONRe,e_OH_Pil10_ONRe,
     +  OH_Pil10_ONa,e_OH_Pil10_ONa,OH_Pil11_NSRe,e_OH_Pil11_NSRe,
     +  OH_Pil11_NSa,e_OH_Pil11_NSa,OH_Cur20_RS32Re,e_OH_Cur20_RS32Re,
     +  OH_Cur20_RS32a,e_OH_Cur20_RS32a,OH_Cur20_R3Re,e_OH_Cur20_R3Re,
     +  OH_Cur20_R3a,e_OH_Cur20_R3a,OH_Cur20_O3O2Re,e_OH_Cur20_O3O2Re,
     +  OH_Cur20_O3O2a,e_OH_Cur20_O3O2a,OH_Cur20_S2Re,e_OH_Cur20_S2Re,
     +  OH_Cur20_S2a,e_OH_Cur20_S2a,OH_Cur20_R2Re,e_OH_Cur20_R2Re,
     +  OH_Cur20_R2a,e_OH_Cur20_R2a,OH_Cur20_N2Re,e_OH_Cur20_N2Re,
     +  OH_Cur20_N2a,e_OH_Cur20_N2a,OH_Cur20_R23Re,e_OH_Cur20_R23Re,
     +  OH_Cur20_R23a,e_OH_Cur20_R23a,OH_Cur20_O3N2Re,
     +  e_OH_Cur20_O3N2Re,OH_Cur20_O3N2a,e_OH_Cur20_O3N2a,
     +  OH_Cur20_O3S2Re,e_OH_Cur20_O3S2Re,OH_Cur20_O3S2a,
     +  e_OH_Cur20_O3S2a,OH_KK04Re,e_OH_KK04Re,OH_KK04a,e_OH_KK04a,
     +  OH_Pil16_RRe,e_OH_Pil16_RRe,OH_Pil16_Ra,e_OH_Pil16_Ra,
     +  OH_Pil16_SRe,e_OH_Pil16_SRe,OH_Pil16_Sa,e_OH_Pil16_Sa,OH_HoRe,
     +  e_OH_HoRe,OH_Hoa,e_OH_Hoa,U_Dors_O32Re,e_U_Dors_O32Re,
     +  U_Dors_O32a,e_U_Dors_O32a,U_Dors_SRe,e_U_Dors_SRe,U_Dors_Sa,
     +  e_U_Dors_Sa,U_Mor16_O23fsRe,e_U_Mor16_O23fsRe,U_Mor16_O23fsa,
     +  e_U_Mor16_O23fsa,U_Mor16_O23tsRe,e_U_Mor16_O23tsRe,
     +  U_Mor16_O23tsa,e_U_Mor16_O23tsa,NH_Pil16RRe,e_NH_Pil16RRe,
     +  NH_Pil16Ra,e_NH_Pil16Ra,NO_Pil16RRe,e_NO_Pil16RRe,NO_Pil16Ra,
     +  e_NO_Pil16Ra,NO_Pil16_HoRRe,e_NO_Pil16_HoRRe,NO_Pil16_HoRa,
     +  e_NO_Pil16_HoRa,NO_Pil16_N2R2Re,e_NO_Pil16_N2R2Re,
     +  NO_Pil16_N2R2a,e_NO_Pil16_N2R2a,Ne_Oster_SRe,e_Ne_Oster_SRe,
     +  Ne_Oster_Sa,e_Ne_Oster_Sa,HdRe,e_HdRe,Hda,e_Hda,HbRe,e_HbRe,
     +  Hba,e_Hba,MgbRe,e_MgbRe,Mgba,e_Mgba,Fe5270Re,e_Fe5270Re,
     +  Fe5270a,e_Fe5270a,Fe5335Re,e_Fe5335Re,Fe5335a,e_Fe5335a,
     +  D4000Re,e_D4000Re,D4000a,e_D4000a,HdmodRe,e_HdmodRe,Hdmoda,
     +  e_Hdmoda,HgRe,e_HgRe,Hga,e_Hga,umag,e_umag,uMag_1,e_uMag_1,
     +  gmag,e_gmag,gMag_1,e_gMag_1,rmag,e_rmag,rMag_1,e_rMag_1,imag,
     +  e_imag,iMag_1,e_iMag_1,Bmag,e_Bmag,BMag_1,e_BMag_1,Vmag,
     +  e_Vmag,VMag_1,e_VMag_1,RJmag,e_RJmag,RJMag_1,e_RJMag_1,R50,
     +  e_R50,R90,e_R90,C,e_C,B_V,e_B_V,B_R,e_B_R,logMph,e_logMph,
     +  SBV1Re,e_SBV1Re,SBVR50,e_SBVR50,nNSA,u_g,g_r,r_i,i_z,P_CD,P_E,
     +  P_S0,P_Sa,P_Sab,P_Sb,P_Sbc,P_Sc,P_Scd,P_Sd,P_Sdm,P_Sm,P_Irr,
     +  Type,Morph,NSAid,w_Vmax,w_Num,Q
        if(ar__(19:23) .EQ. '') Plate = iNULL__
        if(ar__(76:82) .EQ. '') logSFRHa = rNULL__
        if(ar__(126:132) .EQ. '') logSFRs = rNULL__
        if(ar__(134:142) .EQ. '') logNII_Hac = rNULL__
        if(ar__(144:153) .EQ. '') e_logNII_Hac = rNULL__
        if(ar__(155:162) .EQ. '') logOIII_Hbc = rNULL__
        if(ar__(164:170) .EQ. '') e_logOIII_Hbc = rNULL__
        if(ar__(172:179) .EQ. '') logSII_Hac = rNULL__
        if(ar__(181:187) .EQ. '') e_logSII_Hac = rNULL__
        if(ar__(189:196) .EQ. '') logOII_Hbc = rNULL__
        if(ar__(198:204) .EQ. '') e_logOII_Hbc = rNULL__
        if(ar__(206:213) .EQ. '') EWHac = rNULL__
        if(ar__(215:221) .EQ. '') e_EWHac = rNULL__
        if(ar__(231:236) .EQ. '') e_logZReLW = rNULL__
        if(ar__(238:244) .EQ. '') alphaZReLW = rNULL__
        if(ar__(246:251) .EQ. '') e_alphaZReLW = rNULL__
        if(ar__(253:259) .EQ. '') logZReMW = rNULL__
        if(ar__(261:266) .EQ. '') e_logZReMW = rNULL__
        if(ar__(268:274) .EQ. '') alphaZReMW = rNULL__
        if(ar__(276:281) .EQ. '') e_alphaZReMW = rNULL__
        if(ar__(283:289) .EQ. '') logAgeReLW = rNULL__
        if(ar__(291:297) .EQ. '') e_logAgeReLW = rNULL__
        if(ar__(299:306) .EQ. '') alphaAgeReLW = rNULL__
        if(ar__(308:314) .EQ. '') e_alphaAgeReLW = rNULL__
        if(ar__(316:323) .EQ. '') logAgeReMW = rNULL__
        if(ar__(325:331) .EQ. '') e_logAgeReMW = rNULL__
        if(ar__(333:340) .EQ. '') alphaAgeReMW = rNULL__
        if(ar__(342:348) .EQ. '') e_alphaAgeReMW = rNULL__
        if(ar__(350:358) .EQ. '') ReffA = rNULL__
        if(ar__(393:400) .EQ. '') logMassgas = rNULL__
        if(ar__(402:407) .EQ. '') vel_sig = rNULL__
        if(ar__(409:414) .EQ. '') e_vel_sig = rNULL__
        if(ar__(416:422) .EQ. '') logSFRsf = rNULL__
        if(ar__(424:430) .EQ. '') logSFRdc = rNULL__
        if(ar__(432:439) .EQ. '') OH_O3N2c = rNULL__
        if(ar__(441:451) .EQ. '') e_OH_O3N2c = rNULL__
        if(ar__(453:460) .EQ. '') OH_N2c = rNULL__
        if(ar__(462:468) .EQ. '') e_OH_N2c = rNULL__
        if(ar__(470:476) .EQ. '') OH_ONSc = rNULL__
        if(ar__(478:487) .EQ. '') e_OH_ONSc = rNULL__
        if(ar__(489:497) .EQ. '') OH_R23c = rNULL__
        if(ar__(499:506) .EQ. '') e_OH_R23c = rNULL__
        if(ar__(508:514) .EQ. '') OH_pyqzc = rNULL__
        if(ar__(516:522) .EQ. '') e_OH_pyqzc = rNULL__
        if(ar__(524:530) .EQ. '') OH_t2c = rNULL__
        if(ar__(532:538) .EQ. '') e_OH_t2c = rNULL__
        if(ar__(540:547) .EQ. '') OH_M08c = rNULL__
        if(ar__(549:556) .EQ. '') e_OH_M08c = rNULL__
        if(ar__(558:566) .EQ. '') OH_T04c = rNULL__
        if(ar__(568:575) .EQ. '') e_OH_T04c = rNULL__
        if(ar__(577:587) .EQ. '') OH_dopc = rNULL__
        if(ar__(589:599) .EQ. '') e_OH_dopc = rNULL__
        if(ar__(601:608) .EQ. '') OH_O3N2_EPc = rNULL__
        if(ar__(610:616) .EQ. '') e_OH_O3N2_EPc = rNULL__
        if(ar__(618:625) .EQ. '') logOI_Hac = rNULL__
        if(ar__(627:633) .EQ. '') e_logOI_Hac = rNULL__
        if(ar__(635:644) .EQ. '') Ha_Hbc = rNULL__
        if(ar__(646:654) .EQ. '') e_Ha_Hbc = rNULL__
        if(ar__(656:663) .EQ. '') logNII_HaRe = rNULL__
        if(ar__(665:671) .EQ. '') e_logNII_HaRe = rNULL__
        if(ar__(673:680) .EQ. '') logOIII_HbRe = rNULL__
        if(ar__(682:688) .EQ. '') e_logOIII_HbRe = rNULL__
        if(ar__(690:697) .EQ. '') logSII_HaRe = rNULL__
        if(ar__(699:705) .EQ. '') e_logSII_HaRe = rNULL__
        if(ar__(707:714) .EQ. '') logOII_HbRe = rNULL__
        if(ar__(716:722) .EQ. '') e_logOII_HbRe = rNULL__
        if(ar__(724:731) .EQ. '') logOI_HaRe = rNULL__
        if(ar__(733:739) .EQ. '') e_logOI_HaRe = rNULL__
        if(ar__(741:748) .EQ. '') EWHaRe = rNULL__
        if(ar__(750:756) .EQ. '') e_EWHaRe = rNULL__
        if(ar__(758:764) .EQ. '') Ha_HbRe = rNULL__
        if(ar__(766:772) .EQ. '') e_Ha_HbRe = rNULL__
        if(ar__(774:781) .EQ. '') logNII_Ha = rNULL__
        if(ar__(783:789) .EQ. '') e_logNII_Ha = rNULL__
        if(ar__(791:798) .EQ. '') logOIII_Hb = rNULL__
        if(ar__(800:806) .EQ. '') e_logOIII_Hb = rNULL__
        if(ar__(808:815) .EQ. '') logSII_Ha = rNULL__
        if(ar__(817:823) .EQ. '') e_logSII_Ha = rNULL__
        if(ar__(825:832) .EQ. '') logOII_Hb = rNULL__
        if(ar__(834:840) .EQ. '') e_logOII_Hb = rNULL__
        if(ar__(842:849) .EQ. '') logOI_Ha = rNULL__
        if(ar__(851:857) .EQ. '') e_logOI_Ha = rNULL__
        if(ar__(859:866) .EQ. '') EWHa = rNULL__
        if(ar__(868:874) .EQ. '') e_EWHa = rNULL__
        if(ar__(876:881) .EQ. '') Ha_Hb = rNULL__
        if(ar__(883:889) .EQ. '') logsigMc = rNULL__
        if(ar__(891:897) .EQ. '') e_logsigMc = rNULL__
        if(ar__(899:905) .EQ. '') logsigMRe = rNULL__
        if(ar__(907:913) .EQ. '') e_logsigMRe = rNULL__
        if(ar__(915:921) .EQ. '') logsigM = rNULL__
        if(ar__(923:929) .EQ. '') e_logsigM = rNULL__
        if(ar__(947:953) .EQ. '') Z_T30Re = rNULL__
        if(ar__(955:961) .EQ. '') aZ_T30 = rNULL__
        if(ar__(979:985) .EQ. '') Z_T40Re = rNULL__
        if(ar__(987:993) .EQ. '') aZ_T40 = rNULL__
        if(ar__(1011:1017) .EQ. '') Z_T50Re = rNULL__
        if(ar__(1019:1025) .EQ. '') aZ_T50 = rNULL__
        if(ar__(1043:1049) .EQ. '') Z_T60Re = rNULL__
        if(ar__(1051:1057) .EQ. '') aZ_T60 = rNULL__
        if(ar__(1075:1081) .EQ. '') Z_T70Re = rNULL__
        if(ar__(1083:1089) .EQ. '') aZH_T70 = rNULL__
        if(ar__(1107:1113) .EQ. '') Z_T80Re = rNULL__
        if(ar__(1115:1121) .EQ. '') aZ_T80 = rNULL__
        if(ar__(1139:1145) .EQ. '') Z_T90Re = rNULL__
        if(ar__(1147:1153) .EQ. '') aZ_T90 = rNULL__
        if(ar__(1171:1177) .EQ. '') Z_T95Re = rNULL__
        if(ar__(1179:1185) .EQ. '') aZ_T95 = rNULL__
        if(ar__(1202:1208) .EQ. '') Z_T99Re = rNULL__
        if(ar__(1210:1216) .EQ. '') aZ_T99 = rNULL__
        if(ar__(1218:1223) .EQ. '') logMgasAvg = rNULL__
        if(ar__(1225:1230) .EQ. '') logMgasAvs = rNULL__
        if(ar__(1232:1239) .EQ. '') Vels2Re = rNULL__
        if(ar__(1241:1247) .EQ. '') e_Vels2Re = rNULL__
        if(ar__(1249:1256) .EQ. '') VelHa2Re = rNULL__
        if(ar__(1258:1264) .EQ. '') e_VelHa2Re = rNULL__
        if(ar__(1266:1273) .EQ. '') Vels1Re = rNULL__
        if(ar__(1275:1281) .EQ. '') e_Vels1Re = rNULL__
        if(ar__(1283:1290) .EQ. '') VelHa1Re = rNULL__
        if(ar__(1292:1298) .EQ. '') e_VelHa1Re = rNULL__
        if(ar__(1300:1306) .EQ. '') logSFRs100 = rNULL__
        if(ar__(1308:1314) .EQ. '') logSFRs10 = rNULL__
        if(ar__(1316:1322) .EQ. '') VdispHac = rNULL__
        if(ar__(1324:1330) .EQ. '') Vdispsc = rNULL__
        if(ar__(1332:1338) .EQ. '') VdispHa1re = rNULL__
        if(ar__(1340:1346) .EQ. '') Vdisps1Re = rNULL__
        if(ar__(1348:1353) .EQ. '') logM1Re = rNULL__
        if(ar__(1397:1403) .EQ. '') R50V = rNULL__
        if(ar__(1405:1411) .EQ. '') e_R50V = rNULL__
        if(ar__(1413:1419) .EQ. '') R50M = rNULL__
        if(ar__(1421:1427) .EQ. '') e_R50M = rNULL__
        if(ar__(1429:1435) .EQ. '') logMR50V = rNULL__
        if(ar__(1444:1450) .EQ. '') logMgasAvgl = rNULL__
        if(ar__(1452:1461) .EQ. '') Avgas1Re = rNULL__
        if(ar__(1463:1469) .EQ. '') e_Avgas1Re = rNULL__
        if(ar__(1471:1480) .EQ. '') Avs = rNULL__
        if(ar__(1482:1488) .EQ. '') e_Avs = rNULL__
        if(ar__(1490:1496) .EQ. '') Lambda1Re = rNULL__
        if(ar__(1498:1507) .EQ. '') e_Lambda1Re = rNULL__
        if(ar__(1509:1516) .EQ. '') zNSA = rNULL__
        if(ar__(1518:1524) .EQ. '') MNSA = rNULL__
        if(ar__(1526:1532) .EQ. '') IncNSA = rNULL__
        if(ar__(1534:1541) .EQ. '') F_OII_26Re = rNULL__
        if(ar__(1543:1549) .EQ. '') e_F_OII_26Re = rNULL__
        if(ar__(1551:1559) .EQ. '') F_OII_26a = rNULL__
        if(ar__(1561:1567) .EQ. '') e_F_OII_26a = rNULL__
        if(ar__(1569:1576) .EQ. '') F_OII_28Re = rNULL__
        if(ar__(1578:1584) .EQ. '') e_F_OII_28Re = rNULL__
        if(ar__(1586:1594) .EQ. '') F_OII_28a = rNULL__
        if(ar__(1596:1602) .EQ. '') e_F_OII_28a = rNULL__
        if(ar__(1604:1611) .EQ. '') FHI34Re = rNULL__
        if(ar__(1613:1619) .EQ. '') e_FHI34Re = rNULL__
        if(ar__(1621:1629) .EQ. '') FHI34a = rNULL__
        if(ar__(1631:1637) .EQ. '') e_FHI34a = rNULL__
        if(ar__(1639:1646) .EQ. '') FHI97Re = rNULL__
        if(ar__(1648:1654) .EQ. '') e_FHI97Re = rNULL__
        if(ar__(1656:1663) .EQ. '') FHI97a = rNULL__
        if(ar__(1665:1671) .EQ. '') e_FHI97a = rNULL__
        if(ar__(1673:1680) .EQ. '') FHeI38Re = rNULL__
        if(ar__(1682:1688) .EQ. '') e_FHeI38Re = rNULL__
        if(ar__(1690:1697) .EQ. '') FHeI38a = rNULL__
        if(ar__(1699:1705) .EQ. '') e_FHeI38a = rNULL__
        if(ar__(1707:1714) .EQ. '') FHI89Re = rNULL__
        if(ar__(1716:1722) .EQ. '') e_FHI89Re = rNULL__
        if(ar__(1724:1731) .EQ. '') FHI89a = rNULL__
        if(ar__(1733:1739) .EQ. '') e_FHI89a = rNULL__
        if(ar__(1741:1748) .EQ. '') FHeI39Re = rNULL__
        if(ar__(1750:1756) .EQ. '') e_FHeI39Re = rNULL__
        if(ar__(1758:1765) .EQ. '') FHeI39a = rNULL__
        if(ar__(1767:1773) .EQ. '') e_FHeI39a = rNULL__
        if(ar__(1775:1782) .EQ. '') F_NeIII_Re = rNULL__
        if(ar__(1784:1790) .EQ. '') e_F_NeIII_Re = rNULL__
        if(ar__(1792:1799) .EQ. '') F_NeIII_a = rNULL__
        if(ar__(1801:1807) .EQ. '') e_F_NeIII_a = rNULL__
        if(ar__(1809:1816) .EQ. '') FCaIIRe = rNULL__
        if(ar__(1818:1824) .EQ. '') e_FCaIIRe = rNULL__
        if(ar__(1826:1833) .EQ. '') FCaIIa = rNULL__
        if(ar__(1835:1841) .EQ. '') e_FCaIIa = rNULL__
        if(ar__(1843:1850) .EQ. '') FHepsRe = rNULL__
        if(ar__(1852:1858) .EQ. '') e_FHepsRe = rNULL__
        if(ar__(1860:1867) .EQ. '') FHepsa = rNULL__
        if(ar__(1869:1875) .EQ. '') e_FHepsa = rNULL__
        if(ar__(1877:1884) .EQ. '') FHdRe = rNULL__
        if(ar__(1886:1892) .EQ. '') e_FHdRe = rNULL__
        if(ar__(1894:1901) .EQ. '') FHda = rNULL__
        if(ar__(1903:1909) .EQ. '') e_FHda = rNULL__
        if(ar__(1911:1918) .EQ. '') FHgRe = rNULL__
        if(ar__(1920:1926) .EQ. '') e_FHgRe = rNULL__
        if(ar__(1928:1935) .EQ. '') FHga = rNULL__
        if(ar__(1937:1943) .EQ. '') e_FHga = rNULL__
        if(ar__(1945:1952) .EQ. '') FHbRe = rNULL__
        if(ar__(1954:1960) .EQ. '') e_FHbRe = rNULL__
        if(ar__(1962:1969) .EQ. '') FHba = rNULL__
        if(ar__(1971:1977) .EQ. '') e_FHba = rNULL__
        if(ar__(1979:1986) .EQ. '') F_OIII_49Re = rNULL__
        if(ar__(1988:1995) .EQ. '') e_F_OIII_49Re = rNULL__
        if(ar__(1997:2004) .EQ. '') F_OIII_49a = rNULL__
        if(ar__(2006:2013) .EQ. '') e_F_OIII_49a = rNULL__
        if(ar__(2015:2022) .EQ. '') F_OIII_50Re = rNULL__
        if(ar__(2024:2030) .EQ. '') e_F_OIII_50Re = rNULL__
        if(ar__(2032:2040) .EQ. '') F_OIII_50a = rNULL__
        if(ar__(2042:2048) .EQ. '') e_F_OIII_50a = rNULL__
        if(ar__(2050:2057) .EQ. '') FHeI50Re = rNULL__
        if(ar__(2059:2065) .EQ. '') e_FHeI50Re = rNULL__
        if(ar__(2067:2075) .EQ. '') FHeI50a = rNULL__
        if(ar__(2077:2083) .EQ. '') e_FHeI50a = rNULL__
        if(ar__(2085:2092) .EQ. '') F_NI_51Re = rNULL__
        if(ar__(2094:2100) .EQ. '') e_F_NI_51Re = rNULL__
        if(ar__(2102:2109) .EQ. '') F_NI_51a = rNULL__
        if(ar__(2111:2117) .EQ. '') e_F_NI_51a = rNULL__
        if(ar__(2119:2126) .EQ. '') F_NI_52Re = rNULL__
        if(ar__(2128:2134) .EQ. '') e_F_NI_52Re = rNULL__
        if(ar__(2136:2143) .EQ. '') F_NI_52a = rNULL__
        if(ar__(2145:2151) .EQ. '') e_F_NI_52a = rNULL__
        if(ar__(2153:2160) .EQ. '') FHeI58Re = rNULL__
        if(ar__(2162:2168) .EQ. '') e_FHeI58Re = rNULL__
        if(ar__(2170:2177) .EQ. '') FHeI58a = rNULL__
        if(ar__(2179:2185) .EQ. '') e_FHeI58a = rNULL__
        if(ar__(2187:2194) .EQ. '') FNaI89Re = rNULL__
        if(ar__(2196:2202) .EQ. '') e_FNaI89Re = rNULL__
        if(ar__(2204:2211) .EQ. '') FNaI89a = rNULL__
        if(ar__(2213:2219) .EQ. '') e_FNaI89a = rNULL__
        if(ar__(2221:2228) .EQ. '') FNaI95Re = rNULL__
        if(ar__(2230:2236) .EQ. '') e_FNaI95Re = rNULL__
        if(ar__(2238:2245) .EQ. '') FNaI95a = rNULL__
        if(ar__(2247:2253) .EQ. '') e_FNaI95a = rNULL__
        if(ar__(2255:2262) .EQ. '') F_OI_Re = rNULL__
        if(ar__(2264:2270) .EQ. '') e_F_OI_Re = rNULL__
        if(ar__(2272:2279) .EQ. '') F_OI_a = rNULL__
        if(ar__(2281:2287) .EQ. '') e_F_OI_a = rNULL__
        if(ar__(2289:2296) .EQ. '') F_NII_48Re = rNULL__
        if(ar__(2298:2304) .EQ. '') e_F_NII_48Re = rNULL__
        if(ar__(2306:2314) .EQ. '') F_NII_48a = rNULL__
        if(ar__(2316:2322) .EQ. '') e_F_NII_48a = rNULL__
        if(ar__(2324:2331) .EQ. '') FHaRe = rNULL__
        if(ar__(2333:2340) .EQ. '') e_FHaRe = rNULL__
        if(ar__(2342:2350) .EQ. '') FHaa = rNULL__
        if(ar__(2352:2359) .EQ. '') e_FHaa = rNULL__
        if(ar__(2361:2368) .EQ. '') F_NII_83Re = rNULL__
        if(ar__(2370:2376) .EQ. '') e_F_NII_83Re = rNULL__
        if(ar__(2378:2386) .EQ. '') F_NII_83a = rNULL__
        if(ar__(2388:2394) .EQ. '') e_F_NII_83a = rNULL__
        if(ar__(2396:2403) .EQ. '') F_SII_16Re = rNULL__
        if(ar__(2405:2411) .EQ. '') e_F_SII_16Re = rNULL__
        if(ar__(2413:2420) .EQ. '') F_SII_16a = rNULL__
        if(ar__(2422:2428) .EQ. '') e_F_SII_16a = rNULL__
        if(ar__(2430:2437) .EQ. '') F_SII_30Re = rNULL__
        if(ar__(2439:2445) .EQ. '') e_F_SII_30Re = rNULL__
        if(ar__(2447:2454) .EQ. '') F_SII_30a = rNULL__
        if(ar__(2456:2462) .EQ. '') e_F_SII_30a = rNULL__
        if(ar__(2464:2471) .EQ. '') F_ArIII_Re = rNULL__
        if(ar__(2473:2479) .EQ. '') e_F_ArIII_Re = rNULL__
        if(ar__(2481:2488) .EQ. '') F_ArIII_a = rNULL__
        if(ar__(2490:2496) .EQ. '') e_F_ArIII_a = rNULL__
        if(ar__(2498:2505) .EQ. '') FHI14Re = rNULL__
        if(ar__(2507:2513) .EQ. '') e_FHI14Re = rNULL__
        if(ar__(2515:2522) .EQ. '') FHI14a = rNULL__
        if(ar__(2524:2530) .EQ. '') e_FHI14a = rNULL__
        if(ar__(2532:2539) .EQ. '') F_SIII_90Re = rNULL__
        if(ar__(2541:2547) .EQ. '') e_F_SIII_90Re = rNULL__
        if(ar__(2549:2556) .EQ. '') F_SIII_90a = rNULL__
        if(ar__(2558:2567) .EQ. '') e_F_SIII_90a = rNULL__
        if(ar__(2569:2576) .EQ. '') F_FeII_Re = rNULL__
        if(ar__(2578:2584) .EQ. '') e_F_FeII_Re = rNULL__
        if(ar__(2586:2593) .EQ. '') F_FeII_a = rNULL__
        if(ar__(2595:2601) .EQ. '') e_F_FeII_a = rNULL__
        if(ar__(2603:2610) .EQ. '') F_SIII_95Re = rNULL__
        if(ar__(2612:2618) .EQ. '') e_F_SIII_95Re = rNULL__
        if(ar__(2620:2627) .EQ. '') F_SIII_95a = rNULL__
        if(ar__(2629:2635) .EQ. '') e_F_SIII_95a = rNULL__
        if(ar__(2637:2642) .EQ. '') OH_Mar13_N2Re = rNULL__
        if(ar__(2644:2649) .EQ. '') e_OH_Mar13_N2Re = rNULL__
        if(ar__(2651:2658) .EQ. '') OH_Mar13_N2a = rNULL__
        if(ar__(2660:2666) .EQ. '') e_OH_Mar13_N2a = rNULL__
        if(ar__(2668:2674) .EQ. '') OH_Mar13_O3N2Re = rNULL__
        if(ar__(2676:2682) .EQ. '') e_OH_Mar13_O3N2Re = rNULL__
        if(ar__(2684:2691) .EQ. '') OH_Mar13_O3N2a = rNULL__
        if(ar__(2693:2699) .EQ. '') e_OH_Mar13_O3N2a = rNULL__
        if(ar__(2701:2707) .EQ. '') OH_T04Re = rNULL__
        if(ar__(2709:2714) .EQ. '') e_OH_T04Re = rNULL__
        if(ar__(2716:2722) .EQ. '') OH_T04a = rNULL__
        if(ar__(2724:2729) .EQ. '') e_OH_T04a = rNULL__
        if(ar__(2731:2736) .EQ. '') OH_Pet04_N2linRe = rNULL__
        if(ar__(2738:2743) .EQ. '') e_OH_Pet04_N2linRe = rNULL__
        if(ar__(2745:2752) .EQ. '') OH_Pet04_N2lina = rNULL__
        if(ar__(2754:2760) .EQ. '') e_OH_Pet04_N2lina = rNULL__
        if(ar__(2762:2767) .EQ. '') OH_Pet04_N2polRe = rNULL__
        if(ar__(2769:2774) .EQ. '') e_OH_Pet04_N2polRe = rNULL__
        if(ar__(2776:2782) .EQ. '') OH_Pet04_N2pola = rNULL__
        if(ar__(2784:2789) .EQ. '') e_OH_Pet04_N2pola = rNULL__
        if(ar__(2791:2796) .EQ. '') OH_Pet04_O3N2Re = rNULL__
        if(ar__(2798:2803) .EQ. '') e_OH_Pet04_O3N2Re = rNULL__
        if(ar__(2805:2812) .EQ. '') OH_Pet04_O3N2a = rNULL__
        if(ar__(2814:2820) .EQ. '') e_OH_Pet04_O3N2a = rNULL__
        if(ar__(2822:2827) .EQ. '') OH_Kew02_N2O2Re = rNULL__
        if(ar__(2829:2834) .EQ. '') e_OH_Kew02_N2O2Re = rNULL__
        if(ar__(2836:2842) .EQ. '') OH_Kew02_N2O2a = rNULL__
        if(ar__(2844:2849) .EQ. '') e_OH_Kew02_N2O2a = rNULL__
        if(ar__(2851:2857) .EQ. '') OH_Pil10_ONSRe = rNULL__
        if(ar__(2859:2865) .EQ. '') e_OH_Pil10_ONSRe = rNULL__
        if(ar__(2867:2874) .EQ. '') OH_Pil10_ONSa = rNULL__
        if(ar__(2876:2882) .EQ. '') e_OH_Pil10_ONSa = rNULL__
        if(ar__(2884:2889) .EQ. '') OH_Pil10_ONRe = rNULL__
        if(ar__(2891:2896) .EQ. '') e_OH_Pil10_ONRe = rNULL__
        if(ar__(2898:2904) .EQ. '') OH_Pil10_ONa = rNULL__
        if(ar__(2906:2911) .EQ. '') e_OH_Pil10_ONa = rNULL__
        if(ar__(2913:2918) .EQ. '') OH_Pil11_NSRe = rNULL__
        if(ar__(2920:2925) .EQ. '') e_OH_Pil11_NSRe = rNULL__
        if(ar__(2927:2934) .EQ. '') OH_Pil11_NSa = rNULL__
        if(ar__(2936:2942) .EQ. '') e_OH_Pil11_NSa = rNULL__
        if(ar__(2944:2949) .EQ. '') OH_Cur20_RS32Re = rNULL__
        if(ar__(2951:2956) .EQ. '') e_OH_Cur20_RS32Re = rNULL__
        if(ar__(2958:2965) .EQ. '') OH_Cur20_RS32a = rNULL__
        if(ar__(2967:2973) .EQ. '') e_OH_Cur20_RS32a = rNULL__
        if(ar__(2975:2980) .EQ. '') OH_Cur20_R3Re = rNULL__
        if(ar__(2982:2987) .EQ. '') e_OH_Cur20_R3Re = rNULL__
        if(ar__(2989:2996) .EQ. '') OH_Cur20_R3a = rNULL__
        if(ar__(2998:3004) .EQ. '') e_OH_Cur20_R3a = rNULL__
        if(ar__(3006:3011) .EQ. '') OH_Cur20_O3O2Re = rNULL__
        if(ar__(3013:3018) .EQ. '') e_OH_Cur20_O3O2Re = rNULL__
        if(ar__(3020:3026) .EQ. '') OH_Cur20_O3O2a = rNULL__
        if(ar__(3028:3033) .EQ. '') e_OH_Cur20_O3O2a = rNULL__
        if(ar__(3035:3040) .EQ. '') OH_Cur20_S2Re = rNULL__
        if(ar__(3042:3047) .EQ. '') e_OH_Cur20_S2Re = rNULL__
        if(ar__(3049:3055) .EQ. '') OH_Cur20_S2a = rNULL__
        if(ar__(3057:3062) .EQ. '') e_OH_Cur20_S2a = rNULL__
        if(ar__(3064:3069) .EQ. '') OH_Cur20_R2Re = rNULL__
        if(ar__(3071:3076) .EQ. '') e_OH_Cur20_R2Re = rNULL__
        if(ar__(3078:3084) .EQ. '') OH_Cur20_R2a = rNULL__
        if(ar__(3086:3091) .EQ. '') e_OH_Cur20_R2a = rNULL__
        if(ar__(3093:3098) .EQ. '') OH_Cur20_N2Re = rNULL__
        if(ar__(3100:3105) .EQ. '') e_OH_Cur20_N2Re = rNULL__
        if(ar__(3107:3113) .EQ. '') OH_Cur20_N2a = rNULL__
        if(ar__(3115:3120) .EQ. '') e_OH_Cur20_N2a = rNULL__
        if(ar__(3122:3127) .EQ. '') OH_Cur20_R23Re = rNULL__
        if(ar__(3129:3134) .EQ. '') e_OH_Cur20_R23Re = rNULL__
        if(ar__(3136:3142) .EQ. '') OH_Cur20_R23a = rNULL__
        if(ar__(3144:3149) .EQ. '') e_OH_Cur20_R23a = rNULL__
        if(ar__(3151:3157) .EQ. '') OH_Cur20_O3N2Re = rNULL__
        if(ar__(3159:3165) .EQ. '') e_OH_Cur20_O3N2Re = rNULL__
        if(ar__(3167:3174) .EQ. '') OH_Cur20_O3N2a = rNULL__
        if(ar__(3176:3182) .EQ. '') e_OH_Cur20_O3N2a = rNULL__
        if(ar__(3184:3189) .EQ. '') OH_Cur20_O3S2Re = rNULL__
        if(ar__(3191:3196) .EQ. '') e_OH_Cur20_O3S2Re = rNULL__
        if(ar__(3198:3204) .EQ. '') OH_Cur20_O3S2a = rNULL__
        if(ar__(3206:3211) .EQ. '') e_OH_Cur20_O3S2a = rNULL__
        if(ar__(3213:3219) .EQ. '') OH_KK04Re = rNULL__
        if(ar__(3221:3227) .EQ. '') e_OH_KK04Re = rNULL__
        if(ar__(3229:3236) .EQ. '') OH_KK04a = rNULL__
        if(ar__(3238:3244) .EQ. '') e_OH_KK04a = rNULL__
        if(ar__(3246:3251) .EQ. '') OH_Pil16_RRe = rNULL__
        if(ar__(3253:3258) .EQ. '') e_OH_Pil16_RRe = rNULL__
        if(ar__(3260:3266) .EQ. '') OH_Pil16_Ra = rNULL__
        if(ar__(3268:3273) .EQ. '') e_OH_Pil16_Ra = rNULL__
        if(ar__(3275:3280) .EQ. '') OH_Pil16_SRe = rNULL__
        if(ar__(3282:3287) .EQ. '') e_OH_Pil16_SRe = rNULL__
        if(ar__(3289:3295) .EQ. '') OH_Pil16_Sa = rNULL__
        if(ar__(3297:3302) .EQ. '') e_OH_Pil16_Sa = rNULL__
        if(ar__(3304:3309) .EQ. '') OH_HoRe = rNULL__
        if(ar__(3311:3316) .EQ. '') e_OH_HoRe = rNULL__
        if(ar__(3318:3324) .EQ. '') OH_Hoa = rNULL__
        if(ar__(3326:3331) .EQ. '') e_OH_Hoa = rNULL__
        if(ar__(3333:3339) .EQ. '') U_Dors_O32Re = rNULL__
        if(ar__(3341:3346) .EQ. '') e_U_Dors_O32Re = rNULL__
        if(ar__(3348:3354) .EQ. '') U_Dors_O32a = rNULL__
        if(ar__(3356:3361) .EQ. '') e_U_Dors_O32a = rNULL__
        if(ar__(3363:3369) .EQ. '') U_Dors_SRe = rNULL__
        if(ar__(3371:3376) .EQ. '') e_U_Dors_SRe = rNULL__
        if(ar__(3378:3384) .EQ. '') U_Dors_Sa = rNULL__
        if(ar__(3386:3391) .EQ. '') e_U_Dors_Sa = rNULL__
        if(ar__(3393:3399) .EQ. '') U_Mor16_O23fsRe = rNULL__
        if(ar__(3401:3406) .EQ. '') e_U_Mor16_O23fsRe = rNULL__
        if(ar__(3408:3414) .EQ. '') U_Mor16_O23fsa = rNULL__
        if(ar__(3416:3421) .EQ. '') e_U_Mor16_O23fsa = rNULL__
        if(ar__(3423:3429) .EQ. '') U_Mor16_O23tsRe = rNULL__
        if(ar__(3431:3436) .EQ. '') e_U_Mor16_O23tsRe = rNULL__
        if(ar__(3438:3444) .EQ. '') U_Mor16_O23tsa = rNULL__
        if(ar__(3446:3451) .EQ. '') e_U_Mor16_O23tsa = rNULL__
        if(ar__(3453:3458) .EQ. '') NH_Pil16RRe = rNULL__
        if(ar__(3460:3465) .EQ. '') e_NH_Pil16RRe = rNULL__
        if(ar__(3467:3473) .EQ. '') NH_Pil16Ra = rNULL__
        if(ar__(3475:3480) .EQ. '') e_NH_Pil16Ra = rNULL__
        if(ar__(3482:3488) .EQ. '') NO_Pil16RRe = rNULL__
        if(ar__(3490:3495) .EQ. '') e_NO_Pil16RRe = rNULL__
        if(ar__(3497:3503) .EQ. '') NO_Pil16Ra = rNULL__
        if(ar__(3505:3510) .EQ. '') e_NO_Pil16Ra = rNULL__
        if(ar__(3512:3518) .EQ. '') NO_Pil16_HoRRe = rNULL__
        if(ar__(3520:3525) .EQ. '') e_NO_Pil16_HoRRe = rNULL__
        if(ar__(3527:3533) .EQ. '') NO_Pil16_HoRa = rNULL__
        if(ar__(3535:3540) .EQ. '') e_NO_Pil16_HoRa = rNULL__
        if(ar__(3542:3548) .EQ. '') NO_Pil16_N2R2Re = rNULL__
        if(ar__(3550:3555) .EQ. '') e_NO_Pil16_N2R2Re = rNULL__
        if(ar__(3557:3563) .EQ. '') NO_Pil16_N2R2a = rNULL__
        if(ar__(3565:3570) .EQ. '') e_NO_Pil16_N2R2a = rNULL__
        if(ar__(3572:3576) .EQ. '') Ne_Oster_SRe = rNULL__
        if(ar__(3578:3582) .EQ. '') e_Ne_Oster_SRe = rNULL__
        if(ar__(3584:3589) .EQ. '') Ne_Oster_Sa = rNULL__
        if(ar__(3591:3595) .EQ. '') e_Ne_Oster_Sa = rNULL__
        if(ar__(3597:3602) .EQ. '') HdRe = rNULL__
        if(ar__(3604:3608) .EQ. '') e_HdRe = rNULL__
        if(ar__(3610:3616) .EQ. '') Hda = rNULL__
        if(ar__(3618:3622) .EQ. '') e_Hda = rNULL__
        if(ar__(3624:3629) .EQ. '') HbRe = rNULL__
        if(ar__(3631:3635) .EQ. '') e_HbRe = rNULL__
        if(ar__(3637:3643) .EQ. '') Hba = rNULL__
        if(ar__(3645:3649) .EQ. '') e_Hba = rNULL__
        if(ar__(3651:3656) .EQ. '') MgbRe = rNULL__
        if(ar__(3658:3662) .EQ. '') e_MgbRe = rNULL__
        if(ar__(3664:3669) .EQ. '') Mgba = rNULL__
        if(ar__(3671:3675) .EQ. '') e_Mgba = rNULL__
        if(ar__(3677:3682) .EQ. '') Fe5270Re = rNULL__
        if(ar__(3684:3688) .EQ. '') e_Fe5270Re = rNULL__
        if(ar__(3690:3695) .EQ. '') Fe5270a = rNULL__
        if(ar__(3697:3701) .EQ. '') e_Fe5270a = rNULL__
        if(ar__(3703:3708) .EQ. '') Fe5335Re = rNULL__
        if(ar__(3710:3715) .EQ. '') e_Fe5335Re = rNULL__
        if(ar__(3717:3723) .EQ. '') Fe5335a = rNULL__
        if(ar__(3725:3730) .EQ. '') e_Fe5335a = rNULL__
        if(ar__(3732:3737) .EQ. '') D4000Re = rNULL__
        if(ar__(3739:3744) .EQ. '') e_D4000Re = rNULL__
        if(ar__(3746:3752) .EQ. '') D4000a = rNULL__
        if(ar__(3754:3759) .EQ. '') e_D4000a = rNULL__
        if(ar__(3761:3766) .EQ. '') HdmodRe = rNULL__
        if(ar__(3768:3772) .EQ. '') e_HdmodRe = rNULL__
        if(ar__(3774:3780) .EQ. '') Hdmoda = rNULL__
        if(ar__(3782:3786) .EQ. '') e_Hdmoda = rNULL__
        if(ar__(3788:3793) .EQ. '') HgRe = rNULL__
        if(ar__(3795:3799) .EQ. '') e_HgRe = rNULL__
        if(ar__(3801:3806) .EQ. '') Hga = rNULL__
        if(ar__(3808:3812) .EQ. '') e_Hga = rNULL__
        if(ar__(4293:4294) .EQ. '') Type = iNULL__
        if(ar__(4300:4305) .EQ. '') NSAid = iNULL__
        if(ar__(4307:4314) .EQ. '') w_Vmax = rNULL__
        if(ar__(4316:4323) .EQ. '') w_Num = rNULL__
        if(ar__(4325:4325) .EQ. '') Q = iNULL__
c    ..............Just test output...........
        write(6,1)
     +  Name,Plate,IFUdsgn,PlateIFU,MaNGAId,RAdeg,DEdeg,logSFRHa,FoV,
     +  Reff,e_logMass,e_logSFRHa,logMass,logSFRs,logNII_Hac,
     +  e_logNII_Hac,logOIII_Hbc,e_logOIII_Hbc,logSII_Hac,
     +  e_logSII_Hac,logOII_Hbc,e_logOII_Hbc,EWHac,e_EWHac,logZReLW,
     +  e_logZReLW,alphaZReLW,e_alphaZReLW,logZReMW,e_logZReMW,
     +  alphaZReMW,e_alphaZReMW,logAgeReLW,e_logAgeReLW,alphaAgeReLW,
     +  e_alphaAgeReLW,logAgeReMW,e_logAgeReMW,alphaAgeReMW,
     +  e_alphaAgeReMW,ReffA,DL,DA,PA,ell,logMassgas,vel_sig,
     +  e_vel_sig,logSFRsf,logSFRdc,OH_O3N2c,e_OH_O3N2c,OH_N2c,
     +  e_OH_N2c,OH_ONSc,e_OH_ONSc,OH_R23c,e_OH_R23c,OH_pyqzc,
     +  e_OH_pyqzc,OH_t2c,e_OH_t2c,OH_M08c,e_OH_M08c,OH_T04c,
     +  e_OH_T04c,OH_dopc,e_OH_dopc,OH_O3N2_EPc,e_OH_O3N2_EPc,
     +  logOI_Hac,e_logOI_Hac,Ha_Hbc,e_Ha_Hbc,logNII_HaRe,
     +  e_logNII_HaRe,logOIII_HbRe,e_logOIII_HbRe,logSII_HaRe,
     +  e_logSII_HaRe,logOII_HbRe,e_logOII_HbRe,logOI_HaRe,
     +  e_logOI_HaRe,EWHaRe,e_EWHaRe,Ha_HbRe,e_Ha_HbRe,logNII_Ha,
     +  e_logNII_Ha,logOIII_Hb,e_logOIII_Hb,logSII_Ha,e_logSII_Ha,
     +  logOII_Hb,e_logOII_Hb,logOI_Ha,e_logOI_Ha,EWHa,e_EWHa,Ha_Hb,
     +  logsigMc,e_logsigMc,logsigMRe,e_logsigMRe,logsigM,e_logsigM,
     +  T30,Z_T30,Z_T30Re,aZ_T30,T40,Z_T40,Z_T40Re,aZ_T40,T50,Z_T50,
     +  Z_T50Re,aZ_T50,T60,Z_T60,Z_T60Re,aZ_T60,T70,Z_T70,Z_T70Re,
     +  aZH_T70,T80,Z_T80,Z_T80Re,aZ_T80,T90,Z_T90,Z_T90Re,aZ_T90,T95,
     +  Z_T95,Z_T95Re,aZ_T95,T99,Z_T99,Z_T99Re,aZ_T99,logMgasAvg,
     +  logMgasAvs,Vels2Re,e_Vels2Re,VelHa2Re,e_VelHa2Re,Vels1Re,
     +  e_Vels1Re,VelHa1Re,e_VelHa1Re,logSFRs100,logSFRs10,VdispHac,
     +  Vdispsc,VdispHa1re,Vdisps1Re,logM1Re,MLint,MLavg,FHac,e_FHac,
     +  R50V,e_R50V,R50M,e_R50M,logMR50V,e_logMR50V,logMgasAvgl,
     +  Avgas1Re,e_Avgas1Re,Avs,e_Avs,Lambda1Re,e_Lambda1Re,zNSA,MNSA,
     +  IncNSA,F_OII_26Re,e_F_OII_26Re,F_OII_26a,e_F_OII_26a,
     +  F_OII_28Re,e_F_OII_28Re,F_OII_28a,e_F_OII_28a,FHI34Re,
     +  e_FHI34Re,FHI34a,e_FHI34a,FHI97Re,e_FHI97Re,FHI97a,e_FHI97a,
     +  FHeI38Re,e_FHeI38Re,FHeI38a,e_FHeI38a,FHI89Re,e_FHI89Re,
     +  FHI89a,e_FHI89a,FHeI39Re,e_FHeI39Re,FHeI39a,e_FHeI39a,
     +  F_NeIII_Re,e_F_NeIII_Re,F_NeIII_a,e_F_NeIII_a,FCaIIRe,
     +  e_FCaIIRe,FCaIIa,e_FCaIIa,FHepsRe,e_FHepsRe,FHepsa,e_FHepsa,
     +  FHdRe,e_FHdRe,FHda,e_FHda,FHgRe,e_FHgRe,FHga,e_FHga,FHbRe,
     +  e_FHbRe,FHba,e_FHba,F_OIII_49Re,e_F_OIII_49Re,F_OIII_49a,
     +  e_F_OIII_49a,F_OIII_50Re,e_F_OIII_50Re,F_OIII_50a,
     +  e_F_OIII_50a,FHeI50Re,e_FHeI50Re,FHeI50a,e_FHeI50a,F_NI_51Re,
     +  e_F_NI_51Re,F_NI_51a,e_F_NI_51a,F_NI_52Re,e_F_NI_52Re,
     +  F_NI_52a,e_F_NI_52a,FHeI58Re,e_FHeI58Re,FHeI58a,e_FHeI58a,
     +  FNaI89Re,e_FNaI89Re,FNaI89a,e_FNaI89a,FNaI95Re,e_FNaI95Re,
     +  FNaI95a,e_FNaI95a,F_OI_Re,e_F_OI_Re,F_OI_a,e_F_OI_a,
     +  F_NII_48Re,e_F_NII_48Re,F_NII_48a,e_F_NII_48a,FHaRe,e_FHaRe,
     +  FHaa,e_FHaa,F_NII_83Re,e_F_NII_83Re,F_NII_83a,e_F_NII_83a,
     +  F_SII_16Re,e_F_SII_16Re,F_SII_16a,e_F_SII_16a,F_SII_30Re,
     +  e_F_SII_30Re,F_SII_30a,e_F_SII_30a,F_ArIII_Re,e_F_ArIII_Re,
     +  F_ArIII_a,e_F_ArIII_a,FHI14Re,e_FHI14Re,FHI14a,e_FHI14a,
     +  F_SIII_90Re,e_F_SIII_90Re,F_SIII_90a,e_F_SIII_90a,F_FeII_Re,
     +  e_F_FeII_Re,F_FeII_a,e_F_FeII_a,F_SIII_95Re,e_F_SIII_95Re,
     +  F_SIII_95a,e_F_SIII_95a,OH_Mar13_N2Re,e_OH_Mar13_N2Re,
     +  OH_Mar13_N2a,e_OH_Mar13_N2a,OH_Mar13_O3N2Re,e_OH_Mar13_O3N2Re,
     +  OH_Mar13_O3N2a,e_OH_Mar13_O3N2a,OH_T04Re,e_OH_T04Re,OH_T04a,
     +  e_OH_T04a,OH_Pet04_N2linRe,e_OH_Pet04_N2linRe,OH_Pet04_N2lina,
     +  e_OH_Pet04_N2lina,OH_Pet04_N2polRe,e_OH_Pet04_N2polRe,
     +  OH_Pet04_N2pola,e_OH_Pet04_N2pola,OH_Pet04_O3N2Re,
     +  e_OH_Pet04_O3N2Re,OH_Pet04_O3N2a,e_OH_Pet04_O3N2a,
     +  OH_Kew02_N2O2Re,e_OH_Kew02_N2O2Re,OH_Kew02_N2O2a,
     +  e_OH_Kew02_N2O2a,OH_Pil10_ONSRe,e_OH_Pil10_ONSRe,
     +  OH_Pil10_ONSa,e_OH_Pil10_ONSa,OH_Pil10_ONRe,e_OH_Pil10_ONRe,
     +  OH_Pil10_ONa,e_OH_Pil10_ONa,OH_Pil11_NSRe,e_OH_Pil11_NSRe,
     +  OH_Pil11_NSa,e_OH_Pil11_NSa,OH_Cur20_RS32Re,e_OH_Cur20_RS32Re,
     +  OH_Cur20_RS32a,e_OH_Cur20_RS32a,OH_Cur20_R3Re,e_OH_Cur20_R3Re,
     +  OH_Cur20_R3a,e_OH_Cur20_R3a,OH_Cur20_O3O2Re,e_OH_Cur20_O3O2Re,
     +  OH_Cur20_O3O2a,e_OH_Cur20_O3O2a,OH_Cur20_S2Re,e_OH_Cur20_S2Re,
     +  OH_Cur20_S2a,e_OH_Cur20_S2a,OH_Cur20_R2Re,e_OH_Cur20_R2Re,
     +  OH_Cur20_R2a,e_OH_Cur20_R2a,OH_Cur20_N2Re,e_OH_Cur20_N2Re,
     +  OH_Cur20_N2a,e_OH_Cur20_N2a,OH_Cur20_R23Re,e_OH_Cur20_R23Re,
     +  OH_Cur20_R23a,e_OH_Cur20_R23a,OH_Cur20_O3N2Re,
     +  e_OH_Cur20_O3N2Re,OH_Cur20_O3N2a,e_OH_Cur20_O3N2a,
     +  OH_Cur20_O3S2Re,e_OH_Cur20_O3S2Re,OH_Cur20_O3S2a,
     +  e_OH_Cur20_O3S2a,OH_KK04Re,e_OH_KK04Re,OH_KK04a,e_OH_KK04a,
     +  OH_Pil16_RRe,e_OH_Pil16_RRe,OH_Pil16_Ra,e_OH_Pil16_Ra,
     +  OH_Pil16_SRe,e_OH_Pil16_SRe,OH_Pil16_Sa,e_OH_Pil16_Sa,OH_HoRe,
     +  e_OH_HoRe,OH_Hoa,e_OH_Hoa,U_Dors_O32Re,e_U_Dors_O32Re,
     +  U_Dors_O32a,e_U_Dors_O32a,U_Dors_SRe,e_U_Dors_SRe,U_Dors_Sa,
     +  e_U_Dors_Sa,U_Mor16_O23fsRe,e_U_Mor16_O23fsRe,U_Mor16_O23fsa,
     +  e_U_Mor16_O23fsa,U_Mor16_O23tsRe,e_U_Mor16_O23tsRe,
     +  U_Mor16_O23tsa,e_U_Mor16_O23tsa,NH_Pil16RRe,e_NH_Pil16RRe,
     +  NH_Pil16Ra,e_NH_Pil16Ra,NO_Pil16RRe,e_NO_Pil16RRe,NO_Pil16Ra,
     +  e_NO_Pil16Ra,NO_Pil16_HoRRe,e_NO_Pil16_HoRRe,NO_Pil16_HoRa,
     +  e_NO_Pil16_HoRa,NO_Pil16_N2R2Re,e_NO_Pil16_N2R2Re,
     +  NO_Pil16_N2R2a,e_NO_Pil16_N2R2a,Ne_Oster_SRe,e_Ne_Oster_SRe,
     +  Ne_Oster_Sa,e_Ne_Oster_Sa,HdRe,e_HdRe,Hda,e_Hda,HbRe,e_HbRe,
     +  Hba,e_Hba,MgbRe,e_MgbRe,Mgba,e_Mgba,Fe5270Re,e_Fe5270Re,
     +  Fe5270a,e_Fe5270a,Fe5335Re,e_Fe5335Re,Fe5335a,e_Fe5335a,
     +  D4000Re,e_D4000Re,D4000a,e_D4000a,HdmodRe,e_HdmodRe,Hdmoda,
     +  e_Hdmoda,HgRe,e_HgRe,Hga,e_Hga,umag,e_umag,uMag_1,e_uMag_1,
     +  gmag,e_gmag,gMag_1,e_gMag_1,rmag,e_rmag,rMag_1,e_rMag_1,imag,
     +  e_imag,iMag_1,e_iMag_1,Bmag,e_Bmag,BMag_1,e_BMag_1,Vmag,
     +  e_Vmag,VMag_1,e_VMag_1,RJmag,e_RJmag,RJMag_1,e_RJMag_1,R50,
     +  e_R50,R90,e_R90,C,e_C,B_V,e_B_V,B_R,e_B_R,logMph,e_logMph,
     +  SBV1Re,e_SBV1Re,SBVR50,e_SBVR50,nNSA,u_g,g_r,r_i,i_z,P_CD,P_E,
     +  P_S0,P_Sa,P_Sab,P_Sb,P_Sbc,P_Sc,P_Scd,P_Sd,P_Sdm,P_Sm,P_Irr,
     +  Type,Morph,NSAid,w_Vmax,w_Num,Q
c    .......End.of.Just test output...........
      end do
      close(1)

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