output_streams=streams_inbasin,
output_basin=basins_inbasin,
x_outlet=Settings.outlet_point_x,
y_outlet=Settings.outlet_point_y,
output_pour_point=pour_point,
overwrite=True)
# GSFLOW segments: sections of stream that define subbasins
v.gsflow_segments(input=streams_inbasin,
output=segments,
icalc=Settings.icalc,
roughch_value=Settings.channel_Mannings_n,
roughch_raster=Settings.channel_Mannings_n_grid,
roughch_points=Settings.channel_Mannings_n_vector,
roughch_pt_col=Settings.channel_Mannings_n_vector_col,
width1=Settings.channel_width,
width2=Settings.channel_width,
width_points=Settings.channel_width_vector,
width_points_col=Settings.channel_width_vector_col,
fp_width_value=Settings.floodplain_width,
fp_width_pts=Settings.floodplain_width_vector,
fp_width_pts_col=Settings.floodplain_width_vector_col,
overwrite=True)
# MODFLOW grid & basin mask (1s where basin exists and 0 where it doesn't)
# Fill nulls in case of ocean
# Any error-related NULL cells will not be part of the basin, and all cells
# should have elevation > 0, so this hopefully will not cause any problems
r.null(map=DEM, null=0)
v.gsflow_grid(basin=basins_inbasin,
pour_point=pour_point,
# Build stream network
v.stream_network(map=streams_all)
# Restrict to a single basin
v.stream_inbasin(input_streams=streams_all,
input_basins=basins_all,
output_streams=streams_inbasin,
output_basin=basins_inbasin,
x_outlet=outlet_point_x,
y_outlet=outlet_point_y,
output_pour_point=pour_point,
overwrite=True)
# GSFLOW segments: sections of stream that define subbasins
v.gsflow_segments(input=streams_inbasin,
output=segments,
icalc=icalc,
overwrite=True)
# MODFLOW grid & basin mask (1s where basin exists and 0 where it doesn't)
v.gsflow_grid(basin=basins_inbasin,
pour_point=pour_point,
raster_input=DEM,
dx=MODFLOW_grid_resolution,
dy=MODFLOW_grid_resolution,
output=MODFLOW_grid,
mask_output=basin_mask,
bc_cell=bc_cell,
overwrite=True)
# Hydrologically-correct DEM for MODFLOW
r.gsflow_hydrodem(dem=DEM,