Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/893/96 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-Oct-09
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 ! (deg) Right Ascension J2000
real*8 DEdeg ! (deg) Declination J2000
C ---------------------------------- ! (position vector(s) in degrees)
character*19 Reg ! NGC 3184 HII region slit ID (1)
integer*4 RAh ! (h) [10] Hour of Right Ascension (J2000)
integer*4 RAm ! (min) [18] Minute of Right Ascension (J2000)
real*4 RAs ! (s) [1.5/27.8] Second of Right Ascension (J2000)
integer*4 DEd ! (deg) [41] Degree of Declination (J2000)
integer*4 DEm ! (arcmin) [22/27] Arcminute of Declination (J2000)
real*4 DEs ! (arcsec) Arcsecond of Declination (J2000)
real*4 Rgarc ! (arcsec) [2/182.7] De-projected HII region distance
* galaxy center
real*4 R_R25 ! [0.009/0.9] Fraction of R25 radius
real*4 Rgkpc ! (kpc) [0.1/10.4] HII region distance in kpc
real*4 oRA ! (arcsec) [-172.5/123.8] RA offset from the central position
* (RAJ2000=10:18:16.86,DEJ2000=+41:25:26.59)
real*4 oDE ! (arcsec) [-148.2/148.8] DE offset from the central position
* (RAJ2000=10:18:16.86,DEJ2000=+41:25:26.59)
integer*4 dOIII ! [1]? Detection, [OIII] w4363 auroral line (2)
integer*4 dNII ! [1]? Detection, [NII] w5755 auroral line (2)
integer*4 dSIII ! [1]? Detection, [SIII] w6312 auroral line (2)
integer*4 dOII ! [1]? Detection, [OII] w7320,w7330 auroral lines (2)
integer*4 dSII ! [1]? Detection, [SII] w4068,w4076 auroral lines (2)
integer*4 dWR ! [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 ! NGC 3184 HII region slit ID (1)
real*4 H14_w3721 ! [0.009/0.025] Line intensity, H14
* w3721 (2)
real*4 e_H14_w3721 ! [0.001/0.02] Uncertainty in H14-w3721
real*4 OII_w3727 ! [0.3/4.2] Line intensity, [OII]
* w3727 (2)
real*4 e_OII_w3727 ! [0.007/0.8] Uncertainty in OII-w3727
real*4 H13_w3734 ! [0.01/0.031] Line intensity, H13
* w3734 (2)
real*4 e_H13_w3734 ! [0.001/0.02] Uncertainty in H13-w3734
real*4 H12_w3750 ! [0.003/0.06] Line intensity, H12
* w3750 (2)
real*4 e_H12_w3750 ! [0.001/0.015] Uncertainty in H12-w3750
real*4 H11_w3770 ! [0.019/0.09]? Line intensity, H11
* w3770 (2)
real*4 e_H11_w3770 ! [0.001/0.04]? Uncertainty in H11-w3770
real*4 H10_w3797 ! [0.02/0.07] Line intensity, H10
* w3797 (2)
real*4 e_H10_w3797 ! [0.001/0.04] Uncertainty in H10-w3797
real*4 HeI_w3819 ! [0.001/0.05]? Line intensity, HeI
* w3819 (2)
real*4 e_HeI_w3819 ! [0.001/0.01]? Uncertainty in HeI-w3819
real*4 H9_w3835 ! [0.02/0.11] Line intensity, H9 w3835
* (2)
real*4 e_H9_w3835 ! [0.001/0.07] Uncertainty in H9-w3835
real*4 NeIII_w3868 ! [0.003/0.2]? Line intensity, [NeIII]
* w3868 (2)
real*4 e_NeIII_w3868 ! [0.001/0.01]? Uncertainty in
* NeIII-w3868
real*4 HeI_w3888 ! [0.01/0.2] Line intensity, HeI w3888
* (2)
real*4 e_HeI_w3888 ! [0.001/0.07] Uncertainty in HeI-w3888
real*4 H8_w3889 ! [0.04/0.2] Line intensity, H8 w3889
* (2)
real*4 e_H8_w3889 ! [0.001/0.08] Uncertainty in H8-w3889
real*4 HeI_w3964 ! [0.001/0.03]? Line intensity, HeI
* w3964 (2)
real*4 e_HeI_w3964 ! [0.001/0.02]? Uncertainty in HeI-w3964
real*4 NeIII_w3967 ! [0.003/0.2]? Line intensity, [NeIII]
* w3967 (2)
real*4 e_NeIII_w3967 ! [0.001/0.04]? Uncertainty in
* NeIII-w3967
real*4 H7_w3970 ! [0.07/0.2] Line intensity, H7 w3970
* (2)
real*4 e_H7_w3970 ! [0.002/0.11] Uncertainty in H7-w3970
real*4 NeIII_w4011 ! [0.001/0.03]? Line intensity,
* [NeIII] w4011 (2)
real*4 e_NeIII_w4011 ! [0.001/0.006]? Uncertainty in
* NeIII-w4011
real*4 HeI_w4026 ! [0.001/0.04]? Line intensity, HeI
* w4026 (2)
real*4 e_HeI_w4026 ! [0.001/0.004]? Uncertainty in
* HeI-w4026
real*4 SII_w4068 ! [0.001/0.041]? Line intensity, [SII]
* w4068 (2)
real*4 e_SII_w4068 ! [0.001/0.01]? Uncertainty in SII-w4068
real*4 SII_w4076 ! [0.002/0.03]? Line intensity, [SII]
* w4076 (2)
real*4 e_SII_w4076 ! [0.001/0.007]? Uncertainty in
* SII-w4076
real*4 Hd_w4101 ! [0.16/0.4] Line intensity, H{delta}
* w4101 (2)
real*4 e_Hd_w4101 ! [0.001/0.3] Uncertainty in Hd-w4101
real*4 HeI_w4120 ! [0.001/0.018]? Line intensity, HeI
* w4120 (2)
real*4 e_HeI_w4120 ! [0.001/0.009]? Uncertainty in
* HeI-w4120
real*4 HeI_w4143 ! [0.001/0.013]? Line intensity, HeI
* w4143 (2)
real*4 e_HeI_w4143 ! [0.001/0.004]? Uncertainty in
* HeI-w4143
real*4 Hg_w4340 ! [0.39/0.5] Line intensity, H{gamma}
* w4340 (2)
real*4 e_Hg_w4340 ! [0.002/0.5] Uncertainty in Hg-w4340
real*4 OIII_w4363 ! [0.001/0.014]? Line intensity,
* [OIII] w4363 (2)
real*4 e_OIII_w4363 ! [0.001/0.011]? Uncertainty in
* OIII-w4363
real*4 HeI_w4387 ! [0.001/0.02]? Line intensity, HeI
* w4387 (2)
real*4 e_HeI_w4387 ! [0.001/0.004]? Uncertainty in
* HeI-w4387
real*4 HeI_w4471 ! [0.001/0.041]? Line intensity, HeI
* w4471 (2)
real*4 e_HeI_w4471 ! [0.001/0.006]? Uncertainty in
* HeI-w4471
real*4 FeIII_w4658 ! [0.001/0.03]? Line intensity,
* [FeIII] w4658 (2)
real*4 e_FeIII_w4658 ! [0.001/0.03]? Uncertainty in
* FeIII-w4658
real*4 HeII_w4686 ! [0.001/0.03]? Line intensity, HeII
* w4686 (2)
real*4 e_HeII_w4686 ! [0.001/0.03]? Uncertainty in
* HeII-w4686
real*4 Hb_w4861 ! [1] Line intensity, H{beta} w4861 (2)
real*4 e_Hb_w4861 ! [0.003/1] Uncertainty in Hb-w4861
real*4 HeI_w4921 ! [0.001/0.02]? Line intensity, HeI
* w4921 (2)
real*4 e_HeI_w4921 ! [0.001/0.03]? Uncertainty in HeI-w4921
real*4 OIII_w4959 ! [0.008/1.1] Line intensity, [OIII]
* w4959 (2)
real*4 e_OIII_w4959 ! [0.001/0.05] Uncertainty in OIII-w4959
real*4 OIII_w5007 ! [0.03/3] Line intensity, [OIII]
* w5007 (2)
real*4 e_OIII_w5007 ! [0.001/0.12] Uncertainty in OIII-w5007
real*4 HeI_w5015 ! [0.001/0.06] Line intensity, HeI
* w5015 (2)
real*4 e_HeI_w5015 ! [0.001/0.03] Uncertainty in HeI-w5015
real*4 NII_w5755 ! [0.001/0.02]? Line intensity, [NII]
* w5755 (2)
real*4 e_NII_w5755 ! [0.001/0.015]? Uncertainty in
* NII-w5755
real*4 HeI_w5876 ! [0.03/0.4] Line intensity, HeI w5876
* (2)
real*4 e_HeI_w5876 ! [0.001/0.07] Uncertainty in HeI-w5876
real*4 OI_w6300 ! [0.003/0.06]? Line intensity, [OI]
* w6300 (2)
real*4 e_OI_w6300 ! [0.001/0.03]? Uncertainty in OI-w6300
real*4 SIII_w6312 ! [0.001/0.02]? Line intensity, [SIII]
* w6312 (2)
real*4 e_SIII_w6312 ! [0.001/0.008]? Uncertainty in
* SIII-w6312
real*4 OI_w6363 ! [0.001/0.04]? Line intensity, [OI]
* w6363 (2)
real*4 e_OI_w6363 ! [0.001/0.01]? Uncertainty in OI-w6363
real*4 NII_w6548 ! [0.13/0.41] Line intensity, [NII]
* w6548 (2)
real*4 e_NII_w6548 ! [0.004/0.31] Uncertainty in NII-w6548
real*4 Ha_w6563 ! [2.6/3.2] Line intensity, H{alpha}
* w6563 (2)
real*4 e_Ha_w6563 ! [0.02/3] Uncertainty in Ha-w6563
real*4 NII_w6584 ! [0.4/1.3] Line intensity, [NII]
* w6584 (2)
real*4 e_NII_w6584 ! [0.005/1] Uncertainty in NII-w6584
real*4 HeI_w6678 ! [0.006/0.04] Line intensity, HeI
* w6678 (2)
real*4 e_HeI_w6678 ! [0.001/0.03] Uncertainty in HeI-w6678
real*4 SII_w6717 ! [0.1/0.51] Line intensity, [SII]
* w6717 (2)
real*4 e_SII_w6717 ! [0.002/0.4] Uncertainty in SII-w6717
real*4 SII_w6731 ! [0.09/0.4] Line intensity, [SII]
* w6731 (2)
real*4 e_SII_w6731 ! [0.001/0.3] Uncertainty in SII-w6731
real*4 HeI_w7065 ! [0.002/0.033]? Line intensity, HeI
* w7065 (2)
real*4 e_HeI_w7065 ! [0.001/0.009]? Uncertainty in
* HeI-w7065
real*4 ArIII_w7135 ! [0.004/0.11]? Line intensity,
* [ArIII] w7135 (2)
real*4 e_ArIII_w7135 ! [0.001/0.03]? Uncertainty in
* ArIII-w7135
real*4 OII_w7320 ! [0.001/0.03] Line intensity, [OII]
* w7320 (2)
real*4 e_OII_w7320 ! [0.001/0.04] Uncertainty in OII-w7320
real*4 OII_w7330 ! [0.001/0.04]? Line intensity, [OII]
* w7330 (2)
real*4 e_OII_w7330 ! [0.001/0.04]? Uncertainty in OII-w7330
real*4 ArIII_w7751 ! [0.003/0.03]? Line intensity,
* [ArIII] w7751 (2)
real*4 e_ArIII_w7751 ! [0.001/0.008]? Uncertainty in
* ArIII-w7751
real*4 P13_w8665 ! [0.001/0.02]? Line intensity, P13
* w8665 (2)
real*4 e_P13_w8665 ! [0.001/0.03]? Uncertainty in P13-w8665
real*4 P12_w8750 ! [0.004/2.6]? Line intensity, P12
* w8750 (2)
real*4 e_P12_w8750 ! [0.002/0.05]? Uncertainty in P12-w8750
real*4 P11_w8862 ! [0.001/0.08]? Line intensity, P11
* w8862 (2)
real*4 e_P11_w8862 ! [0.002/0.04]? Uncertainty in P11-w8862
real*4 P10_w9015 ! [0.002/0.05]? Line intensity, P10
* w9015 (2)
real*4 e_P10_w9015 ! [0.001/0.03]? Uncertainty in
* P10-w9015
real*4 SIII_w9069 ! [0.05/0.4]? Line intensity, [SIII]
* w9069 (2)
real*4 e_SIII_w9069 ! [0.001/0.2]? Uncertainty in SIII-w9069
real*4 P9_w9229 ! [0.002/0.04]? Line intensity, P9
* w9229 (2)
real*4 e_P9_w9229 ! [0.001/0.03]? Uncertainty in P9-w9229
real*4 SIII_w9532 ! [0.18/1] Line intensity, [SIII]
* w9532 (2)
real*4 e_SIII_w9532 ! [0.003/0.6] Uncertainty in SIII-w9532
real*4 P8_w9546 ! [0.01/0.2]? Line intensity, P8 w9546
* (2)
real*4 e_P8_w9546 ! [0.002/0.03]? Uncertainty in P8-w9546
real*4 E_B_V ! (mag) [0.09/0.7] E(B-V) reddening (3)
real*4 e_E_B_V ! (mag) [0.004/0.03] Uncertainty in E (B-V)
real*4 FHb ! (10-16mW/m2) [9.6/314] FH{beta} line flux,
* 10^-16^erg/s/cm^2^
real*4 e_FHb ! (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 ! NGC 3184 HII region slit ID (1)
integer*4 Te_SII_ ! (K) [5700/23900]? Measured electron temperature
* from the [SII] ion
integer*4 e_Te_SII_ ! (K) [400]? Uncertainty in Te[SII]
integer*4 Te_NII_ ! (K) [6000/9900]? Measured electron temperature
* from the [NII] ion
integer*4 e_Te_NII_ ! (K) [1700]? Uncertainty in Te[NII]
integer*4 Te_OII_ ! (K) [6800/12500]? Measured electron temperature
* from the [OII] ion
integer*4 e_Te_OII_ ! (K) [400]? Uncertainty in Te[OII]
integer*4 Te_SIII_ ! (K) [6400/9400]? Measured electron temperature
* from the [SIII] ion
integer*4 e_Te_SIII_ ! (K) [600/600]? Uncertainty in Te[SIII]
integer*4 Te_OIII_ ! (K) [8600/9600]? Measured electron temperature
* from the [OIII] ion
integer*4 e_Te_OIII_ ! (K) [100]? Uncertainty in Te[OIII]
integer*4 nem ! (cm-3) [10/29500] Measured electron density
integer*4 e_nem ! (cm-3) [0/10] Uncertainty in nem
integer*4 Te_low ! (K) [6000/9900] Adopted ionization-based
* temperature, lower ionization zone (2)
integer*4 e_Te_low ! (K) [100/1900] Uncertainty in Te-low
integer*4 Te_int ! (K) [4700/9900] Adopted ionization-based
* temperature, intermediate ionization zone (2)
integer*4 e_Te_int ! (K) [100/2600] Uncertainty in Te-int
integer*4 Te_high ! (K) [4800/10300] Adopted ionization-based
* temperature, upper ionization zone (2)
integer*4 e_Te_high ! (K) [100/2700] Uncertainty in Te-high
integer*4 neu ! (cm-3) [10/100] Used electron density
integer*4 e_neu ! (cm-3) [10] Uncertainty in neu
real*4 O__H_ ! (10-5) [8.1/62.8] Abundance ratio, O+/H+
real*4 e_O__H_ ! (10-5) [1.4/32] Uncertainty in O+/H+
real*4 O_2_H_ ! (10-5) [1.8/57.7] Abundance ratio, O+2/H+
real*4 e_O_2_H_ ! (10-5) [0.2/67.3] Uncertainty in O+2/H+
real*4 log_O_H ! ([-]) [8.2/9] Oxygen abundance, 12+log (O/H)
real*4 e_log_O_H ! ([-]) [0.02/0.3] Uncertainty in log (O/H)
real*4 N__H_ ! (10-6) [10.2/110.4] Abundance ratio, N+/H+
real*4 e_N__H_ ! (10-6) [1.6/52.5] Uncertainty in N+/H+
real*4 ICF_N ! [1/2.7] Nitrogen ionization correction factor
real*4 e_ICF_N ! [0.09/2.91] Uncertainty in ICF (N)
real*4 log_N_H ! ([-]) [7.3/8.5] Nitrogen-to-hydrogen abundance
* ratio, log
real*4 e_log_N_H ! ([-]) [0.038/0.4] Uncertainty in log (N/H)
real*4 log_N_O ! ([-]) [-1.2/-0.4] Nitrogen-to-oxygen abundance
* ratio, log
real*4 e_log_N_O ! ([-]) [0.03/0.3] Uncertainty in log (N/O)
real*4 S__H_ ! (10-7) [6.2/86.6] Abundance ratio, S+/H+
real*4 e_S__H_ ! (10-7) [1.2/40] Uncertainty in S+/H+
real*4 S2_H_ ! (10-7) [41/227.1] Abundance ratio, S2/H+
real*4 e_S2_H_ ! (10-7) [3.3/173.5] Uncertainty in S2/H+
real*4 ICF_S ! [1/1.8] Sulfur ionization correction factor
real*4 e_ICF_S ! [0.09/1.5] Uncertainty in ICF (S)
real*4 log_S_H ! ([-]) [6.8/7.7] Sulfur-to-hydrogen abundance
* ratio, log
real*4 e_log_S_H ! ([-]) [0.037/0.4] Uncertainty in log (S/H)
real*4 log_S_O ! ([-]) [-1.6/-1.1] Sulfur-to-oxygen abundance
* ratio, log
real*4 e_log_S_O ! ([-]) [0.04/0.4] Uncertainty in log (S/O)
real*4 Ne_2_H_ ! (10-6) [3.9/1185]? Abundance ratio, Ne+2/H+
real*4 e_Ne_2_H_ ! (10-6) [0.7/1519]? Uncertainty in Ne+2/H+
real*4 ICF_Ne ! [1.4/3.1] Neon ionization correction factor
real*4 e_ICF_Ne ! [0.4/9.2] Uncertainty in ICF (Ne)
real*4 log_Ne_H ! ([-]) [7.4/9.3]? Neon-to-hydrogen abundance
* ratio, log
real*4 e_log_Ne_H ! ([-]) [0.07/0.5]? Uncertainty in log (Ne/H)
real*4 log_Ne_O ! ([-]) [-0.92/0.4]? Neon-to-oxygen abundance
* ratio, log
real*4 e_log_Ne_O ! ([-]) [0.04/0.4]? Uncertainty in log (Ne/O)
real*4 Ar_2_H_ ! (10-7) [4.0/20.8] Abundance ratio, Ar+2/H+
real*4 e_Ar_2_H_ ! (10-7) [0.5/18.6] Uncertainty in Ar+2/H+
real*4 ICF_Ar ! [1.4/3.1] Argon ionization correction factor
real*4 e_ICF_Ar ! [0.14/0.31] Uncertainty in ICF (Ar)
real*4 log_Ar_H ! ([-]) [6/6.7] Argon-to-hydrogen abundance ratio,
* log
real*4 e_log_Ar_H ! ([-]) [0.04/0.3] Uncertainty in log (Ar/H)
real*4 log_Ar_O ! ([-]) [-2.61/-1.8] Argon-to-oxygen abundance
* ratio, log
real*4 e_log_Ar_O ! ([-]) [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_ ! NGC 628 HII region
real*4 log_O_H_1 ! ([-]) [8.1/8.8] Oxygen abundance, 12+log(O/H)
real*4 e_log_O_H_1 ! ([-]) [0.01/0.14] Uncertainty in log(O/H)
real*4 log_N_O_1 ! ([-]) [-1.4/-0.6] Nitrogen-to-oxygen abundance
* ratio, log
real*4 e_log_N_O_1 ! ([-]) [0.01/0.2] Uncertainty in log(N/O)
real*4 log_S_O_1 ! ([-]) [-1.5/-0.97] Sulfur-to-oxygen abundance
* ratio, log
real*4 e_log_S_O_1 ! ([-]) [0.02/0.3] Uncertainty in log(S/O)
real*4 log_Ar_O_1 ! ([-]) [-2.5/-1.8] Argon-to-oxygen abundance
* ratio, log
real*4 e_log_Ar_O_1 ! ([-]) [0.04/0.2] Uncertainty in log(Ar/O)
real*4 log_Ne_O_1 ! ([-]) [-1.1/-0.4]? Neon-to-oxygen abundance
* ratio, log
real*4 e_log_Ne_O_1 ! ([-]) [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_ ! NGC 5194 HII region
real*4 log_O_H_1 ! ([-]) [8.46/8.9] Oxygen abundance, 12+log(O/H)
real*4 e_log_O_H_1 ! ([-]) [0.03/0.16] Uncertainty in log(O/H)
real*4 log_N_O_1 ! ([-]) [-0.8/-0.42] Nitrogen-to-oxygen abundance
* ratio, log
real*4 e_log_N_O_1 ! ([-]) [0.04/0.2] Uncertainty in log(N/O)
real*4 log_S_O_1 ! ([-]) [-1.6/-1] Sulfur-to-oxygen abundance ratio,
* log
real*4 e_log_S_O_1 ! ([-]) [0.05/0.3] Uncertainty in log(S/O)
real*4 log_Ar_O_1 ! ([-]) [-2.7/-1.8] Argon-to-oxygen abundance ratio,
* log
real*4 e_log_Ar_O_1 ! ([-]) [0.05/0.24] Uncertainty in log(Ar/O)
real*4 log_Ne_O_1 ! ([-]) [-0.9/-0.6]? Neon-to-oxygen abundance ratio,
* log
real*4 e_log_Ne_O_1 ! ([-]) [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_ ! NGC 5457 HII region
* (<[CPB2016] NGC 5457+xxx.x+yyy.y> in Simbad)
real*4 log_O_H_1 ! ([-]) [7.8/8.9] Oxygen abundance, 12+log(O/H)
real*4 e_log_O_H_1 ! ([-]) [0.01/0.16] Uncertainty in log(O/H)
real*4 log_N_O_1 ! ([-]) [-1.5/-0.6] Nitrogen-to-oxygen abundance
* ratio, log
real*4 e_log_N_O_1 ! ([-]) [0.02/0.3] Uncertainty in log(N/O)
real*4 log_S_O_1 ! ([-]) [-1.7/-1] Sulfur-to-oxygen abundance ratio,
* log
real*4 e_log_S_O_1 ! ([-]) [0.02/0.26] Uncertainty in log(S/O)
real*4 log_Ar_O_1 ! ([-]) [-2.8/-2]? Argon-to-oxygen abundance ratio,
* log
real*4 e_log_Ar_O_1 ! ([-]) [0.04/0.2]? Uncertainty in log(Ar/O)
real*4 log_Ne_O_1 ! ([-]) [-1.1/-0.28]? Neon-to-oxygen abundance ratio,
* log
real*4 e_log_Ne_O_1 ! ([-]) [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,RAh,RAm,RAs,DEd,DEm,DEs,Rgarc,R_R25,Rgkpc,oRA,oDE,dOIII,
+ dNII,dSIII,dOII,dSII,dWR
if(ar__(77:77) .EQ. '') dOIII = iNULL__
if(ar__(79:79) .EQ. '') dNII = iNULL__
if(ar__(81:81) .EQ. '') dSIII = iNULL__
if(ar__(83:83) .EQ. '') dOII = iNULL__
if(ar__(85:85) .EQ. '') dSII = iNULL__
if(ar__(87:88) .EQ. '') dWR = iNULL__
RAdeg = rNULL__
DEdeg = rNULL__
c Derive coordinates RAdeg and DEdeg from input data
c (RAdeg and DEdeg are set to rNULL__ when unknown)
if(RAh .GT. -180) RAdeg=RAh*15.
if(RAm .GT. -180) RAdeg=RAdeg+RAm/4.
if(RAs .GT. -180) RAdeg=RAdeg+RAs/240.
if(DEd .GE. 0) DEdeg=DEd
if(DEm .GE. 0) DEdeg=DEdeg+DEm/60.
if(DEs .GE. 0) DEdeg=DEdeg+DEs/3600.
c ..............Just test output...........
write(6,1)
+ Reg,RAh,RAm,RAs,DEd,DEm,DEs,Rgarc,R_R25,Rgkpc,oRA,oDE,dOIII,
+ dNII,dSIII,dOII,dSII,dWR
write(6,'(6H Pos: 2F8.4)') RAdeg,DEdeg
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,H14_w3721,e_H14_w3721,OII_w3727,e_OII_w3727,H13_w3734,
+ e_H13_w3734,H12_w3750,e_H12_w3750,H11_w3770,e_H11_w3770,
+ H10_w3797,e_H10_w3797,HeI_w3819,e_HeI_w3819,H9_w3835,
+ e_H9_w3835,NeIII_w3868,e_NeIII_w3868,HeI_w3888,e_HeI_w3888,
+ H8_w3889,e_H8_w3889,HeI_w3964,e_HeI_w3964,NeIII_w3967,
+ e_NeIII_w3967,H7_w3970,e_H7_w3970,NeIII_w4011,e_NeIII_w4011,
+ HeI_w4026,e_HeI_w4026,SII_w4068,e_SII_w4068,SII_w4076,
+ e_SII_w4076,Hd_w4101,e_Hd_w4101,HeI_w4120,e_HeI_w4120,
+ HeI_w4143,e_HeI_w4143,Hg_w4340,e_Hg_w4340,OIII_w4363,
+ e_OIII_w4363,HeI_w4387,e_HeI_w4387,HeI_w4471,e_HeI_w4471,
+ FeIII_w4658,e_FeIII_w4658,HeII_w4686,e_HeII_w4686,Hb_w4861,
+ e_Hb_w4861,HeI_w4921,e_HeI_w4921,OIII_w4959,e_OIII_w4959,
+ OIII_w5007,e_OIII_w5007,HeI_w5015,e_HeI_w5015,NII_w5755,
+ e_NII_w5755,HeI_w5876,e_HeI_w5876,OI_w6300,e_OI_w6300,
+ SIII_w6312,e_SIII_w6312,OI_w6363,e_OI_w6363,NII_w6548,
+ e_NII_w6548,Ha_w6563,e_Ha_w6563,NII_w6584,e_NII_w6584,
+ HeI_w6678,e_HeI_w6678,SII_w6717,e_SII_w6717,SII_w6731,
+ e_SII_w6731,HeI_w7065,e_HeI_w7065,ArIII_w7135,e_ArIII_w7135,
+ OII_w7320,e_OII_w7320,OII_w7330,e_OII_w7330,ArIII_w7751,
+ e_ArIII_w7751,P13_w8665,e_P13_w8665,P12_w8750,e_P12_w8750,
+ P11_w8862,e_P11_w8862,P10_w9015,e_P10_w9015,SIII_w9069,
+ e_SIII_w9069,P9_w9229,e_P9_w9229,SIII_w9532,e_SIII_w9532,
+ P8_w9546,e_P8_w9546,E_B_V,e_E_B_V,FHb,e_FHb
if(ar__1(69:73) .EQ. '') H11_w3770 = rNULL__
if(ar__1(75:79) .EQ. '') e_H11_w3770 = rNULL__
if(ar__1(93:97) .EQ. '') HeI_w3819 = rNULL__
if(ar__1(99:103) .EQ. '') e_HeI_w3819 = rNULL__
if(ar__1(117:121) .EQ. '') NeIII_w3868 = rNULL__
if(ar__1(123:127) .EQ. '') e_NeIII_w3868 = rNULL__
if(ar__1(153:157) .EQ. '') HeI_w3964 = rNULL__
if(ar__1(159:163) .EQ. '') e_HeI_w3964 = rNULL__
if(ar__1(165:169) .EQ. '') NeIII_w3967 = rNULL__
if(ar__1(171:175) .EQ. '') e_NeIII_w3967 = rNULL__
if(ar__1(189:193) .EQ. '') NeIII_w4011 = rNULL__
if(ar__1(195:199) .EQ. '') e_NeIII_w4011 = rNULL__
if(ar__1(201:205) .EQ. '') HeI_w4026 = rNULL__
if(ar__1(207:211) .EQ. '') e_HeI_w4026 = rNULL__
if(ar__1(213:217) .EQ. '') SII_w4068 = rNULL__
if(ar__1(219:223) .EQ. '') e_SII_w4068 = rNULL__
if(ar__1(225:229) .EQ. '') SII_w4076 = rNULL__
if(ar__1(231:235) .EQ. '') e_SII_w4076 = rNULL__
if(ar__1(249:253) .EQ. '') HeI_w4120 = rNULL__
if(ar__1(255:259) .EQ. '') e_HeI_w4120 = rNULL__
if(ar__1(261:265) .EQ. '') HeI_w4143 = rNULL__
if(ar__1(267:271) .EQ. '') e_HeI_w4143 = rNULL__
if(ar__1(285:289) .EQ. '') OIII_w4363 = rNULL__
if(ar__1(291:295) .EQ. '') e_OIII_w4363 = rNULL__
if(ar__1(297:301) .EQ. '') HeI_w4387 = rNULL__
if(ar__1(303:307) .EQ. '') e_HeI_w4387 = rNULL__
if(ar__1(309:313) .EQ. '') HeI_w4471 = rNULL__
if(ar__1(315:319) .EQ. '') e_HeI_w4471 = rNULL__
if(ar__1(321:325) .EQ. '') FeIII_w4658 = rNULL__
if(ar__1(327:331) .EQ. '') e_FeIII_w4658 = rNULL__
if(ar__1(333:337) .EQ. '') HeII_w4686 = rNULL__
if(ar__1(339:343) .EQ. '') e_HeII_w4686 = rNULL__
if(ar__1(357:361) .EQ. '') HeI_w4921 = rNULL__
if(ar__1(363:367) .EQ. '') e_HeI_w4921 = rNULL__
if(ar__1(405:409) .EQ. '') NII_w5755 = rNULL__
if(ar__1(411:415) .EQ. '') e_NII_w5755 = rNULL__
if(ar__1(429:433) .EQ. '') OI_w6300 = rNULL__
if(ar__1(435:439) .EQ. '') e_OI_w6300 = rNULL__
if(ar__1(441:445) .EQ. '') SIII_w6312 = rNULL__
if(ar__1(447:451) .EQ. '') e_SIII_w6312 = rNULL__
if(ar__1(453:457) .EQ. '') OI_w6363 = rNULL__
if(ar__1(459:463) .EQ. '') e_OI_w6363 = rNULL__
if(ar__1(537:541) .EQ. '') HeI_w7065 = rNULL__
if(ar__1(543:547) .EQ. '') e_HeI_w7065 = rNULL__
if(ar__1(549:553) .EQ. '') ArIII_w7135 = rNULL__
if(ar__1(555:559) .EQ. '') e_ArIII_w7135 = rNULL__
if(ar__1(573:577) .EQ. '') OII_w7330 = rNULL__
if(ar__1(579:583) .EQ. '') e_OII_w7330 = rNULL__
if(ar__1(585:589) .EQ. '') ArIII_w7751 = rNULL__
if(ar__1(591:595) .EQ. '') e_ArIII_w7751 = rNULL__
if(ar__1(597:601) .EQ. '') P13_w8665 = rNULL__
if(ar__1(603:607) .EQ. '') e_P13_w8665 = rNULL__
if(ar__1(609:613) .EQ. '') P12_w8750 = rNULL__
if(ar__1(615:619) .EQ. '') e_P12_w8750 = rNULL__
if(ar__1(621:625) .EQ. '') P11_w8862 = rNULL__
if(ar__1(627:631) .EQ. '') e_P11_w8862 = rNULL__
if(ar__1(633:637) .EQ. '') P10_w9015 = rNULL__
if(ar__1(639:643) .EQ. '') e_P10_w9015 = rNULL__
if(ar__1(645:649) .EQ. '') SIII_w9069 = rNULL__
if(ar__1(651:655) .EQ. '') e_SIII_w9069 = rNULL__
if(ar__1(657:661) .EQ. '') P9_w9229 = rNULL__
if(ar__1(663:667) .EQ. '') e_P9_w9229 = rNULL__
if(ar__1(681:685) .EQ. '') P8_w9546 = rNULL__
if(ar__1(687:691) .EQ. '') e_P8_w9546 = rNULL__
c ..............Just test output...........
write(6,2)
+ Reg_1,H14_w3721,e_H14_w3721,OII_w3727,e_OII_w3727,H13_w3734,
+ e_H13_w3734,H12_w3750,e_H12_w3750,H11_w3770,e_H11_w3770,
+ H10_w3797,e_H10_w3797,HeI_w3819,e_HeI_w3819,H9_w3835,
+ e_H9_w3835,NeIII_w3868,e_NeIII_w3868,HeI_w3888,e_HeI_w3888,
+ H8_w3889,e_H8_w3889,HeI_w3964,e_HeI_w3964,NeIII_w3967,
+ e_NeIII_w3967,H7_w3970,e_H7_w3970,NeIII_w4011,e_NeIII_w4011,
+ HeI_w4026,e_HeI_w4026,SII_w4068,e_SII_w4068,SII_w4076,
+ e_SII_w4076,Hd_w4101,e_Hd_w4101,HeI_w4120,e_HeI_w4120,
+ HeI_w4143,e_HeI_w4143,Hg_w4340,e_Hg_w4340,OIII_w4363,
+ e_OIII_w4363,HeI_w4387,e_HeI_w4387,HeI_w4471,e_HeI_w4471,
+ FeIII_w4658,e_FeIII_w4658,HeII_w4686,e_HeII_w4686,Hb_w4861,
+ e_Hb_w4861,HeI_w4921,e_HeI_w4921,OIII_w4959,e_OIII_w4959,
+ OIII_w5007,e_OIII_w5007,HeI_w5015,e_HeI_w5015,NII_w5755,
+ e_NII_w5755,HeI_w5876,e_HeI_w5876,OI_w6300,e_OI_w6300,
+ SIII_w6312,e_SIII_w6312,OI_w6363,e_OI_w6363,NII_w6548,
+ e_NII_w6548,Ha_w6563,e_Ha_w6563,NII_w6584,e_NII_w6584,
+ HeI_w6678,e_HeI_w6678,SII_w6717,e_SII_w6717,SII_w6731,
+ e_SII_w6731,HeI_w7065,e_HeI_w7065,ArIII_w7135,e_ArIII_w7135,
+ OII_w7320,e_OII_w7320,OII_w7330,e_OII_w7330,ArIII_w7751,
+ e_ArIII_w7751,P13_w8665,e_P13_w8665,P12_w8750,e_P12_w8750,
+ P11_w8862,e_P11_w8862,P10_w9015,e_P10_w9015,SIII_w9069,
+ e_SIII_w9069,P9_w9229,e_P9_w9229,SIII_w9532,e_SIII_w9532,
+ P8_w9546,e_P8_w9546,E_B_V,e_E_B_V,FHb,e_FHb
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,Te_SII_,e_Te_SII_,Te_NII_,e_Te_NII_,Te_OII_,e_Te_OII_,
+ Te_SIII_,e_Te_SIII_,Te_OIII_,e_Te_OIII_,nem,e_nem,Te_low,
+ e_Te_low,Te_int,e_Te_int,Te_high,e_Te_high,neu,e_neu,O__H_,
+ e_O__H_,O_2_H_,e_O_2_H_,log_O_H,e_log_O_H,N__H_,e_N__H_,ICF_N,
+ e_ICF_N,log_N_H,e_log_N_H,log_N_O,e_log_N_O,S__H_,e_S__H_,
+ S2_H_,e_S2_H_,ICF_S,e_ICF_S,log_S_H,e_log_S_H,log_S_O,
+ e_log_S_O,Ne_2_H_,e_Ne_2_H_,ICF_Ne,e_ICF_Ne,log_Ne_H,
+ e_log_Ne_H,log_Ne_O,e_log_Ne_O,Ar_2_H_,e_Ar_2_H_,ICF_Ar,
+ e_ICF_Ar,log_Ar_H,e_log_Ar_H,log_Ar_O,e_log_Ar_O
if(ar__2(21:25) .EQ. '') Te_SII_ = iNULL__
if(ar__2(27:29) .EQ. '') e_Te_SII_ = iNULL__
if(ar__2(31:34) .EQ. '') Te_NII_ = iNULL__
if(ar__2(36:39) .EQ. '') e_Te_NII_ = iNULL__
if(ar__2(41:45) .EQ. '') Te_OII_ = iNULL__
if(ar__2(47:49) .EQ. '') e_Te_OII_ = iNULL__
if(ar__2(51:54) .EQ. '') Te_SIII_ = iNULL__
if(ar__2(56:58) .EQ. '') e_Te_SIII_ = iNULL__
if(ar__2(60:63) .EQ. '') Te_OIII_ = iNULL__
if(ar__2(65:67) .EQ. '') e_Te_OIII_ = iNULL__
if(ar__2(255:260) .EQ. '') Ne_2_H_ = rNULL__
if(ar__2(262:267) .EQ. '') e_Ne_2_H_ = rNULL__
if(ar__2(281:285) .EQ. '') log_Ne_H = rNULL__
if(ar__2(287:291) .EQ. '') e_log_Ne_H = rNULL__
if(ar__2(293:298) .EQ. '') log_Ne_O = rNULL__
if(ar__2(300:304) .EQ. '') e_log_Ne_O = rNULL__
c ..............Just test output...........
write(6,3)
+ Reg_2,Te_SII_,e_Te_SII_,Te_NII_,e_Te_NII_,Te_OII_,e_Te_OII_,
+ Te_SIII_,e_Te_SIII_,Te_OIII_,e_Te_OIII_,nem,e_nem,Te_low,
+ e_Te_low,Te_int,e_Te_int,Te_high,e_Te_high,neu,e_neu,O__H_,
+ e_O__H_,O_2_H_,e_O_2_H_,log_O_H,e_log_O_H,N__H_,e_N__H_,ICF_N,
+ e_ICF_N,log_N_H,e_log_N_H,log_N_O,e_log_N_O,S__H_,e_S__H_,
+ S2_H_,e_S2_H_,ICF_S,e_ICF_S,log_S_H,e_log_S_H,log_S_O,
+ e_log_S_O,Ne_2_H_,e_Ne_2_H_,ICF_Ne,e_ICF_Ne,log_Ne_H,
+ e_log_Ne_H,log_Ne_O,e_log_Ne_O,Ar_2_H_,e_Ar_2_H_,ICF_Ar,
+ e_ICF_Ar,log_Ar_H,e_log_Ar_H,log_Ar_O,e_log_Ar_O
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_,log_O_H_1,e_log_O_H_1,log_N_O_1,e_log_N_O_1,
+ log_S_O_1,e_log_S_O_1,log_Ar_O_1,e_log_Ar_O_1,log_Ne_O_1,
+ e_log_Ne_O_1
if(ar__3(63:67) .EQ. '') log_Ne_O_1 = rNULL__
if(ar__3(69:72) .EQ. '') e_log_Ne_O_1 = rNULL__
c ..............Just test output...........
write(6,4)
+ v_BSC2015_,log_O_H_1,e_log_O_H_1,log_N_O_1,e_log_N_O_1,
+ log_S_O_1,e_log_S_O_1,log_Ar_O_1,e_log_Ar_O_1,log_Ne_O_1,
+ e_log_Ne_O_1
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_,log_O_H_1,e_log_O_H_1,log_N_O_1,e_log_N_O_1,
+ log_S_O_1,e_log_S_O_1,log_Ar_O_1,e_log_Ar_O_1,log_Ne_O_1,
+ e_log_Ne_O_1
if(ar__4(64:68) .EQ. '') log_Ne_O_1 = rNULL__
if(ar__4(70:73) .EQ. '') e_log_Ne_O_1 = rNULL__
c ..............Just test output...........
write(6,5)
+ v_CPB2015_,log_O_H_1,e_log_O_H_1,log_N_O_1,e_log_N_O_1,
+ log_S_O_1,e_log_S_O_1,log_Ar_O_1,e_log_Ar_O_1,log_Ne_O_1,
+ e_log_Ne_O_1
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_,log_O_H_1,e_log_O_H_1,log_N_O_1,e_log_N_O_1,
+ log_S_O_1,e_log_S_O_1,log_Ar_O_1,e_log_Ar_O_1,log_Ne_O_1,
+ e_log_Ne_O_1
if(ar__5(53:57) .EQ. '') log_Ar_O_1 = rNULL__
if(ar__5(59:62) .EQ. '') e_log_Ar_O_1 = rNULL__
if(ar__5(64:68) .EQ. '') log_Ne_O_1 = rNULL__
if(ar__5(70:73) .EQ. '') e_log_Ne_O_1 = rNULL__
c ..............Just test output...........
write(6,6)
+ v_CPB2016_,log_O_H_1,e_log_O_H_1,log_N_O_1,e_log_N_O_1,
+ log_S_O_1,e_log_S_O_1,log_Ar_O_1,e_log_Ar_O_1,log_Ne_O_1,
+ e_log_Ne_O_1
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end