FORTRAN Generation
(/./ftp/cats/J/A_A/530/A18)

Conversion of standardized ReadMe file for file /./ftp/cats/J/A_A/530/A18 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 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/A+A/530/A18  Nebular NII lines effective recombination coeff.    (Fang+, 2013)
*================================================================================
*New effective recombination coefficients for nebular NII lines (Corrigedum).
*    Fang X., Storey P.J., Liu X.-W.
*   <Astron. Astrophys. 530, A18 (2011) and 550, C2 (2013)>
*   =2011A&A...530A..18F
*   +2013A&A...550C...2F
C=============================================================================

C  Internal variables

      integer*4 i__

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table3.dat'	! Case B effective recombination coefficients
                              for electron density Ne=10^2^cm^-3^

      integer*4 nr__
      parameter (nr__=217)	! Number of records
      character*122 ar__   	! Full-size record

      character*21  Tr          ! General transition term (multiplet)
      character*31  Trans       ! Detailed transition
      real*8        lambda      ! (Angstrom) Transition wavelength {lambda}
      character*1   n_lambda    ! [*] experimentally known energies (1)
      real*4        ERC1        ! (10-15cm3/s) Effective recombination coefficient
*                                       {alpha}_eff_({lambda}) for Te=125K (2)
      real*4        ERC2        ! (10-15cm3/s) Effective recombination coefficient
*                                       {alpha}_eff_({lambda}) for Te=500K (2)
      real*4        ERC3        ! (10-15cm3/s) Effective recombination coefficient
*                                       {alpha}_eff_({lambda}) for Te=1000K (2)
      real*4        ERC4        ! (10-15cm3/s) Effective recombination coefficient
*                                       {alpha}_eff_({lambda}) for Te=5000K (2)
      real*4        ERC5        ! (10-15cm3/s) Effective recombination coefficient
*                                       {alpha}_eff_({lambda}) for Te=10000K (2)
      real*4        ERC6        ! (10-15cm3/s) Effective recombination coefficient
*                                       {alpha}_eff_({lambda}) for Te=15000K (2)
      real*4        ERC7        ! (10-15cm3/s) Effective recombination coefficient
*                                       {alpha}_eff_({lambda}) for Te=20000K (2)
*Note (1): * denotes this transition wavelength was derived from the
*           experimentally known energies of the upper and lower states.
*Note (2): the Case B effective recombination coefficient
*          {alpha}_eff_({lambda}) is defined by emissivity e({lambda}):
*          e({lambda}) = N_e N_+ {alpha}_eff_({lambda}) hc/{lambda}
*          [e({lambda}) in erg/s/cm^-3^ = 100W/m^3^]

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table4.dat'	! Case B effective recombination coefficients
                              for electron density Ne=10^3^cm^-3^

      integer*4 nr__1
      parameter (nr__1=217)	! Number of records
      character*122 ar__1  	! Full-size record

      character*21  Tr_1        ! General transition term (multiplet)
      character*31  Trans_1     ! Detailed transition
      real*8        lambda_1    ! (Angstrom) Transition wavelength {lambda}
      character*1   n_lambda_1  ! [*] experimentally known energies (1)
      real*4        ERC1_1      ! (10-15cm3/s) Effective recombination coefficient
*                                       {alpha}_eff_({lambda}) for Te=125K (2)
      real*4        ERC2_1      ! (10-15cm3/s) Effective recombination coefficient
*                                       {alpha}_eff_({lambda}) for Te=500K (2)
      real*4        ERC3_1      ! (10-15cm3/s) Effective recombination coefficient
*                                       {alpha}_eff_({lambda}) for Te=1000K (2)
      real*4        ERC4_1      ! (10-15cm3/s) Effective recombination coefficient
*                                       {alpha}_eff_({lambda}) for Te=5000K (2)
      real*4        ERC5_1      ! (10-15cm3/s) Effective recombination coefficient
*                                       {alpha}_eff_({lambda}) for Te=10000K (2)
      real*4        ERC6_1      ! (10-15cm3/s) Effective recombination coefficient
*                                       {alpha}_eff_({lambda}) for Te=15000K (2)
      real*4        ERC7_1      ! (10-15cm3/s) Effective recombination coefficient
