Raven

rchlumsk wrote:Hi Elizabeth

I believe you just need to pull the data from within the hydrograph object returned by hyd.read to pass it to the function. The error you get is recreated when a NULL object is passed to the function. Here is an example using the function from the sample data in RavenR, hopefully you can use that to help you with your own data.

Code: Select all

library(RavenR)
data(hydrograph.data) # sample dataset
myhyd <- hydrograph.data # same format as hyd.read function

# manually specify the sim and obs data, or use hyd.extract function
sim <- myhyd$hyd$Sub43
obs <- myhyd$hyd$Sub43_obs
# ?hyd.extract
# dd <- hyd.extract(subs="Sub43",myhyd)

annual.peak(sim,obs)


Cheers,
Rob

Konhee Lee wrote:Hello Elizabeth,

The GAWSER-style snow balance was an attempt to emulate the SNOMLT process described in
acwc.sdp.sirsi.net/client/en_US/search/asset/1001420;jsessionid=22BBBEC081240AB73639002B9F4AE46B.enterprise-15000


Unfortunately, we did not have access to the source code and data, so we were not able to emulate it exactly the same.
We had matching snow water equivalent (SWC) and liquid water content (LWC), but our depth was off by 1 cm compared to the example in the documentation (pg 26).

An example rvi file to use the process can be something like:

#########################################################################
:FileType rvi ASCII Raven 2.5
:WrittenBy Konhee Lee
:CreationDate Nov 2016
# Simple snowmelt model to test GAWSER
#------------------------------------------------------------------------
:RunName run1
:StartDate 2000-01-01 00:00:00
:Duration 21
:Timestep 1
:Method ORDERED_SERIES
:Routing ROUTE_NONE
:SoilModel SOIL_ONE_LAYER
:RainSnowFraction RAINSNOW_DATA

#------------------------------------------------------------------------
#
:HydrologicProcesses
:Precipitation PRECIP_RAVEN ATMOS_PRECIP MULTIPLE
:SnowBalance SNOBAL_GAWSER MULTIPLE MULTIPLE
:Flush RAVEN_DEFAULT PONDED_WATER SURFACE_WATER
:EndHydrologicProcesses
#------------------------------------------------------------------------

#---------------------------------------------------------
# Output Options
#
:WriteForcingFunctions
:DebugMode yes

Let me know if you have any questions.

Regards,
Konhee

ephilip wrote:

rchlumsk wrote:Hi Elizabeth

I believe you just need to pull the data from within the hydrograph object returned by hyd.read to pass it to the function. The error you get is recreated when a NULL object is passed to the function. Here is an example using the function from the sample data in RavenR, hopefully you can use that to help you with your own data.

Code: Select all

library(RavenR)
data(hydrograph.data) # sample dataset
myhyd <- hydrograph.data # same format as hyd.read function

# manually specify the sim and obs data, or use hyd.extract function
sim <- myhyd$hyd$Sub43
obs <- myhyd$hyd$Sub43_obs
# ?hyd.extract
# dd <- hyd.extract(subs="Sub43",myhyd)

annual.peak(sim,obs)


Cheers,
Rob

Good Morning Rob,

I tried the code you provided. It works with the data from RavenR for annual.peak(sim,obs). The function annual.peak(sim,obs) do work with my data as well. My post is not for annual.peak(sim,obs) but for the following functions:

annual.peak.event(sim,obs,add.r2=F)
Error in as.POSIXlt.numeric(zoo::index(x)) : 'origin' must be supplied

annual.peak.event.error(sim,obs)
Error in as.POSIXlt.numeric(zoo::index(x)) : 'origin' must be supplied

annual.peak.timing.error(sim,obs)
Error in as.POSIXlt.numeric(zoo::index(x)) : 'origin' must be supplied

flowdurcurve.plot(myhyd$sim, prd='1999-10-01/2007-09-30',seasonal=F)
Error in as.POSIXlt.default(x, tz = tz(x)) :
do not know how to convert 'x' to class “POSIXlt”

I get the same error with the data from the RavenR as well as from my data.

Your help is highly appreciated.

Thank you
Elizabeth

rchlumsk wrote:Hi Elizabeth

Related to this thread quoted here, the package has been updated and the mentioned functions should not be working. Let me know if you come across any more bugs.

Cheers,
Rob

ephilip wrote:

rchlumsk wrote:Hi Elizabeth

I believe you just need to pull the data from within the hydrograph object returned by hyd.read to pass it to the function. The error you get is recreated when a NULL object is passed to the function. Here is an example using the function from the sample data in RavenR, hopefully you can use that to help you with your own data.

Code: Select all

library(RavenR)
data(hydrograph.data) # sample dataset
myhyd <- hydrograph.data # same format as hyd.read function

# manually specify the sim and obs data, or use hyd.extract function
sim <- myhyd$hyd$Sub43
obs <- myhyd$hyd$Sub43_obs
# ?hyd.extract
# dd <- hyd.extract(subs="Sub43",myhyd)

annual.peak(sim,obs)


Cheers,
Rob

Good Morning Rob,

I tried the code you provided. It works with the data from RavenR for annual.peak(sim,obs). The function annual.peak(sim,obs) do work with my data as well. My post is not for annual.peak(sim,obs) but for the following functions:

annual.peak.event(sim,obs,add.r2=F)
Error in as.POSIXlt.numeric(zoo::index(x)) : 'origin' must be supplied

annual.peak.event.error(sim,obs)
Error in as.POSIXlt.numeric(zoo::index(x)) : 'origin' must be supplied

annual.peak.timing.error(sim,obs)
Error in as.POSIXlt.numeric(zoo::index(x)) : 'origin' must be supplied

flowdurcurve.plot(myhyd$sim, prd='1999-10-01/2007-09-30',seasonal=F)
Error in as.POSIXlt.default(x, tz = tz(x)) :
do not know how to convert 'x' to class “POSIXlt”

I get the same error with the data from the RavenR as well as from my data.

Your help is highly appreciated.

Thank you
Elizabeth