'bedrock__elevation'] + mg.at_node['soil__depth']
#set second boundary layer to keep hillslope from going into baselevel node
second_boundary_layer = mg.nodes_at_bottom_edge + nc
boundary_remove = np.where(second_boundary_layer == channel_nodes[0])
second_boundary_layer = np.delete(second_boundary_layer, boundary_remove)
mg.status_at_node[second_boundary_layer] = CLOSED_BOUNDARY
#inititalize resistant layer
starting_point = max(mg.node_x) - (2 * dx)
#route flow
fr = FlowDirectorD8(mg, 'topographic__elevation')
#instantiate BRaKE, the collection of functions governing river evolution
blocks = bl.BrakeLandlab()
#instantiate hillslope diffusion model component, from Landlab
hillslopeflux = DepthDependentDiffuser(
mg,
linear_diffusivity=linear_diffusivity,
soil_transport_decay_depth=soil_transport_decay_depth)
#instantiate exponential weathering component, from Landlab
weatheringrate = ExponentialWeatherer(
mg,
soil_production__maximum_rate=max_soil_production_rate,
soil_production__decay_depth=soil_production_decay_depth)
#set boundary elevations to enable easier data visualization in Paraview
mg.at_node['topographic__elevation'][mg.boundary_nodes] -= (bl_drop *
'bedrock__elevation'] + mg.at_node['soil__depth']
#set second boundary layer to keep hillslope from going into baselevel node
second_boundary_layer = mg.nodes_at_bottom_edge + nc
boundary_remove = np.where(second_boundary_layer == channel_nodes[0])
second_boundary_layer = np.delete(second_boundary_layer, boundary_remove)
mg.status_at_node[second_boundary_layer] = CLOSED_BOUNDARY
#inititalize resistant layer
starting_point = max(mg.node_x) - (2 * dx)
#route flow
fr = FlowDirectorD8(mg, 'topographic__elevation')
#instantiate BRaKE, the collection of functions governing river evolution
blocks = bl.BrakeLandlab() #create an instance of the class
#instantiate hillslope diffusion model component, from Landlab
hillslopeflux = DepthDependentDiffuser(
mg,
linear_diffusivity=linear_diffusivity,
soil_transport_decay_depth=soil_transport_decay_depth)
#instantiate exponential weathering component, from Landlab
weatheringrate = ExponentialWeatherer(
mg,
soil_production__maximum_rate=max_soil_production_rate,
soil_production__decay_depth=soil_production_decay_depth)
#set boundary elevations to enable easier data visualization in Paraview
mg.at_node['topographic__elevation'][mg.boundary_nodes] -= (bl_drop *