chem/KPP/kpp/kpp-2.1/test/wet.def

   1 #include wet.spc
   2 #include wet.eqn
   3
   4 {#LANGUAGE Fortran77
   5 #INTEGRATOR rodas3
   6 #DRIVER general}
   7
   8 #LOOKATALL
   9
  10 #MONITOR O3;
  11
  12 #INITVALUES
  13   CFACTOR = 2.55E+10; {ppb-to-mcm}
  14   ALL_SPEC= 1.0E-10;
  15 {Variable species}
  16   NO2        = 0.5 ;
  17   NO         = 0.2 ;
  18   O3         = 60.0 ;
  19   HNO3       = 1.6594E-06 ;
  20   HNO2       = 8.4898E-03 ;
  21   H2O2       = 1.5 ;
  22   CO         = 200 ;
  23   HCHO       = 1.0 ;
  24   ALD2       = 0.1 {3.2067E-02} ;
  25   PAN        = 0.1 {1.5022E-02} ;
  26   MGLY       = 7.3500E-03 ;
  27   PAR        = 1.2 ;
  28   OLE        = 6.6095E-03 ;
  29   ETH        = 1.0 {2.4286E-02} ;
  30   TOL        = 3.0184E-02 ;
  31   CRES       = 6.4171E-03 ;
  32   XYL        = 2.0E-02 ;
  33   ROOH       = 7.0E-02 ;
  34   SO2        = 3.3371E-02 ;
  35   PNA        = 1.6736E-02 ;
  36   OPEN       = 1.0E-03 ;
  37   ISOP       = 1.0 ;
  38   HONO       = 1.0E-3 ;
  39   NO3        = 5.0E-05 ;
  40   OH         = 6.0E-04 ;
  41   HO2        = 2.8900E-02 ;
  42   N2O5       = 4.6513E-05 ;
  43   C2O3       = 6.1457E-05 ;
  44   XO2        = 6.2020E-03 ;
  45   XO2N       = 4.0E-3 ;
  46   ROR        = 4.2461E-09 ;
  47   TO2        = 8.6120E-07 ;
  48   NH3        = 1.3637E-02 ;
  49   CRO        = 3.1320E-06 ;
  50   H2O        = 1.35e+7 ;
  51   O2         = 2.1E+8 ;
  52   CH4        = 1.7E+3 ;
  53   CO2aq      = 3.0E-1 ;
  54   HO2aq      = 4.2400E-04 ;
  55   H2O2aq     = 9.6683E-01 ;
  56   O3aq       = 4.8688E-06 ;
  57   HCHOaq     = 2.0245E-02 ;
  58   HCOOH      = 0.2 ;
  59   HCOOHaq    = 2.0180E-03 ;
  60   CH3O2Haq   = 1.1101E-04 ;
  61   CH3O2aq    = 9.0994E-05 ;
  62   HNO3aq     = 2.5563E-06 ;
  63   NO3aq      = 5.2533E-09 ;
  64   NH3aq      = 5.8021E-06 ;
  65   SO2aq      = 3.0071E-07 ;
  66   OHaq       = 4.0141E-05 ;
  67   H2Oaq      = 4.3638E+05 ;
  68   HPLUSaq    = 1.9654E-01 ;
  69   OHaq       = 3.1392E-06 ;
  70   HCOOMINaq  = 3.6325E-03 ;
  71   O2MINaq    = 6.7840E-05 ;
  72   HSO3MINaq  = 5.1188E-05 ;
  73   SO32MINaq  = 3.4183E-08 ;
  74   NH4aq      = 9.8636E-01 ;
  75   NO3aq      = 3.9367E-01 ;
  76   FE3PLUSaq  = 7.8548E-03 ;
  77   FE2PLUSaq  = 3.1419E-02 ;
  78   SO42MINaq  = 2.3564E+00 ;
  79   CU2PLUSaq  = 2.1601E-04 ;
  80
  81 #INLINE F77_INIT
  82         TSTART = 12.D0*3600.D0
  83         TEND = TSTART + 24.D0*3600.D0 * 5
  84         DT = 3600.D0
  85         TEMP = 288.15
  86 #ENDINLINE
  87
  88 #INLINE F90_INIT
  89         TSTART = 12.D0*3600.D0
  90         TEND = TSTART + 24.D0*3600.D0 * 5
  91         DT = 3600.D0
  92         TEMP = 288.15
  93 #ENDINLINE
  94
  95 #INLINE MATLAB_INIT
  96   global TSTART TEND DT TEMP
  97   TSTART = 12.0*3600.0;
  98   TEND = TSTART + 24.0*3600.0*5;
  99   DT = 3600.0;
 100   TEMP = 288.15;
 101 #ENDINLINE
 102
 103 #INLINE C_INIT
 104   TSTART = 12.0*3600.0;
 105   TEND = TSTART + 24.0*3600.0*5;
 106   DT = 3600.0;
 107   TEMP = 288.15;
 108 #ENDINLINE
 109