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

   1 #include strato.spc
   2 #include strato.eqn
   3
   4 {#LANGUAGE Fortran77
   5 #INTEGRATOR rodas3
   6 #DRIVER general}
   7
   8 #LOOKATALL
   9
  10 #MONITOR O;
  11
  12 #INITVALUES
  13
  14 CFACTOR =   8.12E+16 ;
  15
  16 O       =    8.158574E-09 ;
  17 O1D     =    2.773328E-15 ;
  18 O3      =    6.560000E-06 ;
  19 NO      =    1.074516E-08 ;
  20 NO2     =    2.758840E-09 ;
  21 NO3     =    5.667175E-14 ;
  22 N2O5    =    6.905712E-12 ;
  23 HNO3    =    3.486344E-10 ;
  24 HNO4    =    3.827485E-12 ;
  25 H2O     =    6.100000E-06 ;
  26 H       =    4.190815E-14 ;
  27 OH      =    2.041165E-10 ;
  28 HO2     =    1.439703E-10 ;
  29 H2O2    =    7.227409E-11 ;
  30 H2      =    3.700000E-07 ;
  31 CH4     =    4.900000E-07 ;
  32 CH3     =    3.984427E-17 ;
  33 CH3O    =    2.496504E-15 ;
  34 CH3O2   =    5.130799E-12 ;
  35 CH3OOH  =    4.587317E-12 ;
  36 HCO     =    3.010774E-19 ;
  37 CH2O    =    9.913988E-11 ;
  38 CO      =    2.000000E-08 ;
  39 Cl      =    8.756027E-12 ;
  40 Cl2     =    8.665679E-14 ;
  41 ClO     =    1.011874E-09 ;
  42 ClOO    =    2.659306E-15 ;
  43 OClO    =    7.193558E-14 ;
  44 Cl2O2   =    7.103997E-14 ;
  45 HCl     =    2.157657E-09 ;
  46 HOCl    =    2.196369E-10 ;
  47 ClONO2  =    2.166934E-11 ;
  48 Br      =    4.025093E-12 ;
  49 BrCl    =    5.843046E-14 ;
  50 BrO     =    1.132088E-11 ;
  51 HBr     =    4.349534E-13 ;
  52 HOBr    =    5.066596E-12 ;
  53 BrONO2  =    9.405408E-14 ;
  54 O2      =    0.209;
  55 N2      =    0.781;
  56
  57 #INLINE F77_INIT
  58         TSTART = 12.D0*3600.D0
  59         TEND = TSTART + 24.D0*3600.D0 * 3
  60         DT = 3600.D0
  61         TEMP = 241.42
  62 #ENDINLINE
  63
  64 #INLINE F90_INIT
  65         TSTART = 12.D0*3600.D0
  66         TEND = TSTART + 24.D0*3600.D0 * 3
  67         DT = 3600.D0
  68         TEMP = 241.42
  69 #ENDINLINE
  70
  71 #INLINE MATLAB_INIT
  72    global TSTART TEND DT TEMP
  73    TSTART = 12.0*3600.0;
  74    TEND = TSTART + 24.0*3600.0 * 3;
  75    DT = 3600.0;
  76    TEMP = 241.42;
  77 #ENDINLINE
  78
  79 #INLINE C_INIT
  80    TSTART = 12.0*3600.0;
  81    TEND = TSTART + 24.0*3600.0 * 3;
  82    DT = 3600.0;
  83    TEMP = 241.42;
  84 #ENDINLINE