Last change
on this file since 4075 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 |
---|
5 | S1 = ignore; {} |
---|
6 | S2 = ignore; {} |
---|
7 | S3 = ignore; {} |
---|
8 | |
---|
9 | #EQUATIONS |
---|
10 | S1 = PROD : (1); {c1} |
---|
11 | S1 + S1 = S2 : (0.004);{2*c2} |
---|
12 | S2 = S1 + S1 : (0.5); {c3} |
---|
13 | S2 = 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.