Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/893/96 into FORTRAN code for loading all data files into arrays.
Note that special values are assigned to unknown or unspecified
numbers (also called NULL numbers);
when necessary, the coordinate components making up the right ascension
and declination are converted into floating-point numbers
representing these angles in degrees.
program load_ReadMe
C=============================================================================
C F77-compliant program generated by readme2f_1.81 (2015-09-23), on 2026-Oct-13
C=============================================================================
* This code was generated from the ReadMe file documenting a catalogue
* according to the "Standard for Documentation of Astronomical Catalogues"
* currently in use by the Astronomical Data Centers (CDS, ADC, A&A)
* (see full documentation at URL http://vizier.u-strasbg.fr/doc/catstd.htx)
* Please report problems or questions to
C=============================================================================
implicit none
* Unspecified or NULL values, generally corresponding to blank columns,
* are assigned one of the following special values:
* rNULL__ for unknown or NULL floating-point values
* iNULL__ for unknown or NULL integer values
real*4 rNULL__
integer*4 iNULL__
parameter (rNULL__=--2147483648.) ! NULL real number
parameter (iNULL__=(-2147483647-1)) ! NULL int number
integer idig ! testing NULL number
C=============================================================================
Cat. J/ApJ/893/96 CHAOS IV: NGC3184 LBT obs. & 3 other gal. abundances (Berg+, 2020)
*================================================================================
*CHAOS IV: gas-phase abundance trends from the first four CHAOS galaxies.
* Berg D.A., Pogge R.W., Skillman E.D., Croxall K.V., Moustakas J.,
* Rogers N.S.J., Sun J.
* <Astrophys. J., 893, 96 (2020)>
* =2020ApJ...893...96B
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea1.dat' ! NGC 3184 MODS/LBT observations taken
on UT 2012-03-24; 2013-01-06
integer*4 nr__
parameter (nr__=52) ! Number of records
character*88 ar__ ! Full-size record
C J2000 position composed of: RAh RAm RAs DEd DEm DEs
real*8 RAdeg (nr__) ! (deg) Right Ascension J2000
real*8 DEdeg (nr__) ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*19 Reg (nr__) ! NGC 3184 HII region slit ID (1)
integer*4 RAh (nr__) ! (h) [10] Hour of Right Ascension (J2000)
integer*4 RAm (nr__) ! (min) [18] Minute of Right Ascension (J2000)
real*4 RAs (nr__) ! (s) [1.5/27.8] Second of Right Ascension (J2000)
integer*4 DEd (nr__) ! (deg) [41] Degree of Declination (J2000)
integer*4 DEm (nr__) ! (arcmin) [22/27] Arcminute of Declination (J2000)
real*4 DEs (nr__) ! (arcsec) Arcsecond of Declination (J2000)
real*4 Rgarc (nr__) ! (arcsec) [2/182.7] De-projected HII region distance
* galaxy center
real*4 R_R25 (nr__) ! [0.009/0.9] Fraction of R25 radius
real*4 Rgkpc (nr__) ! (kpc) [0.1/10.4] HII region distance in kpc
real*4 oRA (nr__) ! (arcsec) [-172.5/123.8] RA offset from the central position
* (RAJ2000=10:18:16.86,DEJ2000=+41:25:26.59)
real*4 oDE (nr__) ! (arcsec) [-148.2/148.8] DE offset from the central position
* (RAJ2000=10:18:16.86,DEJ2000=+41:25:26.59)
integer*4 dOIII (nr__) ! [1]? Detection, [OIII] w4363 auroral line (2)
integer*4 dNII (nr__) ! [1]? Detection, [NII] w5755 auroral line (2)
integer*4 dSIII (nr__) ! [1]? Detection, [SIII] w6312 auroral line (2)
integer*4 dOII (nr__) ! [1]? Detection, [OII] w7320,w7330 auroral lines (2)
integer*4 dSII (nr__) ! [1]? Detection, [SII] w4068,w4076 auroral lines (2)
integer*4 dWR (nr__) ! [1]? HII regions have Wolf-Rayet feature detections
*Note (1): Slit ID, composed of the galaxy name and the offset in RA and DE,
* in arcseconds, from the central position listed for NGC 3184 in Table 1.
*Note (2): Detected auroral line at 3{sigma} significance.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea2.dat' ! Emission-Line intensities and equivalent widths
for LBT observations of HII regions in NGC 3184
integer*4 nr__1
parameter (nr__1=52) ! Number of records
character*715 ar__1 ! Full-size record
character*19 Reg_1 (nr__1) ! NGC 3184 HII region slit ID (1)
real*4 H14_w3721 (nr__1) ! [0.009/0.025] Line intensity, H14
* w3721 (2)
real*4 e_H14_w3721(nr__1) ! [0.001/0.02] Uncertainty in H14-w3721
real*4 OII_w3727 (nr__1) ! [0.3/4.2] Line intensity, [OII]
* w3727 (2)
real*4 e_OII_w3727(nr__1) ! [0.007/0.8] Uncertainty in OII-w3727
real*4 H13_w3734 (nr__1) ! [0.01/0.031] Line intensity, H13
* w3734 (2)
real*4 e_H13_w3734(nr__1) ! [0.001/0.02] Uncertainty in H13-w3734
real*4 H12_w3750 (nr__1) ! [0.003/0.06] Line intensity, H12
* w3750 (2)
real*4 e_H12_w3750(nr__1) ! [0.001/0.015] Uncertainty in H12-w3750
real*4 H11_w3770 (nr__1) ! [0.019/0.09]? Line intensity, H11
* w3770 (2)
real*4 e_H11_w3770(nr__1) ! [0.001/0.04]? Uncertainty in H11-w3770
real*4 H10_w3797 (nr__1) ! [0.02/0.07] Line intensity, H10
* w3797 (2)
real*4 e_H10_w3797(nr__1) ! [0.001/0.04] Uncertainty in H10-w3797
real*4 HeI_w3819 (nr__1) ! [0.001/0.05]? Line intensity, HeI
* w3819 (2)
real*4 e_HeI_w3819(nr__1) ! [0.001/0.01]? Uncertainty in HeI-w3819
real*4 H9_w3835 (nr__1) ! [0.02/0.11] Line intensity, H9 w3835
* (2)
real*4 e_H9_w3835 (nr__1) ! [0.001/0.07] Uncertainty in H9-w3835
real*4 NeIII_w3868(nr__1) ! [0.003/0.2]? Line intensity, [NeIII]
* w3868 (2)
real*4 e_NeIII_w3868(nr__1) ! [0.001/0.01]? Uncertainty in
* NeIII-w3868
real*4 HeI_w3888 (nr__1) ! [0.01/0.2] Line intensity, HeI w3888
* (2)
real*4 e_HeI_w3888(nr__1) ! [0.001/0.07] Uncertainty in HeI-w3888
real*4 H8_w3889 (nr__1) ! [0.04/0.2] Line intensity, H8 w3889
* (2)
real*4 e_H8_w3889 (nr__1) ! [0.001/0.08] Uncertainty in H8-w3889
real*4 HeI_w3964 (nr__1) ! [0.001/0.03]? Line intensity, HeI
* w3964 (2)
real*4 e_HeI_w3964(nr__1) ! [0.001/0.02]? Uncertainty in HeI-w3964
real*4 NeIII_w3967(nr__1) ! [0.003/0.2]? Line intensity, [NeIII]
* w3967 (2)
real*4 e_NeIII_w3967(nr__1) ! [0.001/0.04]? Uncertainty in
* NeIII-w3967
real*4 H7_w3970 (nr__1) ! [0.07/0.2] Line intensity, H7 w3970
* (2)
real*4 e_H7_w3970 (nr__1) ! [0.002/0.11] Uncertainty in H7-w3970
real*4 NeIII_w4011(nr__1) ! [0.001/0.03]? Line intensity,
* [NeIII] w4011 (2)
real*4 e_NeIII_w4011(nr__1) ! [0.001/0.006]? Uncertainty in
* NeIII-w4011
real*4 HeI_w4026 (nr__1) ! [0.001/0.04]? Line intensity, HeI
* w4026 (2)
real*4 e_HeI_w4026(nr__1) ! [0.001/0.004]? Uncertainty in
* HeI-w4026
real*4 SII_w4068 (nr__1) ! [0.001/0.041]? Line intensity, [SII]
* w4068 (2)
real*4 e_SII_w4068(nr__1) ! [0.001/0.01]? Uncertainty in SII-w4068
real*4 SII_w4076 (nr__1) ! [0.002/0.03]? Line intensity, [SII]
* w4076 (2)
real*4 e_SII_w4076(nr__1) ! [0.001/0.007]? Uncertainty in
* SII-w4076
real*4 Hd_w4101 (nr__1) ! [0.16/0.4] Line intensity, H{delta}
* w4101 (2)
real*4 e_Hd_w4101 (nr__1) ! [0.001/0.3] Uncertainty in Hd-w4101
real*4 HeI_w4120 (nr__1) ! [0.001/0.018]? Line intensity, HeI
* w4120 (2)
real*4 e_HeI_w4120(nr__1) ! [0.001/0.009]? Uncertainty in
* HeI-w4120
real*4 HeI_w4143 (nr__1) ! [0.001/0.013]? Line intensity, HeI
* w4143 (2)
real*4 e_HeI_w4143(nr__1) ! [0.001/0.004]? Uncertainty in
* HeI-w4143
real*4 Hg_w4340 (nr__1) ! [0.39/0.5] Line intensity, H{gamma}
* w4340 (2)
real*4 e_Hg_w4340 (nr__1) ! [0.002/0.5] Uncertainty in Hg-w4340
real*4 OIII_w4363 (nr__1) ! [0.001/0.014]? Line intensity,
* [OIII] w4363 (2)
real*4 e_OIII_w4363(nr__1) ! [0.001/0.011]? Uncertainty in
* OIII-w4363
real*4 HeI_w4387 (nr__1) ! [0.001/0.02]? Line intensity, HeI
* w4387 (2)
real*4 e_HeI_w4387(nr__1) ! [0.001/0.004]? Uncertainty in
* HeI-w4387
real*4 HeI_w4471 (nr__1) ! [0.001/0.041]? Line intensity, HeI
* w4471 (2)
real*4 e_HeI_w4471(nr__1) ! [0.001/0.006]? Uncertainty in
* HeI-w4471
real*4 FeIII_w4658(nr__1) ! [0.001/0.03]? Line intensity,
* [FeIII] w4658 (2)
real*4 e_FeIII_w4658(nr__1) ! [0.001/0.03]? Uncertainty in
* FeIII-w4658
real*4 HeII_w4686 (nr__1) ! [0.001/0.03]? Line intensity, HeII
* w4686 (2)
real*4 e_HeII_w4686(nr__1) ! [0.001/0.03]? Uncertainty in
* HeII-w4686
real*4 Hb_w4861 (nr__1) ! [1] Line intensity, H{beta} w4861 (2)
real*4 e_Hb_w4861 (nr__1) ! [0.003/1] Uncertainty in Hb-w4861
real*4 HeI_w4921 (nr__1) ! [0.001/0.02]? Line intensity, HeI
* w4921 (2)
real*4 e_HeI_w4921(nr__1) ! [0.001/0.03]? Uncertainty in HeI-w4921
real*4 OIII_w4959 (nr__1) ! [0.008/1.1] Line intensity, [OIII]
* w4959 (2)
real*4 e_OIII_w4959(nr__1) ! [0.001/0.05] Uncertainty in OIII-w4959
real*4 OIII_w5007 (nr__1) ! [0.03/3] Line intensity, [OIII]
* w5007 (2)
real*4 e_OIII_w5007(nr__1) ! [0.001/0.12] Uncertainty in OIII-w5007
real*4 HeI_w5015 (nr__1) ! [0.001/0.06] Line intensity, HeI
* w5015 (2)
real*4 e_HeI_w5015(nr__1) ! [0.001/0.03] Uncertainty in HeI-w5015
real*4 NII_w5755 (nr__1) ! [0.001/0.02]? Line intensity, [NII]
* w5755 (2)
real*4 e_NII_w5755(nr__1) ! [0.001/0.015]? Uncertainty in
* NII-w5755
real*4 HeI_w5876 (nr__1) ! [0.03/0.4] Line intensity, HeI w5876
* (2)
real*4 e_HeI_w5876(nr__1) ! [0.001/0.07] Uncertainty in HeI-w5876
real*4 OI_w6300 (nr__1) ! [0.003/0.06]? Line intensity, [OI]
* w6300 (2)
real*4 e_OI_w6300 (nr__1) ! [0.001/0.03]? Uncertainty in OI-w6300
real*4 SIII_w6312 (nr__1) ! [0.001/0.02]? Line intensity, [SIII]
* w6312 (2)
real*4 e_SIII_w6312(nr__1) ! [0.001/0.008]? Uncertainty in
* SIII-w6312
real*4 OI_w6363 (nr__1) ! [0.001/0.04]? Line intensity, [OI]
* w6363 (2)
real*4 e_OI_w6363 (nr__1) ! [0.001/0.01]? Uncertainty in OI-w6363
real*4 NII_w6548 (nr__1) ! [0.13/0.41] Line intensity, [NII]
* w6548 (2)
real*4 e_NII_w6548(nr__1) ! [0.004/0.31] Uncertainty in NII-w6548
real*4 Ha_w6563 (nr__1) ! [2.6/3.2] Line intensity, H{alpha}
* w6563 (2)
real*4 e_Ha_w6563 (nr__1) ! [0.02/3] Uncertainty in Ha-w6563
real*4 NII_w6584 (nr__1) ! [0.4/1.3] Line intensity, [NII]
* w6584 (2)
real*4 e_NII_w6584(nr__1) ! [0.005/1] Uncertainty in NII-w6584
real*4 HeI_w6678 (nr__1) ! [0.006/0.04] Line intensity, HeI
* w6678 (2)
real*4 e_HeI_w6678(nr__1) ! [0.001/0.03] Uncertainty in HeI-w6678
real*4 SII_w6717 (nr__1) ! [0.1/0.51] Line intensity, [SII]
* w6717 (2)
real*4 e_SII_w6717(nr__1) ! [0.002/0.4] Uncertainty in SII-w6717
real*4 SII_w6731 (nr__1) ! [0.09/0.4] Line intensity, [SII]
* w6731 (2)
real*4 e_SII_w6731(nr__1) ! [0.001/0.3] Uncertainty in SII-w6731
real*4 HeI_w7065 (nr__1) ! [0.002/0.033]? Line intensity, HeI
* w7065 (2)
real*4 e_HeI_w7065(nr__1) ! [0.001/0.009]? Uncertainty in
* HeI-w7065
real*4 ArIII_w7135(nr__1) ! [0.004/0.11]? Line intensity,
* [ArIII] w7135 (2)
real*4 e_ArIII_w7135(nr__1) ! [0.001/0.03]? Uncertainty in
* ArIII-w7135
real*4 OII_w7320 (nr__1) ! [0.001/0.03] Line intensity, [OII]
* w7320 (2)
real*4 e_OII_w7320(nr__1) ! [0.001/0.04] Uncertainty in OII-w7320
real*4 OII_w7330 (nr__1) ! [0.001/0.04]? Line intensity, [OII]
* w7330 (2)
real*4 e_OII_w7330(nr__1) ! [0.001/0.04]? Uncertainty in OII-w7330
real*4 ArIII_w7751(nr__1) ! [0.003/0.03]? Line intensity,
* [ArIII] w7751 (2)
real*4 e_ArIII_w7751(nr__1) ! [0.001/0.008]? Uncertainty in
* ArIII-w7751
real*4 P13_w8665 (nr__1) ! [0.001/0.02]? Line intensity, P13
* w8665 (2)
real*4 e_P13_w8665(nr__1) ! [0.001/0.03]? Uncertainty in P13-w8665
real*4 P12_w8750 (nr__1) ! [0.004/2.6]? Line intensity, P12
* w8750 (2)
real*4 e_P12_w8750(nr__1) ! [0.002/0.05]? Uncertainty in P12-w8750
real*4 P11_w8862 (nr__1) ! [0.001/0.08]? Line intensity, P11
* w8862 (2)
real*4 e_P11_w8862(nr__1) ! [0.002/0.04]? Uncertainty in P11-w8862
real*4 P10_w9015 (nr__1) ! [0.002/0.05]? Line intensity, P10
* w9015 (2)
real*4 e_P10_w9015(nr__1) ! [0.001/0.03]? Uncertainty in
* P10-w9015
real*4 SIII_w9069 (nr__1) ! [0.05/0.4]? Line intensity, [SIII]
* w9069 (2)
real*4 e_SIII_w9069(nr__1) ! [0.001/0.2]? Uncertainty in SIII-w9069
real*4 P9_w9229 (nr__1) ! [0.002/0.04]? Line intensity, P9
* w9229 (2)
real*4 e_P9_w9229 (nr__1) ! [0.001/0.03]? Uncertainty in P9-w9229
real*4 SIII_w9532 (nr__1) ! [0.18/1] Line intensity, [SIII]
* w9532 (2)
real*4 e_SIII_w9532(nr__1) ! [0.003/0.6] Uncertainty in SIII-w9532
real*4 P8_w9546 (nr__1) ! [0.01/0.2]? Line intensity, P8 w9546
* (2)
real*4 e_P8_w9546 (nr__1) ! [0.002/0.03]? Uncertainty in P8-w9546
real*4 E_B_V (nr__1) ! (mag) [0.09/0.7] E(B-V) reddening (3)
real*4 e_E_B_V (nr__1) ! (mag) [0.004/0.03] Uncertainty in E (B-V)
real*4 FHb (nr__1) ! (10-16mW/m2) [9.6/314] FH{beta} line flux,
* 10^-16^erg/s/cm^2^
real*4 e_FHb (nr__1) ! (10-16mW/m2) [0.1/1.4] Uncertainty in FH{beta}
*Note (1): Slit ID, composed of the galaxy name and the offset in RA and DE,
* in arcseconds, from the central position listed for NGC 3184 in Table 1.
*Note (2): Reddening-corrected emission line intensities relative to H{beta}.
*Note (3): Note that previous CHAOS papers incorrectly reported the E(B-V)
* quantities as CH{beta}, whereas the values listed here are E(B-V) (mag).
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tablea3.dat' ! Ionic and total abundances for LBT observations
of NGC 3184
integer*4 nr__2
parameter (nr__2=30) ! Number of records
character*351 ar__2 ! Full-size record
character*19 Reg_2 (nr__2) ! NGC 3184 HII region slit ID (1)
integer*4 Te_SII_ (nr__2) ! (K) [5700/23900]? Measured electron temperature
* from the [SII] ion
integer*4 e_Te_SII_ (nr__2) ! (K) [400]? Uncertainty in Te[SII]
integer*4 Te_NII_ (nr__2) ! (K) [6000/9900]? Measured electron temperature
* from the [NII] ion
integer*4 e_Te_NII_ (nr__2) ! (K) [1700]? Uncertainty in Te[NII]
integer*4 Te_OII_ (nr__2) ! (K) [6800/12500]? Measured electron temperature
* from the [OII] ion
integer*4 e_Te_OII_ (nr__2) ! (K) [400]? Uncertainty in Te[OII]
integer*4 Te_SIII_ (nr__2) ! (K) [6400/9400]? Measured electron temperature
* from the [SIII] ion
integer*4 e_Te_SIII_ (nr__2) ! (K) [600/600]? Uncertainty in Te[SIII]
integer*4 Te_OIII_ (nr__2) ! (K) [8600/9600]? Measured electron temperature
* from the [OIII] ion
integer*4 e_Te_OIII_ (nr__2) ! (K) [100]? Uncertainty in Te[OIII]
integer*4 nem (nr__2) ! (cm-3) [10/29500] Measured electron density
integer*4 e_nem (nr__2) ! (cm-3) [0/10] Uncertainty in nem
integer*4 Te_low (nr__2) ! (K) [6000/9900] Adopted ionization-based
* temperature, lower ionization zone (2)
integer*4 e_Te_low (nr__2) ! (K) [100/1900] Uncertainty in Te-low
integer*4 Te_int (nr__2) ! (K) [4700/9900] Adopted ionization-based
* temperature, intermediate ionization zone (2)
integer*4 e_Te_int (nr__2) ! (K) [100/2600] Uncertainty in Te-int
integer*4 Te_high (nr__2) ! (K) [4800/10300] Adopted ionization-based
* temperature, upper ionization zone (2)
integer*4 e_Te_high (nr__2) ! (K) [100/2700] Uncertainty in Te-high
integer*4 neu (nr__2) ! (cm-3) [10/100] Used electron density
integer*4 e_neu (nr__2) ! (cm-3) [10] Uncertainty in neu
real*4 O__H_ (nr__2) ! (10-5) [8.1/62.8] Abundance ratio, O+/H+
real*4 e_O__H_ (nr__2) ! (10-5) [1.4/32] Uncertainty in O+/H+
real*4 O_2_H_ (nr__2) ! (10-5) [1.8/57.7] Abundance ratio, O+2/H+
real*4 e_O_2_H_ (nr__2) ! (10-5) [0.2/67.3] Uncertainty in O+2/H+
real*4 log_O_H (nr__2) ! ([-]) [8.2/9] Oxygen abundance, 12+log (O/H)
real*4 e_log_O_H (nr__2) ! ([-]) [0.02/0.3] Uncertainty in log (O/H)
real*4 N__H_ (nr__2) ! (10-6) [10.2/110.4] Abundance ratio, N+/H+
real*4 e_N__H_ (nr__2) ! (10-6) [1.6/52.5] Uncertainty in N+/H+
real*4 ICF_N (nr__2) ! [1/2.7] Nitrogen ionization correction factor
real*4 e_ICF_N (nr__2) ! [0.09/2.91] Uncertainty in ICF (N)
real*4 log_N_H (nr__2) ! ([-]) [7.3/8.5] Nitrogen-to-hydrogen abundance
* ratio, log
real*4 e_log_N_H (nr__2) ! ([-]) [0.038/0.4] Uncertainty in log (N/H)
real*4 log_N_O (nr__2) ! ([-]) [-1.2/-0.4] Nitrogen-to-oxygen abundance
* ratio, log
real*4 e_log_N_O (nr__2) ! ([-]) [0.03/0.3] Uncertainty in log (N/O)
real*4 S__H_ (nr__2) ! (10-7) [6.2/86.6] Abundance ratio, S+/H+
real*4 e_S__H_ (nr__2) ! (10-7) [1.2/40] Uncertainty in S+/H+
real*4 S2_H_ (nr__2) ! (10-7) [41/227.1] Abundance ratio, S2/H+
real*4 e_S2_H_ (nr__2) ! (10-7) [3.3/173.5] Uncertainty in S2/H+
real*4 ICF_S (nr__2) ! [1/1.8] Sulfur ionization correction factor
real*4 e_ICF_S (nr__2) ! [0.09/1.5] Uncertainty in ICF (S)
real*4 log_S_H (nr__2) ! ([-]) [6.8/7.7] Sulfur-to-hydrogen abundance
* ratio, log
real*4 e_log_S_H (nr__2) ! ([-]) [0.037/0.4] Uncertainty in log (S/H)
real*4 log_S_O (nr__2) ! ([-]) [-1.6/-1.1] Sulfur-to-oxygen abundance
* ratio, log
real*4 e_log_S_O (nr__2) ! ([-]) [0.04/0.4] Uncertainty in log (S/O)
real*4 Ne_2_H_ (nr__2) ! (10-6) [3.9/1185]? Abundance ratio, Ne+2/H+
real*4 e_Ne_2_H_ (nr__2) ! (10-6) [0.7/1519]? Uncertainty in Ne+2/H+
real*4 ICF_Ne (nr__2) ! [1.4/3.1] Neon ionization correction factor
real*4 e_ICF_Ne (nr__2) ! [0.4/9.2] Uncertainty in ICF (Ne)
real*4 log_Ne_H (nr__2) ! ([-]) [7.4/9.3]? Neon-to-hydrogen abundance
* ratio, log
real*4 e_log_Ne_H (nr__2) ! ([-]) [0.07/0.5]? Uncertainty in log (Ne/H)
real*4 log_Ne_O (nr__2) ! ([-]) [-0.92/0.4]? Neon-to-oxygen abundance
* ratio, log
real*4 e_log_Ne_O (nr__2) ! ([-]) [0.04/0.4]? Uncertainty in log (Ne/O)
real*4 Ar_2_H_ (nr__2) ! (10-7) [4.0/20.8] Abundance ratio, Ar+2/H+
real*4 e_Ar_2_H_ (nr__2) ! (10-7) [0.5/18.6] Uncertainty in Ar+2/H+
real*4 ICF_Ar (nr__2) ! [1.4/3.1] Argon ionization correction factor
real*4 e_ICF_Ar (nr__2) ! [0.14/0.31] Uncertainty in ICF (Ar)
real*4 log_Ar_H (nr__2) ! ([-]) [6/6.7] Argon-to-hydrogen abundance ratio,
* log
real*4 e_log_Ar_H (nr__2) ! ([-]) [0.04/0.3] Uncertainty in log (Ar/H)
real*4 log_Ar_O (nr__2) ! ([-]) [-2.61/-1.8] Argon-to-oxygen abundance
* ratio, log
real*4 e_log_Ar_O (nr__2) ! ([-]) [0.06/0.4] Uncertainty in log (Ar/O)
*Note (1): Slit ID, composed of the galaxy name and the offset in RA and DE,
* in arcseconds, from the central position listed for NGC 3184 in Table 1.
*Note (2): Adopted ionization-based temperature based on criteria given
* in Figure 5.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tableb1.dat' ! Updated abundances for NGC 628
using ionization-based Te
integer*4 nr__3
parameter (nr__3=45) ! Number of records
character*72 ar__3 ! Full-size record
character*18 v_BSC2015_ (nr__3) ! NGC 628 HII region
real*4 log_O_H_1 (nr__3) ! ([-]) [8.1/8.8] Oxygen abundance, 12+log(O/H)
real*4 e_log_O_H_1(nr__3) ! ([-]) [0.01/0.14] Uncertainty in log(O/H)
real*4 log_N_O_1 (nr__3) ! ([-]) [-1.4/-0.6] Nitrogen-to-oxygen abundance
* ratio, log
real*4 e_log_N_O_1(nr__3) ! ([-]) [0.01/0.2] Uncertainty in log(N/O)
real*4 log_S_O_1 (nr__3) ! ([-]) [-1.5/-0.97] Sulfur-to-oxygen abundance
* ratio, log
real*4 e_log_S_O_1(nr__3) ! ([-]) [0.02/0.3] Uncertainty in log(S/O)
real*4 log_Ar_O_1 (nr__3) ! ([-]) [-2.5/-1.8] Argon-to-oxygen abundance
* ratio, log
real*4 e_log_Ar_O_1(nr__3) ! ([-]) [0.04/0.2] Uncertainty in log(Ar/O)
real*4 log_Ne_O_1 (nr__3) ! ([-]) [-1.1/-0.4]? Neon-to-oxygen abundance
* ratio, log
real*4 e_log_Ne_O_1(nr__3) ! ([-]) [0.01/0.2]? Uncertainty in log(Ne/O)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tableb2.dat' ! Updated abundances for NGC 5194
using ionization-based Te
integer*4 nr__4
parameter (nr__4=28) ! Number of records
character*73 ar__4 ! Full-size record
character*19 v_CPB2015_ (nr__4) ! NGC 5194 HII region
real*4 log_O_H_1 (nr__4) ! ([-]) [8.46/8.9] Oxygen abundance, 12+log(O/H)
real*4 e_log_O_H_1(nr__4) ! ([-]) [0.03/0.16] Uncertainty in log(O/H)
real*4 log_N_O_1 (nr__4) ! ([-]) [-0.8/-0.42] Nitrogen-to-oxygen abundance
* ratio, log
real*4 e_log_N_O_1(nr__4) ! ([-]) [0.04/0.2] Uncertainty in log(N/O)
real*4 log_S_O_1 (nr__4) ! ([-]) [-1.6/-1] Sulfur-to-oxygen abundance ratio,
* log
real*4 e_log_S_O_1(nr__4) ! ([-]) [0.05/0.3] Uncertainty in log(S/O)
real*4 log_Ar_O_1 (nr__4) ! ([-]) [-2.7/-1.8] Argon-to-oxygen abundance ratio,
* log
real*4 e_log_Ar_O_1(nr__4) ! ([-]) [0.05/0.24] Uncertainty in log(Ar/O)
real*4 log_Ne_O_1 (nr__4) ! ([-]) [-0.9/-0.6]? Neon-to-oxygen abundance ratio,
* log
real*4 e_log_Ne_O_1(nr__4) ! ([-]) [0.04/0.12]? Uncertainty in log(Ne/O)
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'tableb3.dat' ! Updated abundances for NGC 5457
using ionization-based Te
integer*4 nr__5
parameter (nr__5=73) ! Number of records
character*73 ar__5 ! Full-size record
character*19 v_CPB2016_ (nr__5) ! NGC 5457 HII region
* (<[CPB2016] NGC 5457+xxx.x+yyy.y> in Simbad)
real*4 log_O_H_1 (nr__5) ! ([-]) [7.8/8.9] Oxygen abundance, 12+log(O/H)
real*4 e_log_O_H_1(nr__5) ! ([-]) [0.01/0.16] Uncertainty in log(O/H)
real*4 log_N_O_1 (nr__5) ! ([-]) [-1.5/-0.6] Nitrogen-to-oxygen abundance
* ratio, log
real*4 e_log_N_O_1(nr__5) ! ([-]) [0.02/0.3] Uncertainty in log(N/O)
real*4 log_S_O_1 (nr__5) ! ([-]) [-1.7/-1] Sulfur-to-oxygen abundance ratio,
* log
real*4 e_log_S_O_1(nr__5) ! ([-]) [0.02/0.26] Uncertainty in log(S/O)
real*4 log_Ar_O_1 (nr__5) ! ([-]) [-2.8/-2]? Argon-to-oxygen abundance ratio,
* log
real*4 e_log_Ar_O_1(nr__5) ! ([-]) [0.04/0.2]? Uncertainty in log(Ar/O)
real*4 log_Ne_O_1 (nr__5) ! ([-]) [-1.1/-0.28]? Neon-to-oxygen abundance ratio,
* log
real*4 e_log_Ne_O_1(nr__5) ! ([-]) [0.01/0.21]? Uncertainty in log(Ne/O)
C=============================================================================
C Loading file 'tablea1.dat' ! NGC 3184 MODS/LBT observations taken
* on UT 2012-03-24; 2013-01-06
C Format for file interpretation
1 format(
+ A19,1X,I2,1X,I2,1X,F4.1,1X,I2,1X,I2,1X,F5.2,1X,F6.2,1X,F5.3,
+ 1X,F5.2,1X,F6.1,1X,F6.1,1X,I1,1X,I1,1X,I1,1X,I1,1X,I1,1X,I2)
C Effective file loading
open(unit=1,status='old',file=
+'tablea1.dat')
write(6,*) '....Loading file: tablea1.dat'
do i__=1,52
read(1,'(A88)')ar__
read(ar__,1)
+ Reg(i__),RAh(i__),RAm(i__),RAs(i__),DEd(i__),DEm(i__),
+ DEs(i__),Rgarc(i__),R_R25(i__),Rgkpc(i__),oRA(i__),oDE(i__),
+ dOIII(i__),dNII(i__),dSIII(i__),dOII(i__),dSII(i__),dWR(i__)
if(ar__(77:77) .EQ. '') dOIII(i__) = iNULL__
if(ar__(79:79) .EQ. '') dNII(i__) = iNULL__
if(ar__(81:81) .EQ. '') dSIII(i__) = iNULL__
if(ar__(83:83) .EQ. '') dOII(i__) = iNULL__
if(ar__(85:85) .EQ. '') dSII(i__) = iNULL__
if(ar__(87:88) .EQ. '') dWR(i__) = iNULL__
RAdeg(i__) = rNULL__
DEdeg(i__) = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh(i__) .GT. -180) RAdeg(i__)=RAh(i__)*15.
if(RAm(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAm(i__)/4.
if(RAs(i__) .GT. -180) RAdeg(i__)=RAdeg(i__)+RAs(i__)/240.
if(DEd(i__) .GE. 0) DEdeg(i__)=DEd(i__)
if(DEm(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEm(i__)/60.
if(DEs(i__) .GE. 0) DEdeg(i__)=DEdeg(i__)+DEs(i__)/3600.
c ..............Just test output...........
write(6,1)
+ Reg(i__),RAh(i__),RAm(i__),RAs(i__),DEd(i__),DEm(i__),
+ DEs(i__),Rgarc(i__),R_R25(i__),Rgkpc(i__),oRA(i__),oDE(i__),
+ dOIII(i__),dNII(i__),dSIII(i__),dOII(i__),dSII(i__),dWR(i__)
write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea2.dat' ! Emission-Line intensities and equivalent widths
* for LBT observations of HII regions in NGC 3184
C Format for file interpretation
2 format(
+ A19,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,
+ 1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,
+ F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,
+ 1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,
+ F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,
+ 1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,
+ F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,
+ 1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,
+ F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,
+ 1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,
+ F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,
+ 1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,
+ F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,
+ 1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,
+ F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,
+ 1X,F5.3,1X,F5.3,1X,F6.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'tablea2.dat')
write(6,*) '....Loading file: tablea2.dat'
do i__=1,52
read(1,'(A715)')ar__1
read(ar__1,2)
+ Reg_1(i__),H14_w3721(i__),e_H14_w3721(i__),OII_w3727(i__),
+ e_OII_w3727(i__),H13_w3734(i__),e_H13_w3734(i__),
+ H12_w3750(i__),e_H12_w3750(i__),H11_w3770(i__),
+ e_H11_w3770(i__),H10_w3797(i__),e_H10_w3797(i__),
+ HeI_w3819(i__),e_HeI_w3819(i__),H9_w3835(i__),e_H9_w3835(i__),
+ NeIII_w3868(i__),e_NeIII_w3868(i__),HeI_w3888(i__),
+ e_HeI_w3888(i__),H8_w3889(i__),e_H8_w3889(i__),HeI_w3964(i__),
+ e_HeI_w3964(i__),NeIII_w3967(i__),e_NeIII_w3967(i__),
+ H7_w3970(i__),e_H7_w3970(i__),NeIII_w4011(i__),
+ e_NeIII_w4011(i__),HeI_w4026(i__),e_HeI_w4026(i__),
+ SII_w4068(i__),e_SII_w4068(i__),SII_w4076(i__),
+ e_SII_w4076(i__),Hd_w4101(i__),e_Hd_w4101(i__),HeI_w4120(i__),
+ e_HeI_w4120(i__),HeI_w4143(i__),e_HeI_w4143(i__),
+ Hg_w4340(i__),e_Hg_w4340(i__),OIII_w4363(i__),
+ e_OIII_w4363(i__),HeI_w4387(i__),e_HeI_w4387(i__),
+ HeI_w4471(i__),e_HeI_w4471(i__),FeIII_w4658(i__),
+ e_FeIII_w4658(i__),HeII_w4686(i__),e_HeII_w4686(i__),
+ Hb_w4861(i__),e_Hb_w4861(i__),HeI_w4921(i__),e_HeI_w4921(i__),
+ OIII_w4959(i__),e_OIII_w4959(i__),OIII_w5007(i__),
+ e_OIII_w5007(i__),HeI_w5015(i__),e_HeI_w5015(i__),
+ NII_w5755(i__),e_NII_w5755(i__),HeI_w5876(i__),
+ e_HeI_w5876(i__),OI_w6300(i__),e_OI_w6300(i__),
+ SIII_w6312(i__),e_SIII_w6312(i__),OI_w6363(i__),
+ e_OI_w6363(i__),NII_w6548(i__),e_NII_w6548(i__),Ha_w6563(i__),
+ e_Ha_w6563(i__),NII_w6584(i__),e_NII_w6584(i__),
+ HeI_w6678(i__),e_HeI_w6678(i__),SII_w6717(i__),
+ e_SII_w6717(i__),SII_w6731(i__),e_SII_w6731(i__),
+ HeI_w7065(i__),e_HeI_w7065(i__),ArIII_w7135(i__),
+ e_ArIII_w7135(i__),OII_w7320(i__),e_OII_w7320(i__),
+ OII_w7330(i__),e_OII_w7330(i__),ArIII_w7751(i__),
+ e_ArIII_w7751(i__),P13_w8665(i__),e_P13_w8665(i__),
+ P12_w8750(i__),e_P12_w8750(i__),P11_w8862(i__),
+ e_P11_w8862(i__),P10_w9015(i__),e_P10_w9015(i__),
+ SIII_w9069(i__),e_SIII_w9069(i__),P9_w9229(i__),
+ e_P9_w9229(i__),SIII_w9532(i__),e_SIII_w9532(i__),
+ P8_w9546(i__),e_P8_w9546(i__),E_B_V(i__),e_E_B_V(i__),
+ FHb(i__),e_FHb(i__)
if(ar__1(69:73) .EQ. '') H11_w3770(i__) = rNULL__
if(ar__1(75:79) .EQ. '') e_H11_w3770(i__) = rNULL__
if(ar__1(93:97) .EQ. '') HeI_w3819(i__) = rNULL__
if(ar__1(99:103) .EQ. '') e_HeI_w3819(i__) = rNULL__
if(ar__1(117:121) .EQ. '') NeIII_w3868(i__) = rNULL__
if(ar__1(123:127) .EQ. '') e_NeIII_w3868(i__) = rNULL__
if(ar__1(153:157) .EQ. '') HeI_w3964(i__) = rNULL__
if(ar__1(159:163) .EQ. '') e_HeI_w3964(i__) = rNULL__
if(ar__1(165:169) .EQ. '') NeIII_w3967(i__) = rNULL__
if(ar__1(171:175) .EQ. '') e_NeIII_w3967(i__) = rNULL__
if(ar__1(189:193) .EQ. '') NeIII_w4011(i__) = rNULL__
if(ar__1(195:199) .EQ. '') e_NeIII_w4011(i__) = rNULL__
if(ar__1(201:205) .EQ. '') HeI_w4026(i__) = rNULL__
if(ar__1(207:211) .EQ. '') e_HeI_w4026(i__) = rNULL__
if(ar__1(213:217) .EQ. '') SII_w4068(i__) = rNULL__
if(ar__1(219:223) .EQ. '') e_SII_w4068(i__) = rNULL__
if(ar__1(225:229) .EQ. '') SII_w4076(i__) = rNULL__
if(ar__1(231:235) .EQ. '') e_SII_w4076(i__) = rNULL__
if(ar__1(249:253) .EQ. '') HeI_w4120(i__) = rNULL__
if(ar__1(255:259) .EQ. '') e_HeI_w4120(i__) = rNULL__
if(ar__1(261:265) .EQ. '') HeI_w4143(i__) = rNULL__
if(ar__1(267:271) .EQ. '') e_HeI_w4143(i__) = rNULL__
if(ar__1(285:289) .EQ. '') OIII_w4363(i__) = rNULL__
if(ar__1(291:295) .EQ. '') e_OIII_w4363(i__) = rNULL__
if(ar__1(297:301) .EQ. '') HeI_w4387(i__) = rNULL__
if(ar__1(303:307) .EQ. '') e_HeI_w4387(i__) = rNULL__
if(ar__1(309:313) .EQ. '') HeI_w4471(i__) = rNULL__
if(ar__1(315:319) .EQ. '') e_HeI_w4471(i__) = rNULL__
if(ar__1(321:325) .EQ. '') FeIII_w4658(i__) = rNULL__
if(ar__1(327:331) .EQ. '') e_FeIII_w4658(i__) = rNULL__
if(ar__1(333:337) .EQ. '') HeII_w4686(i__) = rNULL__
if(ar__1(339:343) .EQ. '') e_HeII_w4686(i__) = rNULL__
if(ar__1(357:361) .EQ. '') HeI_w4921(i__) = rNULL__
if(ar__1(363:367) .EQ. '') e_HeI_w4921(i__) = rNULL__
if(ar__1(405:409) .EQ. '') NII_w5755(i__) = rNULL__
if(ar__1(411:415) .EQ. '') e_NII_w5755(i__) = rNULL__
if(ar__1(429:433) .EQ. '') OI_w6300(i__) = rNULL__
if(ar__1(435:439) .EQ. '') e_OI_w6300(i__) = rNULL__
if(ar__1(441:445) .EQ. '') SIII_w6312(i__) = rNULL__
if(ar__1(447:451) .EQ. '') e_SIII_w6312(i__) = rNULL__
if(ar__1(453:457) .EQ. '') OI_w6363(i__) = rNULL__
if(ar__1(459:463) .EQ. '') e_OI_w6363(i__) = rNULL__
if(ar__1(537:541) .EQ. '') HeI_w7065(i__) = rNULL__
if(ar__1(543:547) .EQ. '') e_HeI_w7065(i__) = rNULL__
if(ar__1(549:553) .EQ. '') ArIII_w7135(i__) = rNULL__
if(ar__1(555:559) .EQ. '') e_ArIII_w7135(i__) = rNULL__
if(ar__1(573:577) .EQ. '') OII_w7330(i__) = rNULL__
if(ar__1(579:583) .EQ. '') e_OII_w7330(i__) = rNULL__
if(ar__1(585:589) .EQ. '') ArIII_w7751(i__) = rNULL__
if(ar__1(591:595) .EQ. '') e_ArIII_w7751(i__) = rNULL__
if(ar__1(597:601) .EQ. '') P13_w8665(i__) = rNULL__
if(ar__1(603:607) .EQ. '') e_P13_w8665(i__) = rNULL__
if(ar__1(609:613) .EQ. '') P12_w8750(i__) = rNULL__
if(ar__1(615:619) .EQ. '') e_P12_w8750(i__) = rNULL__
if(ar__1(621:625) .EQ. '') P11_w8862(i__) = rNULL__
if(ar__1(627:631) .EQ. '') e_P11_w8862(i__) = rNULL__
if(ar__1(633:637) .EQ. '') P10_w9015(i__) = rNULL__
if(ar__1(639:643) .EQ. '') e_P10_w9015(i__) = rNULL__
if(ar__1(645:649) .EQ. '') SIII_w9069(i__) = rNULL__
if(ar__1(651:655) .EQ. '') e_SIII_w9069(i__) = rNULL__
if(ar__1(657:661) .EQ. '') P9_w9229(i__) = rNULL__
if(ar__1(663:667) .EQ. '') e_P9_w9229(i__) = rNULL__
if(ar__1(681:685) .EQ. '') P8_w9546(i__) = rNULL__
if(ar__1(687:691) .EQ. '') e_P8_w9546(i__) = rNULL__
c ..............Just test output...........
write(6,2)
+ Reg_1(i__),H14_w3721(i__),e_H14_w3721(i__),OII_w3727(i__),
+ e_OII_w3727(i__),H13_w3734(i__),e_H13_w3734(i__),
+ H12_w3750(i__),e_H12_w3750(i__),H11_w3770(i__),
+ e_H11_w3770(i__),H10_w3797(i__),e_H10_w3797(i__),
+ HeI_w3819(i__),e_HeI_w3819(i__),H9_w3835(i__),e_H9_w3835(i__),
+ NeIII_w3868(i__),e_NeIII_w3868(i__),HeI_w3888(i__),
+ e_HeI_w3888(i__),H8_w3889(i__),e_H8_w3889(i__),HeI_w3964(i__),
+ e_HeI_w3964(i__),NeIII_w3967(i__),e_NeIII_w3967(i__),
+ H7_w3970(i__),e_H7_w3970(i__),NeIII_w4011(i__),
+ e_NeIII_w4011(i__),HeI_w4026(i__),e_HeI_w4026(i__),
+ SII_w4068(i__),e_SII_w4068(i__),SII_w4076(i__),
+ e_SII_w4076(i__),Hd_w4101(i__),e_Hd_w4101(i__),HeI_w4120(i__),
+ e_HeI_w4120(i__),HeI_w4143(i__),e_HeI_w4143(i__),
+ Hg_w4340(i__),e_Hg_w4340(i__),OIII_w4363(i__),
+ e_OIII_w4363(i__),HeI_w4387(i__),e_HeI_w4387(i__),
+ HeI_w4471(i__),e_HeI_w4471(i__),FeIII_w4658(i__),
+ e_FeIII_w4658(i__),HeII_w4686(i__),e_HeII_w4686(i__),
+ Hb_w4861(i__),e_Hb_w4861(i__),HeI_w4921(i__),e_HeI_w4921(i__),
+ OIII_w4959(i__),e_OIII_w4959(i__),OIII_w5007(i__),
+ e_OIII_w5007(i__),HeI_w5015(i__),e_HeI_w5015(i__),
+ NII_w5755(i__),e_NII_w5755(i__),HeI_w5876(i__),
+ e_HeI_w5876(i__),OI_w6300(i__),e_OI_w6300(i__),
+ SIII_w6312(i__),e_SIII_w6312(i__),OI_w6363(i__),
+ e_OI_w6363(i__),NII_w6548(i__),e_NII_w6548(i__),Ha_w6563(i__),
+ e_Ha_w6563(i__),NII_w6584(i__),e_NII_w6584(i__),
+ HeI_w6678(i__),e_HeI_w6678(i__),SII_w6717(i__),
+ e_SII_w6717(i__),SII_w6731(i__),e_SII_w6731(i__),
+ HeI_w7065(i__),e_HeI_w7065(i__),ArIII_w7135(i__),
+ e_ArIII_w7135(i__),OII_w7320(i__),e_OII_w7320(i__),
+ OII_w7330(i__),e_OII_w7330(i__),ArIII_w7751(i__),
+ e_ArIII_w7751(i__),P13_w8665(i__),e_P13_w8665(i__),
+ P12_w8750(i__),e_P12_w8750(i__),P11_w8862(i__),
+ e_P11_w8862(i__),P10_w9015(i__),e_P10_w9015(i__),
+ SIII_w9069(i__),e_SIII_w9069(i__),P9_w9229(i__),
+ e_P9_w9229(i__),SIII_w9532(i__),e_SIII_w9532(i__),
+ P8_w9546(i__),e_P8_w9546(i__),E_B_V(i__),e_E_B_V(i__),
+ FHb(i__),e_FHb(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tablea3.dat' ! Ionic and total abundances for LBT observations
* of NGC 3184
C Format for file interpretation
3 format(
+ A19,1X,I5,1X,I3,1X,I4,1X,I4,1X,I5,1X,I3,1X,I4,1X,I3,1X,I4,1X,
+ I3,1X,I5,1X,I2,1X,I4,1X,I4,1X,I4,1X,I4,1X,I5,1X,I4,1X,I3,1X,
+ I2,1X,F4.1,1X,F4.1,1X,F4.1,1X,F4.1,1X,F5.3,1X,F5.3,1X,F5.1,1X,
+ F4.1,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,F6.3,1X,F5.3,1X,F4.1,
+ 1X,F4.1,1X,F5.1,1X,F5.1,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,1X,
+ F6.3,1X,F5.3,1X,F6.1,1X,F6.1,1X,F5.3,1X,F5.3,1X,F5.3,1X,F5.3,
+ 1X,F6.3,1X,F5.3,1X,F4.1,1X,F4.1,1X,F5.3,1X,F5.3,1X,F5.3,1X,
+ F5.3,1X,F6.3,1X,F5.3)
C Effective file loading
open(unit=1,status='old',file=
+'tablea3.dat')
write(6,*) '....Loading file: tablea3.dat'
do i__=1,30
read(1,'(A351)')ar__2
read(ar__2,3)
+ Reg_2(i__),Te_SII_(i__),e_Te_SII_(i__),Te_NII_(i__),
+ e_Te_NII_(i__),Te_OII_(i__),e_Te_OII_(i__),Te_SIII_(i__),
+ e_Te_SIII_(i__),Te_OIII_(i__),e_Te_OIII_(i__),nem(i__),
+ e_nem(i__),Te_low(i__),e_Te_low(i__),Te_int(i__),
+ e_Te_int(i__),Te_high(i__),e_Te_high(i__),neu(i__),e_neu(i__),
+ O__H_(i__),e_O__H_(i__),O_2_H_(i__),e_O_2_H_(i__),
+ log_O_H(i__),e_log_O_H(i__),N__H_(i__),e_N__H_(i__),
+ ICF_N(i__),e_ICF_N(i__),log_N_H(i__),e_log_N_H(i__),
+ log_N_O(i__),e_log_N_O(i__),S__H_(i__),e_S__H_(i__),
+ S2_H_(i__),e_S2_H_(i__),ICF_S(i__),e_ICF_S(i__),log_S_H(i__),
+ e_log_S_H(i__),log_S_O(i__),e_log_S_O(i__),Ne_2_H_(i__),
+ e_Ne_2_H_(i__),ICF_Ne(i__),e_ICF_Ne(i__),log_Ne_H(i__),
+ e_log_Ne_H(i__),log_Ne_O(i__),e_log_Ne_O(i__),Ar_2_H_(i__),
+ e_Ar_2_H_(i__),ICF_Ar(i__),e_ICF_Ar(i__),log_Ar_H(i__),
+ e_log_Ar_H(i__),log_Ar_O(i__),e_log_Ar_O(i__)
if(ar__2(21:25) .EQ. '') Te_SII_(i__) = iNULL__
if(ar__2(27:29) .EQ. '') e_Te_SII_(i__) = iNULL__
if(ar__2(31:34) .EQ. '') Te_NII_(i__) = iNULL__
if(ar__2(36:39) .EQ. '') e_Te_NII_(i__) = iNULL__
if(ar__2(41:45) .EQ. '') Te_OII_(i__) = iNULL__
if(ar__2(47:49) .EQ. '') e_Te_OII_(i__) = iNULL__
if(ar__2(51:54) .EQ. '') Te_SIII_(i__) = iNULL__
if(ar__2(56:58) .EQ. '') e_Te_SIII_(i__) = iNULL__
if(ar__2(60:63) .EQ. '') Te_OIII_(i__) = iNULL__
if(ar__2(65:67) .EQ. '') e_Te_OIII_(i__) = iNULL__
if(ar__2(255:260) .EQ. '') Ne_2_H_(i__) = rNULL__
if(ar__2(262:267) .EQ. '') e_Ne_2_H_(i__) = rNULL__
if(ar__2(281:285) .EQ. '') log_Ne_H(i__) = rNULL__
if(ar__2(287:291) .EQ. '') e_log_Ne_H(i__) = rNULL__
if(ar__2(293:298) .EQ. '') log_Ne_O(i__) = rNULL__
if(ar__2(300:304) .EQ. '') e_log_Ne_O(i__) = rNULL__
c ..............Just test output...........
write(6,3)
+ Reg_2(i__),Te_SII_(i__),e_Te_SII_(i__),Te_NII_(i__),
+ e_Te_NII_(i__),Te_OII_(i__),e_Te_OII_(i__),Te_SIII_(i__),
+ e_Te_SIII_(i__),Te_OIII_(i__),e_Te_OIII_(i__),nem(i__),
+ e_nem(i__),Te_low(i__),e_Te_low(i__),Te_int(i__),
+ e_Te_int(i__),Te_high(i__),e_Te_high(i__),neu(i__),e_neu(i__),
+ O__H_(i__),e_O__H_(i__),O_2_H_(i__),e_O_2_H_(i__),
+ log_O_H(i__),e_log_O_H(i__),N__H_(i__),e_N__H_(i__),
+ ICF_N(i__),e_ICF_N(i__),log_N_H(i__),e_log_N_H(i__),
+ log_N_O(i__),e_log_N_O(i__),S__H_(i__),e_S__H_(i__),
+ S2_H_(i__),e_S2_H_(i__),ICF_S(i__),e_ICF_S(i__),log_S_H(i__),
+ e_log_S_H(i__),log_S_O(i__),e_log_S_O(i__),Ne_2_H_(i__),
+ e_Ne_2_H_(i__),ICF_Ne(i__),e_ICF_Ne(i__),log_Ne_H(i__),
+ e_log_Ne_H(i__),log_Ne_O(i__),e_log_Ne_O(i__),Ar_2_H_(i__),
+ e_Ar_2_H_(i__),ICF_Ar(i__),e_ICF_Ar(i__),log_Ar_H(i__),
+ e_log_Ar_H(i__),log_Ar_O(i__),e_log_Ar_O(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tableb1.dat' ! Updated abundances for NGC 628
* using ionization-based Te
C Format for file interpretation
4 format(
+ A18,1X,F4.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,
+ 1X,F4.2,1X,F5.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'tableb1.dat')
write(6,*) '....Loading file: tableb1.dat'
do i__=1,45
read(1,'(A72)')ar__3
read(ar__3,4)
+ v_BSC2015_(i__),log_O_H_1(i__),e_log_O_H_1(i__),
+ log_N_O_1(i__),e_log_N_O_1(i__),log_S_O_1(i__),
+ e_log_S_O_1(i__),log_Ar_O_1(i__),e_log_Ar_O_1(i__),
+ log_Ne_O_1(i__),e_log_Ne_O_1(i__)
if(ar__3(63:67) .EQ. '') log_Ne_O_1(i__) = rNULL__
if(ar__3(69:72) .EQ. '') e_log_Ne_O_1(i__) = rNULL__
c ..............Just test output...........
write(6,4)
+ v_BSC2015_(i__),log_O_H_1(i__),e_log_O_H_1(i__),
+ log_N_O_1(i__),e_log_N_O_1(i__),log_S_O_1(i__),
+ e_log_S_O_1(i__),log_Ar_O_1(i__),e_log_Ar_O_1(i__),
+ log_Ne_O_1(i__),e_log_Ne_O_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tableb2.dat' ! Updated abundances for NGC 5194
* using ionization-based Te
C Format for file interpretation
5 format(
+ A19,1X,F4.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,
+ 1X,F4.2,1X,F5.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'tableb2.dat')
write(6,*) '....Loading file: tableb2.dat'
do i__=1,28
read(1,'(A73)')ar__4
read(ar__4,5)
+ v_CPB2015_(i__),log_O_H_1(i__),e_log_O_H_1(i__),
+ log_N_O_1(i__),e_log_N_O_1(i__),log_S_O_1(i__),
+ e_log_S_O_1(i__),log_Ar_O_1(i__),e_log_Ar_O_1(i__),
+ log_Ne_O_1(i__),e_log_Ne_O_1(i__)
if(ar__4(64:68) .EQ. '') log_Ne_O_1(i__) = rNULL__
if(ar__4(70:73) .EQ. '') e_log_Ne_O_1(i__) = rNULL__
c ..............Just test output...........
write(6,5)
+ v_CPB2015_(i__),log_O_H_1(i__),e_log_O_H_1(i__),
+ log_N_O_1(i__),e_log_N_O_1(i__),log_S_O_1(i__),
+ e_log_S_O_1(i__),log_Ar_O_1(i__),e_log_Ar_O_1(i__),
+ log_Ne_O_1(i__),e_log_Ne_O_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'tableb3.dat' ! Updated abundances for NGC 5457
* using ionization-based Te
C Format for file interpretation
6 format(
+ A19,1X,F4.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,1X,F4.2,1X,F5.2,
+ 1X,F4.2,1X,F5.2,1X,F4.2)
C Effective file loading
open(unit=1,status='old',file=
+'tableb3.dat')
write(6,*) '....Loading file: tableb3.dat'
do i__=1,73
read(1,'(A73)')ar__5
read(ar__5,6)
+ v_CPB2016_(i__),log_O_H_1(i__),e_log_O_H_1(i__),
+ log_N_O_1(i__),e_log_N_O_1(i__),log_S_O_1(i__),
+ e_log_S_O_1(i__),log_Ar_O_1(i__),e_log_Ar_O_1(i__),
+ log_Ne_O_1(i__),e_log_Ne_O_1(i__)
if(ar__5(53:57) .EQ. '') log_Ar_O_1(i__) = rNULL__
if(ar__5(59:62) .EQ. '') e_log_Ar_O_1(i__) = rNULL__
if(ar__5(64:68) .EQ. '') log_Ne_O_1(i__) = rNULL__
if(ar__5(70:73) .EQ. '') e_log_Ne_O_1(i__) = rNULL__
c ..............Just test output...........
write(6,6)
+ v_CPB2016_(i__),log_O_H_1(i__),e_log_O_H_1(i__),
+ log_N_O_1(i__),e_log_N_O_1(i__),log_S_O_1(i__),
+ e_log_S_O_1(i__),log_Ar_O_1(i__),e_log_Ar_O_1(i__),
+ log_Ne_O_1(i__),e_log_Ne_O_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end