LateralSubsurfaceDischarge: Lateral Subsurface Discharge

Description Usage Arguments See Also Examples

Description

Lateral Subsurface Discharge

Lateral Surface Discharge - Runoff

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
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

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
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.