Perparticle sstp deposition

bindings/python/lib.cpp

@@ -269,6 +269,7 @@ BOOST_PYTHON_MODULE(libcloudphxx)
       .def_readwrite("sedi", &lgr::opts_t<real_t>::sedi)
       .def_readwrite("subs", &lgr::opts_t<real_t>::subs)
       .def_readwrite("cond", &lgr::opts_t<real_t>::cond)
+      .def_readwrite("depo", &lgr::opts_t<real_t>::depo)
       .def_readwrite("coal", &lgr::opts_t<real_t>::coal)
       .def_readwrite("src", &lgr::opts_t<real_t>::src)
       .def_readwrite("rcyc", &lgr::opts_t<real_t>::rcyc)

include/libcloudph++/lgrngn/opts.hpp

@@ -20,7 +20,7 @@ namespace libcloudphxx
     struct opts_t 
     {
       // process toggling
-      bool adve, sedi, subs, cond, coal, src, rlx, rcyc, turb_adve, turb_cond, turb_coal, ice_nucl;
+      bool adve, sedi, subs, cond, coal, src, rlx, rcyc, turb_adve, turb_cond, turb_coal, ice_nucl, depo;
 
       // RH limit for drop growth
       real_t RH_max;       
@@ -42,7 +42,7 @@ namespace libcloudphxx
       opts_t() : 
         adve(true), sedi(true), subs(false), cond(true), coal(true), src(false), rlx(false), rcyc(false),
         chem_dsl(false), chem_dsc(false), chem_rct(false),
-        turb_adve(false), turb_cond(false), turb_coal(false), ice_nucl(false),
+        turb_adve(false), turb_cond(false), turb_coal(false), ice_nucl(false), depo(false),
         RH_max(44), // :) (anything greater than 1.1 would be enough
         dt(-1) // negative means that we do not override dt in this step
       {

include/libcloudph++/lgrngn/opts_init.hpp

@@ -81,7 +81,7 @@ namespace libcloudphxx
       // super-droplet advection scheme
       as_t adve_scheme;
 
-      // RH formula
+      // RH formula, for ice, pv_tet and rv_tet work as pv_cc and rv_cc, respectively, because there is no Tetens formula for ice
       RH_formula_t RH_formula;
 //</listing>