Conversion of standardized ReadMe file for
file /./ftp/cats/J/ApJ/655/814 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 2024-Mar-28 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/655/814 HII regions in IC 342/Maffei group (Fingerhut+, 2007) *================================================================================ *The extinction and distance of Maffei 2 and a new view of the *IC 342/Maffei group. * Fingerhut R.L., Lee H., Mccall M.L., Richer M.G. * <Astrophys. J., 655, 814-830 (2007)> * =2007ApJ...655..814F C============================================================================= C Internal variables integer*4 i__ c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table2.dat' ! Properties of HII regions in Maffei 2 and IC 342 integer*4 nr__ parameter (nr__=8) ! Number of records character*124 ar__ ! Full-size record character*8 Gal (nr__) ! Galaxy name character*10 v_FLM2007_ (nr__) ! HII region name (+xxxx+yyyy) (1) integer*4 r (nr__) ! (arcsec) Deprojected galactocentric radius (2) real*4 EWHb (nr__) ! (0.1nm) ?=- H{beta} equivalent width real*4 Ha_Hb (nr__) ! Observed H{alpha}-to-H{beta} flux ratio (3) real*4 e_Ha_Hb (nr__) ! ? rms uncertainty on Ha/Hb real*4 v_O_H_ (nr__) ! ([-]) Oxygen abundance log(O/H) derived from the * forbidden oxygen lines character*1 n__O_H_ (nr__) ! [a] No O line data (4) integer*4 T (nr__) ! (K) Equilibrium temperature derived from [O/H] real*4 Ha_Hb0 (nr__) ! Intrinsic Balmer flux ratio (5) real*4 tau (nr__) ! Total optical depth at 1um derived from * the observed and intrinsic Balmer decrements real*4 e_tau (nr__) ! ? rms uncertainty on tau real*4 NHI (nr__) ! (10+20cm-2) Annular averaged column density of * neutral hydrogen within the annulus at r real*4 e_NHI (nr__) ! (10+20cm-2) ? rms uncertainty on HI real*4 X (nr__) ! (10+20cm-2/km/s) Conversion factor for CO intensity to H2 * column density (X=N(H2)/I(CO)) * (in 10^20^cm^-2^/(km/s)) real*4 NH2 (nr__) ! (10+20cm-2) Annular averaged column density of * molecular hydrogen within the annulus at r real*4 e_NH2 (nr__) ! (10+20cm-2) ? rms uncertainty on NH2 real*4 nH (nr__) ! (10+20cm-2) Half of the column density of hydrogen * in both atomic and molecular form at r real*4 e_nH (nr__) ! (10+20cm-2) ? rms uncertainty on nH real*4 tauloc (nr__) ! Extragalactic optical depth at 1um (6) real*4 e_tauloc (nr__) ! ? rms uncertainty on tauloc character*4 Ref (nr__) ! References (7) *Note (1): Eastward and northward offset of the H II region from the galaxy * nucleus as projected on the sky, in arcseconds. * Identified as [FLM2007] Gal (+xxxx+yyyy) in Simbad. *Note (2): Deprojected galactocentric radius of the HII region using i and * {phi} from Table 3. *Note (3): corrected for underlying stellar absorption but not for Galactic * reddening, assuming an equivalent width for absorption of 1.9{AA} * (see McCall et al., 1985ApJS...57....1M). *Note (4): a = absence of oxygen line data, resulting in the use of eq. (4) * or eq. (5) to determine the temperature, and hence log(O/H). *Note (5): Intrinsic Balmer flux ratio based on the HII region's equilibrium * temperature and an electron density of 100cm^-3^ *Note (6): Extragalactic optical depth at 1um of the HII region, i.e., the * difference between the total optical depth and the Galactic optical * depth. *Note (7): References as follows: * BK82 = Blair et al., 1982ApJ...254...50B * FL07 = This study * MR85 = McCall et al., 1985ApJS...57....1M * KA81 = Kwitter & Aller, 1981MNRAS.195..939K * RP82 = Rayo et al., 1982ApJ...255....1R * SK96 = Skillman et al., 1996ApJ...462..147S * SS91 = Shields & Skillman, 1991ApJ...371...82S * VP88 = Vilchez et al., 1988MNRAS.235..633V c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table4.dat' ! Properties of HII regions sample integer*4 nr__1 parameter (nr__1=73) ! Number of records character*124 ar__1 ! Full-size record character*8 Gal_1 (nr__1) ! Galaxy name character*10 v_FLM2007__1(nr__1) ! HII region name (+xxxx+yyyy) (1) integer*4 r_1 (nr__1) ! (arcsec) Deprojected galactocentric radius (2) real*4 EWHb_1 (nr__1) ! (0.1nm) ?=- H{beta} equivalent width real*4 Ha_Hb_1 (nr__1) ! Observed H{alpha}-to-H{beta} flux ratio (3) real*4 e_Ha_Hb_1 (nr__1) ! ? rms uncertainty on Ha/Hb real*4 v_O_H__1 (nr__1) ! ([-]) Oxygen abundance log(O/H) derived from the * forbidden oxygen lines character*1 n__O_H__1 (nr__1) ! [a] No O line data (4) integer*4 T_1 (nr__1) ! (K) Equilibrium temperature derived from [O/H] real*4 Ha_Hb0_1 (nr__1) ! Intrinsic Balmer flux ratio (5) real*4 tau_1 (nr__1) ! Total optical depth at 1um derived from * the observed and intrinsic Balmer decrements real*4 e_tau_1 (nr__1) ! ? rms uncertainty on tau real*4 NHI_1 (nr__1) ! (10+20cm-2) Annular averaged column density of * neutral hydrogen within the annulus at r real*4 e_NHI_1 (nr__1) ! (10+20cm-2) ? rms uncertainty on HI real*4 X_1 (nr__1) ! (10+20cm-2/km/s) Conversion factor for CO intensity to H2 * column density (X=N(H2)/I(CO)) * (in 10^20^cm^-2^/(km/s)) real*4 NH2_1 (nr__1) ! (10+20cm-2) Annular averaged column density of * molecular hydrogen within the annulus at r real*4 e_NH2_1 (nr__1) ! (10+20cm-2) ? rms uncertainty on NH2 real*4 nH_1 (nr__1) ! (10+20cm-2) Half of the column density of hydrogen * in both atomic and molecular form at r real*4 e_nH_1 (nr__1) ! (10+20cm-2) ? rms uncertainty on nH real*4 tauloc_1 (nr__1) ! Extragalactic optical depth at 1um (6) real*4 e_tauloc_1 (nr__1) ! ? rms uncertainty on tauloc character*4 Ref_1 (nr__1) ! References (7) *Note (1): Eastward and northward offset of the H II region from the galaxy * nucleus as projected on the sky, in arcseconds. * Identified as [FLM2007] Gal (+xxxx+yyyy) in Simbad. *Note (2): Deprojected galactocentric radius of the HII region using i and * {phi} from Table 3. *Note (3): corrected for underlying stellar absorption but not for Galactic * reddening, assuming an equivalent width for absorption of 1.9{AA} * (see McCall et al., 1985ApJS...57....1M). *Note (4): a = absence of oxygen line data, resulting in the use of eq. (4) * or eq. (5) to determine the temperature, and hence log(O/H). *Note (5): Intrinsic Balmer flux ratio based on the HII region's equilibrium * temperature and an electron density of 100cm^-3^ *Note (6): Extragalactic optical depth at 1um of the HII region, i.e., the * difference between the total optical depth and the Galactic optical * depth. *Note (7): References as follows: * BK82 = Blair et al., 1982ApJ...254...50B * FL07 = This study * MR85 = McCall et al., 1985ApJS...57....1M * KA81 = Kwitter & Aller, 1981MNRAS.195..939K * RP82 = Rayo et al., 1982ApJ...255....1R * SK96 = Skillman et al., 1996ApJ...462..147S * SS91 = Shields & Skillman, 1991ApJ...371...82S * VP88 = Vilchez et al., 1988MNRAS.235..633V c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - C Declarations for 'table3.dat' ! Properties of the spiral galaxies hosting HII regions integer*4 nr__2 parameter (nr__2=12) ! Number of records character*75 ar__2 ! Full-size record C J2000 position composed of: RAh RAm RAs DE- DEd DEm DEs real*8 RAdeg (nr__2) ! (deg) Right Ascension J2000 real*8 DEdeg (nr__2) ! (deg) Declination J2000 C ---------------------------------- ! (position vector(s) in degrees) character*8 Gal_2 (nr__2) ! Galaxy integer*4 TType (nr__2) ! Morphological stage from RC3 (VII/155) integer*4 RAh (nr__2) ! (h) Right ascension (J2000) integer*4 RAm (nr__2) ! (min) Right ascension (J2000) real*4 RAs (nr__2) ! (s) Right ascension (J2000) character*1 DE_ (nr__2) ! Declination sign (J2000) integer*4 DEd (nr__2) ! (deg) Declination (J2000) integer*4 DEm (nr__2) ! (arcmin) Declination (J2000) integer*4 DEs (nr__2) ! (arcsec) Declination (J2000) character*1 r_2 (nr__2) ! [a] a: position from Buta & McCall * (1999ApJS..124...33B), else RC3 (VII/155) integer*4 HRV (nr__2) ! (km/s) Heliocentric radial velocity from the HI * reference character*1 r_HRV (nr__2) ! [b] b: HRV from Newton (1980MNRAS.191..169N), * else RC3 (Cat. VII/155) integer*4 i (nr__2) ! (deg) Inclination angle from the HI reference character*1 r_i (nr__2) ! [b] b: i from Newton (1980MNRAS.191..169N), * else RC3 (Cat. VII/155) integer*4 phi (nr__2) ! (deg) Position angle from the HI reference character*1 r_phi (nr__2) ! [b] b: phi from Newton (1980MNRAS.191..169N), * else RC3 (Cat. VII/155) real*4 r0 (nr__2) ! (arcsec) Disk scale length (1) real*4 E_B_V (nr__2) ! (mag) ?=- Galactic reddening from Schlegel et al. * (1998ApJ...500..525S) real*4 tau_2 (nr__2) ! ?=- Optical depth of Galactic dust at 1um * calculated from E(B-V) character*4 HIref (nr__2) ! Source of the HI map of the galaxy (2) character*4 COref (nr__2) ! Source of the CO distribution of the galaxy (2) *Note (1): The radius at which the surface brightness in B drops to 1/e of * its central value from Skillman et al. (1996ApJ...462..147S) (NGC 4501 * and NGC 4654), Buta & McCall (1999ApJS..124...33B) (Maffei 2 and IC 342), * and Zaritsky et al. (1994ApJ...420...87Z) (for all other galaxies). *Note (2): References as follows: * BG81 = Bosma et al., 1981A&A....93..106B * HT96 = Hurt et al., 1996ApJ...466..135H * KS91 = Kenney et al., 1991ApJ...366..432K * KY88 = Kenney & Young, 1988ApJS...66..261K * MW04 = Mason & Wilson, 2004ApJ...612..860M * N80a = Newton, 1980MNRAS.191..169N * N80b = Newton, 1980MNRAS.191..615N * TA91 = Tilanus & Allen, 1991A&A...244....8T * W88 = Warmels, 1988A&AS...72..427W * WV86 = Wevers et al., 1986A&AS...66..505W * YX95 = Young et al., 1995ApJS...98..219Y C============================================================================= C Loading file 'table2.dat' ! Properties of HII regions in Maffei 2 and IC 342 C Format for file interpretation 1 format( + A8,1X,A10,1X,I4,1X,F5.1,1X,F5.2,1X,F4.2,1X,F5.2,A1,1X,I5,1X, + F4.2,1X,F5.2,1X,F4.2,1X,F5.2,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,A4) C Effective file loading open(unit=1,status='old',file= +'table2.dat') write(6,*) '....Loading file: table2.dat' do i__=1,8 read(1,'(A124)')ar__ read(ar__,1) + Gal(i__),v_FLM2007_(i__),r(i__),EWHb(i__),Ha_Hb(i__), + e_Ha_Hb(i__),v_O_H_(i__),n__O_H_(i__),T(i__),Ha_Hb0(i__), + tau(i__),e_tau(i__),NHI(i__),e_NHI(i__),X(i__),NH2(i__), + e_NH2(i__),nH(i__),e_nH(i__),tauloc(i__),e_tauloc(i__), + Ref(i__) if (idig(ar__(26:30)).EQ.0) EWHb(i__) = rNULL__ if(ar__(38:41) .EQ. '') e_Ha_Hb(i__) = rNULL__ if(ar__(67:70) .EQ. '') e_tau(i__) = rNULL__ if(ar__(78:81) .EQ. '') e_NHI(i__) = rNULL__ if(ar__(94:97) .EQ. '') e_NH2(i__) = rNULL__ if(ar__(105:108) .EQ. '') e_nH(i__) = rNULL__ if(ar__(116:119) .EQ. '') e_tauloc(i__) = rNULL__ c ..............Just test output........... write(6,1) + Gal(i__),v_FLM2007_(i__),r(i__),EWHb(i__),Ha_Hb(i__), + e_Ha_Hb(i__),v_O_H_(i__),n__O_H_(i__),T(i__),Ha_Hb0(i__), + tau(i__),e_tau(i__),NHI(i__),e_NHI(i__),X(i__),NH2(i__), + e_NH2(i__),nH(i__),e_nH(i__),tauloc(i__),e_tauloc(i__), + Ref(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table4.dat' ! Properties of HII regions sample C Format for file interpretation 2 format( + A8,1X,A10,1X,I4,1X,F5.1,1X,F5.2,1X,F4.2,1X,F5.2,A1,1X,I5,1X, + F4.2,1X,F5.2,1X,F4.2,1X,F5.2,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,A4) C Effective file loading open(unit=1,status='old',file= +'table4.dat') write(6,*) '....Loading file: table4.dat' do i__=1,73 read(1,'(A124)')ar__1 read(ar__1,2) + Gal_1(i__),v_FLM2007__1(i__),r_1(i__),EWHb_1(i__), + Ha_Hb_1(i__),e_Ha_Hb_1(i__),v_O_H__1(i__),n__O_H__1(i__), + T_1(i__),Ha_Hb0_1(i__),tau_1(i__),e_tau_1(i__),NHI_1(i__), + e_NHI_1(i__),X_1(i__),NH2_1(i__),e_NH2_1(i__),nH_1(i__), + e_nH_1(i__),tauloc_1(i__),e_tauloc_1(i__),Ref_1(i__) if (idig(ar__1(26:30)).EQ.0) EWHb_1(i__) = rNULL__ if(ar__1(38:41) .EQ. '') e_Ha_Hb_1(i__) = rNULL__ if(ar__1(67:70) .EQ. '') e_tau_1(i__) = rNULL__ if(ar__1(78:81) .EQ. '') e_NHI_1(i__) = rNULL__ if(ar__1(94:97) .EQ. '') e_NH2_1(i__) = rNULL__ if(ar__1(105:108) .EQ. '') e_nH_1(i__) = rNULL__ if(ar__1(116:119) .EQ. '') e_tauloc_1(i__) = rNULL__ c ..............Just test output........... write(6,2) + Gal_1(i__),v_FLM2007__1(i__),r_1(i__),EWHb_1(i__), + Ha_Hb_1(i__),e_Ha_Hb_1(i__),v_O_H__1(i__),n__O_H__1(i__), + T_1(i__),Ha_Hb0_1(i__),tau_1(i__),e_tau_1(i__),NHI_1(i__), + e_NHI_1(i__),X_1(i__),NH2_1(i__),e_NH2_1(i__),nH_1(i__), + e_nH_1(i__),tauloc_1(i__),e_tauloc_1(i__),Ref_1(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= C Loading file 'table3.dat' ! Properties of the spiral galaxies hosting * HII regions C Format for file interpretation 3 format( + A8,1X,I1,1X,I2,1X,I2,1X,F4.1,1X,A1,I2,1X,I2,1X,I2,A1,1X,I4,A1, + 1X,I2,A1,1X,I3,A1,1X,F5.1,1X,F5.3,1X,F5.3,1X,A4,1X,A4) C Effective file loading open(unit=1,status='old',file= +'table3.dat') write(6,*) '....Loading file: table3.dat' do i__=1,12 read(1,'(A75)')ar__2 read(ar__2,3) + Gal_2(i__),TType(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__), + DEd(i__),DEm(i__),DEs(i__),r_2(i__),HRV(i__),r_HRV(i__), + i(i__),r_i(i__),phi(i__),r_phi(i__),r0(i__),E_B_V(i__), + tau_2(i__),HIref(i__),COref(i__) if (idig(ar__2(55:59)).EQ.0) E_B_V(i__) = rNULL__ if (idig(ar__2(61:65)).EQ.0) tau_2(i__) = rNULL__ 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. if(DE_(i__).EQ.'-'.AND.DEdeg(i__).GE.0) DEdeg(i__)=-DEdeg(i__) c ..............Just test output........... write(6,3) + Gal_2(i__),TType(i__),RAh(i__),RAm(i__),RAs(i__),DE_(i__), + DEd(i__),DEm(i__),DEs(i__),r_2(i__),HRV(i__),r_HRV(i__), + i(i__),r_i(i__),phi(i__),r_phi(i__),r0(i__),E_B_V(i__), + tau_2(i__),HIref(i__),COref(i__) write(6,'(6H Pos: 2F8.4)') RAdeg(i__),DEdeg(i__) c .......End.of.Just test output........... end do close(1) C============================================================================= stop end C Locate position of first digit in string; or return 0 integer function idig(c) character*(*) c character*1 c1 integer lc,i lc=len(c) idig=0 do i=1,lc if(c(i:i).ne.' ') go to 1 end do 1 if(i.gt.lc) return c1=c(i:i) if(c1.eq.'.'.or.c1.eq.'-'.or.c1.eq.'+') i=i+1 if(i.gt.lc) return c1=c(i:i) if(c1.ge.'0'.and.c1.le.'9') idig=i return end