*                                       {alpha}_eff_({lambda}) for Te=20000K (2)
*Note (1): * denotes this transition wavelength was derived from the
*           experimentally known energies of the upper and lower states.
*Note (2): the Case B effective recombination coefficient
*          {alpha}_eff_({lambda}) is defined by emissivity e({lambda}):
*          e({lambda}) = N_e N_+ {alpha}_eff_({lambda}) hc/{lambda}
*          [e({lambda}) in erg/s/cm^-3^ = 100W/m^3^]

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table5.dat'	! Case B effective recombination coefficients
                              for electron density Ne=10^4^cm^-3^

      integer*4 nr__2
      parameter (nr__2=217)	! Number of records
      character*122 ar__2  	! Full-size record

      character*21  Tr_2        ! General transition term (multiplet)
      character*31  Trans_2     ! Detailed transition
      real*8        lambda_2    ! (Angstrom) Transition wavelength {lambda}
      character*1   n_lambda_2  ! [*] experimentally known energies (1)
      real*4        ERC1_2      ! (10-15cm3/s) Effective recombination coefficient
*                                       {alpha}_eff_({lambda}) for Te=125K (2)
      real*4        ERC2_2      ! (10-15cm3/s) Effective recombination coefficient
*                                       {alpha}_eff_({lambda}) for Te=500K (2)
      real*4        ERC3_2      ! (10-15cm3/s) Effective recombination coefficient
*                                       {alpha}_eff_({lambda}) for Te=1000K (2)
      real*4        ERC4_2      ! (10-15cm3/s) Effective recombination coefficient
*                                       {alpha}_eff_({lambda}) for Te=5000K (2)
      real*4        ERC5_2      ! (10-15cm3/s) Effective recombination coefficient
*                                       {alpha}_eff_({lambda}) for Te=10000K (2)
      real*4        ERC6_2      ! (10-15cm3/s) Effective recombination coefficient
*                                       {alpha}_eff_({lambda}) for Te=15000K (2)
      real*4        ERC7_2      ! (10-15cm3/s) Effective recombination coefficient
*                                       {alpha}_eff_({lambda}) for Te=20000K (2)
*Note (1): * denotes this transition wavelength was derived from the
*           experimentally known energies of the upper and lower states.
*Note (2): the Case B effective recombination coefficient
*          {alpha}_eff_({lambda}) is defined by emissivity e({lambda}):
*          e({lambda}) = N_e N_+ {alpha}_eff_({lambda}) hc/{lambda}
*          [e({lambda}) in erg/s/cm^-3^ = 100W/m^3^]

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table6.dat'	! Case B effective recombination coefficients
                              for electron density Ne=10^5^cm^-3^

      integer*4 nr__3
      parameter (nr__3=217)	! Number of records
      character*122 ar__3  	! Full-size record

      character*21  Tr_3        ! General transition term (multiplet)
      character*31  Trans_3     ! Detailed transition
      real*8        lambda_3    ! (Angstrom) Transition wavelength {lambda}
      character*1   n_lambda_3  ! [*] experimentally known energies (1)
      real*4        ERC1_3      ! (10-15cm3/s) Effective recombination coefficient
*                                       {alpha}_eff_({lambda}) for Te=125K (2)
      real*4        ERC2_3      ! (10-15cm3/s) Effective recombination coefficient
*                                       {alpha}_eff_({lambda}) for Te=500K (2)
      real*4        ERC3_3      ! (10-15cm3/s) Effective recombination coefficient
*                                       {alpha}_eff_({lambda}) for Te=1000K (2)
      real*4        ERC4_3      ! (10-15cm3/s) Effective recombination coefficient
*                                       {alpha}_eff_({lambda}) for Te=5000K (2)
      real*4        ERC5_3      ! (10-15cm3/s) Effective recombination coefficient
*                                       {alpha}_eff_({lambda}) for Te=10000K (2)
      real*4        ERC6_3      ! (10-15cm3/s) Effective recombination coefficient
*                                       {alpha}_eff_({lambda}) for Te=15000K (2)
      real*4        ERC7_3      ! (10-15cm3/s) Effective recombination coefficient
