Opened 3 years ago
Closed 3 years ago
#131 closed incoherences (fixed)
Appel Guidage (leapfrog)
Reported by: | Ehouarn Millour | Owned by: | Ehouarn Millour |
---|---|---|---|
Priority: | major | Milestone: | |
Component: | LMDZ | Keywords: | |
Cc: |
Description
Message de Francis, le 19/03/2021
Cher poihl, Le guidage est appelé si j'ai bien compris à chaque pas matsuno forward, en particulier au tout premier pas de chaque run, et aussi apparemment au tout dernier. Du coup, on appelle 2x au "même" pas de temps à la fin d'un run et au début du suivant... Pas très grave, sauf que quand on sauve les tendances on se retrouve avec une sortie doublée ce qui est un peu pénible. La solution la plus simple semble être de modifier la ligne dans leapfrog_loc : if (ok_guide .AND..NOT.(itau.EQ.itaufin)) then call guide_main(itau,ucov,vcov,teta,q,masse,ps) ...Mais peut-être que je rate quelque chose... Est-ce que ça semble raisonnable ? Francis
Le problème est probablement plus profond... en l'état le guidage n'est appelé que durant la phase "forward" du pas Matsuno,ce qui n'est pas cohérent... il faudrait soit l'appeler deux fois (à chaque étape du Matsuno) soit plus simplement l'appeler seulement après un pas "backward".
La sortie supplémentaire en fin de run est possiblement lié au fait que l'on fait un (demi) pas Matsuno (la partie forward) pour rien en fin de boucle temporelle.
Change History (3)
comment:1 Changed 3 years ago by
comment:3 Changed 3 years ago by
Resolution: | → fixed |
---|---|
Status: | assigned → closed |
Note: See
TracTickets for help on using
tickets.
Avec un peu de recul...
En fait pas de gros soucis avec l'application du guidage; elle se fait bien avant l'intégration dynamique du pas "forward" (et pas entre pas "forward" et "backward", ce qui aurait été un problème).
Le plus simple semble donc de simplement éviter l'appel à guide_main lors du pas forward "superflu" en fin de run, comme préconisé par Francis.
Note au passage: cet appel à un pas forward "superflu" en fin de run est probablement là pour pouvoir alors appeler fluxstokenc dans le cas "offline".