SoilInfiltrationGreenAmpt
(*args, **kwds)[source]¶Bases: landlab.core.model_component.Component
Infiltrate surface water into a soil following the GreenAmpt method.
This code is based on an overland flow model by Francis Rengers and colleagues, after Julien et al., 1995. The infiltration scheme follows the Green and Ampt equation.
It was implemented in Landlab by DEJH, March ‘16. Please cite Rengers et al., 2016, Model Predictions of Water Runoff in Steep Catchments after Wildfire, WRR.
Construction:
SoilInfiltrationGreenAmpt(grid, hydraulic_conductivity=0.005,
soil_bulk_density=1590, rock_density=2650,
initial_soil_moisture_content=0.15,
soil_type='sandy loam',
volume_fraction_coarse_fragments=0.2,
coarse_sed_flag=False,
surface_water_minimum_depth=1.e8,
soil_pore_size_distribution_index=None,
soil_bubbling_pressure=None,
wetting_front_capillary_pressure_head=None)
Parameters:  grid : RasterModelGrid
hydraulic_conductivity : float, array, or field name (m/s)
soil_bulk_density : float (kg/m**3)
rock_density : float (kg/m**3)
initial_soil_moisture_content : float (m**3/m**3, 0. to 1.)
soil_type : {‘sand’, loamy sand’, ‘sandy loam’, ‘loam’, ‘silt loam’,
volume_fraction_coarse_fragments : float (m**3/m**3, 0. to 1.)
coarse_sed_flag : boolean, optional
surface_water_minimum_depth : float (m), optional
soil_pore_size_distribution_index : float, optional
soil_bubbling_pressure : float (m), optional
wetting_front_capillary_pressure_head : float (m), optional


Examples
>>> from landlab import RasterModelGrid
>>> from landlab.components import SoilInfiltrationGreenAmpt
>>> mg = RasterModelGrid((4,3), spacing=10.)
>>> hydraulic_conductivity = mg.ones('node')*1.e6
>>> hydraulic_conductivity.reshape((4,3))[0:2,:] *= 10000.
>>> h = mg.add_ones('node', 'surface_water__depth')
>>> h *= 0.01
>>> d = mg.add_ones('node', 'soil_water_infiltration__depth', dtype=float)
>>> d *= 0.2
>>> SI = SoilInfiltrationGreenAmpt(
... mg,hydraulic_conductivity=hydraulic_conductivity)
>>> for i in range(10): # 100s total
... SI.run_one_step(10.)
>>> mg.at_node['surface_water__depth']
array([ 1.00000000e08, 1.00000000e08, 1.00000000e08,
1.00000000e08, 1.00000000e08, 1.00000000e08,
9.88530416e03, 9.88530416e03, 9.88530416e03,
9.88530416e03, 9.88530416e03, 9.88530416e03])
>>> mg.at_node['soil_water_infiltration__depth']
array([ 0.20999999, 0.20999999, 0.20999999, 0.20999999, 0.20999999,
0.20999999, 0.2001147 , 0.2001147 , 0.2001147 , 0.2001147 ,
0.2001147 , 0.2001147 ])