Conduct NWRFC style forcing adjustments

Usage

fa_nwrfc(
  dt_hours,
  forcing,
  pars,
  climo = NULL,
  dry_run = FALSE,
  return_adj = FALSE,
  return_climo = FALSE,
  adjust = TRUE
)

dt_hours: timestep in hours
forcing: data frame with with columns for forcing inputs
pars: sac parameters
climo: climotology matrix
dry_run: Do a run without any forcing adjustments, only compute pet and etd. Shorthand for adjust = FALSE; when TRUE it overrides adjust.
return_adj: return monthly adjustment factors only
return_climo: return the computed monthly climo
adjust: Which forcings to adjust: TRUE for all of them, FALSE for none, or a character vector naming a subset of "map", "mat", "ptps" and "pet". Unadjusted forcings still pass through the routine, which computes pet and etd, but with neutral adjustment parameters.

Matrix (1 column per zone) of unrouted channel inflow

data(nrkw1_forcing)
data(nrkw1_pars)
dt_hours <- 6
forcing_adj <- fa_nwrfc(dt_hours, nrkw1_forcing, nrkw1_pars)