source: palm/trunk/UTIL/chemistry/gasphase_preproc/kpp/util/UserRateLaws.f @ 3961

Last change on this file since 3961 was 2696, checked in by kanani, 7 years ago

Merge of branch palm4u into trunk

File size: 2.1 KB
Line 
1C  User-defined Rate Law functions
2C  Note: the default argument type for rate laws, as read from the equations file, is single precision
3C        but all the internal calculations are performed in REAL*8
4
5C 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
17C Simplified Arrhenius, with two arguments
18C 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
Note: See TracBrowser for help on using the repository browser.