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