//chem_gasphase_mod.kpp
//
//Former revisions
//----------------
// $Id: chem_gasphase_mod.kpp 2459 2017-09-13 14:10:33Z forkel $
//      initial revision branch chemistry rev 3090
//
#include      cbm4.spc
#include      cbm4.eqn
#INTEGRATOR rosenbrock
#LANGUAGE   Fortran90
#HESSIAN    on
#STOICMAT   on
//
// *******************************************************************************************
// *** adapt the lines below occurding to the photolysis reactions of your mechanism         *
// *** adapt the number of photolysis frequencies NPHO                                       *
// *** adapt/extend the indices in the INTEGER, PARAMETER,PUBLIC statement below             *
// *** adapt/extend PHOT_NAMES: Note that the order of PHOT_NAMES and the indices must match *
// *******************************************************************************************
//
#INLINE F90_DATA
  !   INLINE F90_DATA: Declaration of global variables for photolysis
  !   REAL(kind=dp) :: phot(nphot) must eventually be moved to GLOBAL later for vector version
  INTEGER, PARAMETER :: nphot = 9
  !   phot Photolysis frequencies 
  REAL(kind=dp) :: phot(nphot)

  INTEGER, PARAMETER,PUBLIC :: j_no2 = 1
  INTEGER, PARAMETER,PUBLIC :: j_o33p = 2
  INTEGER, PARAMETER,PUBLIC :: j_o31d = 3
  INTEGER, PARAMETER,PUBLIC :: j_no3o = 4
  INTEGER, PARAMETER,PUBLIC :: j_no3o2 = 5
  INTEGER, PARAMETER,PUBLIC :: j_hono = 6
  INTEGER, PARAMETER,PUBLIC :: j_h2o2 = 7
  INTEGER, PARAMETER,PUBLIC :: j_ch2or = 8
  INTEGER, PARAMETER,PUBLIC :: j_ch2om = 9

  CHARACTER(LEN=15), PARAMETER, DIMENSION(NPHOT) :: phot_names =  (/ &
     'J_NO2          ','J_O33P         ','J_O31D         ',          &
     'J_NO3O         ','J_NO3O2        ','J_HONO         ',          &
     'J_H2O2         ','J_HCHO_B       ','J_HCHO_A       '/)
#ENDINLINE

#INLINE F90_INIT
  fix(indf_h2o) = qvap 
#ENDINLINE