*                                       {alpha}_eff_({lambda}) for Te=20000K (2)
*Note (1): * denotes this transition wavelength was derived from the
*           experimentally known energies of the upper and lower states.
*Note (2): the Case B effective recombination coefficient
*          {alpha}_eff_({lambda}) is defined by emissivity e({lambda}):
*          e({lambda}) = N_e N_+ {alpha}_eff_({lambda}) hc/{lambda}
*          [e({lambda}) in erg/s/cm^-3^ = 100W/m^3^]

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table7.dat'	! Fit parameters and the average and maximum
                              fitting errors for 125<=Te<=20000K and
                              Ne=10^2^cm^-3^ for Case B recombination

      integer*4 nr__4
      parameter (nr__4=55)	! Number of records
      character*129 ar__4  	! Full-size record

      character*27  Trans_4     ! Transition
      character*3   Mult        ! Multiplet
      real*8        lambda_4    ! (Angstrom) Wavelength of transition {lambda}
      real*8        a           ! a coefficient (1)
      real*8        b           ! b coefficient (1)
      real*8        c           ! c coefficient (1)
      real*8        d           ! d coefficient (1)
      real*8        e           ! e coefficient (1)
      real*8        f           ! f coefficient (1)
      real*8        g           ! g coefficient (1)
      real*8        h           ! h coefficient (1)
      real*4        del         ! (%) Average deviation of the fit (delta)
      real*4        Del_1       ! (%) Maximum deviation of the fit {Delta}
*Note (1): The expression of the effective recombination coefficient is
*   log({alpha}_eff_)+15 = a+b*t+c*t^2^+(d+e*t+f*t^2^)log(t)+g*(log(t))^2^+h/t,
*      where t is the reduced electronic temperature Te/10^4^K.
*   table7.dat, table8.dat, table9.dat and table10.dat are fit results for
*      Case B recombination;
*   table11.dat, table12.dat, table13.dat and table14.dat are fit results for
*      Case A recombination.

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table8.dat'	! Fit parameters and the average and maximum
                              fitting errors for 125<=Te<=20000K and
                              Ne=10^3^cm^-3^ for Case B recombination

      integer*4 nr__5
      parameter (nr__5=55)	! Number of records
      character*129 ar__5  	! Full-size record

      character*27  Trans_5     ! Transition
      character*3   Mult_1      ! Multiplet
      real*8        lambda_5    ! (Angstrom) Wavelength of transition {lambda}
      real*8        a_1         ! a coefficient (1)
      real*8        b_1         ! b coefficient (1)
      real*8        c_1         ! c coefficient (1)
      real*8        d_1         ! d coefficient (1)
      real*8        e_1         ! e coefficient (1)
      real*8        f_1         ! f coefficient (1)
      real*8        g_1         ! g coefficient (1)
      real*8        h_1         ! h coefficient (1)
      real*4        del_2       ! (%) Average deviation of the fit (delta)
      real*4        Del_3       ! (%) Maximum deviation of the fit {Delta}
*Note (1): The expression of the effective recombination coefficient is
*   log({alpha}_eff_)+15 = a+b*t+c*t^2^+(d+e*t+f*t^2^)log(t)+g*(log(t))^2^+h/t,
*      where t is the reduced electronic temperature Te/10^4^K.
*   table7.dat, table8.dat, table9.dat and table10.dat are fit results for
*      Case B recombination;
*   table11.dat, table12.dat, table13.dat and table14.dat are fit results for
*      Case A recombination.

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table9.dat'	! Fit parameters and the average and maximum
                              fitting errors for 125<=Te<=20000K and
                              Ne=10^4^cm^-3^ for Case B recombination

      integer*4 nr__6
      parameter (nr__6=55)	! Number of records
      character*129 ar__6  	! Full-size record

      character*27  Trans_6     ! Transition
      character*3   Mult_2      ! Multiplet
      real*8        lambda_6    ! (Angstrom) Wavelength of transition {lambda}
      real*8        a_2         ! a coefficient (1)
      real*8        b_2         ! b coefficient (1)
      real*8        c_2         ! c coefficient (1)
      real*8        d_2         ! d coefficient (1)
      real*8        e_2         ! e coefficient (1)
      real*8        f_2         ! f coefficient (1)
      real*8        g_2         ! g coefficient (1)
      real*8        h_2         ! h coefficient (1)
      real*4        del_4       ! (%) Average deviation of the fit (delta)
      real*4        Del_5       ! (%) Maximum deviation of the fit {Delta}
*Note (1): The expression of the effective recombination coefficient is
*   log({alpha}_eff_)+15 = a+b*t+c*t^2^+(d+e*t+f*t^2^)log(t)+g*(log(t))^2^+h/t,
*      where t is the reduced electronic temperature Te/10^4^K.
*   table7.dat, table8.dat, table9.dat and table10.dat are fit results for
*      Case B recombination;
*   table11.dat, table12.dat, table13.dat and table14.dat are fit results for
*      Case A recombination.

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table10.dat'	! Fit parameters and the average and maximum
                              fitting errors for 125<=Te<=20000K and
                              Ne=10^5^cm^-3^ for Case B recombination

      integer*4 nr__7
      parameter (nr__7=55)	! Number of records
      character*129 ar__7  	! Full-size record

      character*27  Trans_7     ! Transition
      character*3   Mult_3      ! Multiplet
      real*8        lambda_7    ! (Angstrom) Wavelength of transition {lambda}
      real*8        a_3         ! a coefficient (1)
      real*8        b_3         ! b coefficient (1)
      real*8        c_3         ! c coefficient (1)
      real*8        d_3         ! d coefficient (1)
      real*8        e_3         ! e coefficient (1)
      real*8        f_3         ! f coefficient (1)
      real*8        g_3         ! g coefficient (1)
      real*8        h_3         ! h coefficient (1)
      real*4        del_6       ! (%) Average deviation of the fit (delta)
      real*4        Del_7       ! (%) Maximum deviation of the fit {Delta}
*Note (1): The expression of the effective recombination coefficient is
*   log({alpha}_eff_)+15 = a+b*t+c*t^2^+(d+e*t+f*t^2^)log(t)+g*(log(t))^2^+h/t,
*      where t is the reduced electronic temperature Te/10^4^K.
*   table7.dat, table8.dat, table9.dat and table10.dat are fit results for
*      Case B recombination;
*   table11.dat, table12.dat, table13.dat and table14.dat are fit results for
*      Case A recombination.

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table11.dat'	! Fit parameters and the average and maximum
                              fitting errors for 125<=Te<=20000K and
                              Ne=10^2^cm^-3^ for Case A recombination

      integer*4 nr__8
      parameter (nr__8=55)	! Number of records
      character*129 ar__8  	! Full-size record

      character*27  Trans_8     ! Transition
      character*3   Mult_4      ! Multiplet
      real*8        lambda_8    ! (Angstrom) Wavelength of transition {lambda}
      real*8        a_4         ! a coefficient (1)
      real*8        b_4         ! b coefficient (1)
      real*8        c_4         ! c coefficient (1)
      real*8        d_4         ! d coefficient (1)
      real*8        e_4         ! e coefficient (1)
      real*8        f_4         ! f coefficient (1)
      real*8        g_4         ! g coefficient (1)
      real*8        h_4         ! h coefficient (1)
      real*4        del_8       ! (%) Average deviation of the fit (delta)
      real*4        Del_9       ! (%) Maximum deviation of the fit {Delta}
*Note (1): The expression of the effective recombination coefficient is
*   log({alpha}_eff_)+15 = a+b*t+c*t^2^+(d+e*t+f*t^2^)log(t)+g*(log(t))^2^+h/t,
*      where t is the reduced electronic temperature Te/10^4^K.
*   table7.dat, table8.dat, table9.dat and table10.dat are fit results for
*      Case B recombination;
*   table11.dat, table12.dat, table13.dat and table14.dat are fit results for
*      Case A recombination.

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table12.dat'	! Fit parameters and the average and maximum
                              fitting errors for 125<=Te<=20000K and
                              Ne=10^3^cm^-3^ for Case A recombination

      integer*4 nr__9
      parameter (nr__9=55)	! Number of records
      character*129 ar__9  	! Full-size record

      character*27  Trans_9     ! Transition
      character*3   Mult_5      ! Multiplet
      real*8        lambda_9    ! (Angstrom) Wavelength of transition {lambda}
      real*8        a_5         ! a coefficient (1)
      real*8        b_5         ! b coefficient (1)
      real*8        c_5         ! c coefficient (1)
      real*8        d_5         ! d coefficient (1)
      real*8        e_5         ! e coefficient (1)
      real*8        f_5         ! f coefficient (1)
      real*8        g_5         ! g coefficient (1)
      real*8        h_5         ! h coefficient (1)
      real*4        del_10      ! (%) Average deviation of the fit (delta)
      real*4        Del_11      ! (%) Maximum deviation of the fit {Delta}
