v.gsflow.reaches

Build stream "reaches" that link PRMS segments to MODFLOW cells

v.gsflow.reaches segment_input=name grid_input=name elevation=name output=name [s_min=float] [h_stream=float] [upstream_easting_column_seg=string] [upstream_northing_column_seg=string] [downstream_easting_column_seg=string] [downstream_northing_column_seg=string] [tostream_cat_column_seg=string] [strthick=float] [strhc1=float] [thts=float] [thti=float] [eps=float] [uhc=float] [--overwrite] [--verbose] [--quiet] [--qq] [--ui]

Example:

v.gsflow.reaches segment_input=name grid_input=name elevation=name output=name

grass.script.run_command("v.gsflow.reaches", segment_input, grid_input, elevation, output, s_min=0.0001, h_stream=1, upstream_easting_column_seg="x1", upstream_northing_column_seg="y1", downstream_easting_column_seg="x2", downstream_northing_column_seg="y2", tostream_cat_column_seg="tostream", strthick=1, strhc1=5, thts=0.35, thti=0.3, eps=3.5, uhc=0.3, overwrite=False, verbose=False, quiet=False, superquiet=False)

Example:

gs.run_command("v.gsflow.reaches", segment_input="name", grid_input="name", elevation="name", output="name")

Parameters

Command linePython (grass.script)

segment_input=name [required]
    PRMS stream segments
    Or data source for direct OGR access
grid_input=name [required]
    MODFLOW grid
    Or data source for direct OGR access
elevation=name [required]
    DEM for slope along reaches
    Name of input raster map
output=name [required]
    Reaches for GSFLOW
    Name for output vector map
s_min=float
    Minimum reach slope
    Default: 0.0001
h_stream=float
    Stream channel depth (bank height) [m]
    Default: 1
upstream_easting_column_seg=string
    Upstream easting (or x or lon) column name
    Default: x1
upstream_northing_column_seg=string
    Upstream northing (or y or lat) column name
    Default: y1
downstream_easting_column_seg=string
    Downstream easting (or x or lon) column name
    Default: x2
downstream_northing_column_seg=string
    Downstream northing (or y or lat) column name
    Default: y2
tostream_cat_column_seg=string
    Adjacent downstream segment cat (0 = offmap)
    Default: tostream
strthick=float
    Streambed sediment thickness [m]
    Default: 1
strhc1=float
    Streambed hydraulic conductivity [m/day]
    Default: 5
thts=float
    theta_sat: Streambed saturated water content (i.e. porosity) [unitless]
    Default: 0.35
thti=float
    Streambed initial water content [unitless]
    Default: 0.3
eps=float
    Epsilon: streambed Brooks-Corey exponent [unitless]
    Default: 3.5
uhc=float
    Streambed unsaturated zone saturated hydraulic conductivity [m/day]
    Default: 0.3
--overwrite
    Allow output files to overwrite existing files
--help
    Print usage summary
--verbose
    Verbose module output
--quiet
    Quiet module output
--qq
    Very quiet module output
--ui
    Force launching GUI dialog

segment_input : str, required
    PRMS stream segments
    Or data source for direct OGR access
    Used as: input, vector, name
grid_input : str, required
    MODFLOW grid
    Or data source for direct OGR access
    Used as: input, vector, name
elevation : str, required
    DEM for slope along reaches
    Name of input raster map
    Used as: input, raster, name
output : str, required
    Reaches for GSFLOW
    Name for output vector map
    Used as: output, vector, name
s_min : float, optional
    Minimum reach slope
    Default: 0.0001
h_stream : float, optional
    Stream channel depth (bank height) [m]
    Default: 1
upstream_easting_column_seg : str, optional
    Upstream easting (or x or lon) column name
    Default: x1
upstream_northing_column_seg : str, optional
    Upstream northing (or y or lat) column name
    Default: y1
downstream_easting_column_seg : str, optional
    Downstream easting (or x or lon) column name
    Default: x2
downstream_northing_column_seg : str, optional
    Downstream northing (or y or lat) column name
    Default: y2
tostream_cat_column_seg : str, optional
    Adjacent downstream segment cat (0 = offmap)
    Default: tostream
strthick : float, optional
    Streambed sediment thickness [m]
    Default: 1
strhc1 : float, optional
    Streambed hydraulic conductivity [m/day]
    Default: 5
thts : float, optional
    theta_sat: Streambed saturated water content (i.e. porosity) [unitless]
    Default: 0.35
thti : float, optional
    Streambed initial water content [unitless]
    Default: 0.3
eps : float, optional
    Epsilon: streambed Brooks-Corey exponent [unitless]
    Default: 3.5
uhc : float, optional
    Streambed unsaturated zone saturated hydraulic conductivity [m/day]
    Default: 0.3
overwrite: bool, optional
    Allow output files to overwrite existing files
    Default: False
verbose: bool, optional
    Verbose module output
    Default: False
quiet: bool, optional
    Quiet module output
    Default: False
superquiet: bool, optional
    Very quiet module output
    Default: False

DESCRIPTION

v.gsflow.reaches produces channel "reaches", which are the intersection between PRMS stream segments and MODFLOW grid cells, for the USGS combined groundwater and surface-water model GSFLOW.

REFERENCES

Ng, G.-H. C., A. D. Wickert, R. L. McLaughlin, J. La Frenierre, S. Liess, and L. Sabeeri (2016), Modeling the role of groundwater and vegetation in the hydrological response of tropical glaciated watersheds to climate change, in AGU Fall Meeting Abstracts, H13L–1590, San Francisco, CA.

AUTHOR

Andrew D. Wickert

SOURCE CODE

Available at: v.gsflow.reaches source code (history)
Latest change: Thursday Feb 20 13:02:26 2025 in commit 53de819