07/05/2012

Compte Rendu du Point Hebdomadaire LMDZ du 07/05/2012

Durée 1h25

Jean-Yves :
-----------
Quasiment rien fait. A demandé 200000 heures sur curie.
N'a pas commis.
A continué à travailler sur le nouveau déclenchement de Nicolas.

Catherine :
-----------
A remis le nez dans le code des thermiques.
Deux critères de fin des thermiques en haut :
- annulation de w.
- par détrainement.
Avant c'était seulement w. Maintenant souvent le détrainement.
A ajouté un calcul de zmax pour le second cas.
En a profité pour tester le fait d'utiliser zmax-humide plutôt
que zmax-sec pour la fermeture du modèle du thermique.
Gros impact sur les résultats sur le cas AMMA.
Doit resoumettre le papier révisé avant jeudi.

Marie-Alice :
-------------
Nouvelle couplé.
libigcmv2beta1 : possibilité de tourner en créant beaucoup moins
de fichiers. Regroupement et tar.
Simulation historical 150 ans : 120 000 -> 7 000 fichiers.
Cours début juin pour le lancement des simuls.
Avec les configs revues (notamment simplifiées par Josefine).
Déménagement des fichiers du CCRT prochainement.
Période dure : le CCRT passe son temps à dire qu'on bosse mal.
Merci à Sébastien.
Souhaite se remettre sur curie sur l'aqua-planète.

Lionel :
--------
Correctif sur l'option qsave.
Elle n'est pas par défaut avec mpxlf contrairement à ce que
croyait Lionel.
Il y avait aussi un problème de définition de la taille
des stacks dans les scripts.
En fait, on peut tourner jusqu'à 4 procs sur vargas pour la BR.
Commission pour pouvoir choisir entre le calcul des pressions
par exner ou au milieu des couches.
Commission cette semaine d'un flag.

Adriana :
---------
Regarde les biais des vents.
Toujours des problèmes sur le Sahel.
A trouvé une façon de corriger les vents pour les émissions.
Par points de grille :
séries de temps 6 horaires de LMDZ et ECMWF sur 20 ans
Calcul du rapport entre la moyenne du module de chaque série.
La carte de ces coefs est ensuite réappliquée sur le vent
du modèle pour les émissions.

Ehouarn :
---------
Test du lundi svn 1624
A traduit plein de pages du site web LMDZ en anglais.

Josefine :
----------
N'a pas oublié qu'elle avait promis de faire une testing.
A mis à jour les configs LMDZOR
Propose de ne conserver que la dernière config (3 autres obsolètes).
Page wiki à disposition LMDZORv4.
http://forge.ipsl.jussieu.fr/igcmg/wiki/LMDZOR_v4
dont explications pour refaire tourner CMIP5.

Marie-Pierre :
--------------
RRTM :
Sébastien (thésard de Marion) a chargé sur ciclad le modèle.
Tourne en 3D avec RRTM.
Travaille sur le cas de comparaison Euclipse.
Cas AMMA : il faudrait le faire tourner avec un modèle de sol.


Frédéric :
----------
Retour sur le papier ancienne résolution / LMDZ5A : révisions
mineurs. Ouf.
Travaille sur 1+1=2.
1+1 = 2 avec l'ancienne physique mais sans rayonnement.
Commission d'un if (iflag_pbl/=0) dans physiq.F pour ne pas
appeler
A travaillé sur l'enquête LMDZ avec Véronique.



Suivi des actions en cours !
===========================!


Enquête sur l'utilisation de LMDZ
---------------------------------
A lancer.

Commissions pour Alain
----------------------
1. remplacement de la routine
de diagnostics de couche limite hbtm.F
(hbtm2l, physiq, pbl_surface_mod et les deux phys_output_mod et write)
2. ajout d'un argument dans radlwsw.F
et donc de son appel dans physiq.F (le bilan au sol radsol est
le seul argument sorti qui n'a pas sa version clear-sky, ce qui représente une discrimination peu supportable)
Le tout à partir de la LMDZ20110921.trunk
Les sources, README et doc pour HBTM2L sont dans /d1/alain/LMDZ5V5_mods

Problème de convergence entre 2 compilateurs gfotran sur le cas toga
--------------------------------------------------------------------
Faire un test avec un cas plus simple.

Ajout d'une sous surface pour les lacs :
----------------------------------------
A faire.

Problèmes de dates dans les fichiers de sorties
-----------------------------------------------
Problèmes de dates dans les AMIP de Lidia et les aquaplnètes
d'Abderrahmane.

Lecture des PFTs :
------------------
2012/03/05 : procédure.
Historical distribués OK
Regarder le max des PFT 12 et 13
On voit l'évolution des terre agricoles.
En monitoring par défaut dans le couplé.
Pas disponible pour le moment dans les AMIP LMDZ.
Il faut refaire tourner le monitoring (Abderrahmane).
Disposer d'un diagnostic permettant de vérifier qu'on a fait ce
qu'on croyait en terme de forçage des PFTs (Laurent).
Vérification faite sur les historical et RCP IPSLCM5A
A faire sur les AMIP IPSLCM5A et sur tous les IPSLCM5B.
Ajouter El Nino.

Problèmes de convergence en parallèle :
---------------------------------------
Rm  de Marie-Alice : essayer sur vargas la derniere version de la
version 1.12 du compilateur.
Convergence des versions parallèles. Se donner les moyens de
remplacer les routines suspectes de changer quand on passe en
parallèle. Faire un clone sequentiel qu'on puisse activer.
Commencer par décortiquer le problème en séparant dynamique/
physique, en regardant si on a aussi le problème sur quelques jours.
Problème variable s_therm :
Variable non reproductible dans les fichiers en histoire moyenne.

Conservation avec la nouvelle physique :
----------------------------------------
2012/03/05 : travail en cours (Frédéric cf plus haut).
A regarder.
Lien avec les problèmes identifiés en 1D.
Problème identifié dans les sorties, en 1D.

Prochaine testing :
-------------------
Prochaine testing avec 1+1=2
Non conservation numérique (1+1=2 dans la dynamique, et changement
inidissip).

1+1=2 :
-------
2012/03/05 : Problème avec la FFT : en mode "FFT", le premier
appel est fait avec la transformée MxV.
2012/02/27 : Bug dans ferret.
Problème trouvé dans la dynamique, au moment de l'appel à la dissipation.
Le problème existe aussi sur la prod.
On va commettre.
Josefine va reprendre à partir de là pour la physique.


Réintégrer les sources dyn3dpar à mémoire réduite de Yann :
------------------------- ----------------------------------
2012/03/05 : Laurent a récupéré les sources.
Laurent prend la main.
Discuté avec Yann. Il s'en occupe dans 15 jours si on lui envoie un bench. (Josefine s'occupe du bench).
Josefine a envoyé un bench pour le guidage.

create_etat0_limit
------------------
2012/03/05 : Josefine va faire en sorte qu'on puisse tourner avec
un seul proc.
2012/03/05 : Modif permettant de recréer un fichier limit
à partir d'un restart du couplé, permettant d'éviter le o2a.
Ou au moins vérifier la cohérence de restart et limit.
Démarrage en nouvelle physique : Laurent compte s'y remettre.
Sur vargas : mettre une verue pour que create_etat0_limit se lance
automatiquement sur un seul processeur.
Comprendre pourquoi ça plante sous scientific-linux avec gfortran.

Nouveau rayonnement
-------------------
Discussion après le poihl.
Jean-Louis, Lionel, David, Marion, Abderrahmane, Marie-Pierre,
thésard de Chantal Claud.
On essaie de re-faire tourner ce qui existait en 1D et 3D.

Réglage d'une nouvelle version du modèle
----------------------------------------
coef_eva
S'organiser
Reprendre

IOs
---
Transformer phys_output_write.h en une routine.
Travail en cours (Laurent).
Pouvoir tourner sur un an avec un seul histmth.
Réunion sur le nouveau serveur.
A plus court terme : pouvoir tourner sur un an.
Prochain poihl ouvert autour des IOs.
6 février : présentations + discussion

Développement d'une nouvelle version de la physique.
----------------------------------------------------
Groupe de travail : Laurent, Ehouarn, Frédéric, Lionel, Hubert
Question : que fait on de dudyn et pvtheta dans les arguments
de physique ?

Faire tourner par période de 1 an
---------------------------------
En cours

Inclusion des isotopes
----------------------
Identification d'une façon jolie d'inclure les isotopes dans la
dynamique.
Codé chez Camille.