Description Usage Arguments Author(s) References Examples
The default parameter values are setup to perform the typical disaggregation
for CRSS, based on scaling the Upper Basin inflows, and performing a
direct resampling of Lower Basin tributaries, i.e., no scale factor is
applied to the LB tributaries. Sites 1-20 are scaled (sf_sites = 1:20
);
therefore the remaining sites (21:29) are not scaled.
1 2 | paleo_disagg(x, ann_index_flow, mon_flw, nsite = 29, sf_sites = 1:20,
nsim = 1, ofolder = NULL, index_years = NULL, k_weights = NULL)
|
x |
The annual paleo data to disaagregate. |
ann_index_flow |
Observed annual flow data used for picking index year. |
mon_flw |
Intervening monthly natural flow. Used for spatially and temporaly disaggregating paleo data based on the index year. |
nsite |
The number of sites to diaggregate the data too. |
sf_sites |
The site numbers (indeces), that will scale the index year's volume based on the annual flow being disaggregated. The remaining sites will select the index year directly. See Details. |
nsim |
Number of times to repeat the space/time disaggregation. |
ofolder |
Optional. If specified, the disaggregated flow data and the
selected index years are saved to this folder as csv files. There will be
one csv file for each time the disaggregation is repeated ( |
index_years |
Optional. If specified, these index years will be used instead of selecting years based on weights and sampling. |
k_weights |
If |
Ken Nowak
Nowak, K., Prairie, J., Rajagopalan, B., Lall, U. (2010). A nonparametric stochastic approach for multisite disaggregation of annual to daily streamflow. Water Resources Research.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 | ## Not run:
# read in annual synthetic mon_flw for disag
x <- matrix(scan("data-raw/Meko.txt"), ncol = 2, byrow = TRUE)
# intervening natural flow mon_flw - monthly CY text file
mon_flw <- as.matrix(read.table(
"tests/dp/NFfullbasinWY0608intervening.txt",
sep = "\t"
))
# observed annual flow for picking analog disag yr
ann_index_flow <- as.matrix(read.table("tests/dp/LFWYTotal.txt"))
zz <- paleo_disagg(x, ann_index_flow, mon_flw, 29, 1)
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.