FORTRAN Generation
(/./ftp/cats/J/A_A/706/A95)

Conversion of standardized ReadMe file for file /./ftp/cats/J/A_A/706/A95 into FORTRAN code for loading all data files into arrays.

Note that special values are assigned to unknown or unspecified numbers (also called NULL numbers); when necessary, the coordinate components making up the right ascension and declination are converted into floating-point numbers representing these angles in degrees.



      program load_ReadMe
C=============================================================================
C  F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2026-Aug-11
C=============================================================================
*  This code was generated from the ReadMe file documenting a catalogue
*  according to the "Standard for Documentation of Astronomical Catalogues"
*  currently in use by the Astronomical Data Centers (CDS, ADC, A&A)
*  (see full documentation at URL http://vizier.u-strasbg.fr/doc/catstd.htx)
*  Please report problems or questions to   
C=============================================================================

      implicit none
*  Unspecified or NULL values, generally corresponding to blank columns,
*  are assigned one of the following special values:
*     rNULL__    for unknown or NULL floating-point values
*     iNULL__    for unknown or NULL   integer      values
      real*4     rNULL__
      integer*4  iNULL__
      parameter  (rNULL__=--2147483648.)  	! NULL real number
      parameter  (iNULL__=(-2147483647-1))	! NULL int  number
      integer    idig			! testing NULL number

C=============================================================================
Cat. J/A+A/706/A95       PHANGS-HST Halpha Nebulae Catalogue          (Barnes+, 2026)
*================================================================================
*The PHANGS-MUSE/HST-H{alpha} nebulae catalogue.
*Parsec-scale resolved structure, physical conditions, and stellar associations
*across nearby galaxies.
*    Barnes A.T., Chandar R., Kreckel K., Belfiore F., Pathak D., Thilker D.,
*    Leroy A.K., Groves B., Glover S.C.O., McClain R., Amiri A., Bazzi Z.,
*    Boquien M., Congiu E., Dale D.A., Egorov O.V., Emsellem E., Grasha K.,
*    Gonzalez-Lobos J., Henny K., He H., Indebetouw R., Lee J.C., Li J.,
*    Liang F.-H., Larson K., Maschmann D., Meidt S.E., Mendez-Delgado J.E.,
*    Neumann J., Pan H.-A., Querejeta M., Rosolowsky E., Sarbadhicary S.K.,
*    Scheuermann F., Ubeda L., Williams T.G., Weinbeck T.D., Whitmore B.,
*    Wofford A., the PHANGS collaboration
*    <Astron. Astrophys. 706, A95 (2026)>
*    =2026A&A...706A..95B        (SIMBAD/NED BibCode)
C=============================================================================

C  Internal variables

      integer*4 i__

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

C  Declarations for 'neb_cat.dat'	! Main PHANGS-MUSE/HST Halpha Nebulae Catalogu

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

C  J2000 position composed of: RAdeg DEdeg
      character*7   Galaxy     (nr__) ! Galaxy name (galaxy_name)
      integer*4     Region     (nr__) ! Region ID (region_ID)
      integer*4     fEdgeHST   (nr__) ! Flag=1 if region lies on HST image
*                                           edge (flag_edge_hst)
      integer*4     fTouchHST  (nr__) ! Flag=1 if region touches segmentation
*                                           boundary (flag_touch_hst)
      character*7   fManualHST (nr__) ! ?=- Manual classification flags
*                                           (visual checks) (flag_manual_hst)
      integer*4     fEdgeMUSE  (nr__) ! Flag=1 if region lies on MUSE cube
*                                           edge (flag_edge_muse)
      integer*4     fStarMUSE  (nr__) ! Flag=1 if affected by bright star in
*                                           MUSE data (flag_star_muse)
      real*8        RAdeg      (nr__) ! (deg) Right ascension (J2000) of region
*                                           centroid (ra)
      real*8        DEdeg      (nr__) ! (deg) Declination (J2000) of region
*                                           centroid (dec)
      real*8        RApdeg     (nr__) ! (deg) H{alpha} pixel peak right ascension
*                                           (J2000) (ra_peak)
      real*8        DEpdeg     (nr__) ! (deg) H{alpha} pixel peak declination
*                                           (J2000) (dec_peak)
      real*8        RAMdeg     (nr__) ! (deg) MUSE right ascension (J2000) of
*                                           MUSE centroid (ra_muse)
      real*8        DEMdeg     (nr__) ! (deg) MUSE declination (J2000) of MUSE
*                                           centroid (dec_muse)
      real*8        rR25Deproj (nr__) ! Deprojected galactocentric radius,
*                                           in R25 units (r_R25_deproj)
      real*8        rrefFdeproj(nr__) ! Deprojected galactocentric radius,
*                                           in Re units (r_reff_deproj)
      real*8        phiDeproj  (nr__) ! (rad) Deprojected azimuthal angle within
*                                           host galaxy (phi_deproj)
      integer*4     Env        (nr__) ! [1/9] Environment flag (1)
      real*8        Areacirc   (nr__) ! (arcsec+2) Circularized area (area_circ_arcsec2)
      real*8        Radcirc    (nr__) ! (arcsec) Circularized radius (arcsec)
*                                           (radius_circ_arcsec)
      real*8        Radcircpc  (nr__) ! (pc) Circularized radius (pc)
*                                           (radius_circ_pc)
      real*8        Area2      (nr__) ! (arcsec+2) Second-moment area of region
*                                           (area_mom_arcsec2)
      real*8        amaj2      (nr__) ! (arcsec) Second-moment major axis
*                                           (radius_major_mom_arcsec)
      real*8        bmin2      (nr__) ! (arcsec) Second-moment minor axis
*                                           (radius_minor_mom_arcsec)
      real*8        PA2        (nr__) ! (deg) [] Position angle (moments-based)
*                                           (radius_pos_angle_deg)
      real*8        Rad        (nr__) ! (arcsec) Moment-based radius (arcsec)
*                                           (radius_mom_arcsec)
      real*8        Radpc      (nr__) ! (pc) Moment-based radius (pc)
*                                           (radius_mom_pc)
      real*8        Raddeconv  (nr__) ! (pc) ?=- Deconvolved moment-based radius
*                                           (radius_mom_deconv_pc)
      real*8        RadStr     (nr__) ! (pc) ?=- Stromgren radius
*                                           (radius_stromgren_pc)
      real*8        RadcircpcM (nr__) ! (pc) Circularized radius (pc; PHANGS-MUSE)
*                                           (radius_circ_pc_muse)
      integer*4     ComplexScore(nr__) ! Morphological complexity score
*                                           (Sect. 3.2) (complexity_score)
      real*4        FHaHST     (nr__) ! (mW/m2) HST H{alpha} integrated flux (raw,
*                                           uncorrected) (halpha_flux)
      real*4        e_FHaHST   (nr__) ! (mW/m2) Error on HST H{alpha} flux
*                                           (halpha_flux_err)
      real*4        FHaHSTMax  (nr__) ! (mW/m2) Maximum pixel flux in HST H{alpha} map
*                                           (halpha_flux_max)
      real*4        FHaHSTMin  (nr__) ! (mW/m2) Minimum pixel flux in HST H{alpha} map
*                                           (halpha_flux_min)
      real*4        FHameanHST (nr__) ! (mW/m2) Mean pixel flux in HST H{alpha} map
*                                           (halpha_flux_mean)
      real*4        FHacorrHST (nr__) ! (mW/m2) Extinction-corrected HST H{alpha} flux
*                                           (halpha_flux_corr)
      real*4        e_FHacorrHST(nr__) ! (mW/m2) Error on extinction-corrected HST
*                                           H{alpha} flux (halpha_flux_corr_err)
      real*4        FHa        (nr__) ! (mW/m2) MUSE H{alpha} integrated flux
*                                           (PHANGS-MUSE) (halpha_flux_muse)
      real*4        FHb        (nr__) ! (mW/m2) MUSE H{beta} integrated flux
*                                           (PHANGS-MUSE) (hbeta_flux_muse)
      real*4        FOIII5007  (nr__) ! (mW/m2) MUSE [O III] {lambda}5007 flux
*                                           (PHANGS-MUSE) (oiii5007_flux_muse)
      real*4        FNII6584   (nr__) ! (mW/m2) MUSE [N II] {lambda}6584 flux
*                                           (PHANGS-MUSE) (nii6584_flux_muse)
      real*4        FSII6716   (nr__) ! (mW/m2) MUSE [S II] {lambda}6716 flux
*                                           (PHANGS-MUSE) (sii6716_flux_muse)
      real*4        FSII6731   (nr__) ! (mW/m2) MUSE [S II] {lambda}6731 flux
*                                           (PHANGS-MUSE) (sii6731_flux_muse)
      real*4        FSIII9069  (nr__) ! (mW/m2) MUSE [S III] {lambda}9069 flux
*                                           (PHANGS-MUSE) (siii9069_flux_muse)
      real*4        FOI6300    (nr__) ! (mW/m2) MUSE [O I] {lambda}6300 flux
*                                           (PHANGS-MUSE) (oi6300_flux_muse)
      real*4        FNII5755   (nr__) ! (mW/m2) ?=- MUSE [N II] {lambda}5755 flux
*                                           (PHANGS-MUSE) (nii5755_flux_muse)
      real*4        FHeI5876   (nr__) ! (mW/m2) MUSE He I {lambda}5876 flux
*                                           (PHANGS-MUSE) (hei5876_flux_muse)
      real*4        FSIII6312  (nr__) ! (mW/m2) ?=- MUSE [S III] {lambda}6312 flux
*                                           (PHANGS-MUSE) (siii6312_flux_muse)
      real*4        FOI6363    (nr__) ! (mW/m2) MUSE [O I] {lambda}6363 flux
*                                           (PHANGS-MUSE) (oi6363_flux_muse)
      real*4        FOII7319   (nr__) ! (mW/m2) MUSE [O II] {lambda}7319 flux
*                                           (PHANGS-MUSE) (oii7319_flux_muse)
      real*4        FOII7330   (nr__) ! (mW/m2) MUSE [O II] {lambda}7330 flux
*                                           (PHANGS-MUSE) (oii7330_flux_muse)
      real*4        e_FHa      (nr__) ! (mW/m2) Error on H{alpha} flux (PHANGS-MUSE)
*                                           (halpha_flux_err_muse)
      real*4        e_FHb      (nr__) ! (mW/m2) Error on H{beta} flux (PHANGS-MUSE)
*                                           (hbeta_flux_err_muse)
      real*4        e_FOIII5007(nr__) ! (mW/m2) Error on [O III] {lambda}5007 flux
*                                           (PHANGS-MUSE)
*                                           (oiii5007_flux_err_muse)
      real*4        e_FNII6584 (nr__) ! (mW/m2) Error on [N II] {lambda}6584 flux
*                                           (PHANGS-MUSE) (nii6584_flux_err_muse)
      real*4        e_FSII6716 (nr__) ! (mW/m2) Error on [S II] {lambda}6716 flux
*                                           (PHANGS-MUSE) (sii6716_flux_err_muse)
      real*4        e_FSII6731 (nr__) ! (mW/m2) Error on [S II] {lambda}6731 flux
*                                           (PHANGS-MUSE) (sii6731_flux_err_muse)
      real*4        e_FSIII9069(nr__) ! (mW/m2) Error on [S III] {lambda}9069 flux
*                                           (PHANGS-MUSE)
*                                           (siii9069_flux_err_muse)
      real*4        e_FOI6300  (nr__) ! (mW/m2) Error on [O I] {lambda}6300 flux
*                                           (PHANGS-MUSE) (oi6300_flux_err_muse)
      real*4        e_FNII5755 (nr__) ! (mW/m2) ?=- Error on [N II] {lambda}5755 flux
*                                           (PHANGS-MUSE) (nii5755_flux_err_muse)
      real*4        e_FHeI5876 (nr__) ! (mW/m2) ?=- Error on He I {lambda}5876 flux
*                                           (PHANGS-MUSE) (hei5876_flux_err_muse)
      real*4        e_FSIII6312(nr__) ! (mW/m2) ?=- Error on [Si II] {lambda}6312 flux
*                                           (PHANGS-MUSE)
*                                           (siii6312_flux_err_muse)
      real*4        e_FOI6363  (nr__) ! (mW/m2) Error on [O I] {lambda}6363 flux
*                                           (PHANGS-MUSE) (oi6363_flux_err_muse)
      real*4        e_FOII7319 (nr__) ! (mW/m2) Error on [O I] {lambda}7319 flux
*                                           (PHANGS-MUSE) (oii7319_flux_err_muse)
      real*4        e_FOII7330 (nr__) ! (mW/m2) Error on [O I] {lambda}7319 flux
*                                           (PHANGS-MUSE) (oii7330_flux_err_muse)
      real*4        FHacorr    (nr__) ! (mW/m2) Extinction-corrected H{alpha} flux
*                                           (PHANGS-MUSE) (halpha_flux_corr_muse)
      real*4        FHbcorr    (nr__) ! (mW/m2) Extinction-corrected H{beta} flux
*                                           (PHANGS-MUSE) (hbeta_flux_corr_muse)
      real*4        FOIII5007corr(nr__) ! (mW/m2) Extinction-corrected [O III]
*                                           {lambda}5007 flux (PHANGS-MUSE)
*                                           (oiii5007_flux_corr_muse)
      real*4        FNII6584corr(nr__) ! (mW/m2) Extinction-corrected [N II]
*                                           {lambda}6584 flux (PHANGS-MUSE)
*                                           (nii6584_flux_corr_muse)
      real*4        FSII6716corr(nr__) ! (mW/m2) Extinction-corrected [S II]
*                                           {lambda}6716 flux (PHANGS-MUSE)
*                                           (sii6716_flux_corr_muse)
      real*4        FSII6731corr(nr__) ! (mW/m2) Extinction-corrected [S II]
*                                           {lambda}6731 flux (PHANGS-MUSE)
*                                           (sii6731_flux_corr_muse)
      real*4        FSIII9069corr(nr__) ! (mW/m2) Extinction-corrected [S III]
*                                           {lambda}9069 flux (PHANGS-MUSE)
*                                           (siii9069_flux_corr_muse)
      real*4        FOI6300corr(nr__) ! (mW/m2) Extinction-corrected [O I]
*                                           {lambda}6300 flux (PHANGS-MUSE)
*                                           (oi6300_flux_corr_muse)
      real*4        FNII5755corr(nr__) ! (mW/m2) ?=- Extinction-corrected [N II]
*                                           {lambda}5755 flux (PHANGS-MUSE)
*                                           (nii5755_flux_corr_muse)
      real*4        FHeI5876corr(nr__) ! (mW/m2) Extinction-corrected He I
*                                           {lambda}5876 flux (PHANGS-MUSE)
*                                           (hei5876_flux_corr_muse)
      real*4        FSIII6312corr(nr__) ! (mW/m2) ?=- Extinction-corrected [S III]
*                                           {lambda}6312 flux (PHANGS-MUSE)
*                                           (siii6312_flux_corr_muse)
      real*4        FOI6363corr(nr__) ! (mW/m2) Extinction-corrected [O I]
*                                           {lambda}6363 flux (PHANGS-MUSE)
*                                           (oi6363_flux_corr_muse)
      real*4        FOII7319corr(nr__) ! (mW/m2) ?=- Extinction-corrected [O II]
*                                           {lambda}7319 flux (PHANGS-MUSE)
*                                           (oii7319_flux_corr_muse)
      real*4        FOII7330corr(nr__) ! (mW/m2) ?=- Extinction-corrected [O II]
*                                           {lambda}7330 flux (PHANGS-MUSE)
*                                           (oii7330_flux_corr_muse)
      real*4        e_FHacorr  (nr__) ! (mW/m2) Error on corrected H{alpha} flux
*                                           (PHANGS-MUSE)
*                                           (halpha_flux_corr_err_muse)
      real*4        e_FHbcorr  (nr__) ! (mW/m2) Error on corrected H{beta} flux
*                                           (PHANGS-MUSE)
*                                           (hbeta_flux_corr_err_muse)
      real*4        e_FOIII5007corr(nr__) ! (mW/m2) ?=- Error on corrected [O III]
*                                           {lambda}5007 flux (PHANGS-MUSE)
*                                           (oiii5007_fl ux_corr_err_muse)
      real*4        e_FNII6584corr(nr__) ! (mW/m2) ?=- Error on corrected [N II]
*                                           {lambda}5007 flux (PHANGS-MUSE)
*                                           (nii6584_flux_corr_err_muse)
      real*4        e_FSII6716corr(nr__) ! (mW/m2) ?=- Error on corrected [S II]
*                                           {lambda}6716 flux (PHANGS-MUSE)
*                                           (sii6716_flux_corr_err_muse)
      real*4        e_FSII6731corr(nr__) ! (mW/m2) ?=- Error on corrected [S II]
*                                           {lambda}6731 flux (PHANGS-MUSE)
*                                           (sii6731_flux_corr_err_muse)
      real*4        e_FSIII9069corr(nr__) ! (mW/m2) ?=- Error on corrected [S III]
*                                           {lambda}9069 flux (PHANGS-MUSE)
*                                           (siii9069_fl ux_corr_err_muse)
      real*4        e_FOI6300corr(nr__) ! (mW/m2) ?=- Error on corrected [O I]
*                                           {lambda}6300 flux (PHANGS-MUSE)
*                                           (oi6300_flux_corr_err_muse)
      real*4        e_FNII5755corr(nr__) ! (mW/m2) ?=- Error on corrected [N II]
*                                           {lambda}5755 flux (PHANGS-MUSE)
*                                           (nii5755_flux_corr_err_muse)
      real*4        e_FHeI5876corr(nr__) ! (mW/m2) ?=- Error on corrected He I
*                                           {lambda}5876 flux (PHANGS-MUSE)
*                                           (hei5876_flux_corr_err_muse)
      real*4        e_FSIII6312corr(nr__) ! (mW/m2) ?=- Error on corrected [S III]
*                                           {lambda}6312 flux (PHANGS-MUSE)
*                                           (siii6312_fl ux_corr_err_muse)
      real*4        e_FOI6363corr(nr__) ! (mW/m2) ?=- Error on corrected [O I]
*                                           {lambda}6363 flux (PHANGS-MUSE)
*                                           (oi6363_flux_corr_err_muse)
      real*4        e_FOII7319corr(nr__) ! (mW/m2) ?=- Error on corrected [O II]
*                                           {lambda}7319 flux (PHANGS-MUSE)
*                                           (oii7319_flux_corr_err_muse)
      real*4        e_FOII7330corr(nr__) ! (mW/m2) ?=- Error on corrected [O II]
*                                           {lambda}7330 flux (PHANGS-MUSE)
*                                           (oii7330_flux_corr_err_muse)
      real*4        VHa        (nr__) ! (km/s) H{alpha} velocity (PHANGS-MUSE)
*                                           (halpha_vel_muse)
      real*4        e_VHa      (nr__) ! (km/s) Error on H{alpha} velocity
*                                           (PHANGS-MUSE) (halpha_vel_err_muse)
      real*4        sigmaHa    (nr__) ! (km/s) H{alpha} velocity dispersion
*                                           (PHANGS-MUSE) (halpha_sigma_muse)
      real*4        e_sigmaHa  (nr__) ! (km/s) Error on H{alpha} velocity dispersion
*                                           (PHANGS-MUSE) (halpha_sigma_err_muse)
      real*4        VNII6584   (nr__) ! (km/s) Velocity on [N II] {lambda}6584
*                                           (PHANGS-MUSE) (nii6584_vel_muse)
      real*4        e_VNII6584 (nr__) ! (km/s) Error on velocity of [N II]
*                                           {lambda}6584 (PHANGS-MUSE)
*                                           (nii6584_vel_err_muse)
      real*4        sigmaNII6584(nr__) ! (km/s) Velocity dispersion of [N II]
*                                           {lambda}6584 (PHANGS-MUSE)
*                                           (nii6584_sigma_muse)
      real*4        e_sigmaNII6584(nr__) ! (km/s) Error on velocity dispersion of [N II]
*                                           {lambda}6584 (PHANGS-MUSE)
*                                           (nii6584_sigma_err_muse)
      real*4        VOIII5007  (nr__) ! (km/s) Velocity on [O III] {lambda}5007
*                                           (PHANGS-MUSE) (oiii5007_vel_muse)
      real*4        e_VOIII5007(nr__) ! (km/s) Error on velocity on [O III]
*                                           {lambda}5007 (PHANGS-MUSE)
*                                           (oiii5007_vel_err_muse)
      real*4        sigmaOIII5007(nr__) ! (km/s) Velocity dispersion og [O III]
*                                           {lambda}5007 (PHANGS-MUSE)
*                                           (oiii5007_sigma_muse)
      real*4        e_sigmaOIII5007(nr__) ! (km/s) Error on velocity dispersion of
*                                           [O III] {lambda}5007
*                                           (oiii5007_sigma_err_muse)
      integer*4     BPT_NII    (nr__) ! [-1/3]  BPT class using
*                                           [N II]/H{alpha} (bpt_nii_muse)
      integer*4     BPT_SII    (nr__) ! [-1/3]  BPT class using
*                                           [S II]/H{alpha} (bpt_sii_muse)
      integer*4     BPT_OI     (nr__) ! [-1/3]  BPT class using
*                                           [O I]/H{alpha} (bpt_oi_muse)
      real*4        EWHa       (nr__) ! (0.1nm) H{alpha} equivalent width
*                                           (PHANGS-MUSE) (halpha_ew_muse)
      real*4        e_EWHa     (nr__) ! (0.1nm) [] Uncertainty on H{alpha} equivalent
*                                           width (PHANGS-MUSE)
*                                           (halpha_ew_err_muse)
      real*4        EWHb       (nr__) ! (0.1nm) H{beta} equivalent width (PHANGS-MUSE)
*                                           (hbeta_ew_muse)
      real*4        e_EWHb     (nr__) ! (0.1nm) [] Uncertainty on H{beta} equivalent
*                                           width (PHANGS-MUSE)
*                                           (hbeta_ew_err_muse)
      real*4        EWHaraw    (nr__) ! (0.1nm) H{alpha} equivalent width from raw
*                                           spectrum (no continuum model,
*                                           PHANGS-MUSE) (halpha_ew_raw_muse)
      real*4        e_EWHaraw  (nr__) ! (0.1nm) [] Uncertainty on H{alpha} equivalent
*                                           width from raw spectrum (no continuum
*                                           model, PHANGS-MUSE)
*                                           (halpha_ew_raw_err_muse)
      real*4        EWHbraw    (nr__) ! (0.1nm) H{beta} equivalent width from raw
*                                           spectrum (no continuum model,
*                                           PHANGS-MUSE) (hbeta_ew_raw_muse)
      real*4        e_EWHbraw  (nr__) ! (0.1nm) [] Uncertainty on H{beta} equivalent
*                                           width from raw spectrum (no continuum
*                                           model, PHANGS-MUSE)
*                                           (hbeta_ew_raw_err_muse)
      real*4        EWHafit    (nr__) ! (0.1nm) H{alpha} equivalent width fitted model
*                                           (PHANGS-MUSE) (halpha_ew_fit_muse)
      real*4        e_EWHafit  (nr__) ! (0.1nm) [] Uncertainty on H{alpha} equivalent
*                                           width fitted model (PHANGS-MUSE)
*                                           (halpha_ew_fit_err_muse)
      real*4        EWHbfit    (nr__) ! (0.1nm) H{beta} equivalent width fitted model
*                                           (PHANGS-MUSE) (hbeta_ew_fit_muse)
      real*4        e_EWHbfit  (nr__) ! (0.1nm) [] Uncertainty on H{beta} equivalent
*                                           width fitted model (PHANGS-MUSE)
*                                           (hbeta_ew_fit_err_muse)
      real*4        EWHabgcorr (nr__) ! (0.1nm) ?=- Background-corrected H{beta}
*                                           equivalent width (PHANGS-MUSE)
*                                           (halpha_ew_bgcorr_muse)
      real*4        e_EWHabgcorr(nr__) ! (0.1nm) []?=- Uncertainty on
*                                           background-corrected H{beta}
*                                           equivalent width (PHANGS-MUSE)
*                                           (halpha_ew_bgcorr_err_muse)
      real*4        EWHbbgcorr (nr__) ! (0.1nm) ?=- Background-corrected H{beta}
*                                           equivalent width (PHANGS-MUSE)
*                                           (hbeta_ew_bgcorr_muse)
      real*4        e_EWHbbgcorr(nr__) ! (0.1nm) []?=- Uncertainty on
*                                           background-corrected H{beta}
*                                           equivalent width (PHANGS-MUSE)
*                                           (hbeta_ew_bgcorr_err_muse)
      real*4        LHaHST     (nr__) ! (10-7W) H{alpha} luminosity from HST
*                                           (halpha_lum)
      real*4        e_LHaHST   (nr__) ! (10-7W) Uncertainty on H{alpha} luminosity
*                                           from HST (halpha_lum_err)
      real*4        LHa        (nr__) ! (10-7W) H{alpha} luminosity from MUSE
*                                           (halpha_lum_muse)
      integer*4     HIIclassv2 (nr__) ! [0/1]? PHANGS-MUSE H II classification
*                                           (v2; PHANGS-MUSE) (hii_class_v2_muse)
      integer*4     HIIclassv3 (nr__) ! [0/1]? PHANGS-MUSE H II classification
*                                           (v3; PHANGS-MUSE) (hii_class_v3_muse)
      real*4        TeNII      (nr__) ! (K) Electron temperature from [N II]
*                                           (PHANGS-MUSE) (te_nii_muse)
      real*4        e_TeNII    (nr__) ! (K) ?=- Uncertainty on electron
*                                           temperature from [N II] (PHANGS-MUSE)
*                                           (te_nii_err_muse)
      real*4        TeSIII     (nr__) ! (K) ?=- Electron temperature from [S III]
*                                           (PHANGS-MUSE) (te_siii_muse)
      real*4        e_TeSIII   (nr__) ! (K) ?=- Uncertainty on electron
*                                           temperature from [S III]
*                                           (PHANGS-MUSE) (te_siii_err_muse)
      real*4        neSII      (nr__) ! (cm-3) ?=- Electron density from
*                                           [S II]6716/6731 (PHANGS-MUSE)
*                                           (ne_sii_muse)
      real*4        e_neSII    (nr__) ! (cm-3) ?=- Uncertainty on electron density
*                                           from [S II]6716/6731 (PHANGS-MUSE)
*                                           (ne_sii_err_muse)
      real*4        q          (nr__) ! (s-1) Ionizing photon rate Q(H^0) (q)
      real*4        e_q        (nr__) ! (s-1) Uncertainty on ionizing photon rate
*                                           Q(H^0) (q)
      real*4        logqd91    (nr__) ! ?=- log10 ionization parameter from
*                                           Diaz+1991 calibration (PHANGS-MUSE)
*                                           (logq_d91_muse)
      real*4        e_logqd91  (nr__) ! ?=- Uncertainty on log10 ionization
*                                           parameter (PHANGS-MUSE)
*                                           (logq_d91_err_muse)
      real*4        E_B_V      (nr__) ! (mag) ?=- Color excess E(B-V) from Balmer
*                                           decrement (PHANGS-MUSE) (ebv_muse)
      real*4        e_E_B_V    (nr__) ! (mag) ?=- Uncertainty on color excess E(B-V)
*                                           from Balmer decrement (PHANGS-MUSE)
*                                           (ebv_err_muse)
      real*4        AV         (nr__) ! (mag) ?=- Visual extinction A_V
*                                           (PHANGS-MUSE) (av_muse)
      real*4        MetScal    (nr__) ! ?=- Gas-phase metallicity (Scal code;
*                                           PHANGS-MUSE) (metallicity_scal_muse)
      real*4        e_MetScal  (nr__) ! ?=- Uncertainty on gas-phase
*                                           metallicity (Scal code; PHANGS-MUSE)
*                                           (metallicity_scal_err_muse)
      real*4        deltaMetScal(nr__) ! ?=- Offset from local metallicity
*                                           trend (PHANGS-MUSE)
*                                           (metallicity_scal_delta_ muse)
      integer*4     AgeANUV08pc(nr__) ! (Myr) ?=- Stellar association age
*                                           (NUV;8 pc; PHANGS-HST)
*                                           (age_assoc_nuv_08pc)
      real*4        e_AgeANUV08pc(nr__) ! (Myr) ?=- Uncertainty on stellar association
*                                           age (NUV;8 pc; PHANGS-HST)
*                                           (age_err_assoc_nuv_08pc)
      real*8        DEdegNUV08pc(nr__) ! (deg) ?=- Association declination
*                                           (NUV;8 pc; PHANGS-HST)
*                                           (dec_assoc_nuv_08pc)
      real*4        E_B_V_NUV08pc(nr__) ! (mag) ?=- Color excess E(B-V)
*                                           (NUV; 8 pc; PHANGS-HST)
*                                           (ebv_assoc_nuv_08pc)
      real*4        e_E_B_V_NUV08pc(nr__) ! (mag) ?=- Uncertainty on color excess E(B-V)
*                                           (NUV; 8 pc; PHANGS-HST)
*                                           (ebv_err_assoc_nuv_08pc)
      real*4        MassANUV08pc(nr__) ! (Msun) ?=- Stellar mass (NUV; 8 pc;
*                                           PHANGS-HST) (mass_assoc_nuv_08pc)
      real*4        e_MassANUV08pc(nr__) ! (Msun) ?=- Uncertainty on stellar mass
*                                           (NUV; 8 pc; PHANGS-HST)
*                                           (mass_err_assoc_nuv_08pc)
      integer*4     FlagMultNUV08p(nr__) ! [0/1] Flag: multiple associations
*                                           (NUV; 8 pc; PHANGS-HST)
*                                           (flag_multi_assoc_nuv_08pc)
      integer*4     FlagNoNUV08pc(nr__) ! [0/1] Flag: no association
*                                           (NUV; 8 pc; PHANGS-HST)
*                                           (flag_none_assoc_nuv_08pc)
      integer*4     FlagOneNUV08pc(nr__) ! [0/1] Flag: single association
*                                           (NUV; 8 pc; PHANGS-HST)
*                                           (flag_one_assoc_nuv_08pc)
      real*8        RAdegNUV08pc(nr__) ! (deg) ?=- Association right ascension
*                                           (NUV; 8 pc; PHANGS-HST)
*                                           (ra_assoc_nuv_08pc)
      integer*4     RegionIDNUV08p(nr__) ! ?=- Association identifier
*                                           (NUV; 8 pc; PHANGS-HST)
*                                           (region_ID_assoc_nuv_08pc)
      integer*4     AgeANUV16pc(nr__) ! (Myr) ?=- Association age
*                                           (NUV; 16 pc; PHANGS-HST)
*                                           (age_assoc_nuv_16pc)
      integer*4     e_AgeANUV16pc(nr__) ! (Myr) ?=- Uncertainty on association age
*                                           (NUV; 16 pc; PHANGS-HST)
*                                           (age_err_assoc_nuv_16pc)
      real*8        DEdegNUV16pc(nr__) ! (deg) ?=- Association declination
*                                           (NUV; 16 pc; PHANGS-HST)
*                                           (dec_assoc_nuv_16pc)
      real*4        E_B_V_NUV16pc(nr__) ! (mag) ?=- E(B-V) of association
*                                           (NUV; 16 pc; PHANGS-HST)
*                                           (ebv_assoc_nuv_16pc)
      real*4        e_E_B_V_NUV16pc(nr__) ! (mag) ?=- Uncertainty on E(B-V) of
*                                           association (NUV; 16 pc; PHANGS-HST)
*                                           (ebv_err_assoc_nuv_16pc)
      real*4        MassANUV16pc(nr__) ! (Msun) ?=- Stellar mass (NUV; 16 pc;
*                                           PHANGS-HST) (mass_assoc_nuv_16pc)
      real*4        e_MassANUV16pc(nr__) ! (Msun) ?=- Uncertainty on stellar mass
*                                           (NUV; 16 pc; PHANGS-HST)
*                                           (mass_err_assoc_nuv_16pc)
      integer*4     FlagMultNUV16p(nr__) ! [0/1] Flag: multiple association
*                                           (NUV; 16 pc; PHANGS-HST)
*                                           (flag_multi_assoc_nuv_16pc)
      integer*4     FlagNoNUV16pc(nr__) ! [0/1] Flag: no association
*                                           (NUV; 16 pc; PHANGS-HST)
*                                           (flag_none_assoc_nuv_16pc)
      integer*4     FlagOneNUV16pc(nr__) ! [0/1] Flag: single association
*                                           (NUV; 16 pc; PHANGS-HST)
*                                           (flag_one_assoc_nuv_16pc)
      real*8        RAdegNUV16pc(nr__) ! (deg) ?=- Association right ascension
*                                           (NUV; 16 pc; PHANGS-HST)
*                                           (ra_assoc_nuv_16pc)
      integer*4     RegionIDNUV16p(nr__) ! ?=- Association identifier
*                                           (NUV; 16 pc; PHANGS-HST)
*                                           (region_ID_assoc_nuv_16pc)
      integer*4     AgeANUV32pc(nr__) ! (Myr) ?=- Association age (NUV; 32 pc;
*                                           PHANGS-HST) (age_assoc_nuv_32pc)
      integer*4     e_AgeANUV32pc(nr__) ! (Myr) ?=- Uncertainty on association age
*                                           (NUV; 32 pc; PHANGS-HST)
*                                           (age_err_assoc_nuv_32pc)
      real*8        DEdegNUV32pc(nr__) ! (deg) ?=- Association declination
*                                           (NUV;32 pc; PHANGS-HST)
*                                           (dec_assoc_nuv_32pc)
      real*4        E_B_V_NUV32pc(nr__) ! (mag) ?=- E(B-V) of association
*                                           (NUV; 32 pc; PHANGS-HST)
*                                           (ebv_assoc_nuv_32pc)
      real*4        e_E_B_V_NUV32pc(nr__) ! (mag) ?=- Uncertainty on E(B-V) of
*                                           association (NUV; 32 pc; PHANGS-HST)
*                                           (ebv_err_assoc_nuv_32pc)
      real*4        MassANUV32pc(nr__) ! (Msun) ?=- Stellar mass (NUV; 32 pc;
*                                           PHANGS-HST) (mass_assoc_nuv_32pc)
      real*4        e_MassANUV32pc(nr__) ! (Msun) ?=- Uncertainty on stellar mass
*                                           (NUV; 32 pc; PHANGS-HST)
*                                           (mass_err_assoc_nuv_32pc)
      integer*4     FlagMultNUV32p(nr__) ! [0/1] Flag: multiple associations
*                                           (NUV; 32 pc; PHANGS-HST)
*                                           (flag_multi_assoc_nuv_32pc)
      integer*4     FlagNoNUV32pc(nr__) ! [0/1] Flag: no associations
*                                           (NUV; 32 pc; PHANGS-HST)
*                                           (flag_none_assoc_nuv_32pc)
      integer*4     FlagOneNUV32pc(nr__) ! [0/1] Flag: single associations
*                                           (NUV; 32 pc; PHANGS-HST)
*                                           (flag_one_assoc_nuv_32pc)
      real*8        RAdegNUV32pc(nr__) ! (deg) ?=- Association right ascension
*                                           (NUV; 32 pc; PHANGS-HST)
*                                           (ra_assoc_nuv_32pc)
      integer*4     RegionIDNUV32p(nr__) ! ?=- Association identifier
*                                           (NUV; 32 pc; PHANGS-HST)
*                                           (region_ID_assoc_nuv_32pc)
      integer*4     AgeANUV64pc(nr__) ! (Myr) ?=- Association age (NUV; 64 pc;
*                                           PHANGS-HST) (age_assoc_nuv_64pc)
      integer*4     e_AgeANUV64pc(nr__) ! (Myr) ?=- Uncertainty on association age
*                                           (NUV; 64 pc; PHANGS-HST)
*                                           (age_err_assoc_nuv_64pc)
      real*8        DEdegNUV64pc(nr__) ! (deg) ?=- Association declination
*                                           (NUV; 64 pc; PHANGS-HST)
*                                           (dec_assoc_nuv_64pc)
      real*4        E_B_V_NUV64pc(nr__) ! (mag) ?=- E(B-V) of association
*                                           (NUV; 64 pc; PHANGS-HST)
*                                           (ebv_assoc_nuv_64pc)
      real*4        e_E_B_V_NUV64pc(nr__) ! (mag) ?=- Uncertainty on E(B-V) of
*                                           association (NUV; 64 pc; PHANGS-HST)
*                                           (ebv_err_assoc_nuv_64pc)
      real*4        MassANUV64pc(nr__) ! (Msun) ?=- Stellar mass (NUV; 64 pc;
*                                           PHANGS-HST; PHANGS-HST)
*                                           (mass_assoc_nuv_64pc)
      real*4        e_MassANUV64pc(nr__) ! (Msun) ?=- Uncertainty on stellar mass
*                                           (NUV; 64 pc; PHANGS-HST; PHANGS-HST)
*                                           (mass_err_assoc_nuv_64pc)
      integer*4     FlagMultNUV64p(nr__) ! [0/1] Flag: multiple associations
*                                           (NUV; 64 pc; PHANGS-HST)
*                                           (flag_multi_assoc_nuv_64pc)
      integer*4     FlagNoNUV64pc(nr__) ! [0/1] Flag: no associations
*                                           (NUV; 64 pc; PHANGS-HST)
*                                           (flag_none_assoc_nuv_64pc)
      integer*4     FlagOneNUV64pc(nr__) ! [0/1] Flag: single associations
*                                           (NUV; 64 pc; PHANGS-HST)
*                                           (flag_one_assoc_nuv_64pc)
      real*8        RAdegNUV64pc(nr__) ! (deg) ?=- Association right ascension
*                                           (NUV; 64 pc; PHANGS-HST)
*                                           (ra_assoc_nuv_64pc)
      integer*4     RegionIDNUV64p(nr__) ! ?=- Association identifier
*                                           (NUV; 64 pc; PHANGS-HST)
*                                           (region_ID_assoc_nuv_64pc)
      integer*4     AgeAV08pc  (nr__) ! (Myr) ?=- Stellar association age
*                                            (V;8 pc; PHANGS-HST)
*                                            (age_assoc_v_08pc)
      integer*4     e_AgeAV08pc(nr__) ! (Myr) ?=- Uncertainty on stellar association
*                                           age (V;8 pc; PHANGS-HST)
*                                           (age_err_assoc_v_08pc)
      real*8        DEdegV08pc (nr__) ! (deg) ?=- Association declination
*                                           (V;8 pc; PHANGS-HST)
*                                           (dec_assoc_v_08pc)
      real*4        E_B_V_V08pc(nr__) ! (mag) ?=- Color excess E(B-V)
*                                           (V; 8 pc; PHANGS-HST)
*                                           (ebv_assoc_v_08pc)
      real*4        e_E_B_V_V08pc(nr__) ! (mag) ?=- Uncertainty on color excess E(B-V)
*                                           (V; 8 pc; PHANGS-HST)
*                                           (ebv_err_assoc_v_08pc)
      real*4        MassAV08pc (nr__) ! (Msun) ?=- Stellar mass (V; 8 pc; PHANGS-HST)
*                                           (mass_assoc_v_08pc)
      real*4        e_MassAV08pc(nr__) ! (Msun) ?=- Uncertainty on stellar mass
*                                           (V; 8 pc; PHANGS-HST)
*                                           (mass_err_assoc_v_08pc)
      integer*4     FlagMultV08pc(nr__) ! [0/1] Flag: multiple associations
*                                           (V; 8 pc; PHANGS-HST)
*                                           (flag_multi_assoc_v_08pc)
      integer*4     FlagNoV08pc(nr__) ! [0/1] Flag: no association
*                                           (V; 8 pc; PHANGS-HST)
*                                           (flag_none_assoc_v_08pc)
      integer*4     FlagOneV08pc(nr__) ! [0/1] Flag: single association
*                                           (V; 8 pc; PHANGS-HST)
*                                           (flag_one_assoc_v_08pc)
      real*8        RAdegV08pc (nr__) ! (deg) ?=- Association right ascension
*                                           (V; 8 pc; PHANGS-HST)
*                                           (ra_assoc_v_08pc)
      integer*4     RegionIDV08pc(nr__) ! ?=- Association identifier
*                                           (V; 8 pc; PHANGS-HST)
*                                           (region_ID_assoc_v_08pc)
      integer*4     AgeAV16pc  (nr__) ! (Myr) ?=- Association age (V; 16 pc;
*                                           PHANGS-HST) (age_assoc_v_16pc)
      integer*4     e_AgeAV16pc(nr__) ! (Myr) ?=- Uncertainty on association age
*                                           (V; 16 pc; PHANGS-HST)
*                                           (age_err_assoc_v_16pc)
      real*8        DEdegV16pc (nr__) ! (deg) ?=- Association declination
*                                           (V; 16 pc; PHANGS-HST)
*                                           (dec_assoc_v_16pc)
      real*4        E_B_V_V16pc(nr__) ! (mag) ?=- E(B-V) of association
*                                           (V; 16 pc; PHANGS-HST)
*                                           (ebv_assoc_v_16pc)
      real*4        e_E_B_V_V16pc(nr__) ! (mag) ?=- Uncertainty on E(B-V) of
*                                           association (V; 16 pc; PHANGS-HST)
*                                           (ebv_err_assoc_v_16pc)
      real*4        MassAV16pc (nr__) ! (Msun) ?=- Stellar mass (V; 16 pc;
*                                           PHANGS-HST) (mass_assoc_v_16pc)
      real*4        e_MassAV16pc(nr__) ! (Msun) ?=- Uncertainty on stellar mass
*                                           (V; 16 pc; PHANGS-HST)
*                                           (mass_err_assoc_v_16pc)
      integer*4     FlagMultV16pc(nr__) ! [0/1] Flag: multiple association
*                                           (V; 16 pc; PHANGS-HST)
*                                           (flag_multi_assoc_v_16pc)
      integer*4     FlagNoV16pc(nr__) ! [0/1]  Flag: no association
*                                           (V; 16 pc; PHANGS-HST)
*                                           (flag_none_assoc_v_16pc)
      integer*4     FlagOneV16pc(nr__) ! [0/1] Flag: single association
*                                           (V; 16 pc; PHANGS-HST)
*                                           (flag_one_assoc_v_16pc)
      real*8        RAdegV16pc (nr__) ! (deg) ?=- Association right ascension
*                                           (V; 16 pc; PHANGS-HST)
*                                           (ra_assoc_v_16pc)
      integer*4     RegionIDV16pc(nr__) ! ?=- Association identifier
*                                           (V; 16 pc; PHANGS-HST)
*                                           (region_ID_assoc_v_16pc)
      integer*4     AgeAV32pc  (nr__) ! (Myr) ?=- Association age (V; 32 pc;
*                                           PHANGS-HST) (age_assoc_v_32pc)
      integer*4     e_AgeAV32pc(nr__) ! (Myr) ?=- Uncertainty on association age
*                                           (V; 32 pc; PHANGS-HST)
*                                           (age_err_assoc_v_32pc)
      real*8        DEdegV32pc (nr__) ! (deg) ?=- Association declination (V;32 pc;
*                                           PHANGS-HST) (dec_assoc_v_32pc)
      real*4        E_B_V_V32pc(nr__) ! (mag) ?=- E(B-V) of association (V; 32 pc;
*                                           PHANGS-HST) (ebv_assoc_v_32pc)
      real*4        e_E_B_V_V32pc(nr__) ! (mag) ?=- Uncertainty on E(B-V) of
*                                           association (V; 32 pc; PHANGS-HST)
*                                           (ebv_err_assoc_v_32pc)
      real*4        MassAV32pc (nr__) ! (Msun) ?=- Stellar mass (V; 32 pc;
*                                           PHANGS-HST) (mass_assoc_v_32pc)
      real*4        e_MassAV32pc(nr__) ! (Msun) ?=- Uncertainty on stellar mass
*                                           (V; 32 pc; PHANGS-HST)
*                                           (mass_err_assoc_v_32pc)
      integer*4     FlagMultV32pc(nr__) ! [0/1] Flag: multiple associations
*                                           (V; 32 pc; PHANGS-HST)
*                                           (flag_multi_assoc_v_32pc)
      integer*4     FlagNoV32pc(nr__) ! [0/1]  Flag: no associations
*                                           (V; 32 pc; PHANGS-HST)
*                                           (flag_none_assoc_v_32pc)
      integer*4     FlagOneV32pc(nr__) ! [0/1] Flag: single associations
*                                           (V; 32 pc; PHANGS-HST)
*                                           (flag_one_assoc_v_32pc)
      real*8        RAdegV32pc (nr__) ! (deg) ?=- Association right ascension
*                                           (V; 32 pc; PHANGS-HST)
*                                           (ra_assoc_v_32pc)
      integer*4     RegionIDV32pc(nr__) ! ?=- Association identifier
*                                           (V; 32 pc; PHANGS-HST)
*                                           (region_ID_assoc_v_32pc)
      integer*4     AgeAV64pc  (nr__) ! (Myr) ?=- Association age
*                                           (V; 64 pc; PHANGS-HST)
*                                           (age_assoc_v_64pc)
      integer*4     e_AgeAV64pc(nr__) ! (Myr) ?=- Uncertainty on association age
*                                           (V; 64 pc; PHANGS-HST)
*                                           (age_err_assoc_v_64pc)
      real*8        DEdegV64pc (nr__) ! (deg) ?=- Association declination
*                                           (V; 64 pc; PHANGS-HST)
*                                           (dec_assoc_v_64pc)
      real*4        E_B_V_V64pc(nr__) ! (mag) ?=- E(B-V) of association (V; 64 pc;
*                                           PHANGS-HST) (ebv_assoc_v_64pc)
      real*4        e_E_B_V_V64pc(nr__) ! (mag) ?=- Uncertainty on E(B-V) of
*                                           association (V; 64 pc; PHANGS-HST)
*                                           (ebv_err_assoc_v_64pc)
      real*4        MassAV64pc (nr__) ! (Msun) ?=- Stellar mass (V; 64 pc;
*                                           PHANGS-HST; PHANGS-HST)
*                                           (mass_assoc_v_64pc)
      real*4        e_MassAV64pc(nr__) ! (Msun) ?=- Uncertainty on stellar mass
*                                           (V; 64 pc; PHANGS-HST; PHANGS-HST)
*                                           (mass_err_assoc_v_64pc)
      integer*4     FlagMultV64pc(nr__) ! [0/1] Flag: multiple associations
*                                           (V; 64 pc; PHANGS-HST)
*                                           (flag_multi_assoc_v_64pc)
      integer*4     FlagNoV64pc(nr__) ! [0/1]  Flag: no associations
*                                           (V; 64 pc; PHANGS-HST)
*                                           (flag_none_assoc_v_64pc)
      integer*4     FlagOneV64pc(nr__) ! [0/1] Flag: single associations
*                                           (V; 64 pc; PHANGS-HST)
*                                           (flag_one_assoc_v_64pc)
      real*8        RAdegV64pc (nr__) ! (deg) ?=- Association right ascension
*                                           (V; 64 pc; PHANGS-HST)
*                                           (ra_assoc_v_64pc)
      integer*4     RegionIDV64pc(nr__) ! ?=- Association identifier (V; 64 pc;
*                                          PHANGS-HST) (region_ID_assoc_v_64pc)
*Note (1): Environment fl ag following Querejeta et al. (2021A&A...656A.133Q):
*    Centre (1), Bar (2-3), Arm (4-5), Interarm (6-7), Disc (8-9) (environment)

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

C  Declarations for 'list.dat'	! List of fits maps and masks

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

C  J2000 position composed of: RAdeg DEdeg
      real*8        RAdeg_1    (nr__1) ! (deg) Right Ascension of center (J2000)
      real*8        DEdeg_1    (nr__1) ! (deg) Declination of center (J2000)
      integer*4     Nx         (nr__1) ! Number of pixels along X-axis
      integer*4     Ny         (nr__1) ! Number of pixels along Y-axis
      integer*4     size       (nr__1) ! (Kibyte) Size of FITS file
      character*37  FileName   (nr__1) ! Name of FITS file with subdirectory
      character*36  Title      (nr__1) ! Title of the FITS file

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

C  Loading file 'neb_cat.dat'	! Main PHANGS-MUSE/HST Halpha Nebulae Catalogu

C  Format for file interpretation

    1 format(
     +  A7,1X,I4,1X,I1,1X,I1,1X,A7,1X,I1,1X,I1,1X,F15.11,1X,F15.11,1X,
     +  F15.11,1X,F15.11,1X,F15.11,1X,F15.11,1X,F10.8,1X,F10.8,1X,
     +  F10.8,1X,I1,1X,F11.8,1X,F9.7,1X,F8.4,1X,F12.9,1X,F10.8,1X,
     +  F10.8,1X,F8.3,1X,F10.8,1X,F9.5,1X,F10.6,1X,F10.6,1X,F10.6,1X,
     +  I3,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,
     +  E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,
     +  1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,
     +  E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,
     +  1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,
     +  E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,
     +  1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,
     +  E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,
     +  1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,
     +  E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,
     +  1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,
     +  E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,I2,1X,I2,1X,I2,1X,E12.6,
     +  1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,
     +  E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,
     +  1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,I2,1X,I2,1X,
     +  E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,
     +  1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,E12.6,1X,
     +  E12.6,1X,E12.6,1X,E12.6,1X,I3,1X,F6.1,1X,F8.4,1X,F6.4,1X,F6.4,
     +  1X,E12.6,1X,E12.6,1X,I1,1X,I1,1X,I1,1X,F15.11,1X,I3,1X,I3,1X,
     +  I3,1X,F15.11,1X,F4.2,1X,F4.2,1X,E12.6,1X,E12.6,1X,I1,1X,I1,1X,
     +  I1,1X,F15.11,1X,I4,1X,I5,1X,I3,1X,F15.11,1X,F4.2,1X,F4.2,1X,
     +  E12.6,1X,E12.6,1X,I1,1X,I1,1X,I1,1X,F15.11,1X,I4,1X,I3,1X,I3,
     +  1X,F15.11,1X,F4.2,1X,F4.2,1X,E12.6,1X,E12.6,1X,I1,1X,I1,1X,I1,
     +  1X,F15.11,1X,I4,1X,I5,1X,I5,1X,F15.11,1X,F4.2,1X,F4.2,1X,
     +  E12.6,1X,E12.6,1X,I1,1X,I1,1X,I1,1X,F15.11,1X,I4,1X,I5,1X,I5,
     +  1X,F15.11,1X,F4.2,1X,F4.2,1X,E12.6,1X,E12.6,1X,I1,1X,I1,1X,I1,
     +  1X,F15.11,1X,I4,1X,I5,1X,I5,1X,F15.11,1X,F4.2,1X,F4.2,1X,
     +  E12.6,1X,E12.6,1X,I1,1X,I1,1X,I1,1X,F15.11,1X,I4,1X,I5,1X,I4,
     +  1X,F15.11,1X,F4.2,1X,F4.2,1X,E12.6,1X,E12.6,1X,I1,1X,I1,1X,I1,
     +  1X,F15.11,1X,I4)

C  Effective file loading

      open(unit=1,status='old',file=
     +'neb_cat.dat')
      write(6,*) '....Loading file: neb_cat.dat'
      do i__=1,5467
        read(1,'(A2449)')ar__
        read(ar__,1)
     +  Galaxy(i__),Region(i__),fEdgeHST(i__),fTouchHST(i__),
     +  fManualHST(i__),fEdgeMUSE(i__),fStarMUSE(i__),RAdeg(i__),
     +  DEdeg(i__),RApdeg(i__),DEpdeg(i__),RAMdeg(i__),DEMdeg(i__),
     +  rR25Deproj(i__),rrefFdeproj(i__),phiDeproj(i__),Env(i__),
     +  Areacirc(i__),Radcirc(i__),Radcircpc(i__),Area2(i__),
     +  amaj2(i__),bmin2(i__),PA2(i__),Rad(i__),Radpc(i__),
     +  Raddeconv(i__),RadStr(i__),RadcircpcM(i__),ComplexScore(i__),
     +  FHaHST(i__),e_FHaHST(i__),FHaHSTMax(i__),FHaHSTMin(i__),
     +  FHameanHST(i__),FHacorrHST(i__),e_FHacorrHST(i__),FHa(i__),
     +  FHb(i__),FOIII5007(i__),FNII6584(i__),FSII6716(i__),
     +  FSII6731(i__),FSIII9069(i__),FOI6300(i__),FNII5755(i__),
     +  FHeI5876(i__),FSIII6312(i__),FOI6363(i__),FOII7319(i__),
     +  FOII7330(i__),e_FHa(i__),e_FHb(i__),e_FOIII5007(i__),
     +  e_FNII6584(i__),e_FSII6716(i__),e_FSII6731(i__),
     +  e_FSIII9069(i__),e_FOI6300(i__),e_FNII5755(i__),
     +  e_FHeI5876(i__),e_FSIII6312(i__),e_FOI6363(i__),
     +  e_FOII7319(i__),e_FOII7330(i__),FHacorr(i__),FHbcorr(i__),
     +  FOIII5007corr(i__),FNII6584corr(i__),FSII6716corr(i__),
     +  FSII6731corr(i__),FSIII9069corr(i__),FOI6300corr(i__),
     +  FNII5755corr(i__),FHeI5876corr(i__),FSIII6312corr(i__),
     +  FOI6363corr(i__),FOII7319corr(i__),FOII7330corr(i__),
     +  e_FHacorr(i__),e_FHbcorr(i__),e_FOIII5007corr(i__),
     +  e_FNII6584corr(i__),e_FSII6716corr(i__),e_FSII6731corr(i__),
     +  e_FSIII9069corr(i__),e_FOI6300corr(i__),e_FNII5755corr(i__),
     +  e_FHeI5876corr(i__),e_FSIII6312corr(i__),e_FOI6363corr(i__),
     +  e_FOII7319corr(i__),e_FOII7330corr(i__),VHa(i__),e_VHa(i__),
     +  sigmaHa(i__),e_sigmaHa(i__),VNII6584(i__),e_VNII6584(i__),
     +  sigmaNII6584(i__),e_sigmaNII6584(i__),VOIII5007(i__),
     +  e_VOIII5007(i__),sigmaOIII5007(i__),e_sigmaOIII5007(i__),
     +  BPT_NII(i__),BPT_SII(i__),BPT_OI(i__),EWHa(i__),e_EWHa(i__),
     +  EWHb(i__),e_EWHb(i__),EWHaraw(i__),e_EWHaraw(i__),
     +  EWHbraw(i__),e_EWHbraw(i__),EWHafit(i__),e_EWHafit(i__),
     +  EWHbfit(i__),e_EWHbfit(i__),EWHabgcorr(i__),e_EWHabgcorr(i__),
     +  EWHbbgcorr(i__),e_EWHbbgcorr(i__),LHaHST(i__),e_LHaHST(i__),
     +  LHa(i__),HIIclassv2(i__),HIIclassv3(i__),TeNII(i__),
     +  e_TeNII(i__),TeSIII(i__),e_TeSIII(i__),neSII(i__),
     +  e_neSII(i__),q(i__),e_q(i__),logqd91(i__),e_logqd91(i__),
     +  E_B_V(i__),e_E_B_V(i__),AV(i__),MetScal(i__),e_MetScal(i__),
     +  deltaMetScal(i__),AgeANUV08pc(i__),e_AgeANUV08pc(i__),
     +  DEdegNUV08pc(i__),E_B_V_NUV08pc(i__),e_E_B_V_NUV08pc(i__),
     +  MassANUV08pc(i__),e_MassANUV08pc(i__),FlagMultNUV08p(i__),
     +  FlagNoNUV08pc(i__),FlagOneNUV08pc(i__),RAdegNUV08pc(i__),
     +  RegionIDNUV08p(i__),AgeANUV16pc(i__),e_AgeANUV16pc(i__),
     +  DEdegNUV16pc(i__),E_B_V_NUV16pc(i__),e_E_B_V_NUV16pc(i__),
     +  MassANUV16pc(i__),e_MassANUV16pc(i__),FlagMultNUV16p(i__),
     +  FlagNoNUV16pc(i__),FlagOneNUV16pc(i__),RAdegNUV16pc(i__),
     +  RegionIDNUV16p(i__),AgeANUV32pc(i__),e_AgeANUV32pc(i__),
     +  DEdegNUV32pc(i__),E_B_V_NUV32pc(i__),e_E_B_V_NUV32pc(i__),
     +  MassANUV32pc(i__),e_MassANUV32pc(i__),FlagMultNUV32p(i__),
     +  FlagNoNUV32pc(i__),FlagOneNUV32pc(i__),RAdegNUV32pc(i__),
     +  RegionIDNUV32p(i__),AgeANUV64pc(i__),e_AgeANUV64pc(i__),
     +  DEdegNUV64pc(i__),E_B_V_NUV64pc(i__),e_E_B_V_NUV64pc(i__),
     +  MassANUV64pc(i__),e_MassANUV64pc(i__),FlagMultNUV64p(i__),
     +  FlagNoNUV64pc(i__),FlagOneNUV64pc(i__),RAdegNUV64pc(i__),
     +  RegionIDNUV64p(i__),AgeAV08pc(i__),e_AgeAV08pc(i__),
     +  DEdegV08pc(i__),E_B_V_V08pc(i__),e_E_B_V_V08pc(i__),
     +  MassAV08pc(i__),e_MassAV08pc(i__),FlagMultV08pc(i__),
     +  FlagNoV08pc(i__),FlagOneV08pc(i__),RAdegV08pc(i__),
     +  RegionIDV08pc(i__),AgeAV16pc(i__),e_AgeAV16pc(i__),
     +  DEdegV16pc(i__),E_B_V_V16pc(i__),e_E_B_V_V16pc(i__),
     +  MassAV16pc(i__),e_MassAV16pc(i__),FlagMultV16pc(i__),
     +  FlagNoV16pc(i__),FlagOneV16pc(i__),RAdegV16pc(i__),
     +  RegionIDV16pc(i__),AgeAV32pc(i__),e_AgeAV32pc(i__),
     +  DEdegV32pc(i__),E_B_V_V32pc(i__),e_E_B_V_V32pc(i__),
     +  MassAV32pc(i__),e_MassAV32pc(i__),FlagMultV32pc(i__),
     +  FlagNoV32pc(i__),FlagOneV32pc(i__),RAdegV32pc(i__),
     +  RegionIDV32pc(i__),AgeAV64pc(i__),e_AgeAV64pc(i__),
     +  DEdegV64pc(i__),E_B_V_V64pc(i__),e_E_B_V_V64pc(i__),
     +  MassAV64pc(i__),e_MassAV64pc(i__),FlagMultV64pc(i__),
     +  FlagNoV64pc(i__),FlagOneV64pc(i__),RAdegV64pc(i__),
     +  RegionIDV64pc(i__)
        if (idig(ar__(257:266)).EQ.0) Raddeconv(i__) =  rNULL__
        if (idig(ar__(268:277)).EQ.0) RadStr(i__) =  rNULL__
        if (idig(ar__(489:500)).EQ.0) FNII5755(i__) =  rNULL__
        if (idig(ar__(515:526)).EQ.0) FSIII6312(i__) =  rNULL__
        if (idig(ar__(671:682)).EQ.0) e_FNII5755(i__) =  rNULL__
        if (idig(ar__(684:695)).EQ.0) e_FHeI5876(i__) =  rNULL__
        if (idig(ar__(697:708)).EQ.0) e_FSIII6312(i__) =  rNULL__
        if (idig(ar__(853:864)).EQ.0) FNII5755corr(i__) =  rNULL__
        if (idig(ar__(879:890)).EQ.0) FSIII6312corr(i__) =  rNULL__
        if (idig(ar__(905:916)).EQ.0) FOII7319corr(i__) =  rNULL__
        if (idig(ar__(918:929)).EQ.0) FOII7330corr(i__) =  rNULL__
        if (idig(ar__(957:968)).EQ.0) e_FOIII5007corr(i__) =  rNULL__
        if (idig(ar__(970:981)).EQ.0) e_FNII6584corr(i__) =  rNULL__
        if (idig(ar__(983:994)).EQ.0) e_FSII6716corr(i__) =  rNULL__
        if (idig(ar__(996:1007)).EQ.0) e_FSII6731corr(i__) =  rNULL__
        if (idig(ar__(1009:1020)).EQ.0) e_FSIII9069corr(i__) =  rNULL__
        if (idig(ar__(1022:1033)).EQ.0) e_FOI6300corr(i__) =  rNULL__
        if (idig(ar__(1035:1046)).EQ.0) e_FNII5755corr(i__) =  rNULL__
        if (idig(ar__(1048:1059)).EQ.0) e_FHeI5876corr(i__) =  rNULL__
        if (idig(ar__(1061:1072)).EQ.0) e_FSIII6312corr(i__) =  rNULL__
        if (idig(ar__(1074:1085)).EQ.0) e_FOI6363corr(i__) =  rNULL__
        if (idig(ar__(1087:1098)).EQ.0) e_FOII7319corr(i__) =  rNULL__
        if (idig(ar__(1100:1111)).EQ.0) e_FOII7330corr(i__) =  rNULL__
        if (idig(ar__(1434:1445)).EQ.0) EWHabgcorr(i__) =  rNULL__
        if (idig(ar__(1447:1458)).EQ.0) e_EWHabgcorr(i__) =  rNULL__
        if (idig(ar__(1460:1471)).EQ.0) EWHbbgcorr(i__) =  rNULL__
        if (idig(ar__(1473:1484)).EQ.0) e_EWHbbgcorr(i__) =  rNULL__
        if(ar__(1525:1526) .EQ. '') HIIclassv2(i__) = iNULL__
        if(ar__(1528:1529) .EQ. '') HIIclassv3(i__) = iNULL__
        if (idig(ar__(1544:1555)).EQ.0) e_TeNII(i__) =  rNULL__
        if (idig(ar__(1557:1568)).EQ.0) TeSIII(i__) =  rNULL__
        if (idig(ar__(1570:1581)).EQ.0) e_TeSIII(i__) =  rNULL__
        if (idig(ar__(1583:1594)).EQ.0) neSII(i__) =  rNULL__
        if (idig(ar__(1596:1607)).EQ.0) e_neSII(i__) =  rNULL__
        if (idig(ar__(1635:1646)).EQ.0) logqd91(i__) =  rNULL__
        if (idig(ar__(1648:1659)).EQ.0) e_logqd91(i__) =  rNULL__
        if (idig(ar__(1661:1672)).EQ.0) E_B_V(i__) =  rNULL__
        if (idig(ar__(1674:1685)).EQ.0) e_E_B_V(i__) =  rNULL__
        if (idig(ar__(1687:1698)).EQ.0) AV(i__) =  rNULL__
        if (idig(ar__(1700:1711)).EQ.0) MetScal(i__) =  rNULL__
        if (idig(ar__(1713:1724)).EQ.0) e_MetScal(i__) =  rNULL__
        if (idig(ar__(1726:1737)).EQ.0) deltaMetScal(i__) =  rNULL__
        if (AgeANUV08pc(i__) .EQ. 45) AgeANUV08pc(i__) =  iNULL__
        if (idig(ar__(1743:1748)).EQ.0) e_AgeANUV08pc(i__) =  rNULL__
        if (idig(ar__(1750:1757)).EQ.0) DEdegNUV08pc(i__) =  rNULL__
        if (idig(ar__(1759:1764)).EQ.0) E_B_V_NUV08pc(i__) =  rNULL__
        if (idig(ar__(1766:1771)).EQ.0) e_E_B_V_NUV08pc(i__) =  rNULL__
        if (idig(ar__(1773:1784)).EQ.0) MassANUV08pc(i__) =  rNULL__
        if (idig(ar__(1786:1797)).EQ.0) e_MassANUV08pc(i__) =  rNULL__
        if (idig(ar__(1805:1819)).EQ.0) RAdegNUV08pc(i__) =  rNULL__
        if (RegionIDNUV08p(i__) .EQ. 45) RegionIDNUV08p(i__) =  iNULL__
        if (AgeANUV16pc(i__) .EQ. 45) AgeANUV16pc(i__) =  iNULL__
        if (e_AgeANUV16pc(i__) .EQ. 45) e_AgeANUV16pc(i__) =  iNULL__
        if (idig(ar__(1833:1847)).EQ.0) DEdegNUV16pc(i__) =  rNULL__
        if (idig(ar__(1849:1852)).EQ.0) E_B_V_NUV16pc(i__) =  rNULL__
        if (idig(ar__(1854:1857)).EQ.0) e_E_B_V_NUV16pc(i__) =  rNULL__
        if (idig(ar__(1859:1870)).EQ.0) MassANUV16pc(i__) =  rNULL__
        if (idig(ar__(1872:1883)).EQ.0) e_MassANUV16pc(i__) =  rNULL__
        if (idig(ar__(1891:1905)).EQ.0) RAdegNUV16pc(i__) =  rNULL__
        if (RegionIDNUV16p(i__) .EQ. 45) RegionIDNUV16p(i__) =  iNULL__
        if (AgeANUV32pc(i__) .EQ. 45) AgeANUV32pc(i__) =  iNULL__
        if (e_AgeANUV32pc(i__) .EQ. 45) e_AgeANUV32pc(i__) =  iNULL__
        if (idig(ar__(1922:1936)).EQ.0) DEdegNUV32pc(i__) =  rNULL__
        if (idig(ar__(1938:1941)).EQ.0) E_B_V_NUV32pc(i__) =  rNULL__
        if (idig(ar__(1943:1946)).EQ.0) e_E_B_V_NUV32pc(i__) =  rNULL__
        if (idig(ar__(1948:1959)).EQ.0) MassANUV32pc(i__) =  rNULL__
        if (idig(ar__(1961:1972)).EQ.0) e_MassANUV32pc(i__) =  rNULL__
        if (idig(ar__(1980:1994)).EQ.0) RAdegNUV32pc(i__) =  rNULL__
        if (RegionIDNUV32p(i__) .EQ. 45) RegionIDNUV32p(i__) =  iNULL__
        if (AgeANUV64pc(i__) .EQ. 45) AgeANUV64pc(i__) =  iNULL__
        if (e_AgeANUV64pc(i__) .EQ. 45) e_AgeANUV64pc(i__) =  iNULL__
        if (idig(ar__(2009:2023)).EQ.0) DEdegNUV64pc(i__) =  rNULL__
        if (idig(ar__(2025:2028)).EQ.0) E_B_V_NUV64pc(i__) =  rNULL__
        if (idig(ar__(2030:2033)).EQ.0) e_E_B_V_NUV64pc(i__) =  rNULL__
        if (idig(ar__(2035:2046)).EQ.0) MassANUV64pc(i__) =  rNULL__
        if (idig(ar__(2048:2059)).EQ.0) e_MassANUV64pc(i__) =  rNULL__
        if (idig(ar__(2067:2081)).EQ.0) RAdegNUV64pc(i__) =  rNULL__
        if (RegionIDNUV64p(i__) .EQ. 45) RegionIDNUV64p(i__) =  iNULL__
        if (AgeAV08pc(i__) .EQ. 45) AgeAV08pc(i__) =  iNULL__
        if (e_AgeAV08pc(i__) .EQ. 45) e_AgeAV08pc(i__) =  iNULL__
        if (idig(ar__(2100:2114)).EQ.0) DEdegV08pc(i__) =  rNULL__
        if (idig(ar__(2116:2119)).EQ.0) E_B_V_V08pc(i__) =  rNULL__
        if (idig(ar__(2121:2124)).EQ.0) e_E_B_V_V08pc(i__) =  rNULL__
        if (idig(ar__(2126:2137)).EQ.0) MassAV08pc(i__) =  rNULL__
        if (idig(ar__(2139:2150)).EQ.0) e_MassAV08pc(i__) =  rNULL__
        if (idig(ar__(2158:2172)).EQ.0) RAdegV08pc(i__) =  rNULL__
        if (RegionIDV08pc(i__) .EQ. 45) RegionIDV08pc(i__) =  iNULL__
        if (AgeAV16pc(i__) .EQ. 45) AgeAV16pc(i__) =  iNULL__
        if (e_AgeAV16pc(i__) .EQ. 45) e_AgeAV16pc(i__) =  iNULL__
        if (idig(ar__(2191:2205)).EQ.0) DEdegV16pc(i__) =  rNULL__
        if (idig(ar__(2207:2210)).EQ.0) E_B_V_V16pc(i__) =  rNULL__
        if (idig(ar__(2212:2215)).EQ.0) e_E_B_V_V16pc(i__) =  rNULL__
        if (idig(ar__(2217:2228)).EQ.0) MassAV16pc(i__) =  rNULL__
        if (idig(ar__(2230:2241)).EQ.0) e_MassAV16pc(i__) =  rNULL__
        if (idig(ar__(2249:2263)).EQ.0) RAdegV16pc(i__) =  rNULL__
        if (RegionIDV16pc(i__) .EQ. 45) RegionIDV16pc(i__) =  iNULL__
        if (AgeAV32pc(i__) .EQ. 45) AgeAV32pc(i__) =  iNULL__
        if (e_AgeAV32pc(i__) .EQ. 45) e_AgeAV32pc(i__) =  iNULL__
        if (idig(ar__(2282:2296)).EQ.0) DEdegV32pc(i__) =  rNULL__
        if (idig(ar__(2298:2301)).EQ.0) E_B_V_V32pc(i__) =  rNULL__
        if (idig(ar__(2303:2306)).EQ.0) e_E_B_V_V32pc(i__) =  rNULL__
        if (idig(ar__(2308:2319)).EQ.0) MassAV32pc(i__) =  rNULL__
        if (idig(ar__(2321:2332)).EQ.0) e_MassAV32pc(i__) =  rNULL__
        if (idig(ar__(2340:2354)).EQ.0) RAdegV32pc(i__) =  rNULL__
        if (RegionIDV32pc(i__) .EQ. 45) RegionIDV32pc(i__) =  iNULL__
        if (AgeAV64pc(i__) .EQ. 45) AgeAV64pc(i__) =  iNULL__
        if (e_AgeAV64pc(i__) .EQ. 45) e_AgeAV64pc(i__) =  iNULL__
        if (idig(ar__(2372:2386)).EQ.0) DEdegV64pc(i__) =  rNULL__
        if (idig(ar__(2388:2391)).EQ.0) E_B_V_V64pc(i__) =  rNULL__
        if (idig(ar__(2393:2396)).EQ.0) e_E_B_V_V64pc(i__) =  rNULL__
        if (idig(ar__(2398:2409)).EQ.0) MassAV64pc(i__) =  rNULL__
        if (idig(ar__(2411:2422)).EQ.0) e_MassAV64pc(i__) =  rNULL__
        if (idig(ar__(2430:2444)).EQ.0) RAdegV64pc(i__) =  rNULL__
        if (RegionIDV64pc(i__) .EQ. 45) RegionIDV64pc(i__) =  iNULL__
c    ..............Just test output...........
        write(6,1)
     +  Galaxy(i__),Region(i__),fEdgeHST(i__),fTouchHST(i__),
     +  fManualHST(i__),fEdgeMUSE(i__),fStarMUSE(i__),RAdeg(i__),
     +  DEdeg(i__),RApdeg(i__),DEpdeg(i__),RAMdeg(i__),DEMdeg(i__),
     +  rR25Deproj(i__),rrefFdeproj(i__),phiDeproj(i__),Env(i__),
     +  Areacirc(i__),Radcirc(i__),Radcircpc(i__),Area2(i__),
     +  amaj2(i__),bmin2(i__),PA2(i__),Rad(i__),Radpc(i__),
     +  Raddeconv(i__),RadStr(i__),RadcircpcM(i__),ComplexScore(i__),
     +  FHaHST(i__),e_FHaHST(i__),FHaHSTMax(i__),FHaHSTMin(i__),
     +  FHameanHST(i__),FHacorrHST(i__),e_FHacorrHST(i__),FHa(i__),
     +  FHb(i__),FOIII5007(i__),FNII6584(i__),FSII6716(i__),
     +  FSII6731(i__),FSIII9069(i__),FOI6300(i__),FNII5755(i__),
     +  FHeI5876(i__),FSIII6312(i__),FOI6363(i__),FOII7319(i__),
     +  FOII7330(i__),e_FHa(i__),e_FHb(i__),e_FOIII5007(i__),
     +  e_FNII6584(i__),e_FSII6716(i__),e_FSII6731(i__),
     +  e_FSIII9069(i__),e_FOI6300(i__),e_FNII5755(i__),
     +  e_FHeI5876(i__),e_FSIII6312(i__),e_FOI6363(i__),
     +  e_FOII7319(i__),e_FOII7330(i__),FHacorr(i__),FHbcorr(i__),
     +  FOIII5007corr(i__),FNII6584corr(i__),FSII6716corr(i__),
     +  FSII6731corr(i__),FSIII9069corr(i__),FOI6300corr(i__),
     +  FNII5755corr(i__),FHeI5876corr(i__),FSIII6312corr(i__),
     +  FOI6363corr(i__),FOII7319corr(i__),FOII7330corr(i__),
     +  e_FHacorr(i__),e_FHbcorr(i__),e_FOIII5007corr(i__),
     +  e_FNII6584corr(i__),e_FSII6716corr(i__),e_FSII6731corr(i__),
     +  e_FSIII9069corr(i__),e_FOI6300corr(i__),e_FNII5755corr(i__),
     +  e_FHeI5876corr(i__),e_FSIII6312corr(i__),e_FOI6363corr(i__),
     +  e_FOII7319corr(i__),e_FOII7330corr(i__),VHa(i__),e_VHa(i__),
     +  sigmaHa(i__),e_sigmaHa(i__),VNII6584(i__),e_VNII6584(i__),
     +  sigmaNII6584(i__),e_sigmaNII6584(i__),VOIII5007(i__),
     +  e_VOIII5007(i__),sigmaOIII5007(i__),e_sigmaOIII5007(i__),
     +  BPT_NII(i__),BPT_SII(i__),BPT_OI(i__),EWHa(i__),e_EWHa(i__),
     +  EWHb(i__),e_EWHb(i__),EWHaraw(i__),e_EWHaraw(i__),
     +  EWHbraw(i__),e_EWHbraw(i__),EWHafit(i__),e_EWHafit(i__),
     +  EWHbfit(i__),e_EWHbfit(i__),EWHabgcorr(i__),e_EWHabgcorr(i__),
     +  EWHbbgcorr(i__),e_EWHbbgcorr(i__),LHaHST(i__),e_LHaHST(i__),
     +  LHa(i__),HIIclassv2(i__),HIIclassv3(i__),TeNII(i__),
     +  e_TeNII(i__),TeSIII(i__),e_TeSIII(i__),neSII(i__),
     +  e_neSII(i__),q(i__),e_q(i__),logqd91(i__),e_logqd91(i__),
     +  E_B_V(i__),e_E_B_V(i__),AV(i__),MetScal(i__),e_MetScal(i__),
     +  deltaMetScal(i__),AgeANUV08pc(i__),e_AgeANUV08pc(i__),
     +  DEdegNUV08pc(i__),E_B_V_NUV08pc(i__),e_E_B_V_NUV08pc(i__),
     +  MassANUV08pc(i__),e_MassANUV08pc(i__),FlagMultNUV08p(i__),
     +  FlagNoNUV08pc(i__),FlagOneNUV08pc(i__),RAdegNUV08pc(i__),
     +  RegionIDNUV08p(i__),AgeANUV16pc(i__),e_AgeANUV16pc(i__),
     +  DEdegNUV16pc(i__),E_B_V_NUV16pc(i__),e_E_B_V_NUV16pc(i__),
     +  MassANUV16pc(i__),e_MassANUV16pc(i__),FlagMultNUV16p(i__),
     +  FlagNoNUV16pc(i__),FlagOneNUV16pc(i__),RAdegNUV16pc(i__),
     +  RegionIDNUV16p(i__),AgeANUV32pc(i__),e_AgeANUV32pc(i__),
     +  DEdegNUV32pc(i__),E_B_V_NUV32pc(i__),e_E_B_V_NUV32pc(i__),
     +  MassANUV32pc(i__),e_MassANUV32pc(i__),FlagMultNUV32p(i__),
     +  FlagNoNUV32pc(i__),FlagOneNUV32pc(i__),RAdegNUV32pc(i__),
     +  RegionIDNUV32p(i__),AgeANUV64pc(i__),e_AgeANUV64pc(i__),
     +  DEdegNUV64pc(i__),E_B_V_NUV64pc(i__),e_E_B_V_NUV64pc(i__),
     +  MassANUV64pc(i__),e_MassANUV64pc(i__),FlagMultNUV64p(i__),
     +  FlagNoNUV64pc(i__),FlagOneNUV64pc(i__),RAdegNUV64pc(i__),
     +  RegionIDNUV64p(i__),AgeAV08pc(i__),e_AgeAV08pc(i__),
     +  DEdegV08pc(i__),E_B_V_V08pc(i__),e_E_B_V_V08pc(i__),
     +  MassAV08pc(i__),e_MassAV08pc(i__),FlagMultV08pc(i__),
     +  FlagNoV08pc(i__),FlagOneV08pc(i__),RAdegV08pc(i__),
     +  RegionIDV08pc(i__),AgeAV16pc(i__),e_AgeAV16pc(i__),
     +  DEdegV16pc(i__),E_B_V_V16pc(i__),e_E_B_V_V16pc(i__),
     +  MassAV16pc(i__),e_MassAV16pc(i__),FlagMultV16pc(i__),
     +  FlagNoV16pc(i__),FlagOneV16pc(i__),RAdegV16pc(i__),
     +  RegionIDV16pc(i__),AgeAV32pc(i__),e_AgeAV32pc(i__),
     +  DEdegV32pc(i__),E_B_V_V32pc(i__),e_E_B_V_V32pc(i__),
     +  MassAV32pc(i__),e_MassAV32pc(i__),FlagMultV32pc(i__),
     +  FlagNoV32pc(i__),FlagOneV32pc(i__),RAdegV32pc(i__),
     +  RegionIDV32pc(i__),AgeAV64pc(i__),e_AgeAV64pc(i__),
     +  DEdegV64pc(i__),E_B_V_V64pc(i__),e_E_B_V_V64pc(i__),
     +  MassAV64pc(i__),e_MassAV64pc(i__),FlagMultV64pc(i__),
     +  FlagNoV64pc(i__),FlagOneV64pc(i__),RAdegV64pc(i__),
     +  RegionIDV64pc(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C  Loading file 'list.dat'	! List of fits maps and masks

C  Format for file interpretation

    2 format(F9.5,F9.5,1X,I5,1X,I4,1X,I6,1X,A37,1X,A36)

C  Effective file loading

      open(unit=1,status='old',file=
     +'list.dat')
      write(6,*) '....Loading file: list.dat'
      do i__=1,38
        read(1,'(A111)')ar__1
        read(ar__1,2)
     +  RAdeg_1(i__),DEdeg_1(i__),Nx(i__),Ny(i__),size(i__),
     +  FileName(i__),Title(i__)
c    ..............Just test output...........
        write(6,2)
     +  RAdeg_1(i__),DEdeg_1(i__),Nx(i__),Ny(i__),size(i__),
     +  FileName(i__),Title(i__)
c    .......End.of.Just test output...........
      end do
      close(1)

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

C Locate position of first digit in string; or return 0
      integer function idig(c)
      character*(*) c
      character*1 c1
      integer lc,i
      lc=len(c)
      idig=0
      do i=1,lc
         if(c(i:i).ne.' ') go to 1
      end do
    1 if(i.gt.lc) return
      c1=c(i:i)
      if(c1.eq.'.'.or.c1.eq.'-'.or.c1.eq.'+') i=i+1
      if(i.gt.lc) return
      c1=c(i:i)
      if(c1.ge.'0'.and.c1.le.'9') idig=i
      return
      end