Conversion of standardized ReadMe file for
file /./ftp/cats/J/A_A/513/A61 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-Apr-18
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/513/A61 Neutron-induced astrophysical reaction rates (Panov+, 2010)
*================================================================================
*Neutron-induced astrophysical reaction rates for translead nuclei.
* Panov I.V., Korneev I.Yu., Rauscher T., Martinez-Pinedo G., Kelic-Heil A.,
* Zinner N.T., Thielemann F.-K.
* <Astron. Astrophys. 513, A61 (2010)>
* =2010A&A...513A..61P (SIMBAD/NED BibCode)
C=============================================================================
C Internal variables
integer*4 i__
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table3.dat' ! (n,g) and (n,f) rates with ETFSI-predictions
integer*4 nr__
parameter (nr__=64032) ! Number of records
character*81 ar__ ! Full-size record
character*2 El (nr__) ! Target element (corresponding to Z)
integer*4 A (nr__) ! Atomic mass number of the target
integer*4 Z (nr__) ! [84,118] Charge number of the target
real*4 T9 (nr__) ! (GK) [0.10/10.00] Temperature
real*4 pf (nr__) ! Partition function of the target
real*4 nf (nr__) ! (cm3/s/mol) Neutron-induced fission rate for ground state
real*4 nf_ (nr__) ! (cm3/s/mol) Neutron-induced fission rate with thermally
* populated target states
real*4 ng (nr__) ! (cm3/s/mol) Neutron capture rate for the ground state
real*4 ng_ (nr__) ! (cm3/s/mol) Neutron capture rate with thermally
* populated target states
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table4.dat' ! (n,g) and (n,f) rates with TF-predictions
integer*4 nr__1
parameter (nr__1=71760) ! Number of records
character*81 ar__1 ! Full-size record
character*2 El_1 (nr__1) ! Target element (corresponding to Z)
integer*4 A_1 (nr__1) ! Atomic mass number of the target
integer*4 Z_1 (nr__1) ! [84,118] Charge number of the target
real*4 T9_1 (nr__1) ! (GK) [0.10/10.00] Temperature
real*4 pf_1 (nr__1) ! Partition function of the target
real*4 nf_1 (nr__1) ! (cm3/s/mol) Neutron-induced fission rate for ground state
real*4 nf__1 (nr__1) ! (cm3/s/mol) Neutron-induced fission rate with thermally
* populated target states
real*4 ng_1 (nr__1) ! (cm3/s/mol) Neutron capture rate for the ground state
real*4 ng__1 (nr__1) ! (cm3/s/mol) Neutron capture rate with thermally
* populated target states
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table5.dat' ! (n,g) and (n,f) rates with FRDM+TF-predictions
integer*4 nr__2
parameter (nr__2=70056) ! Number of records
character*81 ar__2 ! Full-size record
character*2 El_2 (nr__2) ! Target element (corresponding to Z)
integer*4 A_2 (nr__2) ! Atomic mass number of the target
integer*4 Z_2 (nr__2) ! [84,118] Charge number of the target
real*4 T9_2 (nr__2) ! (GK) [0.10/10.00] Temperature
real*4 pf_2 (nr__2) ! Partition function of the target
real*4 nf_2 (nr__2) ! (cm3/s/mol) Neutron-induced fission rate for ground state
real*4 nf__2 (nr__2) ! (cm3/s/mol) Neutron-induced fission rate with thermally
* populated target states
real*4 ng_2 (nr__2) ! (cm3/s/mol) Neutron capture rate for the ground state
real*4 ng__2 (nr__2) ! (cm3/s/mol) Neutron capture rate with thermally
* populated target states
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table6.dat' ! (n,g) and (n,f) rates with HFB14-predictions
integer*4 nr__3
parameter (nr__3=31728) ! Number of records
character*81 ar__3 ! Full-size record
character*2 El_3 (nr__3) ! Target element (corresponding to Z)
integer*4 A_3 (nr__3) ! Atomic mass number of the target
integer*4 Z_3 (nr__3) ! [84,118] Charge number of the target
real*4 T9_3 (nr__3) ! (GK) [0.10/10.00] Temperature
real*4 pf_3 (nr__3) ! Partition function of the target
real*4 nf_3 (nr__3) ! (cm3/s/mol) Neutron-induced fission rate for ground state
real*4 nf__3 (nr__3) ! (cm3/s/mol) Neutron-induced fission rate with thermally
* populated target states
real*4 ng_3 (nr__3) ! (cm3/s/mol) Neutron capture rate for the ground state
real*4 ng__3 (nr__3) ! (cm3/s/mol) Neutron capture rate with thermally
* populated target states
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table7.dat' ! Coefficients for (n,g)-rate fit, ETFSI
integer*4 nr__4
parameter (nr__4=2662) ! Number of records
character*124 ar__4 ! Full-size record
character*2 El_4 (nr__4) ! Target element (corresponding to Z)
integer*4 A_4 (nr__4) ! Atomic mass number of the target
integer*4 Z_4 (nr__4) ! [84,118] Charge number of the target
integer*4 i (nr__4) ! [0/1] Number of fitting curve (1)
real*4 a0 (nr__4) ! Coefficient a0 of the forward reaction r_i_ (1)
real*4 a1 (nr__4) ! Coefficient a1 of the forward reaction r_i_ (1)
real*4 a2 (nr__4) ! Coefficient a2 of the forward reaction r_i_ (1)
real*4 a3 (nr__4) ! Coefficient a3 of the forward reaction r_i_ (1)
real*4 a4 (nr__4) ! Coefficient a4 of the forward reaction r_i_ (1)
real*4 a5 (nr__4) ! Coefficient a5 of the forward reaction r_i_ (1)
real*4 a6 (nr__4) ! Coefficient a6 of the forward reaction r_i_ (1)
real*4 zeta (nr__4) ! ? Mean square error for direct reactions (Eq.(5))
*Note (1): The rate r is calculated as r={Sum}r_i_, and each r_i_ is computed as
* r_i_ = exp(a0+a1/T9+a2/T9^1/3^+a3*T9^1/3^+a4*T9+a5*T9^5/3^+a6*ln(T9)).
* When i=0, there is a single seven-parameter set to fit the rate.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table8.dat' ! Coefficients for (n,g)-rate fit, TF
integer*4 nr__5
parameter (nr__5=2984) ! Number of records
character*124 ar__5 ! Full-size record
character*2 El_5 (nr__5) ! Target element (corresponding to Z)
integer*4 A_5 (nr__5) ! Atomic mass number of the target
integer*4 Z_5 (nr__5) ! [84,118] Charge number of the target
integer*4 i_1 (nr__5) ! [0/1] Number of fitting curve (1)
real*4 a0_1 (nr__5) ! Coefficient a0 of the forward reaction r_i_ (1)
real*4 a1_1 (nr__5) ! Coefficient a1 of the forward reaction r_i_ (1)
real*4 a2_1 (nr__5) ! Coefficient a2 of the forward reaction r_i_ (1)
real*4 a3_1 (nr__5) ! Coefficient a3 of the forward reaction r_i_ (1)
real*4 a4_1 (nr__5) ! Coefficient a4 of the forward reaction r_i_ (1)
real*4 a5_1 (nr__5) ! Coefficient a5 of the forward reaction r_i_ (1)
real*4 a6_1 (nr__5) ! Coefficient a6 of the forward reaction r_i_ (1)
real*4 zeta_1 (nr__5) ! ? Mean square error for direct reactions (Eq.(5))
*Note (1): The rate r is calculated as r={Sum}r_i_, and each r_i_ is computed as
* r_i_ = exp(a0+a1/T9+a2/T9^1/3^+a3*T9^1/3^+a4*T9+a5*T9^5/3^+a6*ln(T9)).
* When i=0, there is a single seven-parameter set to fit the rate.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table9.dat' ! Coefficients for (n,g)-rate fit, FRDM+TF
integer*4 nr__6
parameter (nr__6=2911) ! Number of records
character*124 ar__6 ! Full-size record
character*2 El_6 (nr__6) ! Target element (corresponding to Z)
integer*4 A_6 (nr__6) ! Atomic mass number of the target
integer*4 Z_6 (nr__6) ! [84,118] Charge number of the target
integer*4 i_2 (nr__6) ! [0/1] Number of fitting curve (1)
real*4 a0_2 (nr__6) ! Coefficient a0 of the forward reaction r_i_ (1)
real*4 a1_2 (nr__6) ! Coefficient a1 of the forward reaction r_i_ (1)
real*4 a2_2 (nr__6) ! Coefficient a2 of the forward reaction r_i_ (1)
real*4 a3_2 (nr__6) ! Coefficient a3 of the forward reaction r_i_ (1)
real*4 a4_2 (nr__6) ! Coefficient a4 of the forward reaction r_i_ (1)
real*4 a5_2 (nr__6) ! Coefficient a5 of the forward reaction r_i_ (1)
real*4 a6_2 (nr__6) ! Coefficient a6 of the forward reaction r_i_ (1)
real*4 zeta_2 (nr__6) ! ? Mean square error for direct reactions (Eq.(5))
*Note (1): The rate r is calculated as r={Sum}r_i_, and each r_i_ is computed as
* r_i_ = exp(a0+a1/T9+a2/T9^1/3^+a3*T9^1/3^+a4*T9+a5*T9^5/3^+a6*ln(T9)).
* When i=0, there is a single seven-parameter set to fit the rate.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table10.dat' ! Coefficients for (n,g)-rate fit, HFB14
integer*4 nr__7
parameter (nr__7=1468) ! Number of records
character*124 ar__7 ! Full-size record
character*2 El_7 (nr__7) ! Target element (corresponding to Z)
integer*4 A_7 (nr__7) ! Atomic mass number of the target
integer*4 Z_7 (nr__7) ! [84,118] Charge number of the target
integer*4 i_3 (nr__7) ! [0/1] Number of fitting curve (1)
real*4 a0_3 (nr__7) ! Coefficient a0 of the forward reaction r_i_ (1)
real*4 a1_3 (nr__7) ! Coefficient a1 of the forward reaction r_i_ (1)
real*4 a2_3 (nr__7) ! Coefficient a2 of the forward reaction r_i_ (1)
real*4 a3_3 (nr__7) ! Coefficient a3 of the forward reaction r_i_ (1)
real*4 a4_3 (nr__7) ! Coefficient a4 of the forward reaction r_i_ (1)
real*4 a5_3 (nr__7) ! Coefficient a5 of the forward reaction r_i_ (1)
real*4 a6_3 (nr__7) ! Coefficient a6 of the forward reaction r_i_ (1)
real*4 zeta_3 (nr__7) ! ? Mean square error for direct reactions (Eq.(5))
*Note (1): The rate r is calculated as r={Sum}r_i_, and each r_i_ is computed as
* r_i_ = exp(a0+a1/T9+a2/T9^1/3^+a3*T9^1/3^+a4*T9+a5*T9^5/3^+a6*ln(T9)).
* When i=0, there is a single seven-parameter set to fit the rate.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table11.dat' ! Coefficients for (g,n)-rate fit, ETFSI
integer*4 nr__8
parameter (nr__8=2414) ! Number of records
character*124 ar__8 ! Full-size record
character*2 El_8 (nr__8) ! Target element (corresponding to Z)
integer*4 A_8 (nr__8) ! Atomic mass number of the target
integer*4 Z_8 (nr__8) ! [84,118] Charge number of the target
integer*4 i_4 (nr__8) ! [0/1] Number of fitting curve (1)
real*4 a0_4 (nr__8) ! Coefficient a0 of the forward reaction r_i_ (1)
real*4 a1_4 (nr__8) ! Coefficient a1 of the forward reaction r_i_ (1)
real*4 a2_4 (nr__8) ! Coefficient a2 of the forward reaction r_i_ (1)
real*4 a3_4 (nr__8) ! Coefficient a3 of the forward reaction r_i_ (1)
real*4 a4_4 (nr__8) ! Coefficient a4 of the forward reaction r_i_ (1)
real*4 a5_4 (nr__8) ! Coefficient a5 of the forward reaction r_i_ (1)
real*4 a6_4 (nr__8) ! Coefficient a6 of the forward reaction r_i_ (1)
real*4 zeta_4 (nr__8) ! ? Mean square error for direct reactions (Eq.(5))
*Note (1): The rate r is calculated as r={Sum}r_i_, and each r_i_ is computed as
* r_i_ = exp(a0+a1/T9+a2/T9^1/3^+a3*T9^1/3^+a4*T9+a5*T9^5/3^+a6*ln(T9)).
* When i=0, there is a single seven-parameter set to fit the rate.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table12.dat' ! Coefficients for (g,n)-rate fit, TF
integer*4 nr__9
parameter (nr__9=2948) ! Number of records
character*124 ar__9 ! Full-size record
character*2 El_9 (nr__9) ! Target element (corresponding to Z)
integer*4 A_9 (nr__9) ! Atomic mass number of the target
integer*4 Z_9 (nr__9) ! [84,118] Charge number of the target
integer*4 i_5 (nr__9) ! [0/1] Number of fitting curve (1)
real*4 a0_5 (nr__9) ! Coefficient a0 of the forward reaction r_i_ (1)
real*4 a1_5 (nr__9) ! Coefficient a1 of the forward reaction r_i_ (1)
real*4 a2_5 (nr__9) ! Coefficient a2 of the forward reaction r_i_ (1)
real*4 a3_5 (nr__9) ! Coefficient a3 of the forward reaction r_i_ (1)
real*4 a4_5 (nr__9) ! Coefficient a4 of the forward reaction r_i_ (1)
real*4 a5_5 (nr__9) ! Coefficient a5 of the forward reaction r_i_ (1)
real*4 a6_5 (nr__9) ! Coefficient a6 of the forward reaction r_i_ (1)
real*4 zeta_5 (nr__9) ! ? Mean square error for direct reactions (Eq.(5))
*Note (1): The rate r is calculated as r={Sum}r_i_, and each r_i_ is computed as
* r_i_ = exp(a0+a1/T9+a2/T9^1/3^+a3*T9^1/3^+a4*T9+a5*T9^5/3^+a6*ln(T9)).
* When i=0, there is a single seven-parameter set to fit the rate.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table13.dat' ! Coefficients for (g,n)-rate fit, FRDM+TF
integer*4 nr__10
parameter (nr__10=2875) ! Number of records
character*124 ar__10 ! Full-size record
character*2 El_10 (nr__10) ! Target element (corresponding to Z)
integer*4 A_10 (nr__10) ! Atomic mass number of the target
integer*4 Z_10 (nr__10) ! [84,118] Charge number of the target
integer*4 i_6 (nr__10) ! [0/1] Number of fitting curve (1)
real*4 a0_6 (nr__10) ! Coefficient a0 of the forward reaction r_i_ (1)
real*4 a1_6 (nr__10) ! Coefficient a1 of the forward reaction r_i_ (1)
real*4 a2_6 (nr__10) ! Coefficient a2 of the forward reaction r_i_ (1)
real*4 a3_6 (nr__10) ! Coefficient a3 of the forward reaction r_i_ (1)
real*4 a4_6 (nr__10) ! Coefficient a4 of the forward reaction r_i_ (1)
real*4 a5_6 (nr__10) ! Coefficient a5 of the forward reaction r_i_ (1)
real*4 a6_6 (nr__10) ! Coefficient a6 of the forward reaction r_i_ (1)
real*4 zeta_6 (nr__10) ! ? Mean square error for direct reactions (Eq.(5))
*Note (1): The rate r is calculated as r={Sum}r_i_, and each r_i_ is computed as
* r_i_ = exp(a0+a1/T9+a2/T9^1/3^+a3*T9^1/3^+a4*T9+a5*T9^5/3^+a6*ln(T9)).
* When i=0, there is a single seven-parameter set to fit the rate.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table14.dat' ! Coefficients for (g,n)-rate fit, HFB14
integer*4 nr__11
parameter (nr__11=1436) ! Number of records
character*124 ar__11 ! Full-size record
character*2 El_11 (nr__11) ! Target element (corresponding to Z)
integer*4 A_11 (nr__11) ! Atomic mass number of the target
integer*4 Z_11 (nr__11) ! [84,118] Charge number of the target
integer*4 i_7 (nr__11) ! [0/1] Number of fitting curve (1)
real*4 a0_7 (nr__11) ! Coefficient a0 of the forward reaction r_i_ (1)
real*4 a1_7 (nr__11) ! Coefficient a1 of the forward reaction r_i_ (1)
real*4 a2_7 (nr__11) ! Coefficient a2 of the forward reaction r_i_ (1)
real*4 a3_7 (nr__11) ! Coefficient a3 of the forward reaction r_i_ (1)
real*4 a4_7 (nr__11) ! Coefficient a4 of the forward reaction r_i_ (1)
real*4 a5_7 (nr__11) ! Coefficient a5 of the forward reaction r_i_ (1)
real*4 a6_7 (nr__11) ! Coefficient a6 of the forward reaction r_i_ (1)
real*4 zeta_7 (nr__11) ! ? Mean square error for direct reactions (Eq.(5))
*Note (1): The rate r is calculated as r={Sum}r_i_, and each r_i_ is computed as
* r_i_ = exp(a0+a1/T9+a2/T9^1/3^+a3*T9^1/3^+a4*T9+a5*T9^5/3^+a6*ln(T9)).
* When i=0, there is a single seven-parameter set to fit the rate.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table15.dat' ! Coefficients for (n,f)-rate fit, ETFSI
integer*4 nr__12
parameter (nr__12=3288) ! Number of records
character*124 ar__12 ! Full-size record
character*2 El_12 (nr__12) ! Target element (corresponding to Z)
integer*4 A_12 (nr__12) ! Atomic mass number of the target
integer*4 Z_12 (nr__12) ! [84,118] Charge number of the target
integer*4 i_8 (nr__12) ! [0/1] Number of fitting curve (1)
real*4 a0_8 (nr__12) ! Coefficient a0 of the forward reaction r_i_ (1)
real*4 a1_8 (nr__12) ! Coefficient a1 of the forward reaction r_i_ (1)
real*4 a2_8 (nr__12) ! Coefficient a2 of the forward reaction r_i_ (1)
real*4 a3_8 (nr__12) ! Coefficient a3 of the forward reaction r_i_ (1)
real*4 a4_8 (nr__12) ! Coefficient a4 of the forward reaction r_i_ (1)
real*4 a5_8 (nr__12) ! Coefficient a5 of the forward reaction r_i_ (1)
real*4 a6_8 (nr__12) ! Coefficient a6 of the forward reaction r_i_ (1)
real*4 zeta_8 (nr__12) ! ? Mean square error for direct reactions (Eq.(5))
*Note (1): The rate r is calculated as r={Sum}r_i_, and each r_i_ is computed as
* r_i_ = exp(a0+a1/T9+a2/T9^1/3^+a3*T9^1/3^+a4*T9+a5*T9^5/3^+a6*ln(T9)).
* When i=0, there is a single seven-parameter set to fit the rate.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table16.dat' ! Coefficients for (n,f)-rate fit, TF
integer*4 nr__13
parameter (nr__13=3928) ! Number of records
character*124 ar__13 ! Full-size record
character*2 El_13 (nr__13) ! Target element (corresponding to Z)
integer*4 A_13 (nr__13) ! Atomic mass number of the target
integer*4 Z_13 (nr__13) ! [84,118] Charge number of the target
integer*4 i_9 (nr__13) ! [0/1] Number of fitting curve (1)
real*4 a0_9 (nr__13) ! Coefficient a0 of the forward reaction r_i_ (1)
real*4 a1_9 (nr__13) ! Coefficient a1 of the forward reaction r_i_ (1)
real*4 a2_9 (nr__13) ! Coefficient a2 of the forward reaction r_i_ (1)
real*4 a3_9 (nr__13) ! Coefficient a3 of the forward reaction r_i_ (1)
real*4 a4_9 (nr__13) ! Coefficient a4 of the forward reaction r_i_ (1)
real*4 a5_9 (nr__13) ! Coefficient a5 of the forward reaction r_i_ (1)
real*4 a6_9 (nr__13) ! Coefficient a6 of the forward reaction r_i_ (1)
real*4 zeta_9 (nr__13) ! ? Mean square error for direct reactions (Eq.(5))
*Note (1): The rate r is calculated as r={Sum}r_i_, and each r_i_ is computed as
* r_i_ = exp(a0+a1/T9+a2/T9^1/3^+a3*T9^1/3^+a4*T9+a5*T9^5/3^+a6*ln(T9)).
* When i=0, there is a single seven-parameter set to fit the rate.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table17.dat' ! Coefficients for (n,f)-rate fit, FRDM+TF
integer*4 nr__14
parameter (nr__14=3846) ! Number of records
character*124 ar__14 ! Full-size record
character*2 El_14 (nr__14) ! Target element (corresponding to Z)
integer*4 A_14 (nr__14) ! Atomic mass number of the target
integer*4 Z_14 (nr__14) ! [84,118] Charge number of the target
integer*4 i_10 (nr__14) ! [0/1] Number of fitting curve (1)
real*4 a0_10 (nr__14) ! Coefficient a0 of the forward reaction r_i_ (1)
real*4 a1_10 (nr__14) ! Coefficient a1 of the forward reaction r_i_ (1)
real*4 a2_10 (nr__14) ! Coefficient a2 of the forward reaction r_i_ (1)
real*4 a3_10 (nr__14) ! Coefficient a3 of the forward reaction r_i_ (1)
real*4 a4_10 (nr__14) ! Coefficient a4 of the forward reaction r_i_ (1)
real*4 a5_10 (nr__14) ! Coefficient a5 of the forward reaction r_i_ (1)
real*4 a6_10 (nr__14) ! Coefficient a6 of the forward reaction r_i_ (1)
real*4 zeta_10 (nr__14) ! ? Mean square error for direct reactions (Eq.(5))
*Note (1): The rate r is calculated as r={Sum}r_i_, and each r_i_ is computed as
* r_i_ = exp(a0+a1/T9+a2/T9^1/3^+a3*T9^1/3^+a4*T9+a5*T9^5/3^+a6*ln(T9)).
* When i=0, there is a single seven-parameter set to fit the rate.
c - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
C Declarations for 'table18.dat' ! Coefficients for (n,f)-rate fit, HFB14
integer*4 nr__15
parameter (nr__15=1448) ! Number of records
character*124 ar__15 ! Full-size record
character*2 El_15 (nr__15) ! Target element (corresponding to Z)
integer*4 A_15 (nr__15) ! Atomic mass number of the target
integer*4 Z_15 (nr__15) ! [84,118] Charge number of the target
integer*4 i_11 (nr__15) ! [0/1] Number of fitting curve (1)
real*4 a0_11 (nr__15) ! Coefficient a0 of the forward reaction r_i_ (1)
real*4 a1_11 (nr__15) ! Coefficient a1 of the forward reaction r_i_ (1)
real*4 a2_11 (nr__15) ! Coefficient a2 of the forward reaction r_i_ (1)
real*4 a3_11 (nr__15) ! Coefficient a3 of the forward reaction r_i_ (1)
real*4 a4_11 (nr__15) ! Coefficient a4 of the forward reaction r_i_ (1)
real*4 a5_11 (nr__15) ! Coefficient a5 of the forward reaction r_i_ (1)
real*4 a6_11 (nr__15) ! Coefficient a6 of the forward reaction r_i_ (1)
real*4 zeta_11 (nr__15) ! ? Mean square error for direct reactions (Eq.(5))
*Note (1): The rate r is calculated as r={Sum}r_i_, and each r_i_ is computed as
* r_i_ = exp(a0+a1/T9+a2/T9^1/3^+a3*T9^1/3^+a4*T9+a5*T9^5/3^+a6*ln(T9)).
* When i=0, there is a single seven-parameter set to fit the rate.
C=============================================================================
C Loading file 'table3.dat' ! (n,g) and (n,f) rates with ETFSI-predictions
C Format for file interpretation
1 format(
+ A2,2X,I3,1X,I3,3X,F5.2,5X,E9.4,3X,E10.4,2X,E10.4,2X,E9.4,3X,
+ E9.4)
C Effective file loading
open(unit=1,status='old',file=
+'table3.dat')
write(6,*) '....Loading file: table3.dat'
do i__=1,64032
read(1,'(A81)')ar__
read(ar__,1)
+ El(i__),A(i__),Z(i__),T9(i__),pf(i__),nf(i__),nf_(i__),
+ ng(i__),ng_(i__)
c ..............Just test output...........
write(6,1)
+ El(i__),A(i__),Z(i__),T9(i__),pf(i__),nf(i__),nf_(i__),
+ ng(i__),ng_(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table4.dat' ! (n,g) and (n,f) rates with TF-predictions
C Format for file interpretation
2 format(
+ A2,2X,I3,1X,I3,3X,F5.2,5X,E9.4,3X,E10.4,2X,E10.4,2X,E9.4,3X,
+ E9.4)
C Effective file loading
open(unit=1,status='old',file=
+'table4.dat')
write(6,*) '....Loading file: table4.dat'
do i__=1,71760
read(1,'(A81)')ar__1
read(ar__1,2)
+ El_1(i__),A_1(i__),Z_1(i__),T9_1(i__),pf_1(i__),nf_1(i__),
+ nf__1(i__),ng_1(i__),ng__1(i__)
c ..............Just test output...........
write(6,2)
+ El_1(i__),A_1(i__),Z_1(i__),T9_1(i__),pf_1(i__),nf_1(i__),
+ nf__1(i__),ng_1(i__),ng__1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table5.dat' ! (n,g) and (n,f) rates with FRDM+TF-predictions
C Format for file interpretation
3 format(
+ A2,2X,I3,1X,I3,3X,F5.2,5X,E9.4,3X,E10.4,2X,E10.4,2X,E9.4,3X,
+ E9.4)
C Effective file loading
open(unit=1,status='old',file=
+'table5.dat')
write(6,*) '....Loading file: table5.dat'
do i__=1,70056
read(1,'(A81)')ar__2
read(ar__2,3)
+ El_2(i__),A_2(i__),Z_2(i__),T9_2(i__),pf_2(i__),nf_2(i__),
+ nf__2(i__),ng_2(i__),ng__2(i__)
c ..............Just test output...........
write(6,3)
+ El_2(i__),A_2(i__),Z_2(i__),T9_2(i__),pf_2(i__),nf_2(i__),
+ nf__2(i__),ng_2(i__),ng__2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table6.dat' ! (n,g) and (n,f) rates with HFB14-predictions
C Format for file interpretation
4 format(
+ A2,2X,I3,1X,I3,3X,F5.2,5X,E9.4,3X,E10.4,2X,E10.4,2X,E9.4,3X,
+ E9.4)
C Effective file loading
open(unit=1,status='old',file=
+'table6.dat')
write(6,*) '....Loading file: table6.dat'
do i__=1,31728
read(1,'(A81)')ar__3
read(ar__3,4)
+ El_3(i__),A_3(i__),Z_3(i__),T9_3(i__),pf_3(i__),nf_3(i__),
+ nf__3(i__),ng_3(i__),ng__3(i__)
c ..............Just test output...........
write(6,4)
+ El_3(i__),A_3(i__),Z_3(i__),T9_3(i__),pf_3(i__),nf_3(i__),
+ nf__3(i__),ng_3(i__),ng__3(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table7.dat' ! Coefficients for (n,g)-rate fit, ETFSI
C Format for file interpretation
5 format(
+ A2,2X,I3,2X,I3,2X,I1,2X,E13.7,1X,E13.7,1X,E13.7,1X,E13.7,1X,
+ E13.7,1X,E13.7,1X,E13.7,3X,E7.2)
C Effective file loading
open(unit=1,status='old',file=
+'table7.dat')
write(6,*) '....Loading file: table7.dat'
do i__=1,2662
read(1,'(A124)')ar__4
read(ar__4,5)
+ El_4(i__),A_4(i__),Z_4(i__),i(i__),a0(i__),a1(i__),a2(i__),
+ a3(i__),a4(i__),a5(i__),a6(i__),zeta(i__)
if(ar__4(118:124) .EQ. '') zeta(i__) = rNULL__
c ..............Just test output...........
write(6,5)
+ El_4(i__),A_4(i__),Z_4(i__),i(i__),a0(i__),a1(i__),a2(i__),
+ a3(i__),a4(i__),a5(i__),a6(i__),zeta(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table8.dat' ! Coefficients for (n,g)-rate fit, TF
C Format for file interpretation
6 format(
+ A2,2X,I3,2X,I3,2X,I1,2X,E13.7,1X,E13.7,1X,E13.7,1X,E13.7,1X,
+ E13.7,1X,E13.7,1X,E13.7,3X,E7.2)
C Effective file loading
open(unit=1,status='old',file=
+'table8.dat')
write(6,*) '....Loading file: table8.dat'
do i__=1,2984
read(1,'(A124)')ar__5
read(ar__5,6)
+ El_5(i__),A_5(i__),Z_5(i__),i_1(i__),a0_1(i__),a1_1(i__),
+ a2_1(i__),a3_1(i__),a4_1(i__),a5_1(i__),a6_1(i__),zeta_1(i__)
if(ar__5(118:124) .EQ. '') zeta_1(i__) = rNULL__
c ..............Just test output...........
write(6,6)
+ El_5(i__),A_5(i__),Z_5(i__),i_1(i__),a0_1(i__),a1_1(i__),
+ a2_1(i__),a3_1(i__),a4_1(i__),a5_1(i__),a6_1(i__),zeta_1(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table9.dat' ! Coefficients for (n,g)-rate fit, FRDM+TF
C Format for file interpretation
7 format(
+ A2,2X,I3,2X,I3,2X,I1,2X,E13.7,1X,E13.7,1X,E13.7,1X,E13.7,1X,
+ E13.7,1X,E13.7,1X,E13.7,3X,E7.2)
C Effective file loading
open(unit=1,status='old',file=
+'table9.dat')
write(6,*) '....Loading file: table9.dat'
do i__=1,2911
read(1,'(A124)')ar__6
read(ar__6,7)
+ El_6(i__),A_6(i__),Z_6(i__),i_2(i__),a0_2(i__),a1_2(i__),
+ a2_2(i__),a3_2(i__),a4_2(i__),a5_2(i__),a6_2(i__),zeta_2(i__)
if(ar__6(118:124) .EQ. '') zeta_2(i__) = rNULL__
c ..............Just test output...........
write(6,7)
+ El_6(i__),A_6(i__),Z_6(i__),i_2(i__),a0_2(i__),a1_2(i__),
+ a2_2(i__),a3_2(i__),a4_2(i__),a5_2(i__),a6_2(i__),zeta_2(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table10.dat' ! Coefficients for (n,g)-rate fit, HFB14
C Format for file interpretation
8 format(
+ A2,2X,I3,2X,I3,2X,I1,2X,E13.7,1X,E13.7,1X,E13.7,1X,E13.7,1X,
+ E13.7,1X,E13.7,1X,E13.7,3X,E7.2)
C Effective file loading
open(unit=1,status='old',file=
+'table10.dat')
write(6,*) '....Loading file: table10.dat'
do i__=1,1468
read(1,'(A124)')ar__7
read(ar__7,8)
+ El_7(i__),A_7(i__),Z_7(i__),i_3(i__),a0_3(i__),a1_3(i__),
+ a2_3(i__),a3_3(i__),a4_3(i__),a5_3(i__),a6_3(i__),zeta_3(i__)
if(ar__7(118:124) .EQ. '') zeta_3(i__) = rNULL__
c ..............Just test output...........
write(6,8)
+ El_7(i__),A_7(i__),Z_7(i__),i_3(i__),a0_3(i__),a1_3(i__),
+ a2_3(i__),a3_3(i__),a4_3(i__),a5_3(i__),a6_3(i__),zeta_3(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table11.dat' ! Coefficients for (g,n)-rate fit, ETFSI
C Format for file interpretation
9 format(
+ A2,2X,I3,2X,I3,2X,I1,2X,E13.7,1X,E13.7,1X,E13.7,1X,E13.7,1X,
+ E13.7,1X,E13.7,1X,E13.7,3X,E7.2)
C Effective file loading
open(unit=1,status='old',file=
+'table11.dat')
write(6,*) '....Loading file: table11.dat'
do i__=1,2414
read(1,'(A124)')ar__8
read(ar__8,9)
+ El_8(i__),A_8(i__),Z_8(i__),i_4(i__),a0_4(i__),a1_4(i__),
+ a2_4(i__),a3_4(i__),a4_4(i__),a5_4(i__),a6_4(i__),zeta_4(i__)
if(ar__8(118:124) .EQ. '') zeta_4(i__) = rNULL__
c ..............Just test output...........
write(6,9)
+ El_8(i__),A_8(i__),Z_8(i__),i_4(i__),a0_4(i__),a1_4(i__),
+ a2_4(i__),a3_4(i__),a4_4(i__),a5_4(i__),a6_4(i__),zeta_4(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table12.dat' ! Coefficients for (g,n)-rate fit, TF
C Format for file interpretation
10 format(
+ A2,2X,I3,2X,I3,2X,I1,2X,E13.7,1X,E13.7,1X,E13.7,1X,E13.7,1X,
+ E13.7,1X,E13.7,1X,E13.7,3X,E7.2)
C Effective file loading
open(unit=1,status='old',file=
+'table12.dat')
write(6,*) '....Loading file: table12.dat'
do i__=1,2948
read(1,'(A124)')ar__9
read(ar__9,10)
+ El_9(i__),A_9(i__),Z_9(i__),i_5(i__),a0_5(i__),a1_5(i__),
+ a2_5(i__),a3_5(i__),a4_5(i__),a5_5(i__),a6_5(i__),zeta_5(i__)
if(ar__9(118:124) .EQ. '') zeta_5(i__) = rNULL__
c ..............Just test output...........
write(6,10)
+ El_9(i__),A_9(i__),Z_9(i__),i_5(i__),a0_5(i__),a1_5(i__),
+ a2_5(i__),a3_5(i__),a4_5(i__),a5_5(i__),a6_5(i__),zeta_5(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table13.dat' ! Coefficients for (g,n)-rate fit, FRDM+TF
C Format for file interpretation
11 format(
+ A2,2X,I3,2X,I3,2X,I1,2X,E13.7,1X,E13.7,1X,E13.7,1X,E13.7,1X,
+ E13.7,1X,E13.7,1X,E13.7,3X,E7.2)
C Effective file loading
open(unit=1,status='old',file=
+'table13.dat')
write(6,*) '....Loading file: table13.dat'
do i__=1,2875
read(1,'(A124)')ar__10
read(ar__10,11)
+ El_10(i__),A_10(i__),Z_10(i__),i_6(i__),a0_6(i__),a1_6(i__),
+ a2_6(i__),a3_6(i__),a4_6(i__),a5_6(i__),a6_6(i__),zeta_6(i__)
if(ar__10(118:124) .EQ. '') zeta_6(i__) = rNULL__
c ..............Just test output...........
write(6,11)
+ El_10(i__),A_10(i__),Z_10(i__),i_6(i__),a0_6(i__),a1_6(i__),
+ a2_6(i__),a3_6(i__),a4_6(i__),a5_6(i__),a6_6(i__),zeta_6(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table14.dat' ! Coefficients for (g,n)-rate fit, HFB14
C Format for file interpretation
12 format(
+ A2,2X,I3,2X,I3,2X,I1,2X,E13.7,1X,E13.7,1X,E13.7,1X,E13.7,1X,
+ E13.7,1X,E13.7,1X,E13.7,3X,E7.2)
C Effective file loading
open(unit=1,status='old',file=
+'table14.dat')
write(6,*) '....Loading file: table14.dat'
do i__=1,1436
read(1,'(A124)')ar__11
read(ar__11,12)
+ El_11(i__),A_11(i__),Z_11(i__),i_7(i__),a0_7(i__),a1_7(i__),
+ a2_7(i__),a3_7(i__),a4_7(i__),a5_7(i__),a6_7(i__),zeta_7(i__)
if(ar__11(118:124) .EQ. '') zeta_7(i__) = rNULL__
c ..............Just test output...........
write(6,12)
+ El_11(i__),A_11(i__),Z_11(i__),i_7(i__),a0_7(i__),a1_7(i__),
+ a2_7(i__),a3_7(i__),a4_7(i__),a5_7(i__),a6_7(i__),zeta_7(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table15.dat' ! Coefficients for (n,f)-rate fit, ETFSI
C Format for file interpretation
13 format(
+ A2,2X,I3,2X,I3,2X,I1,2X,E13.7,1X,E13.7,1X,E13.7,1X,E13.7,1X,
+ E13.7,1X,E13.7,1X,E13.7,3X,E7.2)
C Effective file loading
open(unit=1,status='old',file=
+'table15.dat')
write(6,*) '....Loading file: table15.dat'
do i__=1,3288
read(1,'(A124)')ar__12
read(ar__12,13)
+ El_12(i__),A_12(i__),Z_12(i__),i_8(i__),a0_8(i__),a1_8(i__),
+ a2_8(i__),a3_8(i__),a4_8(i__),a5_8(i__),a6_8(i__),zeta_8(i__)
if(ar__12(118:124) .EQ. '') zeta_8(i__) = rNULL__
c ..............Just test output...........
write(6,13)
+ El_12(i__),A_12(i__),Z_12(i__),i_8(i__),a0_8(i__),a1_8(i__),
+ a2_8(i__),a3_8(i__),a4_8(i__),a5_8(i__),a6_8(i__),zeta_8(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table16.dat' ! Coefficients for (n,f)-rate fit, TF
C Format for file interpretation
14 format(
+ A2,2X,I3,2X,I3,2X,I1,2X,E13.7,1X,E13.7,1X,E13.7,1X,E13.7,1X,
+ E13.7,1X,E13.7,1X,E13.7,3X,E7.2)
C Effective file loading
open(unit=1,status='old',file=
+'table16.dat')
write(6,*) '....Loading file: table16.dat'
do i__=1,3928
read(1,'(A124)')ar__13
read(ar__13,14)
+ El_13(i__),A_13(i__),Z_13(i__),i_9(i__),a0_9(i__),a1_9(i__),
+ a2_9(i__),a3_9(i__),a4_9(i__),a5_9(i__),a6_9(i__),zeta_9(i__)
if(ar__13(118:124) .EQ. '') zeta_9(i__) = rNULL__
c ..............Just test output...........
write(6,14)
+ El_13(i__),A_13(i__),Z_13(i__),i_9(i__),a0_9(i__),a1_9(i__),
+ a2_9(i__),a3_9(i__),a4_9(i__),a5_9(i__),a6_9(i__),zeta_9(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table17.dat' ! Coefficients for (n,f)-rate fit, FRDM+TF
C Format for file interpretation
15 format(
+ A2,2X,I3,2X,I3,2X,I1,2X,E13.7,1X,E13.7,1X,E13.7,1X,E13.7,1X,
+ E13.7,1X,E13.7,1X,E13.7,3X,E7.2)
C Effective file loading
open(unit=1,status='old',file=
+'table17.dat')
write(6,*) '....Loading file: table17.dat'
do i__=1,3846
read(1,'(A124)')ar__14
read(ar__14,15)
+ El_14(i__),A_14(i__),Z_14(i__),i_10(i__),a0_10(i__),
+ a1_10(i__),a2_10(i__),a3_10(i__),a4_10(i__),a5_10(i__),
+ a6_10(i__),zeta_10(i__)
if(ar__14(118:124) .EQ. '') zeta_10(i__) = rNULL__
c ..............Just test output...........
write(6,15)
+ El_14(i__),A_14(i__),Z_14(i__),i_10(i__),a0_10(i__),
+ a1_10(i__),a2_10(i__),a3_10(i__),a4_10(i__),a5_10(i__),
+ a6_10(i__),zeta_10(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
C Loading file 'table18.dat' ! Coefficients for (n,f)-rate fit, HFB14
C Format for file interpretation
16 format(
+ A2,2X,I3,2X,I3,2X,I1,2X,E13.7,1X,E13.7,1X,E13.7,1X,E13.7,1X,
+ E13.7,1X,E13.7,1X,E13.7,3X,E7.2)
C Effective file loading
open(unit=1,status='old',file=
+'table18.dat')
write(6,*) '....Loading file: table18.dat'
do i__=1,1448
read(1,'(A124)')ar__15
read(ar__15,16)
+ El_15(i__),A_15(i__),Z_15(i__),i_11(i__),a0_11(i__),
+ a1_11(i__),a2_11(i__),a3_11(i__),a4_11(i__),a5_11(i__),
+ a6_11(i__),zeta_11(i__)
if(ar__15(118:124) .EQ. '') zeta_11(i__) = rNULL__
c ..............Just test output...........
write(6,16)
+ El_15(i__),A_15(i__),Z_15(i__),i_11(i__),a0_11(i__),
+ a1_11(i__),a2_11(i__),a3_11(i__),a4_11(i__),a5_11(i__),
+ a6_11(i__),zeta_11(i__)
c .......End.of.Just test output...........
end do
close(1)
C=============================================================================
stop
end