LateralSubsurfaceDischarge: Lateral Subsurface Discharge

In ecor/geotopsim: Analysis of GEOtop Output Data

Description

Lateral Subsurface Discharge

Lateral Surface Discharge - Runoff

Usage

LateralSubsurfaceDischarge(psi, elevation = "DemFile",
  slope = "SlopeMapFile", aspect = NULL, neighbors = 8,
  soilmap = "SoilMapFile", soiltypes = "SoilParFile", Dz = "SoilDz",
  unitelev = 1000, use.soilwater.function = list(fun = "khy", args =
  list(ksat = "LateralHydrConductivity", alpha = "Alpha", n = "N", type_swc =
  "VanGenuchten", type_khy = "Mualem")),
  param.use.soilwater.function = c("LateralHydrConductivity", "Alpha", "N"),
  use_XY = FALSE, unit_input = c("millimeter", "centimeter", "meter"),
  output.discharge = FALSE, ...)

LateralSurfaceDischarge(hsup, elevation = "DemFile", slope = "SlopeMapFile",
  aspect = NULL, neighbors = 8, landmap = "LandCoverMapFile",
  flowres = "SurFlowResLand", expflowres = "SurFlowResExp",
  unit_input = c("millimeter", "centimeter", "meter"),
  output.discharge = FALSE, use_XY = FALSE, ...)

Arguments

psi	soil presssure head. It is RasterBrick-class object.
elevation, slope, aspect	elevation maps of elevation,slope,aspect. Slope and aspect must be expressed in degrees.
neighbors	Integer. Indicating how many neighboring cells to use to compute slope for any cell. Either 8 (queen case) or 4 (rook case). Only used for slope and aspect, see terrain.
soilmap	raster map of soil class types or related keyword of GEOtop file geotop.inpts (inpts.file argument in get.geotop.inpts.keyword.value.
soiltypes	list of soil class types or related keyword of GEOtop file geotop.inpts (inpts.file argument in get.geotop.inpts.keyword.value. See metacat arguments in SoilPropertyMap.
Dz	layer thickness or or related GEOtop keyword used in the soiltypes file.
unitelev	ratio between unis of elevation and psi, e.g. if elevetion is exxxpressed in meters and psi in millimeters , it is 1000 (Default).
use.soilwater.function	list object passed to SoilPropertyMap. It is utilized to calculate unsaturated hydraulic conductivity used in Darcy-Buckingham Law (see also khy).
param.use.soilwater.function	GEOtop keywords of the soil parameters utilized to calculate unsaturated hydraulic conductivity (see also khy)
use_XY	logical vaule. It it TRUE, the function returns the flux along the X (eastward) and Y (northward) directions.
...	further arguments for SoilPropertyMap
hsub	surfdace water depth. It is RasterBrick-class object.

See Also

terrain,,SoilPropertyMap

Examples

wpath <- system.file("geotop_simulation",package="geotopsim")

#' ## Use soilwater function 
library(soilwater)
data(PsiTheta)

t <- 5
psit <- psi[[t]]
hsupt <- hsub[[t]]

q <- LateralSubsurfaceDischarge(psi=psit,wpath=wpath)
qxy <- LateralSubsurfaceDischarge(psi=psit,wpath=wpath,,output.discharge=TRUE)





qdischarge_sub <- LateralSubsurfaceDischarge(psi=psit,wpath=wpath,output.discharge=TRUE)

qdischarge_sup <- LateralSurfaceDischarge(hsup=hsupt,wpath=wpath,output.discharge=TRUE)

ecor/geotopsim documentation built on May 15, 2019, 10:05 p.m.