*Note (1): The expression of the effective recombination coefficient is
*   log({alpha}_eff_)+15 = a+b*t+c*t^2^+(d+e*t+f*t^2^)log(t)+g*(log(t))^2^+h/t,
*      where t is the reduced electronic temperature Te/10^4^K.
*   table7.dat, table8.dat, table9.dat and table10.dat are fit results for
*      Case B recombination;
*   table11.dat, table12.dat, table13.dat and table14.dat are fit results for
*      Case A recombination.

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table13.dat'	! Fit parameters and the average and maximum
                              fitting errors for 125<=Te<=20000K and
                              Ne=10^4^cm^-3^ for Case A recombination

      integer*4 nr__10
      parameter (nr__10=55)	! Number of records
      character*129 ar__10 	! Full-size record

      character*27  Trans_10    ! Transition
      character*3   Mult_6      ! Multiplet
      real*8        lambda_10   ! (Angstrom) Wavelength of transition {lambda}
      real*8        a_6         ! a coefficient (1)
      real*8        b_6         ! b coefficient (1)
      real*8        c_6         ! c coefficient (1)
      real*8        d_6         ! d coefficient (1)
      real*8        e_6         ! e coefficient (1)
      real*8        f_6         ! f coefficient (1)
      real*8        g_6         ! g coefficient (1)
      real*8        h_6         ! h coefficient (1)
      real*4        del_12      ! (%) Average deviation of the fit (delta)
      real*4        Del_13      ! (%) Maximum deviation of the fit {Delta}
*Note (1): The expression of the effective recombination coefficient is
*   log({alpha}_eff_)+15 = a+b*t+c*t^2^+(d+e*t+f*t^2^)log(t)+g*(log(t))^2^+h/t,
*      where t is the reduced electronic temperature Te/10^4^K.
*   table7.dat, table8.dat, table9.dat and table10.dat are fit results for
*      Case B recombination;
*   table11.dat, table12.dat, table13.dat and table14.dat are fit results for
*      Case A recombination.

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table14.dat'	! Fit parameters and the average and maximum
                              fitting errors for 125<=Te<=20000K and
                              Ne=10^5^cm^-3^ for Case A recombination

      integer*4 nr__11
      parameter (nr__11=55)	! Number of records
      character*129 ar__11 	! Full-size record

      character*27  Trans_11    ! Transition
      character*3   Mult_7      ! Multiplet
      real*8        lambda_11   ! (Angstrom) Wavelength of transition {lambda}
      real*8        a_7         ! a coefficient (1)
      real*8        b_7         ! b coefficient (1)
      real*8        c_7         ! c coefficient (1)
      real*8        d_7         ! d coefficient (1)
      real*8        e_7         ! e coefficient (1)
      real*8        f_7         ! f coefficient (1)
      real*8        g_7         ! g coefficient (1)
      real*8        h_7         ! h coefficient (1)
      real*4        del_14      ! (%) Average deviation of the fit (delta)
      real*4        Del_15      ! (%) Maximum deviation of the fit {Delta}
*Note (1): The expression of the effective recombination coefficient is
*   log({alpha}_eff_)+15 = a+b*t+c*t^2^+(d+e*t+f*t^2^)log(t)+g*(log(t))^2^+h/t,
*      where t is the reduced electronic temperature Te/10^4^K.
*   table7.dat, table8.dat, table9.dat and table10.dat are fit results for
*      Case B recombination;
*   table11.dat, table12.dat, table13.dat and table14.dat are fit results for
*      Case A recombination.

c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - 

C  Declarations for 'table15.dat'	! *Comparison of our direct recombination
                              coefficients to states of N^+^ with those of
                              Kisielius & Storey (2002, Cat. J/A+A/387/1135)
                              and of Nahar (1995ApJS..101..423N)

      integer*4 nr__12
      parameter (nr__12=26)	! Number of records
      character*70 ar__12 	! Full-size record

      character*15  State       ! State (1)
      real*4        DRCN3       ! (cm+3/s) ?=- Nahar (1995ApJS..101..423N) direct
*                                      recombination coefficient for Te=1000K
      real*4        DRCK3       ! (cm+3/s) ?=- Kisielius & Storey (J/A+A/387/1135) direct
*                                      recombination coefficient for Te=1000K
      real*4        DRCP3       ! (cm+3/s) Present work direct recombination coefficient
