30/04/2012
Compte rendu du Point Hebdomadaire LMDZ du 30/04/2012
Durée 50min
Lidia
-----
* AMIP nouvelle physique CMIP-5
Amip4CO2 et 4K plantaient -> Lancement d'un contrôle de 10 ans.
- Ne diffère que par l'état initial et par le fait qu'on n'active pas stomate.
- Plante si on ne change que stomate.
* A transferré la simulation Cordex sur dzt de vargas.
L'idris a transferré. Mais dzt trop petit.
N'a transferré que le minimum vital (2To).
Le reste (30To = les hf) est encore sur ces.
Car 6To seulement sur dzt.
Cordex relancé. 1 mois est passé.
Camille
-------
Simulations avec isotopes ancienne physique avec ou sans guidage.
Marche pour O18 et Deuterium et pas pour O17 (résultats bizares).
Nouvelle phsyique : 5 mois mais plante pour raisons isotopiques.
LMDZ-Sechiba. Incompatibilité entre orchidee-isotopes et LMDZ.
A commencé à prendre en main LMDZ1D. Ajouté la possibilité d'imposer
omega. A rajouter dans la version standard.
A essayé avec des omega de -90 à +30 hPa/jour.
Problèmes d'eau négative dans la haute tropo. Regarder le schéma
d'advection verticale de l'eau.
A faire : rentrer les isotopes dans les versions standard du code.
Dans la dynamique : on attend la nouvelle version "eco" de la
dynamique parallèle.
Dans la physique : reprendre les discussion. Faire le lien avec
le travail fait avec Romain.
Jean-Yves
---------
A refilé la glace à Catherine. A faire tourner dans TWPice.
Travaille sur le déclenchement de Nicolas.
A tourné en 3D. Mais la convection ne déclenche plus du tout.
Va passer au splitting de la couche limite entre l'intérieur et
l'extérieur des poches.
Priorités : commissions déclenchement.
tester la glace en 3D.
Frédéric
--------
RAS
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.
Lidia
-----
* AMIP nouvelle physique CMIP-5
Amip4CO2 et 4K plantaient -> Lancement d'un contrôle de 10 ans.
- Ne diffère que par l'état initial et par le fait qu'on n'active pas stomate.
- Plante si on ne change que stomate.
* A transferré la simulation Cordex sur dzt de vargas.
L'idris a transferré. Mais dzt trop petit.
N'a transferré que le minimum vital (2To).
Le reste (30To = les hf) est encore sur ces.
Car 6To seulement sur dzt.
Cordex relancé. 1 mois est passé.
Camille
-------
Simulations avec isotopes ancienne physique avec ou sans guidage.
Marche pour O18 et Deuterium et pas pour O17 (résultats bizares).
Nouvelle phsyique : 5 mois mais plante pour raisons isotopiques.
LMDZ-Sechiba. Incompatibilité entre orchidee-isotopes et LMDZ.
A commencé à prendre en main LMDZ1D. Ajouté la possibilité d'imposer
omega. A rajouter dans la version standard.
A essayé avec des omega de -90 à +30 hPa/jour.
Problèmes d'eau négative dans la haute tropo. Regarder le schéma
d'advection verticale de l'eau.
A faire : rentrer les isotopes dans les versions standard du code.
Dans la dynamique : on attend la nouvelle version "eco" de la
dynamique parallèle.
Dans la physique : reprendre les discussion. Faire le lien avec
le travail fait avec Romain.
Jean-Yves
---------
A refilé la glace à Catherine. A faire tourner dans TWPice.
Travaille sur le déclenchement de Nicolas.
A tourné en 3D. Mais la convection ne déclenche plus du tout.
Va passer au splitting de la couche limite entre l'intérieur et
l'extérieur des poches.
Priorités : commissions déclenchement.
tester la glace en 3D.
Frédéric
--------
RAS
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.