pr_a: Find the a parameter with bias correction for precipitation

In veravavan/AdvDeltaChange: Perform Advanced Delta Change on Data

Description

Find the a parameter with bias correction for precipitation

Usage

pr_a(scen_q, ctrl_q, obs_q, b)

Arguments

scen_q	60 and 90 quantiles of the scenario period precipitation data
ctrl_q	60 and 90 quantiles of the control period precipitation data
obs_q	60 and 90 quantiles of the observed precipitation data in the control period timeline
b	Already calculated b parameter with pr_b

Value

The a parameter value

Examples

load('data/pr_sim.rda')
pr_sim[, pr5 := 60*60*24*frollsum(pr, 5, align = 'c')]
pr_sim = pr_sim[!is.na(pr5)]
pr_ctrl = pr_sim[year(DTM) %in% 1981:2010]
pr_scen = pr_sim[year(DTM) %in% 2071:2100]
ctrl_q = quantile(pr_ctrl$pr5, c(.6,.9))
scen_q = quantile(pr_scen$pr5, c(.6,.9))
load('data/pr.rda')
pr[, pr_o5 := frollsum(pr_o, 5)]
pr = pr[year(DTM) %in% 1981:2010]
obs_q = quantile(pr$pr_o5, c(.6,.9))
b <- pr_b(scen_q, ctrl_q, obs_q)
a <- pr_a(scen_q, ctrl_q, obs_q, b)

veravavan/AdvDeltaChange documentation built on Jan. 18, 2022, 7:05 p.m.