| 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 branch chemistry rev 3090 |
|---|
| 7 | // |
|---|
| 8 | #include cbm4.spc |
|---|
| 9 | #include cbm4.eqn |
|---|
| 10 | #INTEGRATOR rosenbrock |
|---|
| 11 | #LANGUAGE Fortran90 |
|---|
| 12 | #HESSIAN on |
|---|
| 13 | #STOICMAT on |
|---|
| 14 | // |
|---|
| 15 | #INLINE F90_GLOBAL |
|---|
| 16 | ! QVAP - Water vapor |
|---|
| 17 | REAL(kind=dp) :: QVAP |
|---|
| 18 | ! FAKT - Conversion factor |
|---|
| 19 | REAL(kind=dp) :: FAKT |
|---|
| 20 | |
|---|
| 21 | ! Declaration of global variable declarations for photolysis will come from INLINE F90_DATA |
|---|
| 22 | #ENDINLINE |
|---|
| 23 | // |
|---|
| 24 | // ******************************************************************************************* |
|---|
| 25 | // *** adapt the lines below occurding to the photolysis reactions of your mechanism * |
|---|
| 26 | // *** adapt the number of photolysis frequencies NPHO * |
|---|
| 27 | // *** adapt/extend the indices in the INTEGER, PARAMETER,PUBLIC statement below * |
|---|
| 28 | // *** adapt/extend PHOT_NAMES: Note that the order of PHOT_NAMES and the indices must match * |
|---|
| 29 | // ******************************************************************************************* |
|---|
| 30 | // |
|---|
| 31 | #INLINE F90_DATA |
|---|
| 32 | ! INLINE F90_DATA: Declaration of global variables for photolysis |
|---|
| 33 | ! REAL(kind=dp) :: phot(nphot) must eventually be moved to GLOBAL later for vector version |
|---|
| 34 | INTEGER, PARAMETER :: nphot = 9 |
|---|
| 35 | ! phot Photolysis frequencies |
|---|
| 36 | REAL(kind=dp) :: phot(nphot) |
|---|
| 37 | |
|---|
| 38 | INTEGER, PARAMETER,PUBLIC :: j_no2 = 1 |
|---|
| 39 | INTEGER, PARAMETER,PUBLIC :: j_o33p = 2 |
|---|
| 40 | INTEGER, PARAMETER,PUBLIC :: j_o31d = 3 |
|---|
| 41 | INTEGER, PARAMETER,PUBLIC :: j_no3o = 4 |
|---|
| 42 | INTEGER, PARAMETER,PUBLIC :: j_no3o2 = 5 |
|---|
| 43 | INTEGER, PARAMETER,PUBLIC :: j_hono = 6 |
|---|
| 44 | INTEGER, PARAMETER,PUBLIC :: j_h2o2 = 7 |
|---|
| 45 | INTEGER, PARAMETER,PUBLIC :: j_ch2or = 8 |
|---|
| 46 | INTEGER, PARAMETER,PUBLIC :: j_ch2om = 9 |
|---|
| 47 | |
|---|
| 48 | CHARACTER(LEN=15), PARAMETER, DIMENSION(NPHOT) :: phot_names = (/ & |
|---|
| 49 | 'J_NO2 ','J_O33P ','J_O31D ', & |
|---|
| 50 | 'J_NO3O ','J_NO3O2 ','J_HONO ', & |
|---|
| 51 | 'J_H2O2 ','J_HCHO_B ','J_HCHO_A '/) |
|---|
| 52 | #ENDINLINE |
|---|
| 53 | |
|---|
| 54 | #INLINE F90_INIT |
|---|
| 55 | fix(indf_h2o) = qvap |
|---|
| 56 | #ENDINLINE |
|---|
| 57 | |
|---|
