18/06/2012

Compte rendu du Point Hebdomadaire LMDZ du 18/06/2012

Durée 1h10

Laurent
-------
- Les données AMIP et transpose-AMIP nouvelle physique publiées.
- A travaillé avec Yann sur la version basse mémoire.
2 bugs trouvés (un lié à un compilo + exigent et l'autre
à l'utilisation d'une FFT sur vargas).
Tourne en MPI sur vargas.
Soucis avec le mode mixte.
Va commencer à tourner en MPI seul.
- réunion à Grenoble sur l'intégration de Mar (Laurent et Frédéric
présent). Deux options : 1) interface dynamique de MAR // (prête
mais sans forçage réaliste) couplé à phylmd. 2) intégration
de bouts de la physique de MAR directement dans phydev.
Dans le cas 1), études tropicales
2), études polaires.

Ehouarn
-------
RAS

Véronique
---------
Attends un retour sur le dernier questionaire + traduction
"lmdz en bref"
- Nouvelle testing.
- lien vers réunion utilisateurs
- Changement de fonction de Josefine.
- Annonce du questionaire
- déclenchement stochastique

Marie-Pierre
------------
A repris des install avec Arpege. Tout marche.
A terminé les runs d'intercomparaison Euclipse (672 niveaux).
Doit faire aussi des runs à résolution plus standard.

Lionel
------
Passe

Camille
-------
Tests ancienne et nouvelle physique avec les isotopes.
- A des résultats bizares en différence ancienne/nouvelle physique
dans des simulations bizares.
- 3 stagiaires sur equilibre radiative/convectif ou TWPice.
Tests sur humidité, isotopes.
- Couplage avec orchidee / isotopes : va sans doute mettre les isotopes
dans la nouvelle version. Une post-doc viendra mettre des isotopes
dans cette verison (mergehydro). CDD russe pour faire des
simulations isotopes zoomées ou globale (à partir d'octobre).

Jean-Yves
---------
Les Toulousains vont faire tourner des LES (ou CRM) radiatif/convectif

Frédéric
--------
A modifié le script de création des versions du modèle pour install.sh
pour ajouter des tests vérifiant : 1) si le bench 48x36x19,
2) convergence avec la version précédente, 3) test 1+1=2,
4) test de l'initialisation. Tourne pour le moment avec g95
Résultat en sortie de http://www.lmd.jussieu.fr/~lmdz/DistribG95

install version         SVN     Bench   Conv.   1+1=2   Init
run     Nnum.
/prev.
20120504.trunk          1624    OK      OK      -       OK
20120515.trunk          1626    OK      -       OK      OK
20120516.trunk          1628    OK      OK      OK      OK
A rencontré un problème sur 1+1=2 si idissip=iperiod.
2 problème : idissip devrait être un diviseur de day_step
indépendamment de ça

Catherine
---------
Meeting Embrace. Chouet' projet.
* Intercomparaison 1D. 4 modèles. NP pas trop mal mais un peu tôt
Mais pas du tout de nuages convectif. Sans doute un pb lié
à la PDF. Problème aussi que les thermiques sont tués trop
tôt. beta=0.01.
* Forcer les modèles dynamiques avec des taux de chauffage (F. Codron
impliqué) pour les moussins africaine et indienne.
* Simulation mini-planètes.
* MJO. Jean-Philippe Duvel. Test de différentes hypothèses.
* Le WP4 est de faire l'évaluation des modèles. Intégrer des codes
qui calcule des diagnostics. Catherine de l'évaluation du
cycle dirune.
Meeting LES grands domaines.

Suivi des questions sur LMDZ-SVP
================================
- Question de Yves Balkanski et Anne Cozic. Modèle MR plante.
- Commentaires sur grid_noro (Alain. Faire relire)
- Remplacement de la routine de diagnostics de couche limite hbtm.F
(hbtm2l, physiq, pbl_surface_mod et les deux phys_output_mod et write)
- Commission de dump2d (Alain)
- ajout d'un argument dans radlwsw.F (Alain)
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
- Commentaire de run.def à inclure dans les fichiers histoire.
Sans doute pas possible par ioipsl (Alain)
- thermcell_old (Anne Cozic)


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

Enquête sur l'utilisation de LMDZ
---------------------------------
Questionnaire prêt.
Responsable d'équipes presque identifiés.
On valide


Problème de convergence entre 2 compilateurs gfotran
----------------------------------------------------
Problème identifié pour des cas 1D de nuages Toga
Problème de plantage en 3D.
Résultat des courses :
les résultats convergent entre g95 et les versions récentes de gfortran
Mais plantage et résultats différents avec gfortran.4.1.6

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 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
A priori : svn1626,7 ou 8 après vérification.

1+1=2 :
-------
OK avec ancienne physique à partir de svn1626
A faire : convection et poches (OK pour MY et thermiques).
A automatiser dans cretion_modipsl.sh (Frédéric)

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/06/11 : tester ce0l.e avec le dernier gfortran
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
----------------------------------------
En cours

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
---------------------------------
Moyennes à faire à la main.

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