def test_instantiate( soil_type='sandy_loam' ):
soil = soil_base.soil_base( soil_type=soil_type )
soil.CCA = False
soil.DEBUG = False
#-----------------------------------------
# This function adjusts for the platform
# and can be changed in "tf_utils.py".
#-----------------------------------------
cfg_directory = tf_utils.TF_Test_Directory()
os.chdir( cfg_directory )
cfg_prefix = tf_utils.TF_Test_Case_Prefix()
soil.site_prefix = 'Treynor' #############
soil.initialize( cfg_prefix=cfg_prefix, REPORT=True )
print 'soil_type =', c.soil_type
print 'Note: t_r < t_H < t_w < t_f < t_i < t_u < t_s'
print ' '
print 'theta_r =', soil.theta_r # ("absolute" min)
print 'theta_H =', soil.theta_H # (achievable min)
print 'theta_w =', soil.theta_w
print 'theta_f =', soil.theta_f
print 'theta_i =', soil.theta_i
print 'theta_u =', soil.theta_u # (max possible)
print 'theta_s =', soil.theta_s
print ' '
print 'K_s =', soil.K_s
print 'K(theta_s) =', soil.K_of_theta(soil.theta_s)
print ' '
print 'K_i =', soil.K_i
print 'K(theta_i) =', soil.K_of_theta(soil.theta_i)
print ' '
def test_instantiate(soil_type='sandy_loam'):
soil = soil_base.soil_base(soil_type=soil_type)
soil.CCA = False
soil.DEBUG = False
#-----------------------------------------
# This function adjusts for the platform
# and can be changed in "tf_utils.py".
#-----------------------------------------
cfg_directory = tf_utils.TF_Test_Directory()
os.chdir(cfg_directory)
cfg_prefix = tf_utils.TF_Test_Case_Prefix()
soil.site_prefix = 'Treynor' #############
soil.initialize(cfg_prefix=cfg_prefix, REPORT=True)
print('soil_type =', c.soil_type)
print('Note: t_r < t_H < t_w < t_f < t_i < t_u < t_s')
print(' ')
print('theta_r =', soil.theta_r) # ("absolute" min)
print('theta_H =', soil.theta_H) # (achievable min)
print('theta_w =', soil.theta_w)
print('theta_f =', soil.theta_f)
print('theta_i =', soil.theta_i)
print('theta_u =', soil.theta_u) # (max possible)
print('theta_s =', soil.theta_s)
print(' ')
print('K_s =', soil.K_s)
print('K(theta_s) =', soil.K_of_theta(soil.theta_s))
print(' ')
print('K_i =', soil.K_i)
print('K(theta_i) =', soil.K_of_theta(soil.theta_i))
print(' ')
sat.initialize(cfg_file="cfg/Treynor_satzone_darcy_layers.cfg") sat.Rg = 1 sat.update() sat.finalize() ################################################################################## from topoflow.components import soil_base # just functions called by infil_richards_1D soil = soil_base.soil_base() soil.initialize() ################################################################################## from topoflow.components import topoflow_driver import numpy as np topo = topoflow_driver.topoflow_driver() topo.Q_outlet = np.array(1) # numpy array topo.initialize(cfg_file = "cfg/Treynor_topoflow.cfg")