Last change
on this file since 4016 was
4016,
checked in by forkel, 6 years ago
|
fix in simple mechanisms: qvap connected to mechanism
|
-
Property svn:executable set to
*
|
File size:
1.5 KB
|
Line | |
---|
1 | //chem_gasphase_mod.kpp |
---|
2 | // |
---|
3 | //Former revisions |
---|
4 | //---------------- |
---|
5 | // $Id: chem_gasphase_mod.kpp 2459 2017-09-13 14:10:33Z forkel $ |
---|
6 | // initial revision from branch chemistry |
---|
7 | // |
---|
8 | // |
---|
9 | #include simplep.spc |
---|
10 | #include simplep.eqn |
---|
11 | #INTEGRATOR rosenbrock |
---|
12 | #LANGUAGE Fortran90 |
---|
13 | #HESSIAN on |
---|
14 | #STOICMAT on |
---|
15 | // |
---|
16 | // 'simple' gas phase chemistry with additional tracer named PM10 |
---|
17 | // ******************************************************************************************* |
---|
18 | // *** adapt the lines below occurding to the photolysis reactions of your mechanism * |
---|
19 | // *** adapt the number of photolysis frequencies NPHO * |
---|
20 | // *** adapt/extend the indices in the INTEGER, PARAMETER,PUBLIC statement below * |
---|
21 | // *** adapt/extend PHOT_NAMES: Note that the order of PHOT_NAMES and the indices must match * |
---|
22 | // ******************************************************************************************* |
---|
23 | // |
---|
24 | #INLINE F90_DATA |
---|
25 | ! INLINE F90_DATA: Declaration of global variables for photolysis |
---|
26 | ! REAL(kind=dp) :: phot(nphot) must eventually be moved to GLOBAL later for vector version |
---|
27 | INTEGER, PARAMETER :: nphot = 2 |
---|
28 | ! phot Photolysis frequencies |
---|
29 | REAL(kind=dp) :: phot(nphot) |
---|
30 | |
---|
31 | INTEGER, PARAMETER,PUBLIC :: j_no2 = 1 |
---|
32 | INTEGER, PARAMETER,PUBLIC :: j_o31d = 2 |
---|
33 | |
---|
34 | CHARACTER(LEN=15), PARAMETER, DIMENSION(NPHOT) :: phot_names = (/ & |
---|
35 | 'J_NO2 ','J_O31D '/) |
---|
36 | #ENDINLINE |
---|
37 | |
---|
38 | #INLINE F90_INIT |
---|
39 | fix(indf_h2o) = qvap |
---|
40 | #ENDINLINE |
---|
Note: See
TracBrowser
for help on using the repository browser.