Line data    Source code

       1             : 
       2             :       module mo_jpl
       3             : 
       4             :       private
       5             :       public :: jpl
       6             : 
       7             :       contains
       8             : 
       9     1489176 :       subroutine jpl( rate, m, factor, ko, kinf, ncol )
      10             : !-----------------------------------------------------------------
      11             : !        ... Calculate JPL troe rate
      12             : !-----------------------------------------------------------------
      13             :       
      14             :       use shr_kind_mod, only : r8 => shr_kind_r8
      15             : 
      16             :       implicit none
      17             : 
      18             : !-----------------------------------------------------------------
      19             : !        ... Dummy args
      20             : !-----------------------------------------------------------------
      21             :       integer, intent(in)   ::   ncol
      22             :       real(r8), intent(in)  ::   factor
      23             :       real(r8), intent(in)  ::   ko(ncol)
      24             :       real(r8), intent(in)  ::   kinf(ncol)
      25             :       real(r8), intent(in)  ::   m(ncol)
      26             :       real(r8), intent(out) ::   rate(ncol)
      27             : 
      28             : !-----------------------------------------------------------------
      29             : !        ... Local variables
      30             : !-----------------------------------------------------------------
      31     1489176 :       real(r8)  ::  xpo(ncol)
      32             : 
      33  2175512976 :       xpo(:)  = ko(:) * m(:) / kinf(:)
      34  2175512976 :       rate(:) = ko(:) / (1._r8 + xpo(:))
      35  2175512976 :       xpo(:)  = log10( xpo(:) )
      36  2175512976 :       xpo(:)  = 1._r8 / (1._r8 + xpo(:)*xpo(:))
      37  2175512976 :       rate(:) = rate(:) * factor**xpo(:)
      38             : 
      39     1489176 :       end subroutine jpl
      40             : 
      41             :       end module mo_jpl