v.gsflow_reaches(segment_input=segments,
grid_input=MODFLOW_grid,
elevation=DEM,
output=reaches,
overwrite=True)
# GSFLOW HRU parameters
r.slope_aspect(elevation=DEM,
slope=slope,
aspect=aspect,
format='percent',
zscale=0.01,
overwrite=True)
v.gsflow_hruparams(input=basins_inbasin,
elevation=DEM,
output=HRUs,
slope=slope,
aspect=aspect,
overwrite=True)
# GSFLOW gravity reservoirs
v.gsflow_gravres(hru_input=HRUs,
grid_input=MODFLOW_grid,
output=gravity_reservoirs,
overwrite=True)
# Export DEM with MODFLOW resolution
# Also export basin mask -- 1s where basin exists and 0 where it doesn't
# And make sure it is in an appropriate folder
if os.getcwd() != Settings.GIS_output_rootdir:
try:
os.makedirs(Settings.GIS_output_rootdir)
elevation=DEM,
output=reaches,
overwrite=True)
# GSFLOW HRU parameters
r.slope_aspect(elevation=DEM,
slope=slope,
aspect=aspect,
format='percent',
zscale=0.01,
overwrite=True)
v.gsflow_hruparams(input=basins_inbasin,
elevation=DEM,
output=HRUs,
slope=slope,
aspect=aspect,
cov_type=land_cover,
soil_type=soil,
overwrite=True)
# GSFLOW gravity reservoirs
v.gsflow_gravres(hru_input=HRUs,
grid_input=MODFLOW_grid,
output=gravity_reservoirs,
overwrite=True)
# Export DEM with MODFLOW resolution
# Also export basin mask -- 1s where basin exists and 0 where it doesn't
# And make sure it is in an appropriate folder
if os.getcwd() != Settings.GIS_output_rootdir:
try: