FORTRAN Generation
(/./ftp/cats/J/ApJ/655/814)

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