Description Usage Arguments Value Author(s) Examples
View source: R/qm.empqm.bulk.fa.R
qm.empqm.bulk.fa
is the actual QM function for
the bulk series. Equivalent to function qm.empqm.doy.fa but
for the bulk series, i.e. disregarding a day-of-year (DOY) dependence
of the correction function.
1 2 3 | qm.empqm.bulk.fa(series, cdf.mod.bulk, cdf.obs.bulk, cdf.obs.wet.bulk,
cor.fun.bulk, fwet.obs, fwet.mod, pr.wet, method = c("linear", "bin"),
var, output.diag = FALSE, minq, maxq, incq)
|
series |
Time series of modelled data, e.g. transient (control+) scenario series. |
cdf.mod.bulk |
CDF of the modelled time series in the calibration period, estimated with function qm.cdf.bulk (vector[nquantiles]). |
cdf.obs.bulk |
CDF of the observed time series in the calibration period, estimated with function qm.cdf.bulk (vector[nquantiles]). |
cdf.obs.wet.bulk |
CDF of the observed time series in the calibration period for wet days only, estimated with function qm.cdf.bulk (vector[nquantiles]). |
fwet.obs |
Observed wet day frequency (output of function qm.cdf.doy). |
fwet.mod |
Modelled wet day frequency output of function qm.cdf.doy). |
pr.wet |
Internal precipitation wet day threshold [mm/day]. |
method |
QM method: "binary" -> correction for closest quantile is used, "linear" -> linear interpolation of correction function between lower and upper quantile. |
var |
Variable to correct. |
output.diag |
Diagnostic QM console output (TRUE) or not (FALSE); default: FALSE. |
minq |
Minimum quantile for correction function [0.01 .. 0.99]; note: should not be 0 (for correct handling of extremes)! |
maxq |
Maximum quantile for correction function [0.01 .. 0.99]; note: should not be 1 (for correct handling of extremes)! |
incq |
Quantile increment for correction function (bin size). |
corr.fun.bulk |
Additive correction function (cdf A - cdf B), obtained by function qm.corfun.bulk. |
List of 3: $qm.input.series: Input series to QM (input argument series). $qm.corrected.series: Corrected series. $quantile.index: A series of quantile indices wrt. simulated series in calibration period.
Sven Kotlarski (MeteoSwiss), Jan Rajczak (ETH Zurich)
1 2 3 4 5 6 7 8 | ## Not run:
# Carry out bulk QM with frequency adaptation. Standard quantiles.
# Variable: precipitation. Observed wet day frequency of 0.9, modelled
# wet-day frequency of 0.95. Wet-day treshold of 0.1 mm/day. Linear QM,
# diagnostic terminal output.
qm.empqm.bulk.fa(series, cdf.modelled.bulk, cdf.observed.bulk, cdf.observed.wet.bulk, correction.function.bulk, 0.9, 0.95, 0.1, 'linear', 'pr', TRUE, 0.01, 0.00, 0.01)
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.