croco + AGRIF : point

(Jan 2023)

Table of content

A faire

Bug AGRIF de repro MPI sur cas particulier

Pas de plantage mais on perd la repro MPI sur certains découpages MPI

Cas où coincides
1. frontière grille mère/grille fille
2. frontière sous domaine MPI
3. frontière entre un point mer et un point terre.

Rachid:
Ce cas n'est pas traité je pense dans AGRIF et ça doit se coder au plus profond d'AGRIF. Seul Laurent Debreu pourrait faire ça ou m'expliquer comment faire sans trop perdre de temps, mais le seul moyen de le motiver c'est de le croiser. Il a plutôt tendance à se décharger après unne grosse surcharge.

Rachid pense que ça n'a pas été codé non plus dans le AGRIF de NEMO mais que comme c'est assez rare de tomber sur ce cas, c'est passé innaperçu.

Demander quand même à Jérôme Chanut?
Rachid essaiera d'en parler à Laurent s'il le croise.

Tant que ce problème n'est pas résolu, il suffit de vérifier la reproductibilité MPI qui garanti que le découpage qu'on utilise ne tombe pas sur ce cas particulier.

BIO_HADV_WENO5 - Bug AGRIF de repro MPI (pb halo de 3)

Schéma ordre 5 => halo de 3

Rachid :
AGRIF est prêt pour ça, mais croco ne l'appelle pas forcément correctement.
Ca devrait être relativement facile d'y remédier.
En attendant les résultats sont corrects à priori.

Mais toujours des sources de confusion quand on rencontre un problème de ne pas avoir de résultats reproductibles!!!

Quid de ?

# define BIO_HADV_WENO5
# define TS_HADV_WENO5
# define TS_VADV_WENO5
# define UV_HADV_WENO5
# define UV_VADV_WENO5
# define W_HADV_WENO5
# define W_VADV_WENO5
# define BEDLOAD_WENO5

Olivier :
On a attendu pendant très longtemps ce schéma WENO5
WENO5 a réglé beaucoup de problèmes. En plus il est très peu diffusif.
Le problème c'est qu'il n'est pas complètement positif.
Du coup PISCES est obligé de mettre les traceurs à 0 s'ils sont négatifs
=> Pas terrible parce qu'on génère de la masse.
Olivier ne comprend pas pourquoi ils n'ont pas mis un schéma TVD 2eme ordre tout simplement. Absolument positif et beaucoup moins coûteux.
L'advection représente une grosse part dans PISCES.
Il faudrait regarder quand on a des valeurs négatives et si ça arrive souvent. Avec un simple write ou en faisant la somme. Peut-être près de la côte quand l'ordre du schéma diminue par manque de valeurs quand on a plus un halo de 3

WETDRY - Bug AGRIF sur h

Rachid :
Ce bug devrait être assez facile à trouver.

Xavier:
Ca me rassure presque parce que ça fait un moment qu'on pense qu'il y a un bug dans WETDRY

IRENE (AMD & Intel) - pb memory

Un problème de Freed memory apparait avec DDT mais l'execution ne pose pas de problème sans DDT.
Mais les tests sur les bugs AGRIF sont fait sans PISCES. On redoute de retrouver une fois PISCES activé, les problème de deallocate dans les routines PISCES.

Remettre la hotline du TGCC sur le coup maintenant que nous sommes sur un cas plus simple (sans PISCES, ni bulk online, ni rivers... ) et que tout fonctionne sur JEANZAY avec quasiment la même version de DDT.

PISCES - bug avec initialisation des traceurs sur BVHR

Ces problèmes rencontrés avec des jeux de forçages différents entre Rachid et Gildas ne sont ils pas liés à des différences de comportement en fonction du nombres de traceurs pisces.
Regarder si JEANZAY passe avec les jeux de Gildas qui ne passaient pas sur IRENE...

A reprendre une fois que tous les bugs AGRIF sont résolus... ou qu'on est dans un découpage où la reproductibilité MPI est assurée...

Manque

S'assurer de la repro MPI en activant tout ensemble!!