qm.empqm.bulk.fa: The actual QM function for the bulk series considering...

Description Usage Arguments Value Author(s) Examples

View source: R/qm.empqm.bulk.fa.R

Description

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.

Usage

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)

Arguments

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.

Value

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.

Author(s)

Sven Kotlarski (MeteoSwiss), Jan Rajczak (ETH Zurich)

Examples

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)

SvenKotlarski/qmCH2018 documentation built on July 14, 2019, 7:39 p.m.