source: palm/trunk/UTIL/chemistry/gasphase_preproc/mechanisms/def_cbm4/chem_gasphase_mod.kpp @ 3395

Last change on this file since 3395 was 3249, checked in by forkel, 6 years ago

Modification in kpp4palm: Fix of fixed species initialisation, unnecessary Fortran variables deleted

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