dvoverlap: Find Sites with Overlapping Complete Years of Daily Values

In wasquith-usgs/akqdecay: Many Distributional Statistics of USGS Daily Streamflow Hydrographs and Powerful Utilities

Description

Find sites with overlapping complete years of daily values from the daily-value environment.

Usage

dvoverlap(dvenv, begyr=1987, endyr=2016, type=c("wyear", "year"),
                 missing.days=7, silent=TRUE, ...)

Arguments

dvenv	An R environment previously populated by fill_dvenv;
begyr	Starting year;
endyr	Ending year;
type	Are calendar or water years desired;
missing.days	The number of permissible missing days in a given year on which to still report it as a complete year;
silent	Suppress informative calls to message(); and
...	Additional arguments to pass, if ever used.

Value

An R vector of the site numbers from the dvenv meeting the criteria as set forth by the arguments to the function.

Author(s)

W.H. Asquith

See Also

fill_dvenv

Examples

# Six USGS streamgages of the Wolf River in Tennessee
wolfriver <- c("07030392", "07030500", "07031650",
               "07031660", "07031700", "07031740")
wolf.env <- new.env() # the standard declaration of an environment
fill_dvenv(wolfriver, envir=wolf.env) # wolf.env now filled with six tables.
sites_overlap <- dvoverlap(wolf.env, begyr=1987, endyr=2016, silent=FALSE)
# No sites overlap, but change the years and we see
sites_overlap <- dvoverlap(wolf.env, begyr=2000, endyr=2010, silent=FALSE)
# dvoverlap() for 07030392 has 11 'full' years: (2000,2010)
# dvoverlap() for 07030500 --- does not meet criteria
# dvoverlap() for 07031650 has 11 'full' years: (2000,2010)
# dvoverlap() for 07031660 --- does not meet criteria
# dvoverlap() for 07031700 --- does not meet criteria
# dvoverlap() for 07031740 has 11 'full' years: (2000,2010)

# Now let us double check by visualization
for(site in ls(wolf.env)) {
    tmp <- get(site, wolf.env)
    plot(tmp$Date, tmp$Flow, type="l", log="y")
    mtext(site)
} # 

wasquith-usgs/akqdecay documentation built on Nov. 9, 2020, 1:13 p.m.