Skip to contents

Function which creates screen plots giving an overview of the model outputs in the GRiwrm network

Usage

# S3 method for class 'GRiwrmOutputsModel'
plot(x, Qobs = NULL, Vobs = NULL, unit = "m3/s", ...)

Arguments

x

[object of class GRiwrmOutputsModel] see RunModel.GRiwrmInputsModel for details

Qobs

(optional) matrix time series of observed flows (for the same time steps than simulated) (mm/time step) with one column by hydrological model output named with the node ID (See CreateGRiwrm for details)

Vobs

(optional) matrix time series of observed or targeted storage (for the same time steps than simulated) [m3] with one column by hydrological model output named with the node ID (See CreateGRiwrm for details)

unit

(optional) character flows unit ("m3/s" or "mm")

...

Further arguments for airGR::plot.OutputsModel and plot

Value

list of plots.

Details

For examples of use see topics RunModel.GRiwrmInputsModel, RunModel_Reservoir, and RunModel.Supervisor.

The simulated discharge of a node integrates the flows of the whole catchment of the node, while the meteorological time series stored in the OutputsModel of the node (precipitation and potential evapotranspiration) come from the meteorological data of its local sub-basin only. Therefore, when these meteorological time series are available in the GRiwrmOutputsModel object (i.e. requested in Outputs_Sim, see CreateRunOptions), they are replaced on the fly, without modifying the object, by their area-weighted mean over the whole catchment of the node (the node and all its upstream nodes), so that the displayed meteorology reflects the whole catchment (See issue inrae/airGRiwrm#178).

For a CemaNeige node, the same adjustment is applied to the per-layer series of CemaNeigeLayers displayed by airGR::plot.OutputsModel: the mean air temperature (Temp), the solid precipitation (Psol) and the snow pack state (SnowPack) of each elevation layer are replaced by their area-weighted mean over the whole catchment of the node, as well as the mean over the layers that airGR derives from them. The other displayed series are not adjusted because they are either already at the scale of the whole catchment (the simulated discharge Qsim) or simulated at this scale (the actual evapotranspiration AE).

When the ID of a node is shared with a Diversion node (two rows of the network with the same id, see CreateGRiwrm), the local sub-basin area of the shared ID cannot be computed and the node is simply excluded from the weighted mean (graceful behaviour).

See also

adjustMeteoToCatchment() to compute statistics from the meteorological data at the whole-catchment scale.