def __init__(self, wrm=None, coef=1, dT=1.):
if wrm is None:
wrm = wrm_vangenuchten.VanGenuchten(0.8, 1.5e-4)
self.smooth_coef = coef
self.dp = capillary_pressure.pc_ice(273.15-dT, 101325., 1)
print self.dp
self.pm = permafrost_model_explicit_fpd.PermafrostModel(wrm)
def __init__(self, wrm=None, coef=1, dT=1.):
if wrm is None:
wrm = wrm_vangenuchten.VanGenuchten(0.8, 1.5e-4)
self.smooth_coef = coef
self.dp = capillary_pressure.pc_ice(273.15 - dT, 101325., 1)
print self.dp
self.pm = permafrost_model_explicit_fpd.PermafrostModel(wrm)
def wc(self, T, p):
pcliq = capillary_pressure.pc_liq(T, p)
pcice = capillary_pressure.pc_ice(T, p)
sg, sl, si = self.pm.saturations_Tp(T, p)
ni = eos.n_ice(T, p)
nl = eos.n_water(T, p)
ng = eos.n_gas(T, p)
om = eos.mol_frac_gas(T)
return sg * om * ng + sl * nl + si * ni
def wc(self, T, p):
pcliq = capillary_pressure.pc_liq(T,p)
pcice = capillary_pressure.pc_ice(T,p)
sg,sl,si = self.pm.saturations_Tp(T,p)
ni = eos.n_ice(T,p)
nl = eos.n_water(T,p)
ng = eos.n_gas(T,p)
om = eos.mol_frac_gas(T)
return sg*om*ng + sl*nl + si*ni
def __init__(self, wrm=None, coef=1., dT=1.):
if wrm is None:
wrm = wrm_vangenuchten.VanGenuchten(0.8, 1.5e-4)
self.wrm = wrm
self.coef = coef
self.dp = capillary_pressure.pc_ice(273.15 - dT, 101325., coef)
def saturations_Tp(self, T, p):
pcliq = capillary_pressure.pc_liq(T, p)
pcice = capillary_pressure.pc_ice(T, p, self.coef)
return self.saturations_pc(pcliq, pcice)
def dsaturations_dpcice_Tp(self, T, p):
pcliq = capillary_pressure.pc_liq(T,p)
pcice = capillary_pressure.pc_ice(T,p)
return self.dsaturations_dpcice_pc(pcliq, pcice)
def saturations_Tp(self, T, p):
pcliq = capillary_pressure.pc_liq(T,p)
pcice = capillary_pressure.pc_ice(T,p)
return self.saturations_pc(pcliq,pcice)
def __init__(self, wrm=None, coef=1., dT=1.):
if wrm is None:
wrm = wrm_vangenuchten.VanGenuchten(0.8, 1.5e-4)
self.wrm = wrm
self.coef = coef
self.dp = capillary_pressure.pc_ice(273.15-dT, 101325., coef)
