The hardware and bandwidth for this mirror is donated by METANET, the Webhosting and Full Service-Cloud Provider.
If you wish to report a bug, or if you are interested in having us mirror your free-software or open-source project, please feel free to contact us at mirror[@]metanet.ch.
CreateCalibOptions() now uses parameter screening for
RunModel_Lag() which are now expressed in the transformed
space instead of the parameter space. (#156)RunModel_CemaNeige*() now takes into account the case
when dG = 0. (#178)Calibration_Michel() now runs faster as the
ProposeCandidatesGrid() was improved to create the propose
candidates grid. (#157)TransfoParamGR5J() now returns the correct error
message when the number of parameters is incorrect. (#168)frun_* Fortran subroutine does not use anymore the
‘DLLEXPORT’ command. (#180).ErrorCrit() and .FeatModels function are
no more exported. (#173)CreateInputsCrit_Lavenne() to define a composite
criterion based on the formula. Added the
CreateErrorCrit_GAPX() function to compute an error
criterion based on the GAPX formula. (#111)OutputsModel, returned by the
RunModel_*GR*() function, gains a RunOptions
element which is a list and contains 2 sub-elements:
WarmUpQsim (vector series of simulated discharge on the
warm-up period) and Param (vector of the model parameter
values). (#123)plot.OutputsModel() gains a AxisTS
argument in order to manage x-axis representing calendar dates and
times. It avoids to display ugly x-axis. (#122)LatRad argument has been removed from
the PEdaily_Oudin() function. (#81)Qobs argument has been removed from the
CreateInputsCrit() function. (#81)Ind_zeroes argument has been removed
from the CreateInputsCrit() function. (#81)verbose argument has been removed from
the CreateInputsCrit() function. (#81)FUN_CRIT argument has been removed from
the ErrorCrit() function. (#81)SeriesAggreg() now correctly reorders regime time
series when the monthly regime is computed from a time series that does
not start in January. It also keeps original data.frame
column names. (#133)DataAltiExtrapolation_Valery() now correctly extract
HypsoData values for each elevation layers. The selected indices were
wrong (one less than expected) (#144)CreateIniStates() does not return anymore an error
message when IntStore is set and RunModel_GR5H
is used. (#144)RunModel_Lag() now handles warm-up period simulation
(set in CreateRunOptions()). (#132)PE_Oudin() can use inconsistent time series. It allows
to mixing time series from different stations. (#134).GetFeatModel() in
CreateCalibOptions() and CreateIniStates()
functions in order to simplified their codes. (#111).FunTransfo in order to manage the parameter
transformations and to simplified the code of the
CreateCalibOptions() function (#111).ArgumentsCheckGR() function in order to
check the arguments of theRunModel_*() functions and
simplified their codes. (#129).GetOutputsModelGR() function in order to
manage the outputs of theRunModel_*() functions and
simplified their codes. (#129)plot.OutputsModel() function has been
slightly simplified. (#122,
#147)CreateInputsModel() gains a QupstrUnit
argument in order to manage the unit of the flow in the
Qupstream argument in case of the use of a semi-distributed
version of a hydrological model. (#110)RunModel_Lag() gains a QcontribDown
argument containing the time series of the runoff contribution of the
downstream sub-basin in case of the use of a semi-distributed version of
a hydrological model. (#109)RunModel. The RunModel_Lag()
can now be passed to the FUN_MOD argument. (#108)RunModel_Lag(). The function no longer
returns two values for a single time step run. (#102)RunModel_Lag(). The StateEnd
value is now correct when there are more than a single upstream basin.
(#103)RunModel_Lag(). The StateEnd
value is now correct when the upstream flow unit is mm/time step. (#104)RunModel_CemaNeigeGR5H(). The solid
precipitation are now taken into account in the GR5H model. (#105)RunModel_CemaNeige() and
CreateInputsModel(). RunModel_CemaNeige() now
runs at the hourly time step. (#106)LengthHydro must now be set in kilometers (not anymore
in meters) in the CreateInputsModel() function. (#112)TransfoParam_GR5H() now use the same transformation as
TransfoParam_GR4H() for the X1 parameter. The previous
transformation set by Ficchì seems unnecessary as it provokes
irrealistically high X1 values. (#50)RunModel*() functions now run faster. The
computation times are significantly shorter for long times series with
many time steps (e.g. hourly times series), due to a better management
of the missing values in and out the Fortran codes. Only simulation
computation times have been improved (it is largely invisible to the
user for calibration computation times). (#113)RunModel_*()
functions used during the parameter estimation process now run faster
because the outputs contain only the simulated flows (see the
Outputs_Sim argument in the CreateRunOptions()
help page)..FeatModels() and .GetFeatModel()
functions in order to repectively store and get model features
(e.g. name, number of parameters, time unit). Therefore the codes of the
CreateInputsModel() and the CreateRunOptions()
functions have been simplified. (#106)Imax() now returns an error message when
IndPeriod_Run doesn’t select 24 hours by day, instead of
numeric(0). (#92)frun_GR5H.f90. (#93)character dates into POSIXlt in
RunModel_GR1A() example and in SeriesAggreg()
tests in order to avoid bad subsetting on time series due to mixing UTC
and local time (error returned on macOS flavors). (#94)SeriesAggreg S3 method with functions for
InputsModel, OutputsModel, list,
data.frame class objects. This new version of the
SeriesAggreg() function also allows to compute regimes. (#25).GetAggregConvertFun() private function in order
to choose automatically the ConvertFun to apply on each
element of objects used in SeriesAggreg.InputsModel() and
SeriesAggreg.OutputsModel()..AggregConvertFunTable data.frame that allows the
user to see what names of list items or data.frame column names are
guessed and eventually customise this correspondence table.PE_Oudin() now presents a RunFortran
argument to run the code in Fortran or in R. The Fortran mode is the
fastest. (#62)RunModel_Lag() which allows to perform a single
run for the Lag model over the test period in order to run
semi-distributed GR models. (#34)[ S3 method for InputsModel class
object in order to extract subsets of it. (#67)TimeFormat argument is now deprecated in
SeriesAggreg(). (#41)NewTimeFormat argument is now deprecated in
SeriesAggreg() and replaced by the Format
argument. (#41)RunSnowModule argument has been removed
from the CreateRunOptions() function. (#23)SeriesAggreg(). The function now runs when
TimeLag >= 3600. (#41)SeriesAggreg(). The function now runs when
the time series contain some columns entirely filled with missing
values. (#43)RunModel_GR1A(). Reversed PotEvap and
Precip outputs are now reordered (in the previous versions PotEvap
contained the precipitation values and Precip contained the
evapotranspiration values, the Qsim values were already correct). (#65)RunModel_GR2M() function (Ps). (#51)PE_Oudin() can now run for several locations
(i.e. several latitudes) in the Fortran mode
(RunFortran = TRUE). In this case Lat must be
of the same length as Temp. (#62)RunModel() now allows to run semi-distributed GR
models. (#34)ConvertFun argument of the
SeriesAggreg() function can now be set to names of
aggregation functions that return value of length 1 (not only
"sum" or "mean", but e.g. "min",
"max", "Q95"). (#82).FortranOutputs() function is no longer exported in
the namespace.RunModel_GR1A() now uses the Fortran version of the
model code. This code is no longer duplicated: the R version which was
used was removed. (#65)PE_Oudin() and SeriesAggreg() functions. (#37)RunModel_*() funcions were cleaned up, with no effect
on their outputs. (#14).ErrorCrit() function now returns a warning message
when a criterion computed on less than 10 time-steps (whatever the unit
of the time step). (#14)RunModel_GR5H()
documentation. (#49)Exch was renames AExch in the
RunModel_GR2M() output. (#87)RunModel_GR2M() help page. (#88)plot.OutputsModel() function does not check anymore
the time step by comparing the calculation of the difference of the last
two time steps because it is already checked by the class of the
OutputsModel object, which is therefore assumed to be
necessarily valid. (#56)anyNA function.Imax()
function in order to run faster.citation("airGR") command.Imax(). The default value of the
TestedValues argument was wrong due to a mistyped argument
name in the seq() function.plot.Outputsmodel() now allows to draw actual
evapotranspiration when which = "ActEvap" or
which = "All" (overlaid to potential evapotranspiration if
already drawn). (#2)RunModel_GR5H() and
RunModel_CemaNeigeGR5H() functions to run the hourly model
GR5H (with or without the CemaNeige module). These models present an
optional additionnal interception store. (#13)Imax() which allows to estimate the maximum
capacity of the GR5H interception store. (#13)TransfoParam_GR1A(). The number of model
parameters was wrong (2 instead of 1) which caused an error during the
GR1A model calibration. (#1)plot.OutputsModel(). The function does not
return any error message when log_scale = TRUE,
Qobs = NULL and user want to draw flows time series.RunModel_*GR*(). The functions do not
return any error message anymore due to slightly negative values
returned by GR4H, GR4J, GR5J or GR6J Fortran codes (the message was
returned by CreateIniStates() when the final states were
created). The RunModel_*GR*() functions now return zero
instead of these slightly negative values, except for the ExpStore where
negatives values are allowed..ErrorCrit() function. The Box-Cox
transformation formula is now corrected when the
ErrorCrit*() functions are used.RunModel_GR4H() function (Pn, Ps,
AExch1, AExch2).RunModel_GR2M()
documentation.https://gitlab.irstea.fr/HYCAR-Hydro/airgr) and issues
(https://gitlab.irstea.fr/HYCAR-Hydro/airgr/issues).RunModel_CemaNeige(). The function now
runs correctly when IndPeriod_WarmUp = 0L in
CreateRunOptions() in order to completely disable the
warm-up period (e.g. to perform a forecast form a given initial
state).CreateIniStates(). The function now
returns the right number of end states when CemaNeige is used without
hysteresis.RunModel_CemaNeige*() functions. G and
Gthr end states are no more inverted in the output values.NA into the BasinObs
time series of the L0123001 dataset.CreateInputsCrit() now allows power (as a numeric or as
a character) and the Box-Cox transformations in the transfo
argument.RunModel_CemaNeigeGR4H() function to run the
hourly model GR4H with the CemaNeige module.PE_Oudin() function to compute Oudin’s potential
evapotranspiration for hourly or daily time steps.plot.OutputsModel() now presents a
LayoutMat argument (and additionnal related argument:
LayoutWidths, LayoutHeights) to specify
complex plot arrangements.PEdaily_Oudin() function is deprecated and his use
has been replaced by the use of PE_Oudin().plot.OutputsModel(). The function now runs
correctly when the which argument contains the
"CorQQ" value without "CumFreq".plot.OutputsModel() can now draw PE or error time
series if the which argument is set to "all"
or "PotEvap" or "Error".plot.OutputsModel() now allows new values for the which
argument: "all" corresponds to all graphs,
"synth" corresponds to the main graphs (default value;
corresponding to "all" in the previous versions of the
package)
(i.e. c("Precip", "Temp", "SnowPack", "Flows", "Regime", "CumFreq", "CorQQ")),
"ts" corresponds to the time series graphs
(i.e. c("Precip", "PotEvap", "Temp", "SnowPack", "Flows"))
and “perf” corresponds to the performance graphs
(i.e. c("Error", "Regime", "CumFreq", "CorQQ"))..ErrorCrit() private function added to check inputs
into ErrorCrit_*() functions. The
ErrorCrit_*() functions were simplified accordingly.CreateInputsCrit() now returns FUN_CRIT as
a character string.plot.OutputsModel() function.CreateInputsCrit() now presents a VarObs
argument in order to allow to prepare an InputsCrit object
in order to run a criterion on other variables than observed discharges
with the ErrorCrit() function (at the moment SCA and
SWE).CreateInputsCrit() can now prepare an
InputsCrit object in order to compute a single criterion
(Single class), multiple criteria (Multi
class) or a composite criterion (Compo class) with the
ErrorCrit() function.CreateInputsCrit() now presents a Weights
argument in order to allow to prepare an InputsCrit object
in order to compute a composite criterion (Compo class)
with ErrorCrit() or Calibration_Michel().CreateInputsCrit() now returns a idLayer
element to indicate which layer(s) to use for SCA or SWE
aggregation.CreateInputsCrit() now presents a warnings
argument to replace the verbose action (the verbose
argument is kept to print messages).CreateInputsCrit(), it is now possible to set the
following arguments as atomic (as before) or as list:
FUN_CRIT, VarObs, Obs,
BoolCrit, transfo, Weights. If
the list format is chosen, all the lists must have the same length.CreateRunOptions(), CreateIniStates() and
CreateCalibOptions() now present a IsHyst
argument to give the possibility to use the Linear Hysteresis with
CemaNeige. The objects returned present an hysteresis
class.CreateRunOptions() now presents a warnings
argument to replace the verbose action (the verbose
argument is kept to print messages).TransfoParam_CemaNeigeHyst() function in order to
take into account transformation of the parameters of the CemaNeige
module when the Linear Hysteresis is used.X0310010 dataset to run the examples using
the Linear Hysteresis with CemaNeige (it contains necessary SCA
data).Qobs argument is now deprecated in
CreateInputsCrit() and has been renamed
Obs.FUN_CRIT argument is now deprecated in
ErrorCrit(). This function now gets this information from
the InputsCrit argument.FUN_CRIT argument is now deprecated in
Calibration_Michel(). This function now gets this
information from the InputsCrit argument.plot_OutputsModel() had been deprecated in ‘airGR’
1.0.4 (it had been replaced by the use of
plot.OutputsModel() or plot()) and is defunct
now.CreateInputsCrit() now return a list of
InputsCrit (each element is of the Single
class) in the cases of multiple or a composite criteria.ErrorCrit_*() functions now return an error message if
the InputsCrit object is of class Multi or
Compo.ErrorCrit() function can now run on a multiple or a
composite InputsCrit. In these cases, it returns a list of
ErrorCrit.ErrorCrit() and ErrorCrit_*() functions
can now assess Q, SCA or SWE simulations.Calibration_Michel() function can now run on a
composite InputsCrit. It returns a composite value of error
and the formula used to calculate it.RunModel_GR4J(), RunModel_GR5J() and
RunModel_GR6J() functions.plot.OutputsModel() documentation with
?plot.FUN_* arguments (in addition to function objects) in the
following functions: Calibration(),
Calibration_Michel(), CreateCalibOptions(),
CreateIniStates(), CreateIniStates(),
CreateInputsCrit(), CreateInputsModel(),
CreateRunOptions(), ErrorCrit(),
RunModel() and TransfoParam().ErrorCrit_*() functions now return objects of class
ErrorCrit and NSE, KGE,
KGE2 or RMSE..FortranOutputs() private function added to manage
Fortran outputs.frun_GR2M Fortran subroutine were
reordered.DataAltiExtrapolation_Valery() now returns named
elements of lists relative to elevation layer.Calibration() function now returns an error message if
FUN_CALIB is not a function.PEdaily_Oudin() are now checked.PEdaily_Oudin() example corrected (the Julian day was
one day too early).plot.OutputsModel() does not return a warning message
anymore when Qobs = NULL.TransfoParam*() functions are now
checked.CreateRunOptions(). The function now
accounts correctly for leap years when no warm-up period is
defined.CreateRunOptions() was cleaned up, with no effect on
its outputs.vignetteParam*.rda datasets moved to the inst
directory. It contains different objects needed for ‘param_optim’ and
‘param_mcmc’ vignettes.PEdaily_Oudin() now presents a LatUnit
argument which allows to choose the unit of the latitude (radians and
degrees).LatRad argument is now deprecated in
PEdaily_Oudin() and replaced by the Lat
argument.Ind_zeroes argument of the
CreateInputsCrit() function is now deprecated.verbose argument is now deprecated in
CreateInputsCrit() and replaced by the
warnings argument.Calibration_Michel() is now faster during the
grid-screening step when a parameter is set using
FixedParam in CreateCalibOptions().CreateCalibOptions() now returns an error when all the
parameters are set in the FixedParam argument and a warning
message when all the parameters are free (NA) in the
FixedParam argument.CreateInputsCrit() now returns an error when
epsilon is not positive.CreateInputsCrit() now returns a warning message in the
following case: there are zeroes values in Qobs,
epsilon = NULL and transfo = log or
inv.ErrorCrit_*() functions now return a warning message in
the following case: there are zeroes values in Qobs or
Qsim, epsilon = NULL and
transfo = log or inv.L0123001 dataset.Vignette_Param. datasets in order to reduce
runtime during the re-building of vignettes. It contains different
objects needed for param_optim and param_mcmc vignettes.RunModel_GR2M(). The function now returns
the total precipitation (P) instead of the net rainfall (P1).RunModel_GR2M() now returns more explicit precipitation
outputs names.CreateInputsCrit() now returns a warning message when
the KGE (or KGE’) is used with a log transformation on flows.citation("airGR").CreateIniStates() function in order to help user
to format the IniStates argument for
CreateRunOptions().Param_Sets_GR4J dataset. It contains
generalist parameter sets for the GR4J model.RunSnowModule argument is now deprecated in
CreateRunOptions().RunModel_GR4H(): in frun_GR4H
Fortran subroutine, St(2) is now set to 0 (instead of
St(1)) when St(2) < 0.plot.OutputsModel() for the regime plot
when the period is less than 1 year.plot.OutputsModel() when there is no
common data to plot the cumulative frequency or the correlation QQ.plot.OutputsModel() for the y-axis
labelling of flows time series when log_scale = TRUE and
BasinArea is used.RunModel_GR4J(), RunModel_GR5J() and
RunModel_GR6J() (and CemaNeige_GR*J()) now
return Ps, Pn and actual exchanges. See the model Fortran codes for more
details about the calculation of these variables.CreateInputsModel() now returns an error when
DatesR contains duplicated values.RunModel_GR5J now returns StateEnd in the
same order as the other models.plot.OutputsModel() now returns a warning message when
the length of Qobs is different from the length of Qsim.RunModel_*() functions now return a warning message in this
case. RunModel_*() functions now return a warning when X4
< 0.5 and its value is set to 0.5.?L0123001, ?L0123002 and
?L0123003 now return the documentation page related to
BasinObs.DataAltiExtrapolation_Valery() and
CreateInputsModel() now present a PrecipScale
argument which allows rescaling precipitation when it is interpolated on
the elevation layers when CemaNeige is used.DataAltiExtrapolation_Valery(). The
elevation gradients for air temperature returned by
CreateInputsModel() are improved.DataAltiExtrapolation_Valery() has been improved.
DataAltiExtrapolation_Valery() now runs faster (and by
consequence CreateInputsModel() too, when CemaNeige is
used).RunModel_CemaNeige(),
RunModel_CemaNeigeGR4J(),
RunModel_CemaNeigeGR5J() and
RunModel_CemaNeigeGR6J() now return air temperature for
each elevation layer.OutputsModel objects. It
means that the plot_OutputsModel() function is deprecated
and his use has been replaced by the use of
plot.OutputsModel() or plot().plot.OutputsModel() the PlotChoice
argument is deprecated and has been renamed which.plot.OutputsModel() displays air temperature time
series for each layer when CemaNeige is used (argument
which = "Temp" or "all").ErrorCrit_*() functions gain a warnings
argument to replace the verbose action and the verbose
argument now prints the criterion value(s).CreateCalibOptions() when
StartParamList or StartParamDistrib arguments
are used.CreateInputsModel() now returns an error if
NLayers <= 0 when CemaNeige is used.plot_OutputsModel() now displays raw values on the
y-axis when the discharge time series is represented with log scale
(formerly, log values of discharges were displayed on the y-axis).SeriesAggreg() gains a TimeLag argument
that corresponds to a numeric value indicating a time lag (in seconds)
for the time series aggregation (useful to aggregate hourly time series
to the daily time step for instance). In addition, the function now
accepts input dates in both POSIXt formats
(POSIXct and POSIXlt). The output is in
POSIXct format.plot_OutputsModel() gains a log_scale
argument in order to plot the flow with a log scale.vignette("airGR") command.CreateCalibOptions() loses the OptimParam
argument that was redundant with the FixedParam argument.
The Calibration_Michel() was modified to take into account
this change by using directly FixedParam, but this is
transparent to the user.CreateCalibOptions() loses the StartParam
argument that was not used.sort for the transfo argument of
CreateInputsCrit() was not taken into account. It is now
fixed.RunModel_GR6J() and
RunModel_CemaNeigeGR6J() models were modified back to
versions previous to 1.0.1 to prevent from unwanted efficiency criteria
deterioration related to the calibration with
Calibration_Michel(). The actual model codes were not
modified but the TransfoParam_GR6J() and
CreateCalibOptions() functions were modified regarding the
X5 parameter. It is strongly advised to use airGR 1.0.2 for the
RunModel_GR6J() and RunModel_CemaNeigeGR6J()
functions if you are using Calibration_Michel(), as they
are much more efficient. In case you were using your own calibration
algorithm, you will not notice any difference.CreateInputsModel() and
DataAltiExtrapolation_Valery() functions now allow both
POSIXt formats (POSIXct and
POSIXlt).Calibration_HBAN() and
DataAltiExtrapolation_HBAN() functions have respectively
been renamed as Calibration_Michel() and
DataAltiExtrapolation_Valery() after the names of their
creators.Calibration_optim() function has been removed from
the package.verbose = TRUE
argument (formerly quiet = FALSE) in the following
functions: Calibration(),
Calibration_Michel(), CreateInputsModel(),
CreateRunOptions(),
DataAltiExtrapolation_Valery(), ErrorCrit(),
ErrorCrit_KGE(), ErrorCrit_KGE2(),
ErrorCrit_NSE(), ErrorCrit_RMSE(),
plot_OutputsModel(), SeriesAggreg().RunModel_GR5J() function) and CemaNeige
(RunModel_CemaNeigeGR5J() function) models will be
different.TransfoParam_*() functions). Since this
modifies the local search, calibration results will be different .CreateCalibOptions() function). Since
these quantiles constitute the starting point of the calibration
algorithm, calibration results will be different.citation("airGR").CreateInputsModel() that was related to
the handling of missing values.CreateRunOptions() that prevented the
correct use of the IniResLevels argument (to manually set
the filling rate of the production and routing stores).IndPeriod_WarmUp = 0.RunModel_GR4H() function for GR4H (hourly),
RunModel_GR2M() function for GR2M (monthly) and
RunModel_GR1A() function for GR1A (yearly).SeriesAggreg() to easily aggreg
timesteps.ErrorCrit_RMSE() which led to incorrect
calibration (the criterion was maximised instead of minimised).CreateRunOptions(),
CreateCalibOptions() and plot_OutputsModel()
to handle the new models.plot_OutputsModel() function to
allow a selection among available plots.ErrorCrit_KGE() and
ErrorCrit_KGE2() to handle case when only one values in not
NA.BasinData objects.quiet = TRUE or
FALSE) to choose if the warnings should be suppressed or
not.CalibrationAlgo_*() functions were renamed into
Calibration_*().CreateCalibOptions() to handle models with
only one parameter.Calibration_HBAN(). The function was not
working properly with models having only one parameter.MeanAnSolidPrecip
for each elevation layer. The CreateRunOptions() function
is impacted.CreateInputsModel()
and DataAltiExtrapolation_HBAN().CreateInputsModel(), CreateRunOptions(),
CreateInputsCrit(), CreateCalibOptions().ErrorCrit() functions to better
account for the cases with constant flow values or local zeros.plot_OutputsModel function (to
handle 0 in Qobs and Qsim).nsga2() (the
‘mco’ package is required).EfficiencyCrit() have been replaced by
ErrorCrit() to avoid misunderstanding (by default, the
algorithms minimise the error criterion).CalibrationAlgo_HBAN()
function was not working in the very rare case when the diagonal search
was activated and lead to a set outside the authorised range.CalibrationAlgo_HBAN() function
was not working properly with models having only one parameter.ModelDefaultIniOptions()
function (this bug was introduced in the RC7 and caused an error when
IndPeriod_WarmUp = NULL.ModelDefaultIniOptions()
function (the automatic selection of one year for warm-up was not
handling properly missing data).BasinData. This will greatly ease the future use of Rdata
files (instead of txt files) as storage format for the time series of
observation.DefineFunctions_Model() or
DefineFunctions_ErrorCrit()).SelectPer arguments are replaced by
IndPeriod to ease understanding.PE arguments are replaced by PotEvap()
to ease understanding.Fsol arguments are replaced by
FracSolidPrecip to ease understanding.SelectPer_Run() is continuous is now
made in the CheckArg() functions.OutputsAlgo().EfficiencyCrit_NSE_sqrtQ() function to compute NSE
criterion on sqrt flows.RunModelAndCrit from
CalibrationAlgo_optim_stats and
CalibrationAlgo_nlminb_stats.CalibrationAlgo_nlminb_stats argument was wrongly
defined in DefineFunctions_CalibrationAlgo()
(optim instead of nlminb).RunOptions was incorrectly made in
CheckArg() function.These binaries (installable software) and packages are in development.
They may not be fully stable and should be used with caution. We make no claims about them.