*                                    for Te=1000K
      real*4        DRCN4       ! (cm+3/s) ?=- Nahar (1995ApJS..101..423N) direct
*                                      recombination coefficient for Te=10000K
      real*4        DRCK4       ! (cm+3/s) ?=- Kisielius & Storey (J/A+A/387/1135) direct
*                                      recombination coefficient for Te=10000K
      real*4        DRCP4       ! (cm+3/s) Present work direct recombination coefficient
*                                    for Te=10000K
*Note (1): Two last lines are for Sum and Total

C=============================================================================

C  Loading file 'table3.dat'	! Case B effective recombination coefficients
*                              for electron density Ne=10^2^cm^-3^

C  Format for file interpretation

    1 format(
     +  A21,2X,A31,2X,F8.2,1X,A1,1X,E7.2,1X,E7.2,1X,E7.2,1X,E7.2,1X,
     +  E7.2,1X,E7.2,1X,E7.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table3.dat')
      write(6,*) '....Loading file: table3.dat'
      do i__=1,217
        read(1,'(A122)')ar__
        read(ar__,1)
     +  Tr,Trans,lambda,n_lambda,ERC1,ERC2,ERC3,ERC4,ERC5,ERC6,ERC7
c    ..............Just test output...........
        write(6,1)
     +  Tr,Trans,lambda,n_lambda,ERC1,ERC2,ERC3,ERC4,ERC5,ERC6,ERC7
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table4.dat'	! Case B effective recombination coefficients
*                              for electron density Ne=10^3^cm^-3^

C  Format for file interpretation

    2 format(
     +  A21,2X,A31,2X,F8.2,1X,A1,1X,E7.2,1X,E7.2,1X,E7.2,1X,E7.2,1X,
     +  E7.2,1X,E7.2,1X,E7.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table4.dat')
      write(6,*) '....Loading file: table4.dat'
      do i__=1,217
        read(1,'(A122)')ar__1
        read(ar__1,2)
     +  Tr_1,Trans_1,lambda_1,n_lambda_1,ERC1_1,ERC2_1,ERC3_1,ERC4_1,
     +  ERC5_1,ERC6_1,ERC7_1
c    ..............Just test output...........
        write(6,2)
     +  Tr_1,Trans_1,lambda_1,n_lambda_1,ERC1_1,ERC2_1,ERC3_1,ERC4_1,
     +  ERC5_1,ERC6_1,ERC7_1
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table5.dat'	! Case B effective recombination coefficients
*                              for electron density Ne=10^4^cm^-3^

C  Format for file interpretation

    3 format(
     +  A21,2X,A31,2X,F8.2,1X,A1,1X,E7.2,1X,E7.2,1X,E7.2,1X,E7.2,1X,
     +  E7.2,1X,E7.2,1X,E7.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table5.dat')
      write(6,*) '....Loading file: table5.dat'
      do i__=1,217
        read(1,'(A122)')ar__2
        read(ar__2,3)
     +  Tr_2,Trans_2,lambda_2,n_lambda_2,ERC1_2,ERC2_2,ERC3_2,ERC4_2,
     +  ERC5_2,ERC6_2,ERC7_2
c    ..............Just test output...........
        write(6,3)
     +  Tr_2,Trans_2,lambda_2,n_lambda_2,ERC1_2,ERC2_2,ERC3_2,ERC4_2,
     +  ERC5_2,ERC6_2,ERC7_2
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table6.dat'	! Case B effective recombination coefficients
*                              for electron density Ne=10^5^cm^-3^

C  Format for file interpretation

    4 format(
     +  A21,2X,A31,2X,F8.2,1X,A1,1X,E7.2,1X,E7.2,1X,E7.2,1X,E7.2,1X,
     +  E7.2,1X,E7.2,1X,E7.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table6.dat')
      write(6,*) '....Loading file: table6.dat'
      do i__=1,217
        read(1,'(A122)')ar__3
        read(ar__3,4)
     +  Tr_3,Trans_3,lambda_3,n_lambda_3,ERC1_3,ERC2_3,ERC3_3,ERC4_3,
     +  ERC5_3,ERC6_3,ERC7_3
c    ..............Just test output...........
        write(6,4)
     +  Tr_3,Trans_3,lambda_3,n_lambda_3,ERC1_3,ERC2_3,ERC3_3,ERC4_3,
     +  ERC5_3,ERC6_3,ERC7_3
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table7.dat'	! Fit parameters and the average and maximum
*                              fitting errors for 125<=Te<=20000K and
*                              Ne=10^2^cm^-3^ for Case B recombination

C  Format for file interpretation

    5 format(
     +  A27,2X,A3,2X,F7.2,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,
     +  F8.4,1X,F8.4,2X,F9.6,2X,F5.3,2X,F5.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table7.dat')
      write(6,*) '....Loading file: table7.dat'
      do i__=1,55
        read(1,'(A129)')ar__4
        read(ar__4,5)Trans_4,Mult,lambda_4,a,b,c,d,e,f,g,h,del,Del_1
c    ..............Just test output...........
        write(6,5)Trans_4,Mult,lambda_4,a,b,c,d,e,f,g,h,del,Del_1
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table8.dat'	! Fit parameters and the average and maximum
*                              fitting errors for 125<=Te<=20000K and
*                              Ne=10^3^cm^-3^ for Case B recombination

C  Format for file interpretation

    6 format(
     +  A27,2X,A3,2X,F7.2,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,
     +  F8.4,1X,F8.4,2X,F9.6,2X,F5.3,2X,F5.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table8.dat')
      write(6,*) '....Loading file: table8.dat'
      do i__=1,55
        read(1,'(A129)')ar__5
        read(ar__5,6)
     +  Trans_5,Mult_1,lambda_5,a_1,b_1,c_1,d_1,e_1,f_1,g_1,h_1,del_2,
     +  Del_3
c    ..............Just test output...........
        write(6,6)
     +  Trans_5,Mult_1,lambda_5,a_1,b_1,c_1,d_1,e_1,f_1,g_1,h_1,del_2,
     +  Del_3
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table9.dat'	! Fit parameters and the average and maximum
*                              fitting errors for 125<=Te<=20000K and
*                              Ne=10^4^cm^-3^ for Case B recombination

C  Format for file interpretation

    7 format(
     +  A27,2X,A3,2X,F7.2,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,
     +  F8.4,1X,F8.4,2X,F9.6,2X,F5.3,2X,F5.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table9.dat')
      write(6,*) '....Loading file: table9.dat'
      do i__=1,55
        read(1,'(A129)')ar__6
        read(ar__6,7)
     +  Trans_6,Mult_2,lambda_6,a_2,b_2,c_2,d_2,e_2,f_2,g_2,h_2,del_4,
     +  Del_5
c    ..............Just test output...........
        write(6,7)
     +  Trans_6,Mult_2,lambda_6,a_2,b_2,c_2,d_2,e_2,f_2,g_2,h_2,del_4,
     +  Del_5
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table10.dat'	! Fit parameters and the average and maximum
*                              fitting errors for 125<=Te<=20000K and
*                              Ne=10^5^cm^-3^ for Case B recombination

C  Format for file interpretation

    8 format(
     +  A27,2X,A3,2X,F7.2,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,
     +  F8.4,1X,F8.4,2X,F9.6,2X,F5.3,2X,F5.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table10.dat')
      write(6,*) '....Loading file: table10.dat'
      do i__=1,55
        read(1,'(A129)')ar__7
        read(ar__7,8)
     +  Trans_7,Mult_3,lambda_7,a_3,b_3,c_3,d_3,e_3,f_3,g_3,h_3,del_6,
     +  Del_7
c    ..............Just test output...........
        write(6,8)
     +  Trans_7,Mult_3,lambda_7,a_3,b_3,c_3,d_3,e_3,f_3,g_3,h_3,del_6,
     +  Del_7
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table11.dat'	! Fit parameters and the average and maximum
*                              fitting errors for 125<=Te<=20000K and
*                              Ne=10^2^cm^-3^ for Case A recombination

C  Format for file interpretation

    9 format(
     +  A27,2X,A3,2X,F7.2,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,
     +  F8.4,1X,F8.4,2X,F9.6,2X,F5.3,2X,F5.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table11.dat')
      write(6,*) '....Loading file: table11.dat'
      do i__=1,55
        read(1,'(A129)')ar__8
        read(ar__8,9)
     +  Trans_8,Mult_4,lambda_8,a_4,b_4,c_4,d_4,e_4,f_4,g_4,h_4,del_8,
     +  Del_9
c    ..............Just test output...........
        write(6,9)
     +  Trans_8,Mult_4,lambda_8,a_4,b_4,c_4,d_4,e_4,f_4,g_4,h_4,del_8,
     +  Del_9
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table12.dat'	! Fit parameters and the average and maximum
*                              fitting errors for 125<=Te<=20000K and
*                              Ne=10^3^cm^-3^ for Case A recombination

C  Format for file interpretation

   10 format(
     +  A27,2X,A3,2X,F7.2,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,
     +  F8.4,1X,F8.4,2X,F9.6,2X,F5.3,2X,F5.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table12.dat')
      write(6,*) '....Loading file: table12.dat'
      do i__=1,55
        read(1,'(A129)')ar__9
        read(ar__9,10)
     +  Trans_9,Mult_5,lambda_9,a_5,b_5,c_5,d_5,e_5,f_5,g_5,h_5,
     +  del_10,Del_11
c    ..............Just test output...........
        write(6,10)
     +  Trans_9,Mult_5,lambda_9,a_5,b_5,c_5,d_5,e_5,f_5,g_5,h_5,
     +  del_10,Del_11
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table13.dat'	! Fit parameters and the average and maximum
*                              fitting errors for 125<=Te<=20000K and
*                              Ne=10^4^cm^-3^ for Case A recombination

C  Format for file interpretation

   11 format(
     +  A27,2X,A3,2X,F7.2,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,
     +  F8.4,1X,F8.4,2X,F9.6,2X,F5.3,2X,F5.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table13.dat')
      write(6,*) '....Loading file: table13.dat'
      do i__=1,55
        read(1,'(A129)')ar__10
        read(ar__10,11)
     +  Trans_10,Mult_6,lambda_10,a_6,b_6,c_6,d_6,e_6,f_6,g_6,h_6,
     +  del_12,Del_13
c    ..............Just test output...........
        write(6,11)
     +  Trans_10,Mult_6,lambda_10,a_6,b_6,c_6,d_6,e_6,f_6,g_6,h_6,
     +  del_12,Del_13
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table14.dat'	! Fit parameters and the average and maximum
*                              fitting errors for 125<=Te<=20000K and
*                              Ne=10^5^cm^-3^ for Case A recombination

C  Format for file interpretation

   12 format(
     +  A27,2X,A3,2X,F7.2,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,F8.4,1X,
     +  F8.4,1X,F8.4,2X,F9.6,2X,F5.3,2X,F5.3)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table14.dat')
      write(6,*) '....Loading file: table14.dat'
      do i__=1,55
        read(1,'(A129)')ar__11
        read(ar__11,12)
     +  Trans_11,Mult_7,lambda_11,a_7,b_7,c_7,d_7,e_7,f_7,g_7,h_7,
     +  del_14,Del_15
c    ..............Just test output...........
        write(6,12)
     +  Trans_11,Mult_7,lambda_11,a_7,b_7,c_7,d_7,e_7,f_7,g_7,h_7,
     +  del_14,Del_15
c    .......End.of.Just test output...........
      end do
      close(1)

C=============================================================================

C  Loading file 'table15.dat'	! *Comparison of our direct recombination
*                              coefficients to states of N^+^ with those of
*                              Kisielius & Storey (2002, Cat. J/A+A/387/1135)
*                              and of Nahar (1995ApJS..101..423N)

C  Format for file interpretation

   13 format(A15,2X,E8.2,1X,E8.2,1X,E8.2,1X,E8.2,1X,E8.2,1X,E8.2)

C  Effective file loading

      open(unit=1,status='old',file=
     +'table15.dat')
      write(6,*) '....Loading file: table15.dat'
      do i__=1,26
        read(1,'(A70)')ar__12
        read(ar__12,13)State,DRCN3,DRCK3,DRCP3,DRCN4,DRCK4,DRCP4
        if (idig(ar__12(18:25)).EQ.0) DRCN3 =  rNULL__
        if (idig(ar__12(27:34)).EQ.0) DRCK3 =  rNULL__
        if (idig(ar__12(45:52)).EQ.0) DRCN4 =  rNULL__
        if (idig(ar__12(54:61)).EQ.0) DRCK4 =  rNULL__
c    ..............Just test output...........
        write(6,13)State,DRCN3,DRCK3,DRCP3,DRCN4,DRCK4,DRCP4
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