------------------------------ CHANGES TO .mmp FILE ------------------------------ $CAPILLARY_PRESSURE ------------------------------ --> OLD FORMAT ; pc_model -- pb (or alpha for van_genuchten) 6 5000.0 !Problems with this: 1) Pc cannot be independent of relative permeability values. Actually, they should be totally independent. --> NEW FORMAT ; pc_model -- pb -- slr -- slm -- max_pc -- alpha_switch (if alpha_switch>0, Pb is taken as "alpha", where Pb = rho*g/alpha) 6 5000.0 0.0 1.0 1.0e12 0 !Notes: Pc is no longer linked to relative permeability. It carries it's own values. $PERMEABILITY_SATURATION ------------------------------ --> OLD FORMAT ; model -- min_sat -- max_sat -- exponent -- min_kr 6 0.0 1.0 2.0 44 0.0 1.0 0.5 1.0e-9 !Problems with this (these aren't big problems, but since we're reconfiguring Pc and k_rel anyway, may as well make the optimization): 1) Each phase should have k_rel entries that are independent of other phases. --> NEW FORMAT ; model -- slr -- sgr -- exponent -- min_kr 6 0.25 0.025 2.0 1.0e-9 4 0.25 0.025 2.0 1.0e-9 ! Notes. 1) Phases are totally independent. Each has it's own slr, sgr, exponent, minimum_k. 2) Exponents are the same for van Genuchten or Brooks-Corey (i.e., always greater than 1.0. thus for van Genuchten, m = 1 - 1/exponent) 3) No need to distinguish phase specific models.