09/02/2015

Compte Rendu du Point Hebdomadaire LMDZ du 9/2/2015

Durée : 1h50

LIONEL :
-----------
Propose 2 modifications dans les scripts de compilation (makelmdz et makelmdz_fcm) :
- 1ère concerne la valeur de iim qui doit être un multiple de 8
- 2ème concerne la compilation avec dyn3dpar qui est actuellement celle par défaut :
Propose de mettre un arrêt de la compilation (stop) si l'option -mem n'est pas présente !

Ionela :
----------
A lancé 2 tests betaclim sur curie :
- avec ancien rayonnement qui plante après 11ans et demi
- avec RRTM qui plante après 361 jours. Dans ce test, elle n'a pas "1+1=2"
A reçu des fichiers d'observations ocean-atmosphère où elle a calculé la moyenne obs² dans le but de les inclure dans les atlas !
A récupéré les scripts NCO pour les post-traitement MGO, les a modifié et les a testé sur 360j.
A besoin pour utiliser ces scripts d'avoir en sorties journalières : u et v à 850 et 200mb
A vérifier s'il n'y a pas de problème sur les champs niveau pressions standards sortis via XIOS !

Lluis :
--------
A tourné la nouvelle physique de LMDZ dans WRF sur une aqua-planete et ça plante au bout d'une quinzaine de jours dans hgardfou apres la convection.
Ça plante sur de l'ocean a 273K avec une couche limite tres mince (3 niveaux pour les thermiques). Il regarde ça avec Jean-Yves


Jean Yves :
--------------

Avec Lluis : Lluis fait tourner la nouvelle physique de LMDZ dans WRF sur une aqua-planete ; ça plante au bout d'une quinzaine de jours dans hgardfou apres la convection. Ça plante sur de l'ocean a 273K avec une couche limite tres mince (3 niveaux pour les thermiques). On soupçonne que le probleme est lié à cette couche limite très mince (aussi mince que la couche d'alimentation de la convection, qui est de 300m). On a mis des prints dans cv3_yield (qui calcule les tendances convectives). Il apparait que le le flux de masse montant atteint 7 10^6 kg/m2/s (alors qu'on attend 0.01 ou 0.1). On pense donc être près de la découverte du bug.
Avec Catherine : on a élaboré un nouvel épluchage des ascendances adiabatiques dans le schéme d'Emanuel. Avec la nouvelle formulation on peut changer la dependance sur w (la vitesse verticale) de l'epluchage sans changer la fermeture. On devrait ainsi pouvoir explorer divers profils verticaux d'entrainement, depuis un entrainement proportionnel à w (cas actuel) à un entrainement proportionnel à 1/w2 (formulation utilisée dans les modeles de cumulus).


Marie-Pierre :
------------------
A fait plusieurs install.sh
A des soucis avec la version 2188 (print klon/2 qui doit être remplacer par klon/2 +1)
A la convergence numérique à partir de la version 2191
A préparer avec Ionela un fichier de forçage 1D à partir de celui 3D

Laurent :
-----------
 2 runs en machine :
- 59 niveau a fait 79 ans et plante et reparti en divisant par 2 les tetag sur 1 mois
- 79 relancé avec iflag_t_glace = 0, s'est planté au bout de 6 ans sur un point en Alaska toujours dans lsc malgrès les dernières modifs de Fréderic.

Eheouarn :
--------------
A corriger un bug sur la définition (double/simple précision) dans le fichier limit.nc des variable pour les aquaplanètes
Ajout d' "implicit none" là où il n'y en avait pas (grâce aux options debug de pgf90 pour savoir où). =
La fonction "erf" utilisée via pgf90 n'est que celle utilisée en simple précision; compiler avec -r8 ne suffit pas pour utiliser la version double précision, il faut alors utiliser "derf" dans le code. Lionel va regarder pour trouver une solution (en important la fonction erf de slatec?)
Pour info a assisté à la conf de l'ANR HET où est discutée l'interface dynamique/physique de LMDZ
A récupéré les routines modifiées par Francis pour le slab ocen !

David :
---------
A évoqué le script simplifié pour lancer ou enchainer des simulations et qui permet de faire des rebuilds !
A ajouter à la boite à outils LMDZ

Adriana :
------------
A reçu les données Bodex du 27-02 au 13-05-2005 avec des champs de vent d'une fréquence d'1 minute !

Abderrahmane :
---------------------
A une simulation en machine qui a fait 20ans et qui continue de tourner : NPv5.110 avec l'ancien rayonnement
Continue les tests sur le cas qui plante (NPv5.11 avec RRTM) :
 - a lancée cette simulation sans Orchidee
 - a fait un test en remplaçant "omp single" par  "OMP MASTER" avant l'ouverture du fichier ozone
Travail sur le portage des scripts de post-traitement sur curie : Une vingtaine de scripts a adapter en fonction de la machine !
A mis à disposition une partie des post-traitements (bilans, moyennes zonales, classifications en régimes de temps) pour 2 simulations NPv5