deposition | R Documentation |
Run the ‘bLSmodelR’ dry deposition post-processing function.
deposition(x, vDep, rn = NULL, Sensor = NULL, Source = NULL, vDepSpatial = NULL,
ncores = 1, memory_limit = NULL, show_progress = TRUE, variables = 'CE')
x |
A |
vDep |
either a numeric vector providing the model surface deposition velocity (at height d + z0) for the rows specified by argument |
rn |
integer. Subsetting by column ‘rn’. Only matching rows will be post-processed. |
Sensor |
character. Subsetting by column ‘Sensor’. Only matching rows will be post-processed. |
Source |
character. Subsetting by column ‘Source’. Only matching rows will be post-processed. |
vDepSpatial |
a |
ncores |
integer. Number of cores used for parallel processing. |
memory_limit |
character. Maximal vector heap size (memory limit) specified via 'R_MAX_VSIZE' (see also |
show_progress |
logical. Show progress when computing output variables in prallel. Defaults to |
variables |
character vector specifying the output variables of line sensors to calculate (any combination of |
A data.frame
(or data.table
) of class bLSresult
and deposition
containing the results of the model run and the deposition post-processing.
The following attribute
is additionally attached to the existing bLSresult
attributes
:
vector of surface deposition velocities
Christoph Haeni
bLSmodelR-package
, genInputList
, coreModel
.
## Not run:
# Example Run with default values:
Sensor <- genSensors(
PointSensor1=list(x=0,y=0,z=2)
,PointSensor2=list(x=10,y=0,z=2)
)
Sources <- genSources("Circle 1" = list("c", M = c(0,50), R = 10))
Ints <- genInterval(MaxFetch=70)
InputList <- genInputList(Sensor,Ints,Sources)
## run model
Run <- runbLS(InputList,Cat.Path=getwd())
## dry deposition post-processing:
RunDep <- deposition(Run,vDep=0.035,Sensor='PointSensor1')
## spatially inhomogeneous deposition velocity (still experimental!)
DepAreas <- genSources(
"Circle 2" = list("c", M = c(0,30), R = 10)
,"Circle 3" = list("c", M = c(0,10), R = 10)
)
DepVel <- list("Circle 2" = 0.005,"Circle 3" = 0.01)
plot(Sensor,join(Sources,DepAreas),SourceTextArgs=list(labels=paste0('vDep = ',c(0,0.005,0.01))))
RunDep2 <- deposition(Run,vDep=0.035,Sensor='PointSensor1',vDepSpatial=list(DepVel,DepAreas))
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.