Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/706/A95 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-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 ! Galaxy name (galaxy_name)
integer*4 Region ! Region ID (region_ID)
integer*4 fEdgeHST ! Flag=1 if region lies on HST image
* edge (flag_edge_hst)
integer*4 fTouchHST ! Flag=1 if region touches segmentation
* boundary (flag_touch_hst)
character*7 fManualHST ! ?=- Manual classification flags
* (visual checks) (flag_manual_hst)
integer*4 fEdgeMUSE ! Flag=1 if region lies on MUSE cube
* edge (flag_edge_muse)
integer*4 fStarMUSE ! Flag=1 if affected by bright star in
* MUSE data (flag_star_muse)
real*8 RAdeg ! (deg) Right ascension (J2000) of region
* centroid (ra)
real*8 DEdeg ! (deg) Declination (J2000) of region
* centroid (dec)
real*8 RApdeg ! (deg) H{alpha} pixel peak right ascension
* (J2000) (ra_peak)
real*8 DEpdeg ! (deg) H{alpha} pixel peak declination
* (J2000) (dec_peak)
real*8 RAMdeg ! (deg) MUSE right ascension (J2000) of
* MUSE centroid (ra_muse)
real*8 DEMdeg ! (deg) MUSE declination (J2000) of MUSE
* centroid (dec_muse)
real*8 rR25Deproj ! Deprojected galactocentric radius,
* in R25 units (r_R25_deproj)
real*8 rrefFdeproj ! Deprojected galactocentric radius,
* in Re units (r_reff_deproj)
real*8 phiDeproj ! (rad) Deprojected azimuthal angle within
* host galaxy (phi_deproj)
integer*4 Env ! [1/9] Environment flag (1)
real*8 Areacirc ! (arcsec+2) Circularized area (area_circ_arcsec2)
real*8 Radcirc ! (arcsec) Circularized radius (arcsec)
* (radius_circ_arcsec)
real*8 Radcircpc ! (pc) Circularized radius (pc)
* (radius_circ_pc)
real*8 Area2 ! (arcsec+2) Second-moment area of region
* (area_mom_arcsec2)
real*8 amaj2 ! (arcsec) Second-moment major axis
* (radius_major_mom_arcsec)
real*8 bmin2 ! (arcsec) Second-moment minor axis
* (radius_minor_mom_arcsec)
real*8 PA2 ! (deg) [] Position angle (moments-based)
* (radius_pos_angle_deg)
real*8 Rad ! (arcsec) Moment-based radius (arcsec)
* (radius_mom_arcsec)
real*8 Radpc ! (pc) Moment-based radius (pc)
* (radius_mom_pc)
real*8 Raddeconv ! (pc) ?=- Deconvolved moment-based radius
* (radius_mom_deconv_pc)
real*8 RadStr ! (pc) ?=- Stromgren radius
* (radius_stromgren_pc)
real*8 RadcircpcM ! (pc) Circularized radius (pc; PHANGS-MUSE)
* (radius_circ_pc_muse)
integer*4 ComplexScore ! Morphological complexity score
* (Sect. 3.2) (complexity_score)
real*4 FHaHST ! (mW/m2) HST H{alpha} integrated flux (raw,
* uncorrected) (halpha_flux)
real*4 e_FHaHST ! (mW/m2) Error on HST H{alpha} flux
* (halpha_flux_err)
real*4 FHaHSTMax ! (mW/m2) Maximum pixel flux in HST H{alpha} map
* (halpha_flux_max)
real*4 FHaHSTMin ! (mW/m2) Minimum pixel flux in HST H{alpha} map
* (halpha_flux_min)
real*4 FHameanHST ! (mW/m2) Mean pixel flux in HST H{alpha} map
* (halpha_flux_mean)
real*4 FHacorrHST ! (mW/m2) Extinction-corrected HST H{alpha} flux
* (halpha_flux_corr)
real*4 e_FHacorrHST ! (mW/m2) Error on extinction-corrected HST
* H{alpha} flux (halpha_flux_corr_err)
real*4 FHa ! (mW/m2) MUSE H{alpha} integrated flux
* (PHANGS-MUSE) (halpha_flux_muse)
real*4 FHb ! (mW/m2) MUSE H{beta} integrated flux
* (PHANGS-MUSE) (hbeta_flux_muse)
real*4 FOIII5007 ! (mW/m2) MUSE [O III] {lambda}5007 flux
* (PHANGS-MUSE) (oiii5007_flux_muse)
real*4 FNII6584 ! (mW/m2) MUSE [N II] {lambda}6584 flux
* (PHANGS-MUSE) (nii6584_flux_muse)
real*4 FSII6716 ! (mW/m2) MUSE [S II] {lambda}6716 flux
* (PHANGS-MUSE) (sii6716_flux_muse)
real*4 FSII6731 ! (mW/m2) MUSE [S II] {lambda}6731 flux
* (PHANGS-MUSE) (sii6731_flux_muse)
real*4 FSIII9069 ! (mW/m2) MUSE [S III] {lambda}9069 flux
* (PHANGS-MUSE) (siii9069_flux_muse)
real*4 FOI6300 ! (mW/m2) MUSE [O I] {lambda}6300 flux
* (PHANGS-MUSE) (oi6300_flux_muse)
real*4 FNII5755 ! (mW/m2) ?=- MUSE [N II] {lambda}5755 flux
* (PHANGS-MUSE) (nii5755_flux_muse)
real*4 FHeI5876 ! (mW/m2) MUSE He I {lambda}5876 flux
* (PHANGS-MUSE) (hei5876_flux_muse)
real*4 FSIII6312 ! (mW/m2) ?=- MUSE [S III] {lambda}6312 flux
* (PHANGS-MUSE) (siii6312_flux_muse)
real*4 FOI6363 ! (mW/m2) MUSE [O I] {lambda}6363 flux
* (PHANGS-MUSE) (oi6363_flux_muse)
real*4 FOII7319 ! (mW/m2) MUSE [O II] {lambda}7319 flux
* (PHANGS-MUSE) (oii7319_flux_muse)
real*4 FOII7330 ! (mW/m2) MUSE [O II] {lambda}7330 flux
* (PHANGS-MUSE) (oii7330_flux_muse)
real*4 e_FHa ! (mW/m2) Error on H{alpha} flux (PHANGS-MUSE)
* (halpha_flux_err_muse)
real*4 e_FHb ! (mW/m2) Error on H{beta} flux (PHANGS-MUSE)
* (hbeta_flux_err_muse)
real*4 e_FOIII5007 ! (mW/m2) Error on [O III] {lambda}5007 flux
* (PHANGS-MUSE)
* (oiii5007_flux_err_muse)
real*4 e_FNII6584 ! (mW/m2) Error on [N II] {lambda}6584 flux
* (PHANGS-MUSE) (nii6584_flux_err_muse)
real*4 e_FSII6716 ! (mW/m2) Error on [S II] {lambda}6716 flux
* (PHANGS-MUSE) (sii6716_flux_err_muse)
real*4 e_FSII6731 ! (mW/m2) Error on [S II] {lambda}6731 flux
* (PHANGS-MUSE) (sii6731_flux_err_muse)
real*4 e_FSIII9069 ! (mW/m2) Error on [S III] {lambda}9069 flux
* (PHANGS-MUSE)
* (siii9069_flux_err_muse)
real*4 e_FOI6300 ! (mW/m2) Error on [O I] {lambda}6300 flux
* (PHANGS-MUSE) (oi6300_flux_err_muse)
real*4 e_FNII5755 ! (mW/m2) ?=- Error on [N II] {lambda}5755 flux
* (PHANGS-MUSE) (nii5755_flux_err_muse)
real*4 e_FHeI5876 ! (mW/m2) ?=- Error on He I {lambda}5876 flux
* (PHANGS-MUSE) (hei5876_flux_err_muse)
real*4 e_FSIII6312 ! (mW/m2) ?=- Error on [Si II] {lambda}6312 flux
* (PHANGS-MUSE)
* (siii6312_flux_err_muse)
real*4 e_FOI6363 ! (mW/m2) Error on [O I] {lambda}6363 flux
* (PHANGS-MUSE) (oi6363_flux_err_muse)
real*4 e_FOII7319 ! (mW/m2) Error on [O I] {lambda}7319 flux
* (PHANGS-MUSE) (oii7319_flux_err_muse)
real*4 e_FOII7330 ! (mW/m2) Error on [O I] {lambda}7319 flux
* (PHANGS-MUSE) (oii7330_flux_err_muse)
real*4 FHacorr ! (mW/m2) Extinction-corrected H{alpha} flux
* (PHANGS-MUSE) (halpha_flux_corr_muse)
real*4 FHbcorr ! (mW/m2) Extinction-corrected H{beta} flux
* (PHANGS-MUSE) (hbeta_flux_corr_muse)
real*4 FOIII5007corr ! (mW/m2) Extinction-corrected [O III]
* {lambda}5007 flux (PHANGS-MUSE)
* (oiii5007_flux_corr_muse)
real*4 FNII6584corr ! (mW/m2) Extinction-corrected [N II]
* {lambda}6584 flux (PHANGS-MUSE)
* (nii6584_flux_corr_muse)
real*4 FSII6716corr ! (mW/m2) Extinction-corrected [S II]
* {lambda}6716 flux (PHANGS-MUSE)
* (sii6716_flux_corr_muse)
real*4 FSII6731corr ! (mW/m2) Extinction-corrected [S II]
* {lambda}6731 flux (PHANGS-MUSE)
* (sii6731_flux_corr_muse)
real*4 FSIII9069corr ! (mW/m2) Extinction-corrected [S III]
* {lambda}9069 flux (PHANGS-MUSE)
* (siii9069_flux_corr_muse)
real*4 FOI6300corr ! (mW/m2) Extinction-corrected [O I]
* {lambda}6300 flux (PHANGS-MUSE)
* (oi6300_flux_corr_muse)
real*4 FNII5755corr ! (mW/m2) ?=- Extinction-corrected [N II]
* {lambda}5755 flux (PHANGS-MUSE)
* (nii5755_flux_corr_muse)
real*4 FHeI5876corr ! (mW/m2) Extinction-corrected He I
* {lambda}5876 flux (PHANGS-MUSE)
* (hei5876_flux_corr_muse)
real*4 FSIII6312corr ! (mW/m2) ?=- Extinction-corrected [S III]
* {lambda}6312 flux (PHANGS-MUSE)
* (siii6312_flux_corr_muse)
real*4 FOI6363corr ! (mW/m2) Extinction-corrected [O I]
* {lambda}6363 flux (PHANGS-MUSE)
* (oi6363_flux_corr_muse)
real*4 FOII7319corr ! (mW/m2) ?=- Extinction-corrected [O II]
* {lambda}7319 flux (PHANGS-MUSE)
* (oii7319_flux_corr_muse)
real*4 FOII7330corr ! (mW/m2) ?=- Extinction-corrected [O II]
* {lambda}7330 flux (PHANGS-MUSE)
* (oii7330_flux_corr_muse)
real*4 e_FHacorr ! (mW/m2) Error on corrected H{alpha} flux
* (PHANGS-MUSE)
* (halpha_flux_corr_err_muse)
real*4 e_FHbcorr ! (mW/m2) Error on corrected H{beta} flux
* (PHANGS-MUSE)
* (hbeta_flux_corr_err_muse)
real*4 e_FOIII5007corr ! (mW/m2) ?=- Error on corrected [O III]
* {lambda}5007 flux (PHANGS-MUSE)
* (oiii5007_fl ux_corr_err_muse)
real*4 e_FNII6584corr ! (mW/m2) ?=- Error on corrected [N II]
* {lambda}5007 flux (PHANGS-MUSE)
* (nii6584_flux_corr_err_muse)
real*4 e_FSII6716corr ! (mW/m2) ?=- Error on corrected [S II]
* {lambda}6716 flux (PHANGS-MUSE)
* (sii6716_flux_corr_err_muse)
real*4 e_FSII6731corr ! (mW/m2) ?=- Error on corrected [S II]
* {lambda}6731 flux (PHANGS-MUSE)
* (sii6731_flux_corr_err_muse)
real*4 e_FSIII9069corr ! (mW/m2) ?=- Error on corrected [S III]
* {lambda}9069 flux (PHANGS-MUSE)
* (siii9069_fl ux_corr_err_muse)
real*4 e_FOI6300corr ! (mW/m2) ?=- Error on corrected [O I]
* {lambda}6300 flux (PHANGS-MUSE)
* (oi6300_flux_corr_err_muse)
real*4 e_FNII5755corr ! (mW/m2) ?=- Error on corrected [N II]
* {lambda}5755 flux (PHANGS-MUSE)
* (nii5755_flux_corr_err_muse)
real*4 e_FHeI5876corr ! (mW/m2) ?=- Error on corrected He I
* {lambda}5876 flux (PHANGS-MUSE)
* (hei5876_flux_corr_err_muse)
real*4 e_FSIII6312corr ! (mW/m2) ?=- Error on corrected [S III]
* {lambda}6312 flux (PHANGS-MUSE)
* (siii6312_fl ux_corr_err_muse)
real*4 e_FOI6363corr ! (mW/m2) ?=- Error on corrected [O I]
* {lambda}6363 flux (PHANGS-MUSE)
* (oi6363_flux_corr_err_muse)
real*4 e_FOII7319corr ! (mW/m2) ?=- Error on corrected [O II]
* {lambda}7319 flux (PHANGS-MUSE)
* (oii7319_flux_corr_err_muse)
real*4 e_FOII7330corr ! (mW/m2) ?=- Error on corrected [O II]
* {lambda}7330 flux (PHANGS-MUSE)
* (oii7330_flux_corr_err_muse)
real*4 VHa ! (km/s) H{alpha} velocity (PHANGS-MUSE)
* (halpha_vel_muse)
real*4 e_VHa ! (km/s) Error on H{alpha} velocity
* (PHANGS-MUSE) (halpha_vel_err_muse)
real*4 sigmaHa ! (km/s) H{alpha} velocity dispersion
* (PHANGS-MUSE) (halpha_sigma_muse)
real*4 e_sigmaHa ! (km/s) Error on H{alpha} velocity dispersion
* (PHANGS-MUSE) (halpha_sigma_err_muse)
real*4 VNII6584 ! (km/s) Velocity on [N II] {lambda}6584
* (PHANGS-MUSE) (nii6584_vel_muse)
real*4 e_VNII6584 ! (km/s) Error on velocity of [N II]
* {lambda}6584 (PHANGS-MUSE)
* (nii6584_vel_err_muse)
real*4 sigmaNII6584 ! (km/s) Velocity dispersion of [N II]
* {lambda}6584 (PHANGS-MUSE)
* (nii6584_sigma_muse)
real*4 e_sigmaNII6584 ! (km/s) Error on velocity dispersion of [N II]
* {lambda}6584 (PHANGS-MUSE)
* (nii6584_sigma_err_muse)
real*4 VOIII5007 ! (km/s) Velocity on [O III] {lambda}5007
* (PHANGS-MUSE) (oiii5007_vel_muse)
real*4 e_VOIII5007 ! (km/s) Error on velocity on [O III]
* {lambda}5007 (PHANGS-MUSE)
* (oiii5007_vel_err_muse)
real*4 sigmaOIII5007 ! (km/s) Velocity dispersion og [O III]
* {lambda}5007 (PHANGS-MUSE)
* (oiii5007_sigma_muse)
real*4 e_sigmaOIII5007 ! (km/s) Error on velocity dispersion of
* [O III] {lambda}5007
* (oiii5007_sigma_err_muse)
integer*4 BPT_NII ! [-1/3] BPT class using
* [N II]/H{alpha} (bpt_nii_muse)
integer*4 BPT_SII ! [-1/3] BPT class using
* [S II]/H{alpha} (bpt_sii_muse)
integer*4 BPT_OI ! [-1/3] BPT class using
* [O I]/H{alpha} (bpt_oi_muse)
real*4 EWHa ! (0.1nm) H{alpha} equivalent width
* (PHANGS-MUSE) (halpha_ew_muse)
real*4 e_EWHa ! (0.1nm) [] Uncertainty on H{alpha} equivalent
* width (PHANGS-MUSE)
* (halpha_ew_err_muse)
real*4 EWHb ! (0.1nm) H{beta} equivalent width (PHANGS-MUSE)
* (hbeta_ew_muse)
real*4 e_EWHb ! (0.1nm) [] Uncertainty on H{beta} equivalent
* width (PHANGS-MUSE)
* (hbeta_ew_err_muse)
real*4 EWHaraw ! (0.1nm) H{alpha} equivalent width from raw
* spectrum (no continuum model,
* PHANGS-MUSE) (halpha_ew_raw_muse)
real*4 e_EWHaraw ! (0.1nm) [] Uncertainty on H{alpha} equivalent
* width from raw spectrum (no continuum
* model, PHANGS-MUSE)
* (halpha_ew_raw_err_muse)
real*4 EWHbraw ! (0.1nm) H{beta} equivalent width from raw
* spectrum (no continuum model,
* PHANGS-MUSE) (hbeta_ew_raw_muse)
real*4 e_EWHbraw ! (0.1nm) [] Uncertainty on H{beta} equivalent
* width from raw spectrum (no continuum
* model, PHANGS-MUSE)
* (hbeta_ew_raw_err_muse)
real*4 EWHafit ! (0.1nm) H{alpha} equivalent width fitted model
* (PHANGS-MUSE) (halpha_ew_fit_muse)
real*4 e_EWHafit ! (0.1nm) [] Uncertainty on H{alpha} equivalent
* width fitted model (PHANGS-MUSE)
* (halpha_ew_fit_err_muse)
real*4 EWHbfit ! (0.1nm) H{beta} equivalent width fitted model
* (PHANGS-MUSE) (hbeta_ew_fit_muse)
real*4 e_EWHbfit ! (0.1nm) [] Uncertainty on H{beta} equivalent
* width fitted model (PHANGS-MUSE)
* (hbeta_ew_fit_err_muse)
real*4 EWHabgcorr ! (0.1nm) ?=- Background-corrected H{beta}
* equivalent width (PHANGS-MUSE)
* (halpha_ew_bgcorr_muse)
real*4 e_EWHabgcorr ! (0.1nm) []?=- Uncertainty on
* background-corrected H{beta}
* equivalent width (PHANGS-MUSE)
* (halpha_ew_bgcorr_err_muse)
real*4 EWHbbgcorr ! (0.1nm) ?=- Background-corrected H{beta}
* equivalent width (PHANGS-MUSE)
* (hbeta_ew_bgcorr_muse)
real*4 e_EWHbbgcorr ! (0.1nm) []?=- Uncertainty on
* background-corrected H{beta}
* equivalent width (PHANGS-MUSE)
* (hbeta_ew_bgcorr_err_muse)
real*4 LHaHST ! (10-7W) H{alpha} luminosity from HST
* (halpha_lum)
real*4 e_LHaHST ! (10-7W) Uncertainty on H{alpha} luminosity
* from HST (halpha_lum_err)
real*4 LHa ! (10-7W) H{alpha} luminosity from MUSE
* (halpha_lum_muse)
integer*4 HIIclassv2 ! [0/1]? PHANGS-MUSE H II classification
* (v2; PHANGS-MUSE) (hii_class_v2_muse)
integer*4 HIIclassv3 ! [0/1]? PHANGS-MUSE H II classification
* (v3; PHANGS-MUSE) (hii_class_v3_muse)
real*4 TeNII ! (K) Electron temperature from [N II]
* (PHANGS-MUSE) (te_nii_muse)
real*4 e_TeNII ! (K) ?=- Uncertainty on electron
* temperature from [N II] (PHANGS-MUSE)
* (te_nii_err_muse)
real*4 TeSIII ! (K) ?=- Electron temperature from [S III]
* (PHANGS-MUSE) (te_siii_muse)
real*4 e_TeSIII ! (K) ?=- Uncertainty on electron
* temperature from [S III]
* (PHANGS-MUSE) (te_siii_err_muse)
real*4 neSII ! (cm-3) ?=- Electron density from
* [S II]6716/6731 (PHANGS-MUSE)
* (ne_sii_muse)
real*4 e_neSII ! (cm-3) ?=- Uncertainty on electron density
* from [S II]6716/6731 (PHANGS-MUSE)
* (ne_sii_err_muse)
real*4 q ! (s-1) Ionizing photon rate Q(H^0) (q)
real*4 e_q ! (s-1) Uncertainty on ionizing photon rate
* Q(H^0) (q)
real*4 logqd91 ! ?=- log10 ionization parameter from
* Diaz+1991 calibration (PHANGS-MUSE)
* (logq_d91_muse)
real*4 e_logqd91 ! ?=- Uncertainty on log10 ionization
* parameter (PHANGS-MUSE)
* (logq_d91_err_muse)
real*4 E_B_V ! (mag) ?=- Color excess E(B-V) from Balmer
* decrement (PHANGS-MUSE) (ebv_muse)
real*4 e_E_B_V ! (mag) ?=- Uncertainty on color excess E(B-V)
* from Balmer decrement (PHANGS-MUSE)
* (ebv_err_muse)
real*4 AV ! (mag) ?=- Visual extinction A_V
* (PHANGS-MUSE) (av_muse)
real*4 MetScal ! ?=- Gas-phase metallicity (Scal code;
* PHANGS-MUSE) (metallicity_scal_muse)
real*4 e_MetScal ! ?=- Uncertainty on gas-phase
* metallicity (Scal code; PHANGS-MUSE)
* (metallicity_scal_err_muse)
real*4 deltaMetScal ! ?=- Offset from local metallicity
* trend (PHANGS-MUSE)
* (metallicity_scal_delta_ muse)
integer*4 AgeANUV08pc ! (Myr) ?=- Stellar association age
* (NUV;8 pc; PHANGS-HST)
* (age_assoc_nuv_08pc)
real*4 e_AgeANUV08pc ! (Myr) ?=- Uncertainty on stellar association
* age (NUV;8 pc; PHANGS-HST)
* (age_err_assoc_nuv_08pc)
real*8 DEdegNUV08pc ! (deg) ?=- Association declination
* (NUV;8 pc; PHANGS-HST)
* (dec_assoc_nuv_08pc)
real*4 E_B_V_NUV08pc ! (mag) ?=- Color excess E(B-V)
* (NUV; 8 pc; PHANGS-HST)
* (ebv_assoc_nuv_08pc)
real*4 e_E_B_V_NUV08pc ! (mag) ?=- Uncertainty on color excess E(B-V)
* (NUV; 8 pc; PHANGS-HST)
* (ebv_err_assoc_nuv_08pc)
real*4 MassANUV08pc ! (Msun) ?=- Stellar mass (NUV; 8 pc;
* PHANGS-HST) (mass_assoc_nuv_08pc)
real*4 e_MassANUV08pc ! (Msun) ?=- Uncertainty on stellar mass
* (NUV; 8 pc; PHANGS-HST)
* (mass_err_assoc_nuv_08pc)
integer*4 FlagMultNUV08p ! [0/1] Flag: multiple associations
* (NUV; 8 pc; PHANGS-HST)
* (flag_multi_assoc_nuv_08pc)
integer*4 FlagNoNUV08pc ! [0/1] Flag: no association
* (NUV; 8 pc; PHANGS-HST)
* (flag_none_assoc_nuv_08pc)
integer*4 FlagOneNUV08pc ! [0/1] Flag: single association
* (NUV; 8 pc; PHANGS-HST)
* (flag_one_assoc_nuv_08pc)
real*8 RAdegNUV08pc ! (deg) ?=- Association right ascension
* (NUV; 8 pc; PHANGS-HST)
* (ra_assoc_nuv_08pc)
integer*4 RegionIDNUV08p ! ?=- Association identifier
* (NUV; 8 pc; PHANGS-HST)
* (region_ID_assoc_nuv_08pc)
integer*4 AgeANUV16pc ! (Myr) ?=- Association age
* (NUV; 16 pc; PHANGS-HST)
* (age_assoc_nuv_16pc)
integer*4 e_AgeANUV16pc ! (Myr) ?=- Uncertainty on association age
* (NUV; 16 pc; PHANGS-HST)
* (age_err_assoc_nuv_16pc)
real*8 DEdegNUV16pc ! (deg) ?=- Association declination
* (NUV; 16 pc; PHANGS-HST)
* (dec_assoc_nuv_16pc)
real*4 E_B_V_NUV16pc ! (mag) ?=- E(B-V) of association
* (NUV; 16 pc; PHANGS-HST)
* (ebv_assoc_nuv_16pc)
real*4 e_E_B_V_NUV16pc ! (mag) ?=- Uncertainty on E(B-V) of
* association (NUV; 16 pc; PHANGS-HST)
* (ebv_err_assoc_nuv_16pc)
real*4 MassANUV16pc ! (Msun) ?=- Stellar mass (NUV; 16 pc;
* PHANGS-HST) (mass_assoc_nuv_16pc)
real*4 e_MassANUV16pc ! (Msun) ?=- Uncertainty on stellar mass
* (NUV; 16 pc; PHANGS-HST)
* (mass_err_assoc_nuv_16pc)
integer*4 FlagMultNUV16p ! [0/1] Flag: multiple association
* (NUV; 16 pc; PHANGS-HST)
* (flag_multi_assoc_nuv_16pc)
integer*4 FlagNoNUV16pc ! [0/1] Flag: no association
* (NUV; 16 pc; PHANGS-HST)
* (flag_none_assoc_nuv_16pc)
integer*4 FlagOneNUV16pc ! [0/1] Flag: single association
* (NUV; 16 pc; PHANGS-HST)
* (flag_one_assoc_nuv_16pc)
real*8 RAdegNUV16pc ! (deg) ?=- Association right ascension
* (NUV; 16 pc; PHANGS-HST)
* (ra_assoc_nuv_16pc)
integer*4 RegionIDNUV16p ! ?=- Association identifier
* (NUV; 16 pc; PHANGS-HST)
* (region_ID_assoc_nuv_16pc)
integer*4 AgeANUV32pc ! (Myr) ?=- Association age (NUV; 32 pc;
* PHANGS-HST) (age_assoc_nuv_32pc)
integer*4 e_AgeANUV32pc ! (Myr) ?=- Uncertainty on association age
* (NUV; 32 pc; PHANGS-HST)
* (age_err_assoc_nuv_32pc)
real*8 DEdegNUV32pc ! (deg) ?=- Association declination
* (NUV;32 pc; PHANGS-HST)
* (dec_assoc_nuv_32pc)
real*4 E_B_V_NUV32pc ! (mag) ?=- E(B-V) of association
* (NUV; 32 pc; PHANGS-HST)
* (ebv_assoc_nuv_32pc)
real*4 e_E_B_V_NUV32pc ! (mag) ?=- Uncertainty on E(B-V) of
* association (NUV; 32 pc; PHANGS-HST)
* (ebv_err_assoc_nuv_32pc)
real*4 MassANUV32pc ! (Msun) ?=- Stellar mass (NUV; 32 pc;
* PHANGS-HST) (mass_assoc_nuv_32pc)
real*4 e_MassANUV32pc ! (Msun) ?=- Uncertainty on stellar mass
* (NUV; 32 pc; PHANGS-HST)
* (mass_err_assoc_nuv_32pc)
integer*4 FlagMultNUV32p ! [0/1] Flag: multiple associations
* (NUV; 32 pc; PHANGS-HST)
* (flag_multi_assoc_nuv_32pc)
integer*4 FlagNoNUV32pc ! [0/1] Flag: no associations
* (NUV; 32 pc; PHANGS-HST)
* (flag_none_assoc_nuv_32pc)
integer*4 FlagOneNUV32pc ! [0/1] Flag: single associations
* (NUV; 32 pc; PHANGS-HST)
* (flag_one_assoc_nuv_32pc)
real*8 RAdegNUV32pc ! (deg) ?=- Association right ascension
* (NUV; 32 pc; PHANGS-HST)
* (ra_assoc_nuv_32pc)
integer*4 RegionIDNUV32p ! ?=- Association identifier
* (NUV; 32 pc; PHANGS-HST)
* (region_ID_assoc_nuv_32pc)
integer*4 AgeANUV64pc ! (Myr) ?=- Association age (NUV; 64 pc;
* PHANGS-HST) (age_assoc_nuv_64pc)
integer*4 e_AgeANUV64pc ! (Myr) ?=- Uncertainty on association age
* (NUV; 64 pc; PHANGS-HST)
* (age_err_assoc_nuv_64pc)
real*8 DEdegNUV64pc ! (deg) ?=- Association declination
* (NUV; 64 pc; PHANGS-HST)
* (dec_assoc_nuv_64pc)
real*4 E_B_V_NUV64pc ! (mag) ?=- E(B-V) of association
* (NUV; 64 pc; PHANGS-HST)
* (ebv_assoc_nuv_64pc)
real*4 e_E_B_V_NUV64pc ! (mag) ?=- Uncertainty on E(B-V) of
* association (NUV; 64 pc; PHANGS-HST)
* (ebv_err_assoc_nuv_64pc)
real*4 MassANUV64pc ! (Msun) ?=- Stellar mass (NUV; 64 pc;
* PHANGS-HST; PHANGS-HST)
* (mass_assoc_nuv_64pc)
real*4 e_MassANUV64pc ! (Msun) ?=- Uncertainty on stellar mass
* (NUV; 64 pc; PHANGS-HST; PHANGS-HST)
* (mass_err_assoc_nuv_64pc)
integer*4 FlagMultNUV64p ! [0/1] Flag: multiple associations
* (NUV; 64 pc; PHANGS-HST)
* (flag_multi_assoc_nuv_64pc)
integer*4 FlagNoNUV64pc ! [0/1] Flag: no associations
* (NUV; 64 pc; PHANGS-HST)
* (flag_none_assoc_nuv_64pc)
integer*4 FlagOneNUV64pc ! [0/1] Flag: single associations
* (NUV; 64 pc; PHANGS-HST)
* (flag_one_assoc_nuv_64pc)
real*8 RAdegNUV64pc ! (deg) ?=- Association right ascension
* (NUV; 64 pc; PHANGS-HST)
* (ra_assoc_nuv_64pc)
integer*4 RegionIDNUV64p ! ?=- Association identifier
* (NUV; 64 pc; PHANGS-HST)
* (region_ID_assoc_nuv_64pc)
integer*4 AgeAV08pc ! (Myr) ?=- Stellar association age
* (V;8 pc; PHANGS-HST)
* (age_assoc_v_08pc)
integer*4 e_AgeAV08pc ! (Myr) ?=- Uncertainty on stellar association
* age (V;8 pc; PHANGS-HST)
* (age_err_assoc_v_08pc)
real*8 DEdegV08pc ! (deg) ?=- Association declination
* (V;8 pc; PHANGS-HST)
* (dec_assoc_v_08pc)
real*4 E_B_V_V08pc ! (mag) ?=- Color excess E(B-V)
* (V; 8 pc; PHANGS-HST)
* (ebv_assoc_v_08pc)
real*4 e_E_B_V_V08pc ! (mag) ?=- Uncertainty on color excess E(B-V)
* (V; 8 pc; PHANGS-HST)
* (ebv_err_assoc_v_08pc)
real*4 MassAV08pc ! (Msun) ?=- Stellar mass (V; 8 pc; PHANGS-HST)
* (mass_assoc_v_08pc)
real*4 e_MassAV08pc ! (Msun) ?=- Uncertainty on stellar mass
* (V; 8 pc; PHANGS-HST)
* (mass_err_assoc_v_08pc)
integer*4 FlagMultV08pc ! [0/1] Flag: multiple associations
* (V; 8 pc; PHANGS-HST)
* (flag_multi_assoc_v_08pc)
integer*4 FlagNoV08pc ! [0/1] Flag: no association
* (V; 8 pc; PHANGS-HST)
* (flag_none_assoc_v_08pc)
integer*4 FlagOneV08pc ! [0/1] Flag: single association
* (V; 8 pc; PHANGS-HST)
* (flag_one_assoc_v_08pc)
real*8 RAdegV08pc ! (deg) ?=- Association right ascension
* (V; 8 pc; PHANGS-HST)
* (ra_assoc_v_08pc)
integer*4 RegionIDV08pc ! ?=- Association identifier
* (V; 8 pc; PHANGS-HST)
* (region_ID_assoc_v_08pc)
integer*4 AgeAV16pc ! (Myr) ?=- Association age (V; 16 pc;
* PHANGS-HST) (age_assoc_v_16pc)
integer*4 e_AgeAV16pc ! (Myr) ?=- Uncertainty on association age
* (V; 16 pc; PHANGS-HST)
* (age_err_assoc_v_16pc)
real*8 DEdegV16pc ! (deg) ?=- Association declination
* (V; 16 pc; PHANGS-HST)
* (dec_assoc_v_16pc)
real*4 E_B_V_V16pc ! (mag) ?=- E(B-V) of association
* (V; 16 pc; PHANGS-HST)
* (ebv_assoc_v_16pc)
real*4 e_E_B_V_V16pc ! (mag) ?=- Uncertainty on E(B-V) of
* association (V; 16 pc; PHANGS-HST)
* (ebv_err_assoc_v_16pc)
real*4 MassAV16pc ! (Msun) ?=- Stellar mass (V; 16 pc;
* PHANGS-HST) (mass_assoc_v_16pc)
real*4 e_MassAV16pc ! (Msun) ?=- Uncertainty on stellar mass
* (V; 16 pc; PHANGS-HST)
* (mass_err_assoc_v_16pc)
integer*4 FlagMultV16pc ! [0/1] Flag: multiple association
* (V; 16 pc; PHANGS-HST)
* (flag_multi_assoc_v_16pc)
integer*4 FlagNoV16pc ! [0/1] Flag: no association
* (V; 16 pc; PHANGS-HST)
* (flag_none_assoc_v_16pc)
integer*4 FlagOneV16pc ! [0/1] Flag: single association
* (V; 16 pc; PHANGS-HST)
* (flag_one_assoc_v_16pc)
real*8 RAdegV16pc ! (deg) ?=- Association right ascension
* (V; 16 pc; PHANGS-HST)
* (ra_assoc_v_16pc)
integer*4 RegionIDV16pc ! ?=- Association identifier
* (V; 16 pc; PHANGS-HST)
* (region_ID_assoc_v_16pc)
integer*4 AgeAV32pc ! (Myr) ?=- Association age (V; 32 pc;
* PHANGS-HST) (age_assoc_v_32pc)
integer*4 e_AgeAV32pc ! (Myr) ?=- Uncertainty on association age
* (V; 32 pc; PHANGS-HST)
* (age_err_assoc_v_32pc)
real*8 DEdegV32pc ! (deg) ?=- Association declination (V;32 pc;
* PHANGS-HST) (dec_assoc_v_32pc)
real*4 E_B_V_V32pc ! (mag) ?=- E(B-V) of association (V; 32 pc;
* PHANGS-HST) (ebv_assoc_v_32pc)
real*4 e_E_B_V_V32pc ! (mag) ?=- Uncertainty on E(B-V) of
* association (V; 32 pc; PHANGS-HST)
* (ebv_err_assoc_v_32pc)
real*4 MassAV32pc ! (Msun) ?=- Stellar mass (V; 32 pc;
* PHANGS-HST) (mass_assoc_v_32pc)
real*4 e_MassAV32pc ! (Msun) ?=- Uncertainty on stellar mass
* (V; 32 pc; PHANGS-HST)
* (mass_err_assoc_v_32pc)
integer*4 FlagMultV32pc ! [0/1] Flag: multiple associations
* (V; 32 pc; PHANGS-HST)
* (flag_multi_assoc_v_32pc)
integer*4 FlagNoV32pc ! [0/1] Flag: no associations
* (V; 32 pc; PHANGS-HST)
* (flag_none_assoc_v_32pc)
integer*4 FlagOneV32pc ! [0/1] Flag: single associations
* (V; 32 pc; PHANGS-HST)
* (flag_one_assoc_v_32pc)
real*8 RAdegV32pc ! (deg) ?=- Association right ascension
* (V; 32 pc; PHANGS-HST)
* (ra_assoc_v_32pc)
integer*4 RegionIDV32pc ! ?=- Association identifier
* (V; 32 pc; PHANGS-HST)
* (region_ID_assoc_v_32pc)
integer*4 AgeAV64pc ! (Myr) ?=- Association age
* (V; 64 pc; PHANGS-HST)
* (age_assoc_v_64pc)
integer*4 e_AgeAV64pc ! (Myr) ?=- Uncertainty on association age
* (V; 64 pc; PHANGS-HST)
* (age_err_assoc_v_64pc)
real*8 DEdegV64pc ! (deg) ?=- Association declination
* (V; 64 pc; PHANGS-HST)
* (dec_assoc_v_64pc)
real*4 E_B_V_V64pc ! (mag) ?=- E(B-V) of association (V; 64 pc;
* PHANGS-HST) (ebv_assoc_v_64pc)
real*4 e_E_B_V_V64pc ! (mag) ?=- Uncertainty on E(B-V) of
* association (V; 64 pc; PHANGS-HST)
* (ebv_err_assoc_v_64pc)
real*4 MassAV64pc ! (Msun) ?=- Stellar mass (V; 64 pc;
* PHANGS-HST; PHANGS-HST)
* (mass_assoc_v_64pc)
real*4 e_MassAV64pc ! (Msun) ?=- Uncertainty on stellar mass
* (V; 64 pc; PHANGS-HST; PHANGS-HST)
* (mass_err_assoc_v_64pc)
integer*4 FlagMultV64pc ! [0/1] Flag: multiple associations
* (V; 64 pc; PHANGS-HST)
* (flag_multi_assoc_v_64pc)
integer*4 FlagNoV64pc ! [0/1] Flag: no associations
* (V; 64 pc; PHANGS-HST)
* (flag_none_assoc_v_64pc)
integer*4 FlagOneV64pc ! [0/1] Flag: single associations
* (V; 64 pc; PHANGS-HST)
* (flag_one_assoc_v_64pc)
real*8 RAdegV64pc ! (deg) ?=- Association right ascension
* (V; 64 pc; PHANGS-HST)
* (ra_assoc_v_64pc)
integer*4 RegionIDV64pc ! ?=- 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 ! (deg) Right Ascension of center (J2000)
real*8 DEdeg_1 ! (deg) Declination of center (J2000)
integer*4 Nx ! Number of pixels along X-axis
integer*4 Ny ! Number of pixels along Y-axis
integer*4 size ! (Kibyte) Size of FITS file
character*37 FileName ! Name of FITS file with subdirectory
character*36 Title ! 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,Region,fEdgeHST,fTouchHST,fManualHST,fEdgeMUSE,
+ fStarMUSE,RAdeg,DEdeg,RApdeg,DEpdeg,RAMdeg,DEMdeg,rR25Deproj,
+ rrefFdeproj,phiDeproj,Env,Areacirc,Radcirc,Radcircpc,Area2,
+ amaj2,bmin2,PA2,Rad,Radpc,Raddeconv,RadStr,RadcircpcM,
+ ComplexScore,FHaHST,e_FHaHST,FHaHSTMax,FHaHSTMin,FHameanHST,
+ FHacorrHST,e_FHacorrHST,FHa,FHb,FOIII5007,FNII6584,FSII6716,
+ FSII6731,FSIII9069,FOI6300,FNII5755,FHeI5876,FSIII6312,
+ FOI6363,FOII7319,FOII7330,e_FHa,e_FHb,e_FOIII5007,e_FNII6584,
+ e_FSII6716,e_FSII6731,e_FSIII9069,e_FOI6300,e_FNII5755,
+ e_FHeI5876,e_FSIII6312,e_FOI6363,e_FOII7319,e_FOII7330,
+ FHacorr,FHbcorr,FOIII5007corr,FNII6584corr,FSII6716corr,
+ FSII6731corr,FSIII9069corr,FOI6300corr,FNII5755corr,
+ FHeI5876corr,FSIII6312corr,FOI6363corr,FOII7319corr,
+ FOII7330corr,e_FHacorr,e_FHbcorr,e_FOIII5007corr,
+ e_FNII6584corr,e_FSII6716corr,e_FSII6731corr,e_FSIII9069corr,
+ e_FOI6300corr,e_FNII5755corr,e_FHeI5876corr,e_FSIII6312corr,
+ e_FOI6363corr,e_FOII7319corr,e_FOII7330corr,VHa,e_VHa,sigmaHa,
+ e_sigmaHa,VNII6584,e_VNII6584,sigmaNII6584,e_sigmaNII6584,
+ VOIII5007,e_VOIII5007,sigmaOIII5007,e_sigmaOIII5007,BPT_NII,
+ BPT_SII,BPT_OI,EWHa,e_EWHa,EWHb,e_EWHb,EWHaraw,e_EWHaraw,
+ EWHbraw,e_EWHbraw,EWHafit,e_EWHafit,EWHbfit,e_EWHbfit,
+ EWHabgcorr,e_EWHabgcorr,EWHbbgcorr,e_EWHbbgcorr,LHaHST,
+ e_LHaHST,LHa,HIIclassv2,HIIclassv3,TeNII,e_TeNII,TeSIII,
+ e_TeSIII,neSII,e_neSII,q,e_q,logqd91,e_logqd91,E_B_V,e_E_B_V,
+ AV,MetScal,e_MetScal,deltaMetScal,AgeANUV08pc,e_AgeANUV08pc,
+ DEdegNUV08pc,E_B_V_NUV08pc,e_E_B_V_NUV08pc,MassANUV08pc,
+ e_MassANUV08pc,FlagMultNUV08p,FlagNoNUV08pc,FlagOneNUV08pc,
+ RAdegNUV08pc,RegionIDNUV08p,AgeANUV16pc,e_AgeANUV16pc,
+ DEdegNUV16pc,E_B_V_NUV16pc,e_E_B_V_NUV16pc,MassANUV16pc,
+ e_MassANUV16pc,FlagMultNUV16p,FlagNoNUV16pc,FlagOneNUV16pc,
+ RAdegNUV16pc,RegionIDNUV16p,AgeANUV32pc,e_AgeANUV32pc,
+ DEdegNUV32pc,E_B_V_NUV32pc,e_E_B_V_NUV32pc,MassANUV32pc,
+ e_MassANUV32pc,FlagMultNUV32p,FlagNoNUV32pc,FlagOneNUV32pc,
+ RAdegNUV32pc,RegionIDNUV32p,AgeANUV64pc,e_AgeANUV64pc,
+ DEdegNUV64pc,E_B_V_NUV64pc,e_E_B_V_NUV64pc,MassANUV64pc,
+ e_MassANUV64pc,FlagMultNUV64p,FlagNoNUV64pc,FlagOneNUV64pc,
+ RAdegNUV64pc,RegionIDNUV64p,AgeAV08pc,e_AgeAV08pc,DEdegV08pc,
+ E_B_V_V08pc,e_E_B_V_V08pc,MassAV08pc,e_MassAV08pc,
+ FlagMultV08pc,FlagNoV08pc,FlagOneV08pc,RAdegV08pc,
+ RegionIDV08pc,AgeAV16pc,e_AgeAV16pc,DEdegV16pc,E_B_V_V16pc,
+ e_E_B_V_V16pc,MassAV16pc,e_MassAV16pc,FlagMultV16pc,
+ FlagNoV16pc,FlagOneV16pc,RAdegV16pc,RegionIDV16pc,AgeAV32pc,
+ e_AgeAV32pc,DEdegV32pc,E_B_V_V32pc,e_E_B_V_V32pc,MassAV32pc,
+ e_MassAV32pc,FlagMultV32pc,FlagNoV32pc,FlagOneV32pc,
+ RAdegV32pc,RegionIDV32pc,AgeAV64pc,e_AgeAV64pc,DEdegV64pc,
+ E_B_V_V64pc,e_E_B_V_V64pc,MassAV64pc,e_MassAV64pc,
+ FlagMultV64pc,FlagNoV64pc,FlagOneV64pc,RAdegV64pc,
+ RegionIDV64pc
if (idig(ar__(257:266)).EQ.0) Raddeconv = rNULL__
if (idig(ar__(268:277)).EQ.0) RadStr = rNULL__
if (idig(ar__(489:500)).EQ.0) FNII5755 = rNULL__
if (idig(ar__(515:526)).EQ.0) FSIII6312 = rNULL__
if (idig(ar__(671:682)).EQ.0) e_FNII5755 = rNULL__
if (idig(ar__(684:695)).EQ.0) e_FHeI5876 = rNULL__
if (idig(ar__(697:708)).EQ.0) e_FSIII6312 = rNULL__
if (idig(ar__(853:864)).EQ.0) FNII5755corr = rNULL__
if (idig(ar__(879:890)).EQ.0) FSIII6312corr = rNULL__
if (idig(ar__(905:916)).EQ.0) FOII7319corr = rNULL__
if (idig(ar__(918:929)).EQ.0) FOII7330corr = rNULL__
if (idig(ar__(957:968)).EQ.0) e_FOIII5007corr = rNULL__
if (idig(ar__(970:981)).EQ.0) e_FNII6584corr = rNULL__
if (idig(ar__(983:994)).EQ.0) e_FSII6716corr = rNULL__
if (idig(ar__(996:1007)).EQ.0) e_FSII6731corr = rNULL__
if (idig(ar__(1009:1020)).EQ.0) e_FSIII9069corr = rNULL__
if (idig(ar__(1022:1033)).EQ.0) e_FOI6300corr = rNULL__
if (idig(ar__(1035:1046)).EQ.0) e_FNII5755corr = rNULL__
if (idig(ar__(1048:1059)).EQ.0) e_FHeI5876corr = rNULL__
if (idig(ar__(1061:1072)).EQ.0) e_FSIII6312corr = rNULL__
if (idig(ar__(1074:1085)).EQ.0) e_FOI6363corr = rNULL__
if (idig(ar__(1087:1098)).EQ.0) e_FOII7319corr = rNULL__
if (idig(ar__(1100:1111)).EQ.0) e_FOII7330corr = rNULL__
if (idig(ar__(1434:1445)).EQ.0) EWHabgcorr = rNULL__
if (idig(ar__(1447:1458)).EQ.0) e_EWHabgcorr = rNULL__
if (idig(ar__(1460:1471)).EQ.0) EWHbbgcorr = rNULL__
if (idig(ar__(1473:1484)).EQ.0) e_EWHbbgcorr = rNULL__
if(ar__(1525:1526) .EQ. '') HIIclassv2 = iNULL__
if(ar__(1528:1529) .EQ. '') HIIclassv3 = iNULL__
if (idig(ar__(1544:1555)).EQ.0) e_TeNII = rNULL__
if (idig(ar__(1557:1568)).EQ.0) TeSIII = rNULL__
if (idig(ar__(1570:1581)).EQ.0) e_TeSIII = rNULL__
if (idig(ar__(1583:1594)).EQ.0) neSII = rNULL__
if (idig(ar__(1596:1607)).EQ.0) e_neSII = rNULL__
if (idig(ar__(1635:1646)).EQ.0) logqd91 = rNULL__
if (idig(ar__(1648:1659)).EQ.0) e_logqd91 = rNULL__
if (idig(ar__(1661:1672)).EQ.0) E_B_V = rNULL__
if (idig(ar__(1674:1685)).EQ.0) e_E_B_V = rNULL__
if (idig(ar__(1687:1698)).EQ.0) AV = rNULL__
if (idig(ar__(1700:1711)).EQ.0) MetScal = rNULL__
if (idig(ar__(1713:1724)).EQ.0) e_MetScal = rNULL__
if (idig(ar__(1726:1737)).EQ.0) deltaMetScal = rNULL__
if (AgeANUV08pc .EQ. 45) AgeANUV08pc = iNULL__
if (idig(ar__(1743:1748)).EQ.0) e_AgeANUV08pc = rNULL__
if (idig(ar__(1750:1757)).EQ.0) DEdegNUV08pc = rNULL__
if (idig(ar__(1759:1764)).EQ.0) E_B_V_NUV08pc = rNULL__
if (idig(ar__(1766:1771)).EQ.0) e_E_B_V_NUV08pc = rNULL__
if (idig(ar__(1773:1784)).EQ.0) MassANUV08pc = rNULL__
if (idig(ar__(1786:1797)).EQ.0) e_MassANUV08pc = rNULL__
if (idig(ar__(1805:1819)).EQ.0) RAdegNUV08pc = rNULL__
if (RegionIDNUV08p .EQ. 45) RegionIDNUV08p = iNULL__
if (AgeANUV16pc .EQ. 45) AgeANUV16pc = iNULL__
if (e_AgeANUV16pc .EQ. 45) e_AgeANUV16pc = iNULL__
if (idig(ar__(1833:1847)).EQ.0) DEdegNUV16pc = rNULL__
if (idig(ar__(1849:1852)).EQ.0) E_B_V_NUV16pc = rNULL__
if (idig(ar__(1854:1857)).EQ.0) e_E_B_V_NUV16pc = rNULL__
if (idig(ar__(1859:1870)).EQ.0) MassANUV16pc = rNULL__
if (idig(ar__(1872:1883)).EQ.0) e_MassANUV16pc = rNULL__
if (idig(ar__(1891:1905)).EQ.0) RAdegNUV16pc = rNULL__
if (RegionIDNUV16p .EQ. 45) RegionIDNUV16p = iNULL__
if (AgeANUV32pc .EQ. 45) AgeANUV32pc = iNULL__
if (e_AgeANUV32pc .EQ. 45) e_AgeANUV32pc = iNULL__
if (idig(ar__(1922:1936)).EQ.0) DEdegNUV32pc = rNULL__
if (idig(ar__(1938:1941)).EQ.0) E_B_V_NUV32pc = rNULL__
if (idig(ar__(1943:1946)).EQ.0) e_E_B_V_NUV32pc = rNULL__
if (idig(ar__(1948:1959)).EQ.0) MassANUV32pc = rNULL__
if (idig(ar__(1961:1972)).EQ.0) e_MassANUV32pc = rNULL__
if (idig(ar__(1980:1994)).EQ.0) RAdegNUV32pc = rNULL__
if (RegionIDNUV32p .EQ. 45) RegionIDNUV32p = iNULL__
if (AgeANUV64pc .EQ. 45) AgeANUV64pc = iNULL__
if (e_AgeANUV64pc .EQ. 45) e_AgeANUV64pc = iNULL__
if (idig(ar__(2009:2023)).EQ.0) DEdegNUV64pc = rNULL__
if (idig(ar__(2025:2028)).EQ.0) E_B_V_NUV64pc = rNULL__
if (idig(ar__(2030:2033)).EQ.0) e_E_B_V_NUV64pc = rNULL__
if (idig(ar__(2035:2046)).EQ.0) MassANUV64pc = rNULL__
if (idig(ar__(2048:2059)).EQ.0) e_MassANUV64pc = rNULL__
if (idig(ar__(2067:2081)).EQ.0) RAdegNUV64pc = rNULL__
if (RegionIDNUV64p .EQ. 45) RegionIDNUV64p = iNULL__
if (AgeAV08pc .EQ. 45) AgeAV08pc = iNULL__
if (e_AgeAV08pc .EQ. 45) e_AgeAV08pc = iNULL__
if (idig(ar__(2100:2114)).EQ.0) DEdegV08pc = rNULL__
if (idig(ar__(2116:2119)).EQ.0) E_B_V_V08pc = rNULL__
if (idig(ar__(2121:2124)).EQ.0) e_E_B_V_V08pc = rNULL__
if (idig(ar__(2126:2137)).EQ.0) MassAV08pc = rNULL__
if (idig(ar__(2139:2150)).EQ.0) e_MassAV08pc = rNULL__
if (idig(ar__(2158:2172)).EQ.0) RAdegV08pc = rNULL__
if (RegionIDV08pc .EQ. 45) RegionIDV08pc = iNULL__
if (AgeAV16pc .EQ. 45) AgeAV16pc = iNULL__
if (e_AgeAV16pc .EQ. 45) e_AgeAV16pc = iNULL__
if (idig(ar__(2191:2205)).EQ.0) DEdegV16pc = rNULL__
if (idig(ar__(2207:2210)).EQ.0) E_B_V_V16pc = rNULL__
if (idig(ar__(2212:2215)).EQ.0) e_E_B_V_V16pc = rNULL__
if (idig(ar__(2217:2228)).EQ.0) MassAV16pc = rNULL__
if (idig(ar__(2230:2241)).EQ.0) e_MassAV16pc = rNULL__
if (idig(ar__(2249:2263)).EQ.0) RAdegV16pc = rNULL__
if (RegionIDV16pc .EQ. 45) RegionIDV16pc = iNULL__
if (AgeAV32pc .EQ. 45) AgeAV32pc = iNULL__
if (e_AgeAV32pc .EQ. 45) e_AgeAV32pc = iNULL__
if (idig(ar__(2282:2296)).EQ.0) DEdegV32pc = rNULL__
if (idig(ar__(2298:2301)).EQ.0) E_B_V_V32pc = rNULL__
if (idig(ar__(2303:2306)).EQ.0) e_E_B_V_V32pc = rNULL__
if (idig(ar__(2308:2319)).EQ.0) MassAV32pc = rNULL__
if (idig(ar__(2321:2332)).EQ.0) e_MassAV32pc = rNULL__
if (idig(ar__(2340:2354)).EQ.0) RAdegV32pc = rNULL__
if (RegionIDV32pc .EQ. 45) RegionIDV32pc = iNULL__
if (AgeAV64pc .EQ. 45) AgeAV64pc = iNULL__
if (e_AgeAV64pc .EQ. 45) e_AgeAV64pc = iNULL__
if (idig(ar__(2372:2386)).EQ.0) DEdegV64pc = rNULL__
if (idig(ar__(2388:2391)).EQ.0) E_B_V_V64pc = rNULL__
if (idig(ar__(2393:2396)).EQ.0) e_E_B_V_V64pc = rNULL__
if (idig(ar__(2398:2409)).EQ.0) MassAV64pc = rNULL__
if (idig(ar__(2411:2422)).EQ.0) e_MassAV64pc = rNULL__
if (idig(ar__(2430:2444)).EQ.0) RAdegV64pc = rNULL__
if (RegionIDV64pc .EQ. 45) RegionIDV64pc = iNULL__
c ..............Just test output...........
write(6,1)
+ Galaxy,Region,fEdgeHST,fTouchHST,fManualHST,fEdgeMUSE,
+ fStarMUSE,RAdeg,DEdeg,RApdeg,DEpdeg,RAMdeg,DEMdeg,rR25Deproj,
+ rrefFdeproj,phiDeproj,Env,Areacirc,Radcirc,Radcircpc,Area2,
+ amaj2,bmin2,PA2,Rad,Radpc,Raddeconv,RadStr,RadcircpcM,
+ ComplexScore,FHaHST,e_FHaHST,FHaHSTMax,FHaHSTMin,FHameanHST,
+ FHacorrHST,e_FHacorrHST,FHa,FHb,FOIII5007,FNII6584,FSII6716,
+ FSII6731,FSIII9069,FOI6300,FNII5755,FHeI5876,FSIII6312,
+ FOI6363,FOII7319,FOII7330,e_FHa,e_FHb,e_FOIII5007,e_FNII6584,
+ e_FSII6716,e_FSII6731,e_FSIII9069,e_FOI6300,e_FNII5755,
+ e_FHeI5876,e_FSIII6312,e_FOI6363,e_FOII7319,e_FOII7330,
+ FHacorr,FHbcorr,FOIII5007corr,FNII6584corr,FSII6716corr,
+ FSII6731corr,FSIII9069corr,FOI6300corr,FNII5755corr,
+ FHeI5876corr,FSIII6312corr,FOI6363corr,FOII7319corr,
+ FOII7330corr,e_FHacorr,e_FHbcorr,e_FOIII5007corr,
+ e_FNII6584corr,e_FSII6716corr,e_FSII6731corr,e_FSIII9069corr,
+ e_FOI6300corr,e_FNII5755corr,e_FHeI5876corr,e_FSIII6312corr,
+ e_FOI6363corr,e_FOII7319corr,e_FOII7330corr,VHa,e_VHa,sigmaHa,
+ e_sigmaHa,VNII6584,e_VNII6584,sigmaNII6584,e_sigmaNII6584,
+ VOIII5007,e_VOIII5007,sigmaOIII5007,e_sigmaOIII5007,BPT_NII,
+ BPT_SII,BPT_OI,EWHa,e_EWHa,EWHb,e_EWHb,EWHaraw,e_EWHaraw,
+ EWHbraw,e_EWHbraw,EWHafit,e_EWHafit,EWHbfit,e_EWHbfit,
+ EWHabgcorr,e_EWHabgcorr,EWHbbgcorr,e_EWHbbgcorr,LHaHST,
+ e_LHaHST,LHa,HIIclassv2,HIIclassv3,TeNII,e_TeNII,TeSIII,
+ e_TeSIII,neSII,e_neSII,q,e_q,logqd91,e_logqd91,E_B_V,e_E_B_V,
+ AV,MetScal,e_MetScal,deltaMetScal,AgeANUV08pc,e_AgeANUV08pc,
+ DEdegNUV08pc,E_B_V_NUV08pc,e_E_B_V_NUV08pc,MassANUV08pc,
+ e_MassANUV08pc,FlagMultNUV08p,FlagNoNUV08pc,FlagOneNUV08pc,
+ RAdegNUV08pc,RegionIDNUV08p,AgeANUV16pc,e_AgeANUV16pc,
+ DEdegNUV16pc,E_B_V_NUV16pc,e_E_B_V_NUV16pc,MassANUV16pc,
+ e_MassANUV16pc,FlagMultNUV16p,FlagNoNUV16pc,FlagOneNUV16pc,
+ RAdegNUV16pc,RegionIDNUV16p,AgeANUV32pc,e_AgeANUV32pc,
+ DEdegNUV32pc,E_B_V_NUV32pc,e_E_B_V_NUV32pc,MassANUV32pc,
+ e_MassANUV32pc,FlagMultNUV32p,FlagNoNUV32pc,FlagOneNUV32pc,
+ RAdegNUV32pc,RegionIDNUV32p,AgeANUV64pc,e_AgeANUV64pc,
+ DEdegNUV64pc,E_B_V_NUV64pc,e_E_B_V_NUV64pc,MassANUV64pc,
+ e_MassANUV64pc,FlagMultNUV64p,FlagNoNUV64pc,FlagOneNUV64pc,
+ RAdegNUV64pc,RegionIDNUV64p,AgeAV08pc,e_AgeAV08pc,DEdegV08pc,
+ E_B_V_V08pc,e_E_B_V_V08pc,MassAV08pc,e_MassAV08pc,
+ FlagMultV08pc,FlagNoV08pc,FlagOneV08pc,RAdegV08pc,
+ RegionIDV08pc,AgeAV16pc,e_AgeAV16pc,DEdegV16pc,E_B_V_V16pc,
+ e_E_B_V_V16pc,MassAV16pc,e_MassAV16pc,FlagMultV16pc,
+ FlagNoV16pc,FlagOneV16pc,RAdegV16pc,RegionIDV16pc,AgeAV32pc,
+ e_AgeAV32pc,DEdegV32pc,E_B_V_V32pc,e_E_B_V_V32pc,MassAV32pc,
+ e_MassAV32pc,FlagMultV32pc,FlagNoV32pc,FlagOneV32pc,
+ RAdegV32pc,RegionIDV32pc,AgeAV64pc,e_AgeAV64pc,DEdegV64pc,
+ E_B_V_V64pc,e_E_B_V_V64pc,MassAV64pc,e_MassAV64pc,
+ FlagMultV64pc,FlagNoV64pc,FlagOneV64pc,RAdegV64pc,
+ RegionIDV64pc
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,DEdeg_1,Nx,Ny,size,FileName,Title
c ..............Just test output...........
write(6,2)RAdeg_1,DEdeg_1,Nx,Ny,size,FileName,Title
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