source: palm/trunk/UTIL/chemistry/gasphase_preproc/kpp/models/stochastic_dimer.def @ 2898

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

Merge of branch palm4u into trunk

File size: 1.0 KB
Line 
1{Decaying-dimerizing reaction set
2 (Gillespie, J. Chem. Phys. 115(4), p. 1716, 2001)}
3
4#DEFVAR
5S1  = ignore; {}
6S2  = ignore; {}
7S3  = ignore; {}
8
9#EQUATIONS
10S1        = PROD        : (1);    {c1}
11S1 + S1   = S2          : (0.004);{2*c2}
12S2        = S1 + S1     : (0.5);  {c3}
13S2        = S3          : (0.04); {c4}
14
15#LOOKATALL
16
17#INITVALUES
18  CFACTOR = 1.0;
19  S1   = 1.0e+5;
20  S2   = 0.0;
21  S3   = 0.0;
22
23#INLINE F77_INIT
24      TSTART = 0.d0
25      TEND   = 30.0d0
26      DT     = 0.5d0;
27      DO i=1,NVAR
28        RTOL(i) = 1.0e-4
29        ATOL(i) = 1.0e-8
30      END DO
31      Volume = 1.0d0
32#ENDINLINE
33
34#INLINE F90_INIT
35  TSTART = 0.d0
36  TEND   = 30.0d0
37  DT     = 0.5d0
38  RTOL(1:NVAR) = 1.0e-4
39  ATOL(1:NVAR) = 1.0e-8
40  Volume = 1.0d0
41#ENDINLINE
42
43#INLINE MATLAB_INIT
44  global TSTART TEND DT TEMP
45  TSTART = 0;
46  TEND   = 30;
47  DT     = 0.5;
48  RTOL(1:NVAR) = 1.0e-4;
49  ATOL(1:NVAR) = 1.0e-8;
50  Volume = 1;
51#ENDINLINE
52
53#INLINE C_INIT
54  TSTART = 0.0;
55  TEND   = 30.0;
56  DT     = 0.5;
57  for(i=0; i<NVAR; i++) {
58      RTOL[i] = 1.0e-4;
59      ATOL[i] = 1.0e-8;
60  }
61  Volume = 1.0;
62#ENDINLINE
Note: See TracBrowser for help on using the repository browser.