'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 *