#41 closed defect (fixed)
restartabilite : 1+1=2
| Reported by: | jghattas | Owned by: | Laurent Fairhead |
|---|---|---|---|
| Priority: | major | Milestone: | |
| Component: | LMDZ | Keywords: | |
| Cc: |
Description
Avec 1+1=2 on entend dire que quand on tourne 1j puis 1j on obtient les memes restart que si on tourne 2j de suite. Les restarts doivent etre identique numeriquement.
LMDZ actuellement n'a pas 1+1=2 quand on tourne avec avec la physique.
Depuis revision 1614 le dynamque a 1+1=2, c.a.d. un run avec iflag_phys=0 donne les meme restart.nc. => => MAIS pour avoir les resultats identique il ne faut pas utilise le filtre FFT. Il faut donc mettre use_filtre_fft=n.
Change History (3)
comment:1 Changed 13 years ago by
comment:2 Changed 9 years ago by
| Resolution: | → fixed |
|---|---|
| Status: | new → closed |
Note: See
TracTickets for help on using
tickets.
Pourquoi le FFT ne donne pas 1+1=2, mail de Yann Meurdesoif le 28/02/2012 :
Bonjour,
Je sais pourquoi. Dans leapfrog_p, avant la boucle en temps, on fait appel une fois à filtreg séquentiel. C'est OK puisque les données sont encore globale. Or le filtre séquentiel utilise obligatoirement le filtre classique. Or le filtre classique n'est pas parfaitement équivalent au filtre FFT (numériquement)... Normalement, ce problème est résolu dans la version à mémoire réduite. Je pense que la correction est assez simple : remplacer filtreg par filtreg_p.
Yann M.