<li><p>An article describing the airGR package has been published. Its reference has been added and will be displayed with <code>citation(“airGR”)</code>.</p></li>
<li><p>Added <code>CreateIniStates()</code> function in order to help user to format <code>IniStates</code> argument for <code>CreateRunOptions()</code>.</p></li>
<li><p>Added (<code>Param_Sets_GR4J</code>) dataset. It contains generalist parameter sets for the GR4J model.</p></li>
<li><p>Added <code>CreateIniStates()</code> function in order to help user to format the <code>IniStates</code> argument for <code>CreateRunOptions()</code>.</p></li>
<li><p>Added the <code>Param_Sets_GR4J</code> dataset. It contains generalist parameter sets for the GR4J model.</p></li>
<li><p>Three vignettes have been added. They are relative to different calibration methods (including the generalist parameters sets of the GR4J model).</p></li>
</ul>
</div>
<div id="bug-fixes" class="section level4">
<div id="bug-fixes-1" class="section level4">
<h4>Bug fixes</h4>
<ul>
<li><p>Fixed bug in <code>RunModel_GR4H()</code>: in <code>frun_GR4H</code> Fortran subroutine, <code>St(2)</code> is now set to 0 (and not <code>St(1)</code>) when <code>St(2) < 0</code>.</p></li>
<li><p>Fixed bug in <code>RunModel_GR4H()</code>: in <code>frun_GR4H</code> Fortran subroutine, <code>St(2)</code> is now set to 0 (instead of <code>St(1)</code>) when <code>St(2) < 0</code>.</p></li>
<li><p>Fixed bug in <code>plot.OutputsModel()</code> for the regime plot when the period is less than 1 year.</p></li>
<li><p>Fixed bug in <code>plot.OutputsModel()</code> when there is no common data to plot the cumulative frequency or the correlation QQ.</p></li>
<li><p>Fixed bug in <code>plot.OutputsModel()</code> for y-axis labelling of flows time series when <code>log_scale = TRUE</code> and <code>BasinArea</code> used.</p></li>
<li><p>Fixed bug in <code>plot.OutputsModel()</code> for the y-axis labelling of flows time series when <code>log_scale = TRUE</code> and <code>BasinArea</code> is used.</p></li>
<li><p><code>RunModel_GR4J()</code>, <code>RunModel_GR5J()</code> and <code>RunModel_GR6J()</code> (and <code>CemaNeige_GR*J()</code>) now return Ps, Pn and actual exchanges. See the model Fortran codes for more details about the calculation of these variables.</p></li>
...
...
@@ -236,18 +264,18 @@ div.tocify {
<li><p><code>RunModel_GR5J</code> now returns <code>StateEnd</code> in the same order as the other models.</p></li>
<li><p><code>plot.OutputsModel()</code> now returns a warning message when the length of Qobs is different from the length of Qsim.</p></li>
<li><p>The X1, X3 (and X6) parameters from GR4H, GR4J, GR2M, GR5J (and GR6J) are now set to 1e-2 when they are fixed to lower values. <code>RunModel*()</code> functions now return a warning message. <code>RunModel*()</code> functions now return a warning when X4 < 0.5 and set it to 0.5.</p></li>
<li><p>The X1 parameter from GR4H, GR4J, GR2M, GR5J and GR6J, the X3 parameter from GR4H, GR4J, GR5J and GR6J and the X6 parameter from GR6J are now set to 1e-2 when they are fixed to lower values. <code>RunModel*()</code> functions now return a warning message in this case. <code>RunModel*()</code> functions now return a warning when X4 < 0.5 and its value is set to 0.5.</p></li>
<li><p>The commands <code>?L0123001</code>, <code>?L0123002</code> and <code>?L0123003</code> now return the documentation page related to <code>BasinObs</code>.</p></li>
<li><p>Many functions of the package were cleant or slightly modified, with no effect on their outputs.</p></li>
<li><p>The documentation and help of several functions were improved.</p></li>
<li><code>DataAltiExtrapolation_Valery()</code> and <code>CreateInputsModel()</code> now present a PrecipScale argument which allows rescaling precipitation when it is interpolated on the elevation layers when CemaNeige is used.</li>
</ul>
</div>
<div id="bug-fixes-1" class="section level4">
<div id="bug-fixes-2" class="section level4">
<h4>Bug fixes</h4>
<ul>
<li>Fixed bug in <code>DataAltiExtrapolation_Valery()</code>. The elevation gradients for air temperature returned by <code>CreateInputsModel()</code> are improved.</li>
...
...
@@ -307,7 +335,7 @@ div.tocify {
<li><code>ErrorCrit*()</code> functions gain a warnings argument to replace the verbose action and the verbose argument now prints the criterion value(s).</li>
</ul>
</div>
<div id="bug-fixes-2" class="section level4">
<div id="bug-fixes-3" class="section level4">
<h4>Bug fixes</h4>
<ul>
<li>Fixed bug in <code>CreateCalibOptions()</code> when <code>StartParamList</code> or <code>StartParamDistrib</code> arguments are used.</li>
...
...
@@ -329,10 +357,10 @@ div.tocify {
<ul>
<li><p><code>SeriesAggreg()</code> 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.</p></li>
<li><p><code>plot_OutputsModel()</code> gains a <code>log_scale</code> argument in order to plot the flow with a log scale.</p></li>
<li><p>A tutorial is available online on the following link: from <a href="http://webgr.irstea.fr/airGR" class="uri">http://webgr.irstea.fr/airGR</a>. It can also be displayed with the <code>vignette(“airGR”)</code> command</p></li>
<li><p>A tutorial is available online on the following link: from <a href="http://webgr.irstea.fr/airGR" class="uri">http://webgr.irstea.fr/airGR</a>. It can also be displayed with the <code>vignette(“airGR”)</code> command.</p></li>
</ul>
</div>
<div id="bug-fixes-3" class="section level4">
<div id="bug-fixes-4" class="section level4">
<h4>Bug fixes</h4>
<ul>
<li>The value <code>sort</code> for the <code>transfo</code> argument of <code>CreateInputsCrit()</code> was not taken into account. It is now fixed.</li>
...
...
@@ -345,13 +373,13 @@ div.tocify {
<li><p><code>CreateCalibOptions()</code> loses the StartParam argument that was not used.</p></li>
<li>The <code>RunModel_GR6J()</code> and <code>RunModel_CemaNeigeGR6J()</code> models were modified back to versions previous to 1.0.1 to prevent from unwanted efficiency criteria deterioration related to the calibration with <code>Calibration_Michel()</code>. The actual model codes were not modified but the <code>TransfoParam_GR6J()</code> and <code>CreateCalibOptions()</code> functions were modified regarding the X5 parameter. It is strongly advised to use airGR 1.0.2 for the <code>RunModel_GR6J()</code> and <code>RunModel_CemaNeigeGR6J()</code> functions if you are using <code>Calibration_Michel()</code>, as they are much more efficient. In case you were using your own calibration algorithm, you will not notice any difference.</li>
<li><code>CreateInputsModel()</code> and <code>DataAltiExtrapolation_Valery()</code> functions now allow both POSIXt formats (POSIXct and POSIXlt).</li>
...
...
@@ -369,7 +397,7 @@ div.tocify {
<li><p>The silent mode is now defined by the <code>verbose = TRUE</code> argument (formerly <code>quiet = FALSE</code>) in the following functions: <code>Calibration()</code>, <code>Calibration_Michel()</code>, <code>CreateInputsModel()</code>, <code>CreateRunOptions()</code>, <code>DataAltiExtrapolation_Valery()</code>, <code>ErrorCrit()</code>, <code>ErrorCrit_KGE()</code>, <code>ErrorCrit_KGE2()</code>, <code>ErrorCrit_NSE()</code>, <code>ErrorCrit_RMSE()</code>, <code>plot_OutputsModel()</code>, <code>SeriesAggreg()</code>.</p></li>
<li><p>The GR5J model has been modified: previously, two unit hydrographs were used, now only one is remaining. As a consequence, simulations from the GR5J (<code>RunModel_GR5J()</code> function) and CemaNeige (<code>RunModel_CemaNeigeGR5J()</code> function) models will be different.</p></li>
...
...
@@ -377,7 +405,7 @@ div.tocify {
<li><p>The quantiles of the parameters have been recalculated with the new transformations (<code>CreateCalibOptions()</code> function). Since these quantiles constitute the starting point of the calibration algorithm, calibration results will be different.</p></li>
<li><p>Bug fxed in <code>CreateInputsModel()</code> related to the handling of missing values.</p></li>
<li><p>Bug fxed in <code>CreateRunOptions()</code> preventing the correct use of the <code>IniResLevels</code> argument (to manually set the filling rate of the production and routing stores).</p></li>
<li><p>Fixed bug in <code>CreateInputsModel()</code> that was related to the handling of missing values.</p></li>
<li><p>Fixed bug in <code>CreateRunOptions()</code> that prevented the correct use of the <code>IniResLevels</code> argument (to manually set the filling rate of the production and routing stores).</p></li>
<li><p>Update of the functions <code>CreateRunOptions()</code>, <code>CreateCalibOptions()</code> and <code>plot_OutputsModel()</code> to handle the new models.</p></li>
<li><p>Modification of CemaNeige Fortran code to add an update of Gratio after the SnowPack update (no impact on snow simulation).</p></li>
<li><p>CemaNeige users must now specify one <code>MeanAnSolidPrecip</code> for each elevation layer. The <code>CreateRunOptions()</code> function is impacted.</p></li>
<li><p>CemaNeige users can now specify the mean elevation of the input series (before it was always considered equal to the catchment median elevation). The impacted functions are <code>CreateInputsModel()</code> and <code>DataAltiExtrapolation_HBAN()</code>.</p></li>
<li><p>RC11 bug correction: the automatic selection of the warm-up period was not working properly when no data was available from warm-up (i.e. when the user had set the run to start at the very first index).</p></li>
<li><p>RC10 bug correction: the <code>CalibrationAlgo_HBAN()</code> function was not working in the very rare case when the diagonal search was activated and lead to a set outside the authorised range.</p></li>
<li><p>RC9 bug correction: the <code>CalibrationAlgo_HBAN()</code> function was not working properly with models having only one parameter.</p></li>
<li><p>RC8 bug correction of the <code>ModelDefaultIniOptions()</code> function (this bug was introduced in the RC7 and caused an error when <code>IndPeriod_WarmUp = NULL</code></p></li>
<li><p>RC8 bug correction of the <code>ModelDefaultIniOptions()</code> function (this bug was introduced in the RC7 and caused an error when <code>IndPeriod_WarmUp = NULL</code>.</p></li>
<li><p>RC7 bug correction of the <code>ModelDefaultIniOptions()</code> function (the automatic selection of one year for warm-up was not handling properly missing data).</p></li>
<li><p>RC6 correction of the help files (the description of CemaNeige parameters were inverted).</p></li>
<li><p>RC5 differs from previous releases in the way the data are read and stored (in a list instead of individual vectors). The package is similar, only the examples of Main and the files in MyScriptBlocks have changed. All basin data are now stored inside a list named <code>BasinData</code>. This will greatly ease the future use of Rdata files (instead of txt files) as storage format for the time series of observation.</p></li>
...
...
@@ -558,13 +586,13 @@ div.tocify {
<li><p>Add multi-objective calibration using <code>nsga2()</code> (the <code>{mco}</code> package is required).</p></li>
<li>The definition of the generic function is now made in a much simpler way (e.g. see <code>DefineFunctions_Model()</code> or <code>DefineFunctions_ErrorCrit()</code>).</li>
<li>Missing values in Fortran are now -999.999 instead of -9.999.</li>
...
...
@@ -622,10 +650,10 @@ div.tocify {
<div id="new-features-9" class="section level4">
<h4>New features</h4>
<ul>
<li>New <code>EfficiencyCrit_NSE_sqrtQ()</code> function to compute NSE criterion on sqrt flows</li>
<li>New <code>EfficiencyCrit_NSE_sqrtQ()</code> function to compute NSE criterion on sqrt flows.</li>
</ul>
</div>
<div id="bug-fixes-9" class="section level4">
<div id="bug-fixes-10" class="section level4">
<h4>Bug fixes</h4>
<ul>
<li><p>Incorrect arguments in the call to <code>RunModelAndCrit</code> from <code>CalibrationAlgo_optim_stats</code> and <code>CalibrationAlgo_nlminb_stats</code>.</p></li>
