chem/KPP/kpp/kpp-2.1/util/UserRateLaws.m

   1 %  User-defined Rate Law functions
   2 %  Note: insert this file at the end of Update_RCONST
   3
   4 %---  Arrhenius
   5    function [rate] =  ARR( A0,B0,C0 )
   6       global TEMP CFACTOR
   7       rate =  (A0) * exp(-(B0)/TEMP) * (TEMP/300.0)^(C0) ;
   8    return %  ARR
   9
  10 %--- Simplified Arrhenius, with two arguments
  11 %--- Note: The argument B0 has a changed sign when compared to ARR
  12    function [rate] =  ARR2( A0,B0 )
  13       global TEMP CFACTOR
  14       rate =  (A0) * exp( (B0)/TEMP ) ;
  15    return %  ARR2
  16
  17    function [rate] =  EP2(A0,C0,A2,C2,A3,C3)
  18       global TEMP CFACTOR
  19       K0 = (A0) * exp(-C0/TEMP);
  20       K2 = (A2) * exp(-C2/TEMP);
  21       K3 = (A3) * exp(-C3/TEMP);
  22       K3 = K3*CFACTOR*1.0e+6;
  23       rate = K0 + K3/(1.0+K3/K2) ;
  24    return %  EP2
  25
  26    function [rate] =  EP3(A1,C1,A2,C2)
  27       global TEMP CFACTOR
  28       K1 = (A1) * exp(-(C1)/TEMP);
  29       K2 = (A2) * exp(-(C2)/TEMP);
  30       rate = K1 + K2*(1.0e+6*CFACTOR);
  31    return %  EP3
  32
  33    function [rate] =  FALL ( A0,B0,C0,A1,B1,C1,CF)
  34       global TEMP CFACTOR
  35       K0 = A0 * exp(-B0/TEMP)* (TEMP/300.0)^(C0);
  36       K1 = A1 * exp(-B1/TEMP)* (TEMP/300.0)^(C1);
  37       K0 = K0*CFACTOR*1.0e+6;
  38       K1 = K0/K1;
  39       rate = (K0/(1.0+K1))*(CF)^(1.0/(1.0+(log(K1))^2));
  40    return %  FALL
  41