LCOV - code coverage report
Current view: top level - chemistry/pp_trop_strat_mam5_ts2 - mo_setrxt.F90 (source / functions) Hit Total Coverage
Test: coverage.info Lines: 917 946 96.9 %
Date: 2025-03-14 01:18:36 Functions: 1 2 50.0 %

          Line data    Source code
       1             : 
       2             :       module mo_setrxt
       3             : 
       4             :       use shr_kind_mod, only : r8 => shr_kind_r8
       5             : 
       6             :       private
       7             :       public :: setrxt
       8             :       public :: setrxt_hrates
       9             : 
      10             :       contains
      11             : 
      12       72960 :       subroutine setrxt( rate, temp, m, ncol )
      13             :  
      14             :       use ppgrid, only : pcols, pver
      15             : 
      16             : 
      17             :       use chem_mods, only : rxntot
      18             :       use mo_jpl,    only : jpl
      19             : 
      20             :       implicit none
      21             : 
      22             : !-------------------------------------------------------
      23             : !       ... dummy arguments
      24             : !-------------------------------------------------------
      25             :       integer, intent(in) :: ncol
      26             :       real(r8), intent(in)    :: temp(pcols,pver)
      27             :       real(r8), intent(in)    :: m(ncol*pver)
      28             :       real(r8), intent(inout) :: rate(ncol*pver,max(1,rxntot))
      29             : 
      30             : !-------------------------------------------------------
      31             : !       ... local variables
      32             : !-------------------------------------------------------
      33             :       integer   ::  n
      34             :       integer   ::  offset
      35      145920 :       real(r8)  :: itemp(ncol*pver)
      36      145920 :       real(r8)  :: exp_fac(ncol*pver)
      37      145920 :       real(r8)  :: ko(ncol*pver)
      38      145920 :       real(r8)  :: kinf(ncol*pver)
      39             : 
      40    58907904 :       rate(:,168) = 1.2e-10_r8
      41    58907904 :       rate(:,172) = 1.2e-10_r8
      42    58907904 :       rate(:,178) = 6.9e-12_r8
      43    58907904 :       rate(:,179) = 7.2e-11_r8
      44    58907904 :       rate(:,180) = 1.6e-12_r8
      45    58907904 :       rate(:,186) = 1.8e-12_r8
      46    58907904 :       rate(:,190) = 1.8e-12_r8
      47    58907904 :       rate(:,202) = 3.5e-12_r8
      48    58907904 :       rate(:,204) = 1.3e-11_r8
      49    58907904 :       rate(:,205) = 2.2e-11_r8
      50    58907904 :       rate(:,206) = 5e-11_r8
      51    58907904 :       rate(:,241) = 1.7e-13_r8
      52    58907904 :       rate(:,243) = 2.607e-10_r8
      53    58907904 :       rate(:,244) = 9.75e-11_r8
      54    58907904 :       rate(:,245) = 2.07e-10_r8
      55    58907904 :       rate(:,246) = 2.088e-10_r8
      56    58907904 :       rate(:,247) = 1.17e-10_r8
      57    58907904 :       rate(:,248) = 4.644e-11_r8
      58    58907904 :       rate(:,249) = 1.204e-10_r8
      59    58907904 :       rate(:,250) = 9.9e-11_r8
      60    58907904 :       rate(:,251) = 3.3e-12_r8
      61    58907904 :       rate(:,270) = 4.5e-11_r8
      62    58907904 :       rate(:,271) = 4.62e-10_r8
      63    58907904 :       rate(:,272) = 1.2e-10_r8
      64    58907904 :       rate(:,273) = 9e-11_r8
      65    58907904 :       rate(:,274) = 3e-11_r8
      66    58907904 :       rate(:,279) = 2.14e-11_r8
      67    58907904 :       rate(:,280) = 1.9e-10_r8
      68    58907904 :       rate(:,293) = 2.57e-10_r8
      69    58907904 :       rate(:,294) = 1.8e-10_r8
      70    58907904 :       rate(:,295) = 1.794e-10_r8
      71    58907904 :       rate(:,296) = 1.3e-10_r8
      72    58907904 :       rate(:,297) = 7.65e-11_r8
      73    58907904 :       rate(:,310) = 4e-13_r8
      74    58907904 :       rate(:,315) = 1.31e-10_r8
      75    58907904 :       rate(:,316) = 3.5e-11_r8
      76    58907904 :       rate(:,317) = 9e-12_r8
      77    58907904 :       rate(:,324) = 6.8e-14_r8
      78    58907904 :       rate(:,325) = 2e-13_r8
      79    58907904 :       rate(:,340) = 1e-12_r8
      80    58907904 :       rate(:,344) = 1e-14_r8
      81    58907904 :       rate(:,345) = 1e-11_r8
      82    58907904 :       rate(:,346) = 1.15e-11_r8
      83    58907904 :       rate(:,347) = 3.3e-11_r8
      84    58907904 :       rate(:,348) = 3.4e-12_r8
      85    58907904 :       rate(:,349) = 4e-14_r8
      86    58907904 :       rate(:,362) = 3e-12_r8
      87    58907904 :       rate(:,363) = 1.2e-11_r8
      88    58907904 :       rate(:,364) = 6.7e-13_r8
      89    58907904 :       rate(:,374) = 3.5e-13_r8
      90    58907904 :       rate(:,375) = 5.4e-11_r8
      91    58907904 :       rate(:,376) = 3.77e-11_r8
      92    58907904 :       rate(:,379) = 2e-12_r8
      93    58907904 :       rate(:,380) = 1.29e-11_r8
      94    58907904 :       rate(:,382) = 4.5e-14_r8
      95    58907904 :       rate(:,387) = 3.77e-11_r8
      96    58907904 :       rate(:,393) = 4e-12_r8
      97    58907904 :       rate(:,399) = 1.78e-12_r8
      98    58907904 :       rate(:,401) = 6.1e-13_r8
      99    58907904 :       rate(:,405) = 4.8e-11_r8
     100    58907904 :       rate(:,408) = 1.6e-12_r8
     101    58907904 :       rate(:,410) = 6.7e-12_r8
     102    58907904 :       rate(:,413) = 3.5e-12_r8
     103    58907904 :       rate(:,418) = 6.42e-11_r8
     104    58907904 :       rate(:,425) = 1.6e-13_r8
     105    58907904 :       rate(:,431) = 1.4e-12_r8
     106    58907904 :       rate(:,436) = 7.5e-13_r8
     107    58907904 :       rate(:,437) = 1.4e-13_r8
     108    58907904 :       rate(:,438) = 7.5e-13_r8
     109    58907904 :       rate(:,439) = 3.6e-13_r8
     110    58907904 :       rate(:,440) = 6.5e-13_r8
     111    58907904 :       rate(:,441) = 2.1e-13_r8
     112    58907904 :       rate(:,442) = 6.5e-13_r8
     113    58907904 :       rate(:,443) = 4.9e-13_r8
     114    58907904 :       rate(:,445) = 1.2e-12_r8
     115    58907904 :       rate(:,449) = 9.8e-13_r8
     116    58907904 :       rate(:,452) = 1.85e-11_r8
     117    58907904 :       rate(:,453) = 1.63e-12_r8
     118    58907904 :       rate(:,454) = 2.5e-11_r8
     119    58907904 :       rate(:,455) = 1.1e-11_r8
     120    58907904 :       rate(:,456) = 3.3e-11_r8
     121    58907904 :       rate(:,459) = 2.8e-17_r8
     122    58907904 :       rate(:,460) = 8e-11_r8
     123    58907904 :       rate(:,463) = 3e-11_r8
     124    58907904 :       rate(:,466) = 4.2e-11_r8
     125    58907904 :       rate(:,469) = 2.8e-17_r8
     126    58907904 :       rate(:,470) = 1.1e-10_r8
     127    58907904 :       rate(:,472) = 3.9e-11_r8
     128    58907904 :       rate(:,475) = 1.3e-12_r8
     129    58907904 :       rate(:,477) = 5e-12_r8
     130    58907904 :       rate(:,478) = 2.3e-12_r8
     131    58907904 :       rate(:,481) = 3.9e-11_r8
     132    58907904 :       rate(:,484) = 2.8e-17_r8
     133    58907904 :       rate(:,485) = 9.2e-11_r8
     134    58907904 :       rate(:,488) = 3.85e-11_r8
     135    58907904 :       rate(:,492) = 1.2e-12_r8
     136    58907904 :       rate(:,496) = 9.8e-13_r8
     137    58907904 :       rate(:,501) = 4.4e-18_r8
     138    58907904 :       rate(:,502) = 3.6e-11_r8
     139    58907904 :       rate(:,554) = 4.7e-11_r8
     140    58907904 :       rate(:,567) = 2.1e-12_r8
     141    58907904 :       rate(:,568) = 2.8e-13_r8
     142    58907904 :       rate(:,576) = 1.7e-11_r8
     143    58907904 :       rate(:,582) = 8.4e-11_r8
     144    58907904 :       rate(:,585) = 5.3e-13_r8
     145    58907904 :       rate(:,587) = 2e-12_r8
     146    58907904 :       rate(:,590) = 2.3e-12_r8
     147    58907904 :       rate(:,595) = 2e-12_r8
     148    58907904 :       rate(:,598) = 2.3e-12_r8
     149    58907904 :       rate(:,604) = 1.9e-11_r8
     150    58907904 :       rate(:,605) = 5.3e-13_r8
     151    58907904 :       rate(:,607) = 2e-12_r8
     152    58907904 :       rate(:,610) = 2.3e-12_r8
     153    58907904 :       rate(:,615) = 2e-12_r8
     154    58907904 :       rate(:,618) = 2.3e-12_r8
     155    58907904 :       rate(:,622) = 1.2e-14_r8
     156    58907904 :       rate(:,623) = 2e-10_r8
     157    58907904 :       rate(:,624) = 2.5e-12_r8
     158    58907904 :       rate(:,625) = 5.3e-13_r8
     159    58907904 :       rate(:,627) = 2e-12_r8
     160    58907904 :       rate(:,630) = 2.3e-12_r8
     161    58907904 :       rate(:,635) = 2e-12_r8
     162    58907904 :       rate(:,638) = 2.3e-12_r8
     163    58907904 :       rate(:,644) = 1.2e-11_r8
     164    58907904 :       rate(:,646) = 2e-12_r8
     165    58907904 :       rate(:,648) = 5.3e-13_r8
     166    58907904 :       rate(:,650) = 2.3e-12_r8
     167    58907904 :       rate(:,655) = 2e-12_r8
     168    58907904 :       rate(:,658) = 2.3e-12_r8
     169    58907904 :       rate(:,664) = 1.1e-11_r8
     170    58907904 :       rate(:,666) = 2e-12_r8
     171    58907904 :       rate(:,668) = 5.3e-13_r8
     172    58907904 :       rate(:,670) = 2.3e-12_r8
     173    58907904 :       rate(:,675) = 2e-12_r8
     174    58907904 :       rate(:,678) = 2.3e-12_r8
     175    58907904 :       rate(:,683) = 2.1e-10_r8
     176    58907904 :       rate(:,689) = 8.9e-11_r8
     177    58907904 :       rate(:,690) = 8.9e-11_r8
     178    58907904 :       rate(:,694) = 2e-12_r8
     179    58907904 :       rate(:,697) = 2.3e-12_r8
     180    58907904 :       rate(:,705) = 4e-12_r8
     181    58907904 :       rate(:,708) = 2e-14_r8
     182    58907904 :       rate(:,710) = 2e-12_r8
     183    58907904 :       rate(:,713) = 2.3e-12_r8
     184    58907904 :       rate(:,718) = 2.52e-11_r8
     185    58907904 :       rate(:,723) = 4e-12_r8
     186    58907904 :       rate(:,727) = 2e-14_r8
     187    58907904 :       rate(:,729) = 2e-12_r8
     188    58907904 :       rate(:,732) = 2.3e-12_r8
     189    58907904 :       rate(:,737) = 1.92e-11_r8
     190    58907904 :       rate(:,739) = 2e-12_r8
     191    58907904 :       rate(:,742) = 2.3e-12_r8
     192    58907904 :       rate(:,746) = 8.8e-12_r8
     193    58907904 :       rate(:,747) = 8.8e-12_r8
     194    58907904 :       rate(:,748) = 8.8e-12_r8
     195    58907904 :       rate(:,753) = 4e-12_r8
     196    58907904 :       rate(:,755) = 2e-14_r8
     197    58907904 :       rate(:,757) = 3.66e-12_r8
     198    58907904 :       rate(:,758) = 2.8e-11_r8
     199    58907904 :       rate(:,759) = 2.6e-13_r8
     200    58907904 :       rate(:,762) = 8.3e-18_r8
     201    58907904 :       rate(:,763) = 1.1e-10_r8
     202    58907904 :       rate(:,767) = 1.1e-16_r8
     203    58907904 :       rate(:,769) = 3.64e-12_r8
     204    58907904 :       rate(:,770) = 2.8e-11_r8
     205    58907904 :       rate(:,771) = 1.7e-11_r8
     206    58907904 :       rate(:,774) = 1.1e-10_r8
     207    58907904 :       rate(:,775) = 9.58e-12_r8
     208    58907904 :       rate(:,778) = 1.1e-10_r8
     209    58907904 :       rate(:,779) = 1.23e-11_r8
     210    58907904 :       rate(:,782) = 1.1e-10_r8
     211    58907904 :       rate(:,783) = 3.64e-12_r8
     212    58907904 :       rate(:,786) = 1.1e-10_r8
     213    58907904 :       rate(:,787) = 5.5e-12_r8
     214    58907904 :       rate(:,788) = 4.65e-11_r8
     215    58907904 :       rate(:,789) = 2.8e-11_r8
     216    58907904 :       rate(:,797) = 2.3e-12_r8
     217    58907904 :       rate(:,799) = 1.2e-11_r8
     218    58907904 :       rate(:,800) = 5.7e-11_r8
     219    58907904 :       rate(:,801) = 2.8e-11_r8
     220    58907904 :       rate(:,802) = 6.6e-11_r8
     221    58907904 :       rate(:,803) = 1.4e-11_r8
     222    58907904 :       rate(:,806) = 1.9e-12_r8
     223    58907904 :       rate(:,829) = 6.34e-08_r8
     224    58907904 :       rate(:,846) = 1.9e-11_r8
     225    58907904 :       rate(:,849) = 1.2e-14_r8
     226    58907904 :       rate(:,850) = 2e-10_r8
     227    58907904 :       rate(:,854) = 2.5e-12_r8
     228    58907904 :       rate(:,866) = 1.34e-11_r8
     229    58907904 :       rate(:,867) = 1.2e-11_r8
     230    58907904 :       rate(:,872) = 1.1e-11_r8
     231    58907904 :       rate(:,876) = 2.1e-10_r8
     232    58907904 :       rate(:,877) = 1.34e-11_r8
     233    58907904 :       rate(:,881) = 1.7e-11_r8
     234    58907904 :       rate(:,901) = 1.29e-07_r8
     235    58907904 :       rate(:,902) = 2.31e-07_r8
     236    58907904 :       rate(:,903) = 2.31e-06_r8
     237    58907904 :       rate(:,904) = 4.63e-07_r8
     238             :  
     239     4158720 :       do n = 1,pver
     240     4085760 :         offset = (n-1)*ncol
     241    62993664 :         itemp(offset+1:offset+ncol) = 1._r8 / temp(:ncol,n)
     242             :       end do
     243             : 
     244    58907904 :       rate(:,169) = 1.63e-10_r8 * exp( 60._r8 * itemp(:) )
     245    58907904 :       rate(:,170) = 2.15e-11_r8 * exp( 110._r8 * itemp(:) )
     246    58907904 :       rate(:,171) = 3.3e-11_r8 * exp( 55._r8 * itemp(:) )
     247    58907904 :       rate(:,173) = 8e-12_r8 * exp( -2060._r8 * itemp(:) )
     248    58907904 :       rate(:,176) = 1.6e-11_r8 * exp( -4570._r8 * itemp(:) )
     249    58907904 :       rate(:,177) = 1.4e-12_r8 * exp( -2000._r8 * itemp(:) )
     250    58907904 :       exp_fac(:) = exp( 200._r8 * itemp(:) )
     251    58907904 :       rate(:,182) = 3e-11_r8 * exp_fac(:)
     252    58907904 :       rate(:,268) = 5.5e-12_r8 * exp_fac(:)
     253    58907904 :       rate(:,307) = 3.8e-12_r8 * exp_fac(:)
     254    58907904 :       rate(:,329) = 3.8e-12_r8 * exp_fac(:)
     255    58907904 :       rate(:,358) = 3.8e-12_r8 * exp_fac(:)
     256    58907904 :       rate(:,367) = 3.8e-12_r8 * exp_fac(:)
     257    58907904 :       rate(:,371) = 3.8e-12_r8 * exp_fac(:)
     258    58907904 :       rate(:,397) = 3.8e-12_r8 * exp_fac(:)
     259    58907904 :       rate(:,412) = 3.8e-12_r8 * exp_fac(:)
     260    58907904 :       rate(:,489) = 5.53e-12_r8 * exp_fac(:)
     261    58907904 :       rate(:,546) = 3.8e-12_r8 * exp_fac(:)
     262    58907904 :       rate(:,549) = 3.8e-12_r8 * exp_fac(:)
     263    58907904 :       rate(:,553) = 3.8e-12_r8 * exp_fac(:)
     264    58907904 :       rate(:,569) = 3.8e-12_r8 * exp_fac(:)
     265    58907904 :       rate(:,573) = 3.8e-12_r8 * exp_fac(:)
     266    58907904 :       rate(:,579) = 3.8e-12_r8 * exp_fac(:)
     267    58907904 :       rate(:,583) = 3.8e-12_r8 * exp_fac(:)
     268    58907904 :       rate(:,183) = 1e-14_r8 * exp( -490._r8 * itemp(:) )
     269    58907904 :       rate(:,184) = 1.4e-10_r8 * exp( -470._r8 * itemp(:) )
     270    58907904 :       rate(:,185) = 2.8e-12_r8 * exp( -1800._r8 * itemp(:) )
     271    58907904 :       exp_fac(:) = exp( 250._r8 * itemp(:) )
     272    58907904 :       rate(:,187) = 4.8e-11_r8 * exp_fac(:)
     273    58907904 :       rate(:,266) = 1.7e-11_r8 * exp_fac(:)
     274    58907904 :       exp_fac(:) = exp( 180._r8 * itemp(:) )
     275    58907904 :       rate(:,188) = 1.8e-11_r8 * exp_fac(:)
     276    58907904 :       rate(:,342) = 4.2e-12_r8 * exp_fac(:)
     277    58907904 :       rate(:,357) = 4.2e-12_r8 * exp_fac(:)
     278    58907904 :       rate(:,366) = 4.2e-12_r8 * exp_fac(:)
     279    58907904 :       rate(:,395) = 4.2e-12_r8 * exp_fac(:)
     280    58907904 :       rate(:,189) = 1.7e-12_r8 * exp( -940._r8 * itemp(:) )
     281    58907904 :       rate(:,193) = 4.5e-13_r8 * exp( 610._r8 * itemp(:) )
     282    58907904 :       rate(:,194) = 2.1e-11_r8 * exp( 100._r8 * itemp(:) )
     283    58907904 :       exp_fac(:) = exp( 220._r8 * itemp(:) )
     284    58907904 :       rate(:,195) = 2.9e-12_r8 * exp_fac(:)
     285    58907904 :       rate(:,196) = 1.45e-12_r8 * exp_fac(:)
     286    58907904 :       rate(:,197) = 1.45e-12_r8 * exp_fac(:)
     287    58907904 :       rate(:,198) = 3.3e-12_r8 * exp( -3150._r8 * itemp(:) )
     288    58907904 :       rate(:,199) = 5.1e-12_r8 * exp( 210._r8 * itemp(:) )
     289    58907904 :       exp_fac(:) = exp( -2450._r8 * itemp(:) )
     290    58907904 :       rate(:,200) = 1.2e-13_r8 * exp_fac(:)
     291    58907904 :       rate(:,226) = 3e-11_r8 * exp_fac(:)
     292    58907904 :       exp_fac(:) = exp( 125._r8 * itemp(:) )
     293    58907904 :       rate(:,203) = 1.7e-11_r8 * exp_fac(:)
     294    58907904 :       rate(:,301) = 5.5e-12_r8 * exp_fac(:)
     295    58907904 :       exp_fac(:) = exp( 260._r8 * itemp(:) )
     296    58907904 :       rate(:,207) = 3.44e-12_r8 * exp_fac(:)
     297    58907904 :       rate(:,259) = 2.3e-12_r8 * exp_fac(:)
     298    58907904 :       rate(:,262) = 8.8e-12_r8 * exp_fac(:)
     299    58907904 :       exp_fac(:) = exp( -1500._r8 * itemp(:) )
     300    58907904 :       rate(:,208) = 3e-12_r8 * exp_fac(:)
     301    58907904 :       rate(:,267) = 5.8e-12_r8 * exp_fac(:)
     302    58907904 :       exp_fac(:) = exp( 20._r8 * itemp(:) )
     303    58907904 :       rate(:,210) = 7.26e-11_r8 * exp_fac(:)
     304    58907904 :       rate(:,211) = 4.64e-11_r8 * exp_fac(:)
     305    58907904 :       rate(:,218) = 8.1e-11_r8 * exp( -30._r8 * itemp(:) )
     306    58907904 :       exp_fac(:) = exp( -1270._r8 * itemp(:) )
     307    58907904 :       rate(:,219) = 7.1e-12_r8 * exp_fac(:)
     308    58907904 :       rate(:,642) = 1.35e-15_r8 * exp_fac(:)
     309    58907904 :       rate(:,857) = 1.35e-15_r8 * exp_fac(:)
     310    58907904 :       rate(:,220) = 3.05e-11_r8 * exp( -2270._r8 * itemp(:) )
     311    58907904 :       rate(:,221) = 1.1e-11_r8 * exp( -980._r8 * itemp(:) )
     312    58907904 :       exp_fac(:) = exp( 270._r8 * itemp(:) )
     313    58907904 :       rate(:,222) = 1.4e-11_r8 * exp_fac(:)
     314    58907904 :       rate(:,236) = 7.4e-12_r8 * exp_fac(:)
     315    58907904 :       rate(:,338) = 8.1e-12_r8 * exp_fac(:)
     316    58907904 :       rate(:,392) = 8.1e-12_r8 * exp_fac(:)
     317    58907904 :       rate(:,704) = 8.1e-12_r8 * exp_fac(:)
     318    58907904 :       rate(:,722) = 8.1e-12_r8 * exp_fac(:)
     319    58907904 :       rate(:,752) = 8.1e-12_r8 * exp_fac(:)
     320    58907904 :       rate(:,223) = 3.6e-11_r8 * exp( -375._r8 * itemp(:) )
     321    58907904 :       rate(:,224) = 2.3e-11_r8 * exp( -200._r8 * itemp(:) )
     322    58907904 :       rate(:,225) = 3.3e-12_r8 * exp( -115._r8 * itemp(:) )
     323    58907904 :       rate(:,227) = 1e-12_r8 * exp( -1590._r8 * itemp(:) )
     324    58907904 :       rate(:,228) = 3.5e-13_r8 * exp( -1370._r8 * itemp(:) )
     325    58907904 :       exp_fac(:) = exp( 290._r8 * itemp(:) )
     326    58907904 :       rate(:,229) = 2.6e-12_r8 * exp_fac(:)
     327    58907904 :       rate(:,230) = 6.4e-12_r8 * exp_fac(:)
     328    58907904 :       rate(:,260) = 4.1e-13_r8 * exp_fac(:)
     329    58907904 :       rate(:,542) = 7.5e-12_r8 * exp_fac(:)
     330    58907904 :       rate(:,556) = 7.5e-12_r8 * exp_fac(:)
     331    58907904 :       rate(:,559) = 7.5e-12_r8 * exp_fac(:)
     332    58907904 :       rate(:,562) = 7.5e-12_r8 * exp_fac(:)
     333    58907904 :       rate(:,231) = 6.5e-12_r8 * exp( 135._r8 * itemp(:) )
     334    58907904 :       exp_fac(:) = exp( -840._r8 * itemp(:) )
     335    58907904 :       rate(:,233) = 3.6e-12_r8 * exp_fac(:)
     336    58907904 :       rate(:,282) = 2e-12_r8 * exp_fac(:)
     337    58907904 :       rate(:,234) = 1.2e-12_r8 * exp( -330._r8 * itemp(:) )
     338    58907904 :       rate(:,235) = 2.8e-11_r8 * exp( 85._r8 * itemp(:) )
     339    58907904 :       exp_fac(:) = exp( 230._r8 * itemp(:) )
     340    58907904 :       rate(:,237) = 6e-13_r8 * exp_fac(:)
     341    58907904 :       rate(:,257) = 1.5e-12_r8 * exp_fac(:)
     342    58907904 :       rate(:,265) = 1.9e-11_r8 * exp_fac(:)
     343    58907904 :       rate(:,238) = 1e-11_r8 * exp( -3300._r8 * itemp(:) )
     344    58907904 :       rate(:,239) = 1.8e-12_r8 * exp( -250._r8 * itemp(:) )
     345    58907904 :       rate(:,240) = 3.4e-12_r8 * exp( -130._r8 * itemp(:) )
     346    58907904 :       exp_fac(:) = exp( -500._r8 * itemp(:) )
     347    58907904 :       rate(:,242) = 3e-12_r8 * exp_fac(:)
     348    58907904 :       rate(:,276) = 1.4e-10_r8 * exp_fac(:)
     349    58907904 :       exp_fac(:) = exp( -800._r8 * itemp(:) )
     350    58907904 :       rate(:,254) = 1.7e-11_r8 * exp_fac(:)
     351    58907904 :       rate(:,281) = 6.3e-12_r8 * exp_fac(:)
     352    58907904 :       rate(:,255) = 4.8e-12_r8 * exp( -310._r8 * itemp(:) )
     353    58907904 :       rate(:,256) = 1.6e-11_r8 * exp( -780._r8 * itemp(:) )
     354    58907904 :       rate(:,258) = 9.5e-13_r8 * exp( 550._r8 * itemp(:) )
     355    58907904 :       exp_fac(:) = exp( 460._r8 * itemp(:) )
     356    58907904 :       rate(:,261) = 4.5e-12_r8 * exp_fac(:)
     357    58907904 :       rate(:,643) = 1.62e-11_r8 * exp_fac(:)
     358    58907904 :       rate(:,858) = 1.62e-11_r8 * exp_fac(:)
     359    58907904 :       rate(:,264) = 1.9e-11_r8 * exp( 215._r8 * itemp(:) )
     360    58907904 :       rate(:,269) = 1.2e-10_r8 * exp( -430._r8 * itemp(:) )
     361    58907904 :       rate(:,275) = 1.6e-10_r8 * exp( -260._r8 * itemp(:) )
     362    58907904 :       exp_fac(:) = exp( 0._r8 * itemp(:) )
     363    58907904 :       rate(:,277) = 1.4e-11_r8 * exp_fac(:)
     364    58907904 :       rate(:,279) = 2.14e-11_r8 * exp_fac(:)
     365    58907904 :       rate(:,280) = 1.9e-10_r8 * exp_fac(:)
     366    58907904 :       rate(:,293) = 2.57e-10_r8 * exp_fac(:)
     367    58907904 :       rate(:,294) = 1.8e-10_r8 * exp_fac(:)
     368    58907904 :       rate(:,295) = 1.794e-10_r8 * exp_fac(:)
     369    58907904 :       rate(:,296) = 1.3e-10_r8 * exp_fac(:)
     370    58907904 :       rate(:,297) = 7.65e-11_r8 * exp_fac(:)
     371    58907904 :       rate(:,310) = 4e-13_r8 * exp_fac(:)
     372    58907904 :       rate(:,315) = 1.31e-10_r8 * exp_fac(:)
     373    58907904 :       rate(:,316) = 3.5e-11_r8 * exp_fac(:)
     374    58907904 :       rate(:,317) = 9e-12_r8 * exp_fac(:)
     375    58907904 :       rate(:,324) = 6.8e-14_r8 * exp_fac(:)
     376    58907904 :       rate(:,325) = 2e-13_r8 * exp_fac(:)
     377    58907904 :       rate(:,340) = 1e-12_r8 * exp_fac(:)
     378    58907904 :       rate(:,344) = 1e-14_r8 * exp_fac(:)
     379    58907904 :       rate(:,345) = 1e-11_r8 * exp_fac(:)
     380    58907904 :       rate(:,346) = 1.15e-11_r8 * exp_fac(:)
     381    58907904 :       rate(:,347) = 3.3e-11_r8 * exp_fac(:)
     382    58907904 :       rate(:,348) = 3.4e-12_r8 * exp_fac(:)
     383    58907904 :       rate(:,349) = 4e-14_r8 * exp_fac(:)
     384    58907904 :       rate(:,362) = 3e-12_r8 * exp_fac(:)
     385    58907904 :       rate(:,363) = 1.2e-11_r8 * exp_fac(:)
     386    58907904 :       rate(:,364) = 6.7e-13_r8 * exp_fac(:)
     387    58907904 :       rate(:,374) = 3.5e-13_r8 * exp_fac(:)
     388    58907904 :       rate(:,375) = 5.4e-11_r8 * exp_fac(:)
     389    58907904 :       rate(:,376) = 3.77e-11_r8 * exp_fac(:)
     390    58907904 :       rate(:,379) = 2e-12_r8 * exp_fac(:)
     391    58907904 :       rate(:,380) = 1.29e-11_r8 * exp_fac(:)
     392    58907904 :       rate(:,382) = 4.5e-14_r8 * exp_fac(:)
     393    58907904 :       rate(:,387) = 3.77e-11_r8 * exp_fac(:)
     394    58907904 :       rate(:,393) = 4e-12_r8 * exp_fac(:)
     395    58907904 :       rate(:,399) = 1.78e-12_r8 * exp_fac(:)
     396    58907904 :       rate(:,401) = 6.1e-13_r8 * exp_fac(:)
     397    58907904 :       rate(:,405) = 4.8e-11_r8 * exp_fac(:)
     398    58907904 :       rate(:,408) = 1.6e-12_r8 * exp_fac(:)
     399    58907904 :       rate(:,410) = 6.7e-12_r8 * exp_fac(:)
     400    58907904 :       rate(:,413) = 3.5e-12_r8 * exp_fac(:)
     401    58907904 :       rate(:,418) = 6.42e-11_r8 * exp_fac(:)
     402    58907904 :       rate(:,425) = 1.6e-13_r8 * exp_fac(:)
     403    58907904 :       rate(:,431) = 1.4e-12_r8 * exp_fac(:)
     404    58907904 :       rate(:,436) = 7.5e-13_r8 * exp_fac(:)
     405    58907904 :       rate(:,437) = 1.4e-13_r8 * exp_fac(:)
     406    58907904 :       rate(:,438) = 7.5e-13_r8 * exp_fac(:)
     407    58907904 :       rate(:,439) = 3.6e-13_r8 * exp_fac(:)
     408    58907904 :       rate(:,440) = 6.5e-13_r8 * exp_fac(:)
     409    58907904 :       rate(:,441) = 2.1e-13_r8 * exp_fac(:)
     410    58907904 :       rate(:,442) = 6.5e-13_r8 * exp_fac(:)
     411    58907904 :       rate(:,443) = 4.9e-13_r8 * exp_fac(:)
     412    58907904 :       rate(:,445) = 1.2e-12_r8 * exp_fac(:)
     413    58907904 :       rate(:,449) = 9.8e-13_r8 * exp_fac(:)
     414    58907904 :       rate(:,452) = 1.85e-11_r8 * exp_fac(:)
     415    58907904 :       rate(:,453) = 1.63e-12_r8 * exp_fac(:)
     416    58907904 :       rate(:,454) = 2.5e-11_r8 * exp_fac(:)
     417    58907904 :       rate(:,455) = 1.1e-11_r8 * exp_fac(:)
     418    58907904 :       rate(:,456) = 3.3e-11_r8 * exp_fac(:)
     419    58907904 :       rate(:,459) = 2.8e-17_r8 * exp_fac(:)
     420    58907904 :       rate(:,460) = 8e-11_r8 * exp_fac(:)
     421    58907904 :       rate(:,463) = 3e-11_r8 * exp_fac(:)
     422    58907904 :       rate(:,466) = 4.2e-11_r8 * exp_fac(:)
     423    58907904 :       rate(:,469) = 2.8e-17_r8 * exp_fac(:)
     424    58907904 :       rate(:,470) = 1.1e-10_r8 * exp_fac(:)
     425    58907904 :       rate(:,472) = 3.9e-11_r8 * exp_fac(:)
     426    58907904 :       rate(:,475) = 1.3e-12_r8 * exp_fac(:)
     427    58907904 :       rate(:,477) = 5e-12_r8 * exp_fac(:)
     428    58907904 :       rate(:,478) = 2.3e-12_r8 * exp_fac(:)
     429    58907904 :       rate(:,481) = 3.9e-11_r8 * exp_fac(:)
     430    58907904 :       rate(:,484) = 2.8e-17_r8 * exp_fac(:)
     431    58907904 :       rate(:,485) = 9.2e-11_r8 * exp_fac(:)
     432    58907904 :       rate(:,488) = 3.85e-11_r8 * exp_fac(:)
     433    58907904 :       rate(:,492) = 1.2e-12_r8 * exp_fac(:)
     434    58907904 :       rate(:,496) = 9.8e-13_r8 * exp_fac(:)
     435    58907904 :       rate(:,501) = 4.4e-18_r8 * exp_fac(:)
     436    58907904 :       rate(:,502) = 3.6e-11_r8 * exp_fac(:)
     437    58907904 :       rate(:,554) = 4.7e-11_r8 * exp_fac(:)
     438    58907904 :       rate(:,567) = 2.1e-12_r8 * exp_fac(:)
     439    58907904 :       rate(:,568) = 2.8e-13_r8 * exp_fac(:)
     440    58907904 :       rate(:,576) = 1.7e-11_r8 * exp_fac(:)
     441    58907904 :       rate(:,582) = 8.4e-11_r8 * exp_fac(:)
     442    58907904 :       rate(:,585) = 5.3e-13_r8 * exp_fac(:)
     443    58907904 :       rate(:,587) = 2e-12_r8 * exp_fac(:)
     444    58907904 :       rate(:,590) = 2.3e-12_r8 * exp_fac(:)
     445    58907904 :       rate(:,595) = 2e-12_r8 * exp_fac(:)
     446    58907904 :       rate(:,598) = 2.3e-12_r8 * exp_fac(:)
     447    58907904 :       rate(:,604) = 1.9e-11_r8 * exp_fac(:)
     448    58907904 :       rate(:,605) = 5.3e-13_r8 * exp_fac(:)
     449    58907904 :       rate(:,607) = 2e-12_r8 * exp_fac(:)
     450    58907904 :       rate(:,610) = 2.3e-12_r8 * exp_fac(:)
     451    58907904 :       rate(:,615) = 2e-12_r8 * exp_fac(:)
     452    58907904 :       rate(:,618) = 2.3e-12_r8 * exp_fac(:)
     453    58907904 :       rate(:,622) = 1.2e-14_r8 * exp_fac(:)
     454    58907904 :       rate(:,623) = 2e-10_r8 * exp_fac(:)
     455    58907904 :       rate(:,624) = 2.5e-12_r8 * exp_fac(:)
     456    58907904 :       rate(:,625) = 5.3e-13_r8 * exp_fac(:)
     457    58907904 :       rate(:,627) = 2e-12_r8 * exp_fac(:)
     458    58907904 :       rate(:,630) = 2.3e-12_r8 * exp_fac(:)
     459    58907904 :       rate(:,635) = 2e-12_r8 * exp_fac(:)
     460    58907904 :       rate(:,638) = 2.3e-12_r8 * exp_fac(:)
     461    58907904 :       rate(:,644) = 1.2e-11_r8 * exp_fac(:)
     462    58907904 :       rate(:,646) = 2e-12_r8 * exp_fac(:)
     463    58907904 :       rate(:,648) = 5.3e-13_r8 * exp_fac(:)
     464    58907904 :       rate(:,650) = 2.3e-12_r8 * exp_fac(:)
     465    58907904 :       rate(:,655) = 2e-12_r8 * exp_fac(:)
     466    58907904 :       rate(:,658) = 2.3e-12_r8 * exp_fac(:)
     467    58907904 :       rate(:,664) = 1.1e-11_r8 * exp_fac(:)
     468    58907904 :       rate(:,666) = 2e-12_r8 * exp_fac(:)
     469    58907904 :       rate(:,668) = 5.3e-13_r8 * exp_fac(:)
     470    58907904 :       rate(:,670) = 2.3e-12_r8 * exp_fac(:)
     471    58907904 :       rate(:,675) = 2e-12_r8 * exp_fac(:)
     472    58907904 :       rate(:,678) = 2.3e-12_r8 * exp_fac(:)
     473    58907904 :       rate(:,683) = 2.1e-10_r8 * exp_fac(:)
     474    58907904 :       rate(:,689) = 8.9e-11_r8 * exp_fac(:)
     475    58907904 :       rate(:,690) = 8.9e-11_r8 * exp_fac(:)
     476    58907904 :       rate(:,694) = 2e-12_r8 * exp_fac(:)
     477    58907904 :       rate(:,697) = 2.3e-12_r8 * exp_fac(:)
     478    58907904 :       rate(:,705) = 4e-12_r8 * exp_fac(:)
     479    58907904 :       rate(:,708) = 2e-14_r8 * exp_fac(:)
     480    58907904 :       rate(:,710) = 2e-12_r8 * exp_fac(:)
     481    58907904 :       rate(:,713) = 2.3e-12_r8 * exp_fac(:)
     482    58907904 :       rate(:,718) = 2.52e-11_r8 * exp_fac(:)
     483    58907904 :       rate(:,723) = 4e-12_r8 * exp_fac(:)
     484    58907904 :       rate(:,727) = 2e-14_r8 * exp_fac(:)
     485    58907904 :       rate(:,729) = 2e-12_r8 * exp_fac(:)
     486    58907904 :       rate(:,732) = 2.3e-12_r8 * exp_fac(:)
     487    58907904 :       rate(:,737) = 1.92e-11_r8 * exp_fac(:)
     488    58907904 :       rate(:,739) = 2e-12_r8 * exp_fac(:)
     489    58907904 :       rate(:,742) = 2.3e-12_r8 * exp_fac(:)
     490    58907904 :       rate(:,746) = 8.8e-12_r8 * exp_fac(:)
     491    58907904 :       rate(:,747) = 8.8e-12_r8 * exp_fac(:)
     492    58907904 :       rate(:,748) = 8.8e-12_r8 * exp_fac(:)
     493    58907904 :       rate(:,753) = 4e-12_r8 * exp_fac(:)
     494    58907904 :       rate(:,755) = 2e-14_r8 * exp_fac(:)
     495    58907904 :       rate(:,757) = 3.66e-12_r8 * exp_fac(:)
     496    58907904 :       rate(:,758) = 2.8e-11_r8 * exp_fac(:)
     497    58907904 :       rate(:,759) = 2.6e-13_r8 * exp_fac(:)
     498    58907904 :       rate(:,762) = 8.3e-18_r8 * exp_fac(:)
     499    58907904 :       rate(:,763) = 1.1e-10_r8 * exp_fac(:)
     500    58907904 :       rate(:,767) = 1.1e-16_r8 * exp_fac(:)
     501    58907904 :       rate(:,769) = 3.64e-12_r8 * exp_fac(:)
     502    58907904 :       rate(:,770) = 2.8e-11_r8 * exp_fac(:)
     503    58907904 :       rate(:,771) = 1.7e-11_r8 * exp_fac(:)
     504    58907904 :       rate(:,774) = 1.1e-10_r8 * exp_fac(:)
     505    58907904 :       rate(:,775) = 9.58e-12_r8 * exp_fac(:)
     506    58907904 :       rate(:,778) = 1.1e-10_r8 * exp_fac(:)
     507    58907904 :       rate(:,779) = 1.23e-11_r8 * exp_fac(:)
     508    58907904 :       rate(:,782) = 1.1e-10_r8 * exp_fac(:)
     509    58907904 :       rate(:,783) = 3.64e-12_r8 * exp_fac(:)
     510    58907904 :       rate(:,786) = 1.1e-10_r8 * exp_fac(:)
     511    58907904 :       rate(:,787) = 5.5e-12_r8 * exp_fac(:)
     512    58907904 :       rate(:,788) = 4.65e-11_r8 * exp_fac(:)
     513    58907904 :       rate(:,789) = 2.8e-11_r8 * exp_fac(:)
     514    58907904 :       rate(:,797) = 2.3e-12_r8 * exp_fac(:)
     515    58907904 :       rate(:,799) = 1.2e-11_r8 * exp_fac(:)
     516    58907904 :       rate(:,800) = 5.7e-11_r8 * exp_fac(:)
     517    58907904 :       rate(:,801) = 2.8e-11_r8 * exp_fac(:)
     518    58907904 :       rate(:,802) = 6.6e-11_r8 * exp_fac(:)
     519    58907904 :       rate(:,803) = 1.4e-11_r8 * exp_fac(:)
     520    58907904 :       rate(:,806) = 1.9e-12_r8 * exp_fac(:)
     521    58907904 :       rate(:,829) = 6.34e-08_r8 * exp_fac(:)
     522    58907904 :       rate(:,846) = 1.9e-11_r8 * exp_fac(:)
     523    58907904 :       rate(:,849) = 1.2e-14_r8 * exp_fac(:)
     524    58907904 :       rate(:,850) = 2e-10_r8 * exp_fac(:)
     525    58907904 :       rate(:,854) = 2.5e-12_r8 * exp_fac(:)
     526    58907904 :       rate(:,866) = 1.34e-11_r8 * exp_fac(:)
     527    58907904 :       rate(:,867) = 1.2e-11_r8 * exp_fac(:)
     528    58907904 :       rate(:,872) = 1.1e-11_r8 * exp_fac(:)
     529    58907904 :       rate(:,876) = 2.1e-10_r8 * exp_fac(:)
     530    58907904 :       rate(:,877) = 1.34e-11_r8 * exp_fac(:)
     531    58907904 :       rate(:,881) = 1.7e-11_r8 * exp_fac(:)
     532    58907904 :       rate(:,901) = 1.29e-07_r8 * exp_fac(:)
     533    58907904 :       rate(:,902) = 2.31e-07_r8 * exp_fac(:)
     534    58907904 :       rate(:,903) = 2.31e-06_r8 * exp_fac(:)
     535    58907904 :       rate(:,904) = 4.63e-07_r8 * exp_fac(:)
     536    58907904 :       rate(:,278) = 6e-12_r8 * exp( 400._r8 * itemp(:) )
     537    58907904 :       rate(:,283) = 1.46e-11_r8 * exp( -1040._r8 * itemp(:) )
     538    58907904 :       rate(:,284) = 1.42e-12_r8 * exp( -1150._r8 * itemp(:) )
     539    58907904 :       exp_fac(:) = exp( -1520._r8 * itemp(:) )
     540    58907904 :       rate(:,285) = 1.64e-12_r8 * exp_fac(:)
     541    58907904 :       rate(:,403) = 8.5e-16_r8 * exp_fac(:)
     542    58907904 :       exp_fac(:) = exp( -1100._r8 * itemp(:) )
     543    58907904 :       rate(:,286) = 2.03e-11_r8 * exp_fac(:)
     544    58907904 :       rate(:,805) = 3.4e-12_r8 * exp_fac(:)
     545    58907904 :       rate(:,287) = 1.96e-12_r8 * exp( -1200._r8 * itemp(:) )
     546    58907904 :       rate(:,288) = 4.85e-12_r8 * exp( -850._r8 * itemp(:) )
     547    58907904 :       rate(:,289) = 9e-13_r8 * exp( -360._r8 * itemp(:) )
     548    58907904 :       exp_fac(:) = exp( -1600._r8 * itemp(:) )
     549    58907904 :       rate(:,290) = 1.25e-12_r8 * exp_fac(:)
     550    58907904 :       rate(:,300) = 3.4e-11_r8 * exp_fac(:)
     551    58907904 :       rate(:,291) = 1.3e-12_r8 * exp( -1770._r8 * itemp(:) )
     552    58907904 :       rate(:,292) = 9.2e-13_r8 * exp( -1560._r8 * itemp(:) )
     553    58907904 :       rate(:,298) = 9.7e-15_r8 * exp( 625._r8 * itemp(:) )
     554    58907904 :       rate(:,299) = 6e-13_r8 * exp( -2058._r8 * itemp(:) )
     555    58907904 :       rate(:,302) = 5e-13_r8 * exp( -424._r8 * itemp(:) )
     556    58907904 :       rate(:,303) = 1.9e-14_r8 * exp( 706._r8 * itemp(:) )
     557    58907904 :       rate(:,304) = 4.1e-13_r8 * exp( 750._r8 * itemp(:) )
     558    58907904 :       exp_fac(:) = exp( 300._r8 * itemp(:) )
     559    58907904 :       rate(:,305) = 2.8e-12_r8 * exp_fac(:)
     560    58907904 :       rate(:,370) = 2.9e-12_r8 * exp_fac(:)
     561    58907904 :       rate(:,306) = 2.9e-12_r8 * exp( -345._r8 * itemp(:) )
     562    58907904 :       rate(:,308) = 2.45e-12_r8 * exp( -1775._r8 * itemp(:) )
     563    58907904 :       exp_fac(:) = exp( 500._r8 * itemp(:) )
     564    58907904 :       rate(:,311) = 1.3e-12_r8 * exp_fac(:)
     565    58907904 :       rate(:,335) = 2.9e-12_r8 * exp_fac(:)
     566    58907904 :       rate(:,336) = 2e-12_r8 * exp_fac(:)
     567    58907904 :       rate(:,368) = 7.1e-13_r8 * exp_fac(:)
     568    58907904 :       rate(:,381) = 2e-12_r8 * exp_fac(:)
     569    58907904 :       rate(:,388) = 2.9e-12_r8 * exp_fac(:)
     570    58907904 :       rate(:,389) = 2e-12_r8 * exp_fac(:)
     571    58907904 :       rate(:,391) = 2.9e-12_r8 * exp_fac(:)
     572    58907904 :       rate(:,400) = 2e-12_r8 * exp_fac(:)
     573    58907904 :       rate(:,424) = 2e-12_r8 * exp_fac(:)
     574    58907904 :       rate(:,430) = 2e-12_r8 * exp_fac(:)
     575    58907904 :       rate(:,444) = 2e-12_r8 * exp_fac(:)
     576    58907904 :       rate(:,448) = 2e-12_r8 * exp_fac(:)
     577    58907904 :       rate(:,474) = 2e-12_r8 * exp_fac(:)
     578    58907904 :       rate(:,491) = 2e-12_r8 * exp_fac(:)
     579    58907904 :       rate(:,495) = 2e-12_r8 * exp_fac(:)
     580    58907904 :       rate(:,586) = 2e-12_r8 * exp_fac(:)
     581    58907904 :       rate(:,591) = 2e-12_r8 * exp_fac(:)
     582    58907904 :       rate(:,592) = 2e-12_r8 * exp_fac(:)
     583    58907904 :       rate(:,593) = 2e-12_r8 * exp_fac(:)
     584    58907904 :       rate(:,594) = 2e-12_r8 * exp_fac(:)
     585    58907904 :       rate(:,599) = 2e-12_r8 * exp_fac(:)
     586    58907904 :       rate(:,600) = 2e-12_r8 * exp_fac(:)
     587    58907904 :       rate(:,601) = 2e-12_r8 * exp_fac(:)
     588    58907904 :       rate(:,606) = 2e-12_r8 * exp_fac(:)
     589    58907904 :       rate(:,611) = 2e-12_r8 * exp_fac(:)
     590    58907904 :       rate(:,612) = 2e-12_r8 * exp_fac(:)
     591    58907904 :       rate(:,613) = 2e-12_r8 * exp_fac(:)
     592    58907904 :       rate(:,614) = 2e-12_r8 * exp_fac(:)
     593    58907904 :       rate(:,619) = 2e-12_r8 * exp_fac(:)
     594    58907904 :       rate(:,620) = 2e-12_r8 * exp_fac(:)
     595    58907904 :       rate(:,621) = 2e-12_r8 * exp_fac(:)
     596    58907904 :       rate(:,626) = 2e-12_r8 * exp_fac(:)
     597    58907904 :       rate(:,631) = 2e-12_r8 * exp_fac(:)
     598    58907904 :       rate(:,632) = 2e-12_r8 * exp_fac(:)
     599    58907904 :       rate(:,633) = 2e-12_r8 * exp_fac(:)
     600    58907904 :       rate(:,634) = 2e-12_r8 * exp_fac(:)
     601    58907904 :       rate(:,639) = 2e-12_r8 * exp_fac(:)
     602    58907904 :       rate(:,640) = 2e-12_r8 * exp_fac(:)
     603    58907904 :       rate(:,641) = 2e-12_r8 * exp_fac(:)
     604    58907904 :       rate(:,645) = 2e-12_r8 * exp_fac(:)
     605    58907904 :       rate(:,651) = 2e-12_r8 * exp_fac(:)
     606    58907904 :       rate(:,652) = 2e-12_r8 * exp_fac(:)
     607    58907904 :       rate(:,653) = 2e-12_r8 * exp_fac(:)
     608    58907904 :       rate(:,654) = 2e-12_r8 * exp_fac(:)
     609    58907904 :       rate(:,659) = 2e-12_r8 * exp_fac(:)
     610    58907904 :       rate(:,660) = 2e-12_r8 * exp_fac(:)
     611    58907904 :       rate(:,661) = 2e-12_r8 * exp_fac(:)
     612    58907904 :       rate(:,665) = 2e-12_r8 * exp_fac(:)
     613    58907904 :       rate(:,671) = 2e-12_r8 * exp_fac(:)
     614    58907904 :       rate(:,672) = 2e-12_r8 * exp_fac(:)
     615    58907904 :       rate(:,673) = 2e-12_r8 * exp_fac(:)
     616    58907904 :       rate(:,674) = 2e-12_r8 * exp_fac(:)
     617    58907904 :       rate(:,679) = 2e-12_r8 * exp_fac(:)
     618    58907904 :       rate(:,680) = 2e-12_r8 * exp_fac(:)
     619    58907904 :       rate(:,681) = 2e-12_r8 * exp_fac(:)
     620    58907904 :       rate(:,693) = 2e-12_r8 * exp_fac(:)
     621    58907904 :       rate(:,698) = 2e-12_r8 * exp_fac(:)
     622    58907904 :       rate(:,699) = 2e-12_r8 * exp_fac(:)
     623    58907904 :       rate(:,700) = 2e-12_r8 * exp_fac(:)
     624    58907904 :       rate(:,701) = 2.9e-12_r8 * exp_fac(:)
     625    58907904 :       rate(:,702) = 2e-12_r8 * exp_fac(:)
     626    58907904 :       rate(:,706) = 2.9e-12_r8 * exp_fac(:)
     627    58907904 :       rate(:,707) = 2.9e-12_r8 * exp_fac(:)
     628    58907904 :       rate(:,709) = 2e-12_r8 * exp_fac(:)
     629    58907904 :       rate(:,714) = 2e-12_r8 * exp_fac(:)
     630    58907904 :       rate(:,715) = 2e-12_r8 * exp_fac(:)
     631    58907904 :       rate(:,716) = 2e-12_r8 * exp_fac(:)
     632    58907904 :       rate(:,719) = 2.9e-12_r8 * exp_fac(:)
     633    58907904 :       rate(:,720) = 2e-12_r8 * exp_fac(:)
     634    58907904 :       rate(:,724) = 2.9e-12_r8 * exp_fac(:)
     635    58907904 :       rate(:,725) = 2.9e-12_r8 * exp_fac(:)
     636    58907904 :       rate(:,726) = 2.9e-12_r8 * exp_fac(:)
     637    58907904 :       rate(:,728) = 2e-12_r8 * exp_fac(:)
     638    58907904 :       rate(:,733) = 2e-12_r8 * exp_fac(:)
     639    58907904 :       rate(:,734) = 2e-12_r8 * exp_fac(:)
     640    58907904 :       rate(:,735) = 2e-12_r8 * exp_fac(:)
     641    58907904 :       rate(:,738) = 2e-12_r8 * exp_fac(:)
     642    58907904 :       rate(:,743) = 2e-12_r8 * exp_fac(:)
     643    58907904 :       rate(:,744) = 2e-12_r8 * exp_fac(:)
     644    58907904 :       rate(:,745) = 2e-12_r8 * exp_fac(:)
     645    58907904 :       rate(:,749) = 2.9e-12_r8 * exp_fac(:)
     646    58907904 :       rate(:,750) = 2e-12_r8 * exp_fac(:)
     647    58907904 :       rate(:,754) = 2.9e-12_r8 * exp_fac(:)
     648    58907904 :       rate(:,312) = 5.6e-15_r8 * exp( 2300._r8 * itemp(:) )
     649    58907904 :       rate(:,313) = 2.4e+12_r8 * exp( -7000._r8 * itemp(:) )
     650    58907904 :       rate(:,314) = 2.6e-12_r8 * exp( 265._r8 * itemp(:) )
     651    58907904 :       rate(:,318) = 1.08e-10_r8 * exp( 105._r8 * itemp(:) )
     652    58907904 :       rate(:,323) = 1.2e-14_r8 * exp( -2630._r8 * itemp(:) )
     653    58907904 :       exp_fac(:) = exp( 700._r8 * itemp(:) )
     654    58907904 :       rate(:,326) = 7.5e-13_r8 * exp_fac(:)
     655    58907904 :       rate(:,341) = 7.5e-13_r8 * exp_fac(:)
     656    58907904 :       rate(:,356) = 7.5e-13_r8 * exp_fac(:)
     657    58907904 :       rate(:,365) = 7.5e-13_r8 * exp_fac(:)
     658    58907904 :       rate(:,369) = 8.6e-13_r8 * exp_fac(:)
     659    58907904 :       rate(:,394) = 7.5e-13_r8 * exp_fac(:)
     660    58907904 :       rate(:,409) = 7.5e-13_r8 * exp_fac(:)
     661    58907904 :       rate(:,544) = 7.5e-13_r8 * exp_fac(:)
     662    58907904 :       rate(:,548) = 7.5e-13_r8 * exp_fac(:)
     663    58907904 :       rate(:,551) = 7.5e-13_r8 * exp_fac(:)
     664    58907904 :       rate(:,564) = 7.5e-13_r8 * exp_fac(:)
     665    58907904 :       rate(:,571) = 7.5e-13_r8 * exp_fac(:)
     666    58907904 :       rate(:,577) = 7.5e-13_r8 * exp_fac(:)
     667    58907904 :       rate(:,580) = 7.5e-13_r8 * exp_fac(:)
     668    58907904 :       rate(:,852) = 7.5e-13_r8 * exp_fac(:)
     669    58907904 :       rate(:,864) = 7.5e-13_r8 * exp_fac(:)
     670    58907904 :       rate(:,879) = 7.5e-13_r8 * exp_fac(:)
     671    58907904 :       rate(:,882) = 7.5e-13_r8 * exp_fac(:)
     672    58907904 :       exp_fac(:) = exp( 365._r8 * itemp(:) )
     673    58907904 :       rate(:,327) = 2.6e-12_r8 * exp_fac(:)
     674    58907904 :       rate(:,545) = 2.6e-12_r8 * exp_fac(:)
     675    58907904 :       rate(:,550) = 2.6e-12_r8 * exp_fac(:)
     676    58907904 :       rate(:,552) = 2.6e-12_r8 * exp_fac(:)
     677    58907904 :       rate(:,565) = 2.6e-12_r8 * exp_fac(:)
     678    58907904 :       rate(:,572) = 2.6e-12_r8 * exp_fac(:)
     679    58907904 :       rate(:,578) = 2.6e-12_r8 * exp_fac(:)
     680    58907904 :       rate(:,581) = 2.6e-12_r8 * exp_fac(:)
     681    58907904 :       rate(:,853) = 2.6e-12_r8 * exp_fac(:)
     682    58907904 :       rate(:,865) = 2.6e-12_r8 * exp_fac(:)
     683    58907904 :       rate(:,880) = 2.6e-12_r8 * exp_fac(:)
     684    58907904 :       rate(:,883) = 2.6e-12_r8 * exp_fac(:)
     685    58907904 :       rate(:,328) = 6.9e-12_r8 * exp( -230._r8 * itemp(:) )
     686    58907904 :       rate(:,330) = 7.2e-11_r8 * exp( -70._r8 * itemp(:) )
     687    58907904 :       rate(:,331) = 7.66e-12_r8 * exp( -1020._r8 * itemp(:) )
     688    58907904 :       exp_fac(:) = exp( -1900._r8 * itemp(:) )
     689    58907904 :       rate(:,332) = 1.4e-12_r8 * exp_fac(:)
     690    58907904 :       rate(:,354) = 6.5e-15_r8 * exp_fac(:)
     691    58907904 :       rate(:,333) = 4.63e-12_r8 * exp( 350._r8 * itemp(:) )
     692    58907904 :       rate(:,334) = 7.8e-13_r8 * exp( -1050._r8 * itemp(:) )
     693    58907904 :       exp_fac(:) = exp( 1040._r8 * itemp(:) )
     694    58907904 :       rate(:,337) = 4.3e-13_r8 * exp_fac(:)
     695    58907904 :       rate(:,390) = 4.3e-13_r8 * exp_fac(:)
     696    58907904 :       rate(:,541) = 4.3e-13_r8 * exp_fac(:)
     697    58907904 :       rate(:,555) = 4.3e-13_r8 * exp_fac(:)
     698    58907904 :       rate(:,558) = 4.3e-13_r8 * exp_fac(:)
     699    58907904 :       rate(:,561) = 4.3e-13_r8 * exp_fac(:)
     700    58907904 :       rate(:,703) = 4.3e-13_r8 * exp_fac(:)
     701    58907904 :       rate(:,721) = 4.3e-13_r8 * exp_fac(:)
     702    58907904 :       rate(:,751) = 4.3e-13_r8 * exp_fac(:)
     703    58907904 :       rate(:,339) = 3.15e-14_r8 * exp( 920._r8 * itemp(:) )
     704    58907904 :       rate(:,343) = 1.6e+11_r8 * exp( -4150._r8 * itemp(:) )
     705    58907904 :       rate(:,353) = 4.6e-13_r8 * exp( -1156._r8 * itemp(:) )
     706    58907904 :       rate(:,355) = 3.75e-13_r8 * exp( -40._r8 * itemp(:) )
     707    58907904 :       rate(:,359) = 9.19e-12_r8 * exp( -630._r8 * itemp(:) )
     708    58907904 :       rate(:,360) = 1.4e-12_r8 * exp( -1860._r8 * itemp(:) )
     709    58907904 :       rate(:,361) = 8.4e-13_r8 * exp( 830._r8 * itemp(:) )
     710    58907904 :       rate(:,377) = 4.8e-12_r8 * exp( 120._r8 * itemp(:) )
     711    58907904 :       rate(:,378) = 5.1e-14_r8 * exp( 693._r8 * itemp(:) )
     712    58907904 :       exp_fac(:) = exp( 1300._r8 * itemp(:) )
     713    58907904 :       rate(:,383) = 2.11e-13_r8 * exp_fac(:)
     714    58907904 :       rate(:,402) = 2.11e-13_r8 * exp_fac(:)
     715    58907904 :       rate(:,421) = 2.38e-13_r8 * exp_fac(:)
     716    58907904 :       rate(:,426) = 2.12e-13_r8 * exp_fac(:)
     717    58907904 :       rate(:,432) = 2.12e-13_r8 * exp_fac(:)
     718    58907904 :       rate(:,446) = 2.12e-13_r8 * exp_fac(:)
     719    58907904 :       rate(:,450) = 2.12e-13_r8 * exp_fac(:)
     720    58907904 :       rate(:,457) = 2.6e-13_r8 * exp_fac(:)
     721    58907904 :       rate(:,461) = 2.6e-13_r8 * exp_fac(:)
     722    58907904 :       rate(:,464) = 2.6e-13_r8 * exp_fac(:)
     723    58907904 :       rate(:,467) = 2.6e-13_r8 * exp_fac(:)
     724    58907904 :       rate(:,471) = 2.6e-13_r8 * exp_fac(:)
     725    58907904 :       rate(:,476) = 2.47e-13_r8 * exp_fac(:)
     726    58907904 :       rate(:,479) = 2.64e-13_r8 * exp_fac(:)
     727    58907904 :       rate(:,482) = 2.64e-13_r8 * exp_fac(:)
     728    58907904 :       rate(:,493) = 2.12e-13_r8 * exp_fac(:)
     729    58907904 :       rate(:,497) = 2.12e-13_r8 * exp_fac(:)
     730    58907904 :       rate(:,499) = 2.6e-13_r8 * exp_fac(:)
     731    58907904 :       rate(:,588) = 2.71e-13_r8 * exp_fac(:)
     732    58907904 :       rate(:,596) = 2.6e-13_r8 * exp_fac(:)
     733    58907904 :       rate(:,608) = 2.78e-13_r8 * exp_fac(:)
     734    58907904 :       rate(:,616) = 2.75e-13_r8 * exp_fac(:)
     735    58907904 :       rate(:,628) = 2.71e-13_r8 * exp_fac(:)
     736    58907904 :       rate(:,636) = 2.6e-13_r8 * exp_fac(:)
     737    58907904 :       rate(:,647) = 2.71e-13_r8 * exp_fac(:)
     738    58907904 :       rate(:,656) = 2.6e-13_r8 * exp_fac(:)
     739    58907904 :       rate(:,667) = 2.71e-13_r8 * exp_fac(:)
     740    58907904 :       rate(:,676) = 2.6e-13_r8 * exp_fac(:)
     741    58907904 :       rate(:,687) = 2.71e-13_r8 * exp_fac(:)
     742    58907904 :       rate(:,691) = 2.71e-13_r8 * exp_fac(:)
     743    58907904 :       rate(:,695) = 2.54e-13_r8 * exp_fac(:)
     744    58907904 :       rate(:,711) = 2.62e-13_r8 * exp_fac(:)
     745    58907904 :       rate(:,730) = 2.66e-13_r8 * exp_fac(:)
     746    58907904 :       rate(:,740) = 2.51e-13_r8 * exp_fac(:)
     747    58907904 :       rate(:,760) = 2.68e-13_r8 * exp_fac(:)
     748    58907904 :       rate(:,765) = 2.47e-13_r8 * exp_fac(:)
     749    58907904 :       rate(:,772) = 2.76e-13_r8 * exp_fac(:)
     750    58907904 :       rate(:,776) = 2.76e-13_r8 * exp_fac(:)
     751    58907904 :       rate(:,780) = 2.75e-13_r8 * exp_fac(:)
     752    58907904 :       rate(:,784) = 2.75e-13_r8 * exp_fac(:)
     753    58907904 :       rate(:,842) = 2.6e-13_r8 * exp_fac(:)
     754    58907904 :       rate(:,847) = 2.75e-13_r8 * exp_fac(:)
     755    58907904 :       rate(:,855) = 2.6e-13_r8 * exp_fac(:)
     756    58907904 :       rate(:,860) = 2.12e-13_r8 * exp_fac(:)
     757    58907904 :       rate(:,868) = 2.6e-13_r8 * exp_fac(:)
     758    58907904 :       rate(:,873) = 2.6e-13_r8 * exp_fac(:)
     759    58907904 :       rate(:,384) = 2.9e+07_r8 * exp( -5297._r8 * itemp(:) )
     760    58907904 :       rate(:,385) = 1.5e-15_r8 * exp( -2100._r8 * itemp(:) )
     761    58907904 :       exp_fac(:) = exp( 360._r8 * itemp(:) )
     762    58907904 :       rate(:,386) = 9.6e-12_r8 * exp_fac(:)
     763    58907904 :       rate(:,589) = 2.7e-12_r8 * exp_fac(:)
     764    58907904 :       rate(:,597) = 2.7e-12_r8 * exp_fac(:)
     765    58907904 :       rate(:,609) = 2.7e-12_r8 * exp_fac(:)
     766    58907904 :       rate(:,617) = 2.7e-12_r8 * exp_fac(:)
     767    58907904 :       rate(:,629) = 2.7e-12_r8 * exp_fac(:)
     768    58907904 :       rate(:,637) = 2.7e-12_r8 * exp_fac(:)
     769    58907904 :       rate(:,649) = 2.7e-12_r8 * exp_fac(:)
     770    58907904 :       rate(:,657) = 2.7e-12_r8 * exp_fac(:)
     771    58907904 :       rate(:,669) = 2.7e-12_r8 * exp_fac(:)
     772    58907904 :       rate(:,677) = 2.7e-12_r8 * exp_fac(:)
     773    58907904 :       rate(:,688) = 2.7e-12_r8 * exp_fac(:)
     774    58907904 :       rate(:,692) = 2.7e-12_r8 * exp_fac(:)
     775    58907904 :       rate(:,696) = 2.7e-12_r8 * exp_fac(:)
     776    58907904 :       rate(:,712) = 2.7e-12_r8 * exp_fac(:)
     777    58907904 :       rate(:,731) = 2.7e-12_r8 * exp_fac(:)
     778    58907904 :       rate(:,741) = 2.7e-12_r8 * exp_fac(:)
     779    58907904 :       rate(:,761) = 2.7e-12_r8 * exp_fac(:)
     780    58907904 :       rate(:,766) = 2.7e-12_r8 * exp_fac(:)
     781    58907904 :       rate(:,773) = 2.7e-12_r8 * exp_fac(:)
     782    58907904 :       rate(:,777) = 2.7e-12_r8 * exp_fac(:)
     783    58907904 :       rate(:,781) = 2.7e-12_r8 * exp_fac(:)
     784    58907904 :       rate(:,785) = 2.7e-12_r8 * exp_fac(:)
     785    58907904 :       rate(:,843) = 2.7e-12_r8 * exp_fac(:)
     786    58907904 :       rate(:,848) = 2.7e-12_r8 * exp_fac(:)
     787    58907904 :       rate(:,856) = 2.7e-12_r8 * exp_fac(:)
     788    58907904 :       rate(:,861) = 2.7e-12_r8 * exp_fac(:)
     789    58907904 :       rate(:,869) = 2.7e-12_r8 * exp_fac(:)
     790    58907904 :       rate(:,874) = 2.7e-12_r8 * exp_fac(:)
     791    58907904 :       rate(:,396) = 2.3e-12_r8 * exp( -170._r8 * itemp(:) )
     792    58907904 :       rate(:,404) = 2.7e-12_r8 * exp( 580._r8 * itemp(:) )
     793    58907904 :       rate(:,411) = 5.4e-14_r8 * exp( 870._r8 * itemp(:) )
     794    58907904 :       exp_fac(:) = exp( 450._r8 * itemp(:) )
     795    58907904 :       rate(:,414) = 1.17e-11_r8 * exp_fac(:)
     796    58907904 :       rate(:,415) = 1.17e-11_r8 * exp_fac(:)
     797    58907904 :       exp_fac(:) = exp( 390._r8 * itemp(:) )
     798    58907904 :       rate(:,416) = 2.2e-11_r8 * exp_fac(:)
     799    58907904 :       rate(:,417) = 3.5e-11_r8 * exp_fac(:)
     800    58907904 :       rate(:,487) = 2.7e-11_r8 * exp_fac(:)
     801    58907904 :       rate(:,490) = 2.08e-11_r8 * exp_fac(:)
     802    58907904 :       rate(:,768) = 2.7e-11_r8 * exp_fac(:)
     803    58907904 :       rate(:,863) = 2.7e-11_r8 * exp_fac(:)
     804    58907904 :       exp_fac(:) = exp( 410._r8 * itemp(:) )
     805    58907904 :       rate(:,419) = 9.85e-12_r8 * exp_fac(:)
     806    58907904 :       rate(:,603) = 1.34e-11_r8 * exp_fac(:)
     807    58907904 :       rate(:,845) = 1.34e-11_r8 * exp_fac(:)
     808    58907904 :       exp_fac(:) = exp( -400._r8 * itemp(:) )
     809    58907904 :       rate(:,420) = 4.43e-11_r8 * exp_fac(:)
     810    58907904 :       rate(:,422) = 4.43e-11_r8 * exp_fac(:)
     811    58907904 :       rate(:,423) = 3.22e-11_r8 * exp_fac(:)
     812    58907904 :       rate(:,427) = 1.04e+11_r8 * exp( -9746._r8 * itemp(:) )
     813    58907904 :       rate(:,428) = 2.24e+15_r8 * exp( -10865._r8 * itemp(:) )
     814    58907904 :       rate(:,429) = 2.22e+15_r8 * exp( -10355._r8 * itemp(:) )
     815    58907904 :       rate(:,433) = 1.88e+11_r8 * exp( -9752._r8 * itemp(:) )
     816    58907904 :       rate(:,434) = 2.49e+15_r8 * exp( -11112._r8 * itemp(:) )
     817    58907904 :       rate(:,435) = 2.49e+15_r8 * exp( -10890._r8 * itemp(:) )
     818    58907904 :       rate(:,447) = 1.83e+14_r8 * exp( -8930._r8 * itemp(:) )
     819    58907904 :       rate(:,451) = 2.08e+14_r8 * exp( -9400._r8 * itemp(:) )
     820    58907904 :       exp_fac(:) = exp( -10000._r8 * itemp(:) )
     821    58907904 :       rate(:,458) = 1.256e+13_r8 * exp_fac(:)
     822    58907904 :       rate(:,462) = 1.875e+13_r8 * exp_fac(:)
     823    58907904 :       rate(:,465) = 1.875e+13_r8 * exp_fac(:)
     824    58907904 :       rate(:,468) = 5.092e+12_r8 * exp_fac(:)
     825    58907904 :       rate(:,480) = 8.72e+12_r8 * exp_fac(:)
     826    58907904 :       rate(:,483) = 6.55e+12_r8 * exp_fac(:)
     827    58907904 :       exp_fac(:) = exp( -450._r8 * itemp(:) )
     828    58907904 :       rate(:,473) = 2.95e-12_r8 * exp_fac(:)
     829    58907904 :       rate(:,764) = 2.95e-12_r8 * exp_fac(:)
     830    58907904 :       rate(:,859) = 2.95e-12_r8 * exp_fac(:)
     831    58907904 :       exp_fac(:) = exp( -1995._r8 * itemp(:) )
     832    58907904 :       rate(:,486) = 1.03e-14_r8 * exp_fac(:)
     833    58907904 :       rate(:,862) = 1.03e-14_r8 * exp_fac(:)
     834    58907904 :       rate(:,494) = 1.79e+14_r8 * exp( -8830._r8 * itemp(:) )
     835    58907904 :       rate(:,498) = 1.75e+14_r8 * exp( -9054._r8 * itemp(:) )
     836    58907904 :       rate(:,500) = 1e+07_r8 * exp( -5000._r8 * itemp(:) )
     837    58907904 :       exp_fac(:) = exp( -193._r8 * itemp(:) )
     838    58907904 :       rate(:,543) = 2.3e-12_r8 * exp_fac(:)
     839    58907904 :       rate(:,851) = 2.3e-12_r8 * exp_fac(:)
     840    58907904 :       rate(:,547) = 5.9e-12_r8 * exp( 225._r8 * itemp(:) )
     841    58907904 :       rate(:,566) = 4.7e-13_r8 * exp( 1220._r8 * itemp(:) )
     842    58907904 :       exp_fac(:) = exp( 352._r8 * itemp(:) )
     843    58907904 :       rate(:,574) = 1.7e-12_r8 * exp_fac(:)
     844    58907904 :       rate(:,878) = 1.7e-12_r8 * exp_fac(:)
     845    58907904 :       exp_fac(:) = exp( 490._r8 * itemp(:) )
     846    58907904 :       rate(:,584) = 1.2e-12_r8 * exp_fac(:)
     847    58907904 :       rate(:,841) = 1.2e-12_r8 * exp_fac(:)
     848    58907904 :       exp_fac(:) = exp( -640._r8 * itemp(:) )
     849    58907904 :       rate(:,602) = 8.05e-16_r8 * exp_fac(:)
     850    58907904 :       rate(:,844) = 8.05e-16_r8 * exp_fac(:)
     851    58907904 :       exp_fac(:) = exp( -770._r8 * itemp(:) )
     852    58907904 :       rate(:,662) = 2.8e-15_r8 * exp_fac(:)
     853    58907904 :       rate(:,870) = 2.8e-15_r8 * exp_fac(:)
     854    58907904 :       exp_fac(:) = exp( 470._r8 * itemp(:) )
     855    58907904 :       rate(:,663) = 3.41e-11_r8 * exp_fac(:)
     856    58907904 :       rate(:,871) = 3.41e-11_r8 * exp_fac(:)
     857    58907904 :       exp_fac(:) = exp( -520._r8 * itemp(:) )
     858    58907904 :       rate(:,682) = 2.65e-15_r8 * exp_fac(:)
     859    58907904 :       rate(:,875) = 2.65e-15_r8 * exp_fac(:)
     860    58907904 :       exp_fac(:) = exp( 600._r8 * itemp(:) )
     861    58907904 :       rate(:,717) = 5.2e-12_r8 * exp_fac(:)
     862    58907904 :       rate(:,736) = 5.2e-12_r8 * exp_fac(:)
     863    58907904 :       rate(:,756) = 5.2e-12_r8 * exp_fac(:)
     864    58907904 :       rate(:,793) = 1.9e-13_r8 * exp( 520._r8 * itemp(:) )
     865    58907904 :       rate(:,794) = 1.1e-11_r8 * exp( -280._r8 * itemp(:) )
     866    58907904 :       rate(:,795) = 2.1e-11_r8 * exp( -2200._r8 * itemp(:) )
     867    58907904 :       rate(:,796) = 7.2e-14_r8 * exp( -1070._r8 * itemp(:) )
     868    58907904 :       rate(:,804) = 1.6e-13_r8 * exp( -2280._r8 * itemp(:) )
     869    58907904 :       rate(:,807) = 2.6e-11_r8 * exp( 330._r8 * itemp(:) )
     870    58907904 :       rate(:,810) = 1.7e-12_r8 * exp( -710._r8 * itemp(:) )
     871             : 
     872    58907904 :       itemp(:) = 300._r8 * itemp(:)
     873             :  
     874       72960 :       n = ncol*pver
     875             : 
     876    58907904 :       ko(:) = 5.3e-32_r8 * itemp(:)**1.8_r8
     877    58907904 :       kinf(:) = 9.5e-11_r8 * itemp(:)**(-0.4_r8)
     878       72960 :       call jpl( rate(:,181), m, 0.6_r8, ko, kinf, n )
     879             : 
     880    58907904 :       ko(:) = 6.9e-31_r8 * itemp(:)**1._r8
     881    58907904 :       kinf(:) = 2.6e-11_r8
     882       72960 :       call jpl( rate(:,191), m, 0.6_r8, ko, kinf, n )
     883             : 
     884    58907904 :       ko(:) = 2.5e-31_r8 * itemp(:)**1.8_r8
     885    58907904 :       kinf(:) = 2.2e-11_r8 * itemp(:)**0.7_r8
     886       72960 :       call jpl( rate(:,201), m, 0.6_r8, ko, kinf, n )
     887             : 
     888    58907904 :       ko(:) = 9e-32_r8 * itemp(:)**1.5_r8
     889    58907904 :       kinf(:) = 3e-11_r8
     890       72960 :       call jpl( rate(:,209), m, 0.6_r8, ko, kinf, n )
     891             : 
     892    58907904 :       ko(:) = 1.9e-31_r8 * itemp(:)**3.4_r8
     893    58907904 :       kinf(:) = 4e-12_r8 * itemp(:)**0.3_r8
     894       72960 :       call jpl( rate(:,212), m, 0.6_r8, ko, kinf, n )
     895             : 
     896    58907904 :       ko(:) = 2.4e-30_r8 * itemp(:)**3._r8
     897    58907904 :       kinf(:) = 1.6e-12_r8 * itemp(:)**(-0.1_r8)
     898       72960 :       call jpl( rate(:,213), m, 0.6_r8, ko, kinf, n )
     899             : 
     900    58907904 :       ko(:) = 1.8e-30_r8 * itemp(:)**3._r8
     901    58907904 :       kinf(:) = 2.8e-11_r8
     902       72960 :       call jpl( rate(:,214), m, 0.6_r8, ko, kinf, n )
     903             : 
     904    58907904 :       ko(:) = 1.8e-31_r8 * itemp(:)**3.4_r8
     905    58907904 :       kinf(:) = 1.5e-11_r8 * itemp(:)**1.9_r8
     906       72960 :       call jpl( rate(:,232), m, 0.6_r8, ko, kinf, n )
     907             : 
     908    58907904 :       ko(:) = 1.9e-32_r8 * itemp(:)**3.6_r8
     909    58907904 :       kinf(:) = 3.7e-12_r8 * itemp(:)**1.6_r8
     910       72960 :       call jpl( rate(:,252), m, 0.6_r8, ko, kinf, n )
     911             : 
     912    58907904 :       ko(:) = 5.2e-31_r8 * itemp(:)**3.2_r8
     913    58907904 :       kinf(:) = 6.9e-12_r8 * itemp(:)**2.9_r8
     914       72960 :       call jpl( rate(:,263), m, 0.6_r8, ko, kinf, n )
     915             : 
     916    58907904 :       ko(:) = 6.1e-33_r8 * itemp(:)**1.5_r8
     917    58907904 :       kinf(:) = 9.8e-15_r8 * itemp(:)**(-4.6_r8)
     918       72960 :       call jpl( rate(:,309), m, 0.8_r8, ko, kinf, n )
     919             : 
     920    58907904 :       ko(:) = 5.2e-30_r8 * itemp(:)**2.4_r8
     921    58907904 :       kinf(:) = 2.2e-10_r8 * itemp(:)**0.7_r8
     922       72960 :       call jpl( rate(:,320), m, 0.6_r8, ko, kinf, n )
     923             : 
     924    58907904 :       ko(:) = 5.5e-30_r8
     925    58907904 :       kinf(:) = 8.3e-13_r8 * itemp(:)**(-2._r8)
     926       72960 :       call jpl( rate(:,321), m, 0.6_r8, ko, kinf, n )
     927             : 
     928    58907904 :       ko(:) = 1.6e-29_r8 * itemp(:)**3.3_r8
     929    58907904 :       kinf(:) = 3.1e-10_r8 * itemp(:)
     930       72960 :       call jpl( rate(:,322), m, 0.6_r8, ko, kinf, n )
     931             : 
     932    58907904 :       ko(:) = 8.6e-29_r8 * itemp(:)**3.1_r8
     933    58907904 :       kinf(:) = 9e-12_r8 * itemp(:)**0.85_r8
     934       72960 :       call jpl( rate(:,350), m, 0.48_r8, ko, kinf, n )
     935             : 
     936    58907904 :       ko(:) = 7.3e-29_r8 * itemp(:)**4.1_r8
     937    58907904 :       kinf(:) = 9.5e-12_r8 * itemp(:)**1.6_r8
     938       72960 :       call jpl( rate(:,351), m, 0.6_r8, ko, kinf, n )
     939             : 
     940    58907904 :       ko(:) = 8e-27_r8 * itemp(:)**3.5_r8
     941    58907904 :       kinf(:) = 3e-11_r8
     942       72960 :       call jpl( rate(:,372), m, 0.5_r8, ko, kinf, n )
     943             : 
     944    58907904 :       ko(:) = 8e-27_r8 * itemp(:)**3.5_r8
     945    58907904 :       kinf(:) = 3e-11_r8
     946       72960 :       call jpl( rate(:,398), m, 0.5_r8, ko, kinf, n )
     947             : 
     948    58907904 :       ko(:) = 9.7e-29_r8 * itemp(:)**5.6_r8
     949    58907904 :       kinf(:) = 9.3e-12_r8 * itemp(:)**1.5_r8
     950       72960 :       call jpl( rate(:,406), m, 0.6_r8, ko, kinf, n )
     951             : 
     952    58907904 :       ko(:) = 9.7e-29_r8 * itemp(:)**5.6_r8
     953    58907904 :       kinf(:) = 9.3e-12_r8 * itemp(:)**1.5_r8
     954       72960 :       call jpl( rate(:,557), m, 0.6_r8, ko, kinf, n )
     955             : 
     956    58907904 :       ko(:) = 9.7e-29_r8 * itemp(:)**5.6_r8
     957    58907904 :       kinf(:) = 9.3e-12_r8 * itemp(:)**1.5_r8
     958       72960 :       call jpl( rate(:,560), m, 0.6_r8, ko, kinf, n )
     959             : 
     960    58907904 :       ko(:) = 9.7e-29_r8 * itemp(:)**5.6_r8
     961    58907904 :       kinf(:) = 9.3e-12_r8 * itemp(:)**1.5_r8
     962       72960 :       call jpl( rate(:,563), m, 0.6_r8, ko, kinf, n )
     963             : 
     964    58907904 :       ko(:) = 9.7e-29_r8 * itemp(:)**5.6_r8
     965    58907904 :       kinf(:) = 9.3e-12_r8 * itemp(:)**1.5_r8
     966       72960 :       call jpl( rate(:,570), m, 0.6_r8, ko, kinf, n )
     967             : 
     968    58907904 :       ko(:) = 9.7e-29_r8 * itemp(:)**5.6_r8
     969    58907904 :       kinf(:) = 9.3e-12_r8 * itemp(:)**1.5_r8
     970       72960 :       call jpl( rate(:,684), m, 0.6_r8, ko, kinf, n )
     971             : 
     972    58907904 :       ko(:) = 9.7e-29_r8 * itemp(:)**5.6_r8
     973    58907904 :       kinf(:) = 9.3e-12_r8 * itemp(:)**1.5_r8
     974       72960 :       call jpl( rate(:,685), m, 0.6_r8, ko, kinf, n )
     975             : 
     976    58907904 :       ko(:) = 9.7e-29_r8 * itemp(:)**5.6_r8
     977    58907904 :       kinf(:) = 9.3e-12_r8 * itemp(:)**1.5_r8
     978       72960 :       call jpl( rate(:,686), m, 0.6_r8, ko, kinf, n )
     979             : 
     980    58907904 :       ko(:) = 2.9e-31_r8 * itemp(:)**4.1_r8
     981    58907904 :       kinf(:) = 1.7e-12_r8 * itemp(:)**(-0.2_r8)
     982       72960 :       call jpl( rate(:,798), m, 0.6_r8, ko, kinf, n )
     983             : 
     984       72960 :       end subroutine setrxt
     985             : 
     986             : 
     987           0 :       subroutine setrxt_hrates( rate, temp, m, ncol, kbot )
     988             :  
     989             :       use ppgrid, only : pcols, pver
     990             : 
     991             : 
     992             :       use chem_mods, only : rxntot
     993             :       use mo_jpl,    only : jpl
     994             : 
     995             :       implicit none
     996             : 
     997             : !-------------------------------------------------------
     998             : !       ... dummy arguments
     999             : !-------------------------------------------------------
    1000             :       integer, intent(in) :: ncol
    1001             :       integer, intent(in) :: kbot
    1002             :       real(r8), intent(in)    :: temp(pcols,pver)
    1003             :       real(r8), intent(in)    :: m(ncol*pver)
    1004             :       real(r8), intent(inout) :: rate(ncol*pver,max(1,rxntot))
    1005             : 
    1006             : !-------------------------------------------------------
    1007             : !       ... local variables
    1008             : !-------------------------------------------------------
    1009             :       integer   ::  n
    1010             :       integer   ::  offset
    1011             :       integer   ::  k
    1012           0 :       real(r8)  :: itemp(ncol*kbot)
    1013             :       real(r8)  :: exp_fac(ncol*kbot)
    1014           0 :       real(r8)  :: ko(ncol*kbot)
    1015           0 :       real(r8)  :: kinf(ncol*kbot)
    1016           0 :       real(r8)  :: wrk(ncol*kbot)
    1017             :  
    1018           0 :       n = ncol*kbot
    1019             : 
    1020           0 :       rate(:n,178) = 6.9e-12_r8
    1021             :  
    1022           0 :       do k = 1,kbot
    1023           0 :         offset = (k-1)*ncol
    1024           0 :         itemp(offset+1:offset+ncol) = 1._r8 / temp(:ncol,k)
    1025             :       end do
    1026             : 
    1027           0 :       rate(:n,170) = 2.15e-11_r8 * exp( 110._r8 * itemp(:) )
    1028           0 :       rate(:n,173) = 8e-12_r8 * exp( -2060._r8 * itemp(:) )
    1029           0 :       rate(:n,182) = 3e-11_r8 * exp( 200._r8 * itemp(:) )
    1030           0 :       rate(:n,183) = 1e-14_r8 * exp( -490._r8 * itemp(:) )
    1031           0 :       rate(:n,184) = 1.4e-10_r8 * exp( -470._r8 * itemp(:) )
    1032           0 :       rate(:n,187) = 4.8e-11_r8 * exp( 250._r8 * itemp(:) )
    1033           0 :       rate(:n,188) = 1.8e-11_r8 * exp( 180._r8 * itemp(:) )
    1034           0 :       rate(:n,189) = 1.7e-12_r8 * exp( -940._r8 * itemp(:) )
    1035           0 :       rate(:n,194) = 2.1e-11_r8 * exp( 100._r8 * itemp(:) )
    1036           0 :       rate(:n,198) = 3.3e-12_r8 * exp( -3150._r8 * itemp(:) )
    1037           0 :       rate(:n,199) = 5.1e-12_r8 * exp( 210._r8 * itemp(:) )
    1038           0 :       rate(:n,207) = 3.44e-12_r8 * exp( 260._r8 * itemp(:) )
    1039           0 :       rate(:n,208) = 3e-12_r8 * exp( -1500._r8 * itemp(:) )
    1040             : 
    1041           0 :       itemp(:) = 300._r8 * itemp(:)
    1042             : 
    1043           0 :       ko(:) = 5.3e-32_r8 * itemp(:)**1.8_r8
    1044           0 :       kinf(:) = 9.5e-11_r8 * itemp(:)**(-0.4_r8)
    1045           0 :       call jpl( wrk, m, 0.6_r8, ko, kinf, n )
    1046           0 :       rate(:n,181) = wrk(:)
    1047             : 
    1048             : 
    1049             : 
    1050             : 
    1051             : 
    1052             : 
    1053             : 
    1054             : 
    1055             : 
    1056             : 
    1057             : 
    1058             : 
    1059             : 
    1060             : 
    1061             : 
    1062             : 
    1063             : 
    1064             : 
    1065             : 
    1066             : 
    1067             : 
    1068             : 
    1069             : 
    1070             : 
    1071             : 
    1072             : 
    1073             : 
    1074           0 :       end subroutine setrxt_hrates
    1075             : 
    1076             :       end module mo_setrxt

Generated by: LCOV version 1.14