|
Last change
on this file since 3944 was
3585,
checked in by forkel, 6 years ago
|
|
Comments in all chem_gasphase_mod.kpp, add def_salsagas/*.f90, removed executables
|
-
Property svn:executable set to
*
|
|
File size:
1.5 KB
|
| Rev | Line | |
|---|
| [2696] | 1 | //chem_gasphase_mod.kpp |
|---|
| 2 | // |
|---|
| [3298] | 3 | // |
|---|
| [2696] | 4 | //Former revisions |
|---|
| 5 | //---------------- |
|---|
| 6 | // $Id: chem_gasphase_mod.kpp 2459 2017-09-13 14:10:33Z forkel $ |
|---|
| [3585] | 7 | // initial revision branch chemistry rev 2459 |
|---|
| [2696] | 8 | // |
|---|
| 9 | #include smog.spc |
|---|
| 10 | #include smog.eqn |
|---|
| 11 | #INTEGRATOR rosenbrock |
|---|
| 12 | #LANGUAGE Fortran90 |
|---|
| 13 | #HESSIAN on |
|---|
| 14 | #STOICMAT on |
|---|
| 15 | // |
|---|
| 16 | // ******************************************************************************************* |
|---|
| 17 | // *** adapt the lines below occurding to the photolysis reactions of your mechanism * |
|---|
| 18 | // *** adapt the number of photolysis frequencies NPHO * |
|---|
| 19 | // *** adapt/extend the indices in the INTEGER, PARAMETER,PUBLIC statement below * |
|---|
| 20 | // *** adapt/extend PHOT_NAMES: Note that the order of PHOT_NAMES and the indices must match * |
|---|
| 21 | // ******************************************************************************************* |
|---|
| 22 | // |
|---|
| 23 | #INLINE F90_DATA |
|---|
| 24 | ! INLINE F90_DATA: Declaration of global variables for photolysis |
|---|
| 25 | ! REAL(kind=dp) :: phot(nphot) must eventually be moved to GLOBAL later for vector version |
|---|
| 26 | INTEGER, PARAMETER :: nphot = 2 |
|---|
| 27 | ! phot Photolysis frequencies |
|---|
| 28 | REAL(kind=dp) :: phot(nphot) |
|---|
| 29 | |
|---|
| 30 | INTEGER, PARAMETER,PUBLIC :: j_no2 = 1 |
|---|
| 31 | INTEGER, PARAMETER,PUBLIC :: j_rcho = 2 |
|---|
| 32 | |
|---|
| 33 | CHARACTER(LEN=15), PARAMETER, DIMENSION(NPHOT) :: phot_names = (/ & |
|---|
| 34 | 'J_NO2 ','J_RCHO '/) |
|---|
| 35 | #ENDINLINE |
|---|
| 36 | |
|---|
| [3298] | 37 | #INLINE F90_INIT |
|---|
| 38 | fix(indf_h2o) = qvap |
|---|
| 39 | fix(indf_o2) = 0.2e+6_dp * fakt |
|---|
| 40 | fix(indf_co2) = 400.0_dp * fakt |
|---|
| 41 | #ENDINLINE |
|---|
Note: See
TracBrowser
for help on using the repository browser.
