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

   1 C  User-defined Rate Law functions
   2 C  Note: the default argument type for rate laws, as read from the equations file, is single precision
   3 C        but all the internal calculations are performed in REAL*8
   4
   5 C Arrhenius
   6       KPP_REAL FUNCTION ARR( A0,B0,C0 )
   7       INCLUDE 'KPP_ROOT_Parameters.h'
   8       INCLUDE 'KPP_ROOT_Global.h'
   9
  10       REAL A0,B0,C0
  11       ARR =  DBLE(A0) * EXP(-DBLE(B0)/TEMP) * (TEMP/300.0D0)**DBLE(C0)
  12
  13       RETURN
  14       END
  15
  16
  17 C Simplified Arrhenius, with two arguments
  18 C Note: The argument B0 has a changed sign when compared to ARR
  19       KPP_REAL FUNCTION ARR2( A0,B0 )
  20       INCLUDE 'KPP_ROOT_Parameters.h'
  21       INCLUDE 'KPP_ROOT_Global.h'
  22
  23       REAL A0,B0
  24       ARR2 =  DBLE(A0) * EXP( DBLE(B0)/TEMP )
  25
  26       RETURN
  27       END
  28
  29       KPP_REAL FUNCTION EP2(A0,C0,A2,C2,A3,C3)
  30       INCLUDE 'KPP_ROOT_Parameters.h'
  31       INCLUDE 'KPP_ROOT_Global.h'
  32
  33       REAL A0,C0,A2,C2,A3,C3
  34       REAL*8 K0,K2,K3
  35
  36       K0 = DBLE(A0) * EXP(-DBLE(C0)/TEMP)
  37       K2 = DBLE(A2) * EXP(-DBLE(C2)/TEMP)
  38       K3 = DBLE(A3) * EXP(-DBLE(C3)/TEMP)
  39       K3 = K3*CFACTOR*1.0d6
  40       EP2 = K0 + K3/(1.0d0+K3/K2 )
  41
  42       RETURN
  43       END
  44
  45
  46       KPP_REAL FUNCTION EP3(A1,C1,A2,C2)
  47       INCLUDE 'KPP_ROOT_Parameters.h'
  48       INCLUDE 'KPP_ROOT_Global.h'
  49
  50       REAL A1, C1, A2, C2
  51       REAL*8 K1, K2
  52
  53       K1 = DBLE(A1) * EXP(-DBLE(C1)/TEMP)
  54       K2 = DBLE(A2) * EXP(-DBLE(C2)/TEMP)
  55       EP3 = K1 + K2*(1.0d6*CFACTOR)
  56
  57       RETURN
  58       END
  59
  60
  61       KPP_REAL FUNCTION FALL ( A0,B0,C0,A1,B1,C1,CF)
  62       INCLUDE 'KPP_ROOT_Parameters.h'
  63       INCLUDE 'KPP_ROOT_Global.h'
  64
  65       REAL A0,B0,C0,A1,B1,C1,CF
  66       REAL*8 K0, K1
  67
  68       K0 = DBLE(A0) * EXP(-DBLE(B0)/TEMP)* (TEMP/300.0D0)**DBLE(C0)
  69       K1 = DBLE(A1) * EXP(-DBLE(B1)/TEMP)* (TEMP/300.0D0)**DBLE(C1)
  70       K0 = K0*CFACTOR*1.0D6
  71       K1 = K0/K1
  72       FALL = (K0/(1.0d0+K1))*
  73      *      DBLE(CF)**(1.0d0/(1.0d0+(DLOG10(K1))**2))
  74
  75       RETURN
  76       END