Correct small issue in the vignette
See below (in French):
Salut Guillaume,
En fait en lisant la vignette je n'ai pas compris ce qu'était QsimEns : je croyais que c'était le débit d'ensemble juste après mise à jour des états à partir des débits mesurés à J+0, avant de passer au pas de temps suivant, alors qu'en fait c'est la prévision à J+0 à partir des débits mis à jour à J-1 avec forçage d'ensemble, ce qui est je trouve un peu complexe pour ne pas être explicité. Je pense qu'il serait utile qu'il y ait une phrase expliquant que la performance des modèles est comparée en prévision à J+1 avec prévisions météo parfaites perturbées.
Effectivement tu as raison, ce n'est pas du tout la même chose que l'open-loop. Et je suis d'accord, celle-ci devrait être calculée avec RunModel_DA avec DaMethod="none" plutôt qu'avec RunModel pour avoir les perturbations.
Merci !
A bientôt, Paul
De: "Thirel Guillaume" Envoyé: Vendredi 17 Décembre 2021 13:58:29 Objet: Re: Commentaire sur vignette airGRdatassim
Salut Paul,
Je ne suis pas certain de comprendre ta remarque à 100 %. En effet, il s'agit des débits simulés à partir de la mise à jour du pas de temps précédent qu'on retrouve dans QsimEns. Quand tu dis " on compare la même chose avec ou sans filtre ", ce n'est du coup pas la même chose, car l'un a l'assimilation de données, l'autre non. L'open loop c'est juste mettre de la météo perturbée et laisser le modèle faire, du début à la fin de la période. Donc à mon avis, la chose incorrecte c'est que dans le code l'open loop ça devrait être du RunModel_DA avec DaMethod = "none", non ?
A+ Guillaume.
De: "Royer Gaspard Paul" Objet: Commentaire sur vignette airGRdatassim
Bonjour Olivier, Guillaume,
Je souhaitais vous signaler une petite imprécision dans la vignette de airGRdatassim https://gitlab.irstea.fr/HYCAR-Hydro/airgrdatassim/-/blob/master/vignettes/get_started.Rmd
A la section "Comparative performance Assessment", il est écrit qu'appliquer les codes capturés en PJ permet de comparer les débits d'ensemble open-loop et avec application d'un filtre. En réalité, la variable QsimEns renvoyée par la fonction RunModel_DA correspond aux débits simulés en perturbant les entrées météo avant mise à jour des états. Donc on compare la même chose avec ou sans filtre.
Il faudrait peut-être préciser que ce ne sont pas les débits mis à jour, quand vous aurez le temps.
A bientôt, Paul