Last change
on this file since 4127 was
2696,
checked in by kanani, 7 years ago
|
Merge of branch palm4u into trunk
|
File size:
1.0 KB
|
Rev | Line | |
---|
[2696] | 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.