25/06/2012
Compte rendu du Point Hebdomadaire LMDZ du 25/06/2012
Durée 1h15
Marie-Pierre
------------
1. S'est renseignée sur RRTM auprès du thésard de Marion.
A fait tourner sur 3 ans avec 2, 4 ou 6 bandes dans le SW.
Demande des obs. Pas évident que ça ait un sens.
Suggestion : déjà bien documenter la différence entre les versions.
2. Termine l'intercomparaison rayonnement d'Euclipse en 1D basse resol.
3. Relance du bilan Dephy.
Lionel
------
A envoyé un message sur la sortie des flux de masse convectifs.
Va ajouter prw dans tous les cas.
Où en est-on de la confiance en OpenMP ?
Pas disponible avec Orchidee
Jean-Yves
---------
A surtout travaillé sur la thèse de Romain.
Abderrahmane
------------
RAS
Catherine
---------
Papier sur la fermeture : minor revisions !!!!!!!!!!!!!!!!!!!
A regardé sur TWPice les simulations de réglage nouvelle physique.
Sur TWPice, ce qui joue le plus c'est coefevap et le profile du
ratqs changent beaucoup.
Camille
-------
Arrivée d'une stagiaire L3 sur la nouvelle physique.
Laurent (par mail)
------------------
A fond sur l'intégration de la
basse mémoire (pas évident, il faut jongler avec 3 fichiers source en
même temps pour la moindre modification)
Frédéric
--------
A rédigé une réponse à l'éditeur sur le papier de
Question Marie-Alice sur curie
Passage à g95 dans install.sh si compilaeur trop ancien.
Romain a écrit une bonne partie de l'interface de Arpege avec
LMDZ et a installé avec succès (après une petite itération)
le 1D LMDZ chez lui.
Conclusions réunion du comité de pilotage du pôle
=================================================
1. Préparation des configurations futures
On va changer d'océan. Passage au 1/4 de degré a priori 75 niveaux
On vise deux configuration :
a. LMDZ5A-HR (HR=280x280) pour la variabilité
b. LMDZ5B-MR pour le cycle diurne et les tropiques
Sans doute avec Orchidee 11 couches.
2. Nécessite dyn3dmem et IOserveur.
3. Réfléchir à la reproductibilité d'une machine à l'autre à partir
de simulations couplées longues.
4. Discussion du passage aux traceurs. Idée d'utiliser le 1D avec
les traceurs de l'échantillonage. Réunion avec LSCE.
5. Ecrire au moins les bilans MissTerre et Dephy pour le 7 septembre
accompagnés d'une lettre d'intention.
6. Soulevé : la question de l'intégration des isotopes dans le
couplé.
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
---------------------------------
Questionaire 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.
Marie-Pierre
------------
1. S'est renseignée sur RRTM auprès du thésard de Marion.
A fait tourner sur 3 ans avec 2, 4 ou 6 bandes dans le SW.
Demande des obs. Pas évident que ça ait un sens.
Suggestion : déjà bien documenter la différence entre les versions.
2. Termine l'intercomparaison rayonnement d'Euclipse en 1D basse resol.
3. Relance du bilan Dephy.
Lionel
------
A envoyé un message sur la sortie des flux de masse convectifs.
Va ajouter prw dans tous les cas.
Où en est-on de la confiance en OpenMP ?
Pas disponible avec Orchidee
Jean-Yves
---------
A surtout travaillé sur la thèse de Romain.
Abderrahmane
------------
RAS
Catherine
---------
Papier sur la fermeture : minor revisions !!!!!!!!!!!!!!!!!!!
A regardé sur TWPice les simulations de réglage nouvelle physique.
Sur TWPice, ce qui joue le plus c'est coefevap et le profile du
ratqs changent beaucoup.
Camille
-------
Arrivée d'une stagiaire L3 sur la nouvelle physique.
Laurent (par mail)
------------------
A fond sur l'intégration de la
basse mémoire (pas évident, il faut jongler avec 3 fichiers source en
même temps pour la moindre modification)
Frédéric
--------
A rédigé une réponse à l'éditeur sur le papier de
Question Marie-Alice sur curie
Passage à g95 dans install.sh si compilaeur trop ancien.
Romain a écrit une bonne partie de l'interface de Arpege avec
LMDZ et a installé avec succès (après une petite itération)
le 1D LMDZ chez lui.
Conclusions réunion du comité de pilotage du pôle
=================================================
1. Préparation des configurations futures
On va changer d'océan. Passage au 1/4 de degré a priori 75 niveaux
On vise deux configuration :
a. LMDZ5A-HR (HR=280x280) pour la variabilité
b. LMDZ5B-MR pour le cycle diurne et les tropiques
Sans doute avec Orchidee 11 couches.
2. Nécessite dyn3dmem et IOserveur.
3. Réfléchir à la reproductibilité d'une machine à l'autre à partir
de simulations couplées longues.
4. Discussion du passage aux traceurs. Idée d'utiliser le 1D avec
les traceurs de l'échantillonage. Réunion avec LSCE.
5. Ecrire au moins les bilans MissTerre et Dephy pour le 7 septembre
accompagnés d'une lettre d'intention.
6. Soulevé : la question de l'intégration des isotopes dans le
couplé.
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
---------------------------------
Questionaire 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.