23/11/2015
Compte Rendu du Point Hebdomadaire LMDZ du 23/11/2015
Durée : 1h25 minutes
Formation LMDZ :
================
Essayer de passer les scripts d'install parallele en OpenMP seul.
Marie-Pierre va tester.
Discussion par axe :
===================
Axe 1 : Surfaces continentales
------------------------------
- Principal soucis : biais froid
- Test en cours
* Intégré nouvelle discrétisation verticale et code de Fuxing.
L'ancienne discrétisation n'est plus disponible.
En cours d'analyse.
* Nouveau module de neige du CNRM. Commis dans orchidee.
* Albedo : hors albedo de la neige. Travail d'optimisation de Philippe Peylin
non intégré.
* z0 : Fuxing et Nicolas (Vuichard) on travaillé dessus.
Réduit l'évaporation l'été.
- analyse du biais froid en hiver
* Il y avait déjà un biais froid en hiver.
* Aussi bien avec Orchidee 2 que 11.
* Très fortes inversions en surface. Sans doute un pbolème de couche limite.
Présente même en moyenne mensuelle.
* L'été : rayonnement de surface beaucoup trop faible. Sans doute trop de nuages.
Axe 2 : variabilité des précipitations tropicales
-------------------------------------------------
- Beaucoup trop de pluies océaniques avec pas assez de variabilité
- evaluation
* daignostics à mettre dans les atlas.
* systématiser les diagnostics MJO.
* Réactiver les diagnostics en régimes et les potentiels de vitesse.
- Amélioration de la convection.
* Réactiver l'ancien mélange. Jean-Yves a activé et commis la thermo de la glace
avec l'ancien mélange. Mais bug soupçoné.
* Catherine a des développements à commettre sur l'épluchage de l'ascendance
adiabatique.
* Répartition pluc/plul : la disparition de plul peut avoir un lien
avec la disparition de la variabilité.
Idée : coupler la convection à fisrtilp. En exploration.
* Fermeture du schéma de convection : travailler sur le lien entre la densité
de cumulus et la densité
* Frédéric : effets de brises.
- tests de sensibilité.
* Tests des développements récents à faire.
* Eventuellement mettre des paramètres différents sur océan et continent.
Axe 3 : nuages / rayonnement
----------------------------
- Questions identifiées
* hautes latitudes trop froides avec des gros biais radiatifs.
Biais opposés sur le ciel clair et les nuages.
* Biais d'albedo sur la banquise. Albedo beaucoup trop faible et
beaucoup plus faible qu'avant.
* Trop de nuages bas l'été sur les continents
* Non prise en compte des précipitations glacées dans le rayonnement.
* Va regarder les effets d'angle aux hautes latitudes.
- A commencé à modifier les atlas pour ajouter des diagnostics
Travaille sur les graphiques latitude/altitude des nuages.
Va utiliser ce qu'a fait Lionel sur l'interpollation verticale.
- A fait des cycles saisonniers haut/moyens bas.
Axe 4 : climats polaire
-----------------------
- problème principal : températures des basses couches beaucoup trop froides.
Biais de 10K sur océans.
Coeff de mélange à tuner.
LW descendant trop faible en hiver.
- Biais autour de l'Antartique du forçage des nuages.
SW descendant trop fort malgré des nuages trop présents.
Axe 5 : Moyennes latitudes
--------------------------
- problème essentiel : position des jets et storms-track
Couplages avec les billans d'énergie
- Un moyen de séparer : faire de temps en temps des simulations guidées.
- Circulation méridienne.
- une énigme : la thermohaline s'effondre dans des simulations très
guidées en vent.
- Besoin de bilans d'eau par bassins et régions océaniques.
Axe 6 : Haute atmosphère
------------------------
- Marion travaille sur la connexion LMDZ/retrobus.
Veut ajouter des sorties spécifiques.
- Travail en prévision sur la lecture des forçages d'ozone.
Avant : lecture de fichiers de INCA+Reprobus
Maintenant : Ozone "communautaire".
Problème d'interpollation verticale/lissage.
Nouvelles interpolations horizontales.
Vertcialement : le but est de recaler les tropopauses
par rapport aux données. Nécessite un étirement de la grille
verticale.
François premier utilisateur.
- A faire : prise en charge des choix de .def et évaluation sur
les simulations longues. Notamment diagnostics QBO.
Nouveaux diagnostics
--------------------
* Diagnostics en régimes, y compris les taux de chauffage (3D)
* Potentiel de vitesse : existe dans ncl. A mettre en place dans
climaf. Demander à Jérome Servonnat ?
* dynzon : marche sans doute mais avec des sorties un peu bidon.
Francis doit faire le point.
* Précipitations par bassins.
* Rapatrier des fichiers jour à jour sur ciclad sur toute la période.
pluc, plul, topl, topl0 + u,v,T,q à 850hPa
Idées transverses à tester :
----------------------------
- Prise en compte des chutes de glace dans le transfert radiatif
- Prise en compte d'une fraction de nuages effective pour les angles
rasants.
- Contrôle des couches limites stables, pour les biais froids mais
également les plantages.
Plantages :
-----------
Laurent continue à tourner sans planter avec la version de l'été.
Abderrahmane a extrait un plantage qui selon Jean-Yves semble concerner
une couche limite avec une tr_s forte inversion.
svn upgrade (pour mémoire ...)
------------------------------
Besoin de faire un svn upgrade avant svn update avec des svn trop récent.
mettre un set -e avant svn upgrade
et set +e
Consolidation (Inchangé depuis les dernier poihl)
-------------
Séparer gcm.def en gcm.def et vert.def
Faire le point au prochain poihl sur ces différents aspects :
* Problème d'évaluation des SLP : Marie
* Problème de fichiers qui ne sortent pas normalement sur
ada (Ionela et Marie-Alice)
* Calendrier.
Mettre au carré. Faire le lien avec le forçage d'Olivier.
* Vérifier qu'il n'y a pas un problème avec les variables SAVE
(notamment l0) dans yamada4.
* A FAIRE : Test XIOS versus IOIPSL
* AFAIRE : Protection de la lecture des fichiers .data n'co.
Reprendre la routine getin_p des planetos.
Le infotrac n'est pas protégé mais il est en dehors du
domaine openMP.
FAIT : Ehouarn a mis la routine getin_p dans les sources.
* 1+1=2 avec les dernières versions de la physique (Laurent)
Tourner avec plusieurs physiques.
* Sensibilité au pas de temps de la fermeture stochastique
(récréation pour Jean-Yves).
* Traceurs et isotopes.
Pour le moment les isotopes ne marchent pas.
* Nettoyer le calcul de ALE et ALP dans physiq.F90 avec le problème
particulier de zw2.
* Faire que le adjust tourne en mpi-omp. Vérifier que ca tournerait avec 1 omp.
* trusting : monter en puissance sur la fréquence et le nombre des configuration
testées.
* Conservation de l'énergie avec l'eau.
* Choses bizares sur le filtrage de la vorticité.
* Question d'Olivier : rétro-compatibilité des start avec le nouvel albedo.
Laurent transmet la question à Olivier Boucher
* Problème des valeurs de la gravité.
* Sur les états initiaux :
- Généraliser l'utilisation des pondérations barycentriques
(plutôt que des interpollation de points de grille)
* grossismx : si on passe de 1. à 1.+epsilon, on décalle d'un quart
de point de grille.
Nouveaux développements
-----------------------
* Mélange convectif dans Emanuel (pour repartir à iflag_mix=1)
* Efficacité de précipitation.
* Splitting des poches
* Entrainement convectif
* Drags de surface sur océan (Frédéric, Traore)
* Pdfs verticales pour les nuages bas (à tester)
* Tube dentifrice (à tester)
* Richardson minimum en surface.
* Eau avec source méthane.
Optimisation
------------
En suspend
Formation LMDZ :
================
Essayer de passer les scripts d'install parallele en OpenMP seul.
Marie-Pierre va tester.
Discussion par axe :
===================
Axe 1 : Surfaces continentales
------------------------------
- Principal soucis : biais froid
- Test en cours
* Intégré nouvelle discrétisation verticale et code de Fuxing.
L'ancienne discrétisation n'est plus disponible.
En cours d'analyse.
* Nouveau module de neige du CNRM. Commis dans orchidee.
* Albedo : hors albedo de la neige. Travail d'optimisation de Philippe Peylin
non intégré.
* z0 : Fuxing et Nicolas (Vuichard) on travaillé dessus.
Réduit l'évaporation l'été.
- analyse du biais froid en hiver
* Il y avait déjà un biais froid en hiver.
* Aussi bien avec Orchidee 2 que 11.
* Très fortes inversions en surface. Sans doute un pbolème de couche limite.
Présente même en moyenne mensuelle.
* L'été : rayonnement de surface beaucoup trop faible. Sans doute trop de nuages.
Axe 2 : variabilité des précipitations tropicales
-------------------------------------------------
- Beaucoup trop de pluies océaniques avec pas assez de variabilité
- evaluation
* daignostics à mettre dans les atlas.
* systématiser les diagnostics MJO.
* Réactiver les diagnostics en régimes et les potentiels de vitesse.
- Amélioration de la convection.
* Réactiver l'ancien mélange. Jean-Yves a activé et commis la thermo de la glace
avec l'ancien mélange. Mais bug soupçoné.
* Catherine a des développements à commettre sur l'épluchage de l'ascendance
adiabatique.
* Répartition pluc/plul : la disparition de plul peut avoir un lien
avec la disparition de la variabilité.
Idée : coupler la convection à fisrtilp. En exploration.
* Fermeture du schéma de convection : travailler sur le lien entre la densité
de cumulus et la densité
* Frédéric : effets de brises.
- tests de sensibilité.
* Tests des développements récents à faire.
* Eventuellement mettre des paramètres différents sur océan et continent.
Axe 3 : nuages / rayonnement
----------------------------
- Questions identifiées
* hautes latitudes trop froides avec des gros biais radiatifs.
Biais opposés sur le ciel clair et les nuages.
* Biais d'albedo sur la banquise. Albedo beaucoup trop faible et
beaucoup plus faible qu'avant.
* Trop de nuages bas l'été sur les continents
* Non prise en compte des précipitations glacées dans le rayonnement.
* Va regarder les effets d'angle aux hautes latitudes.
- A commencé à modifier les atlas pour ajouter des diagnostics
Travaille sur les graphiques latitude/altitude des nuages.
Va utiliser ce qu'a fait Lionel sur l'interpollation verticale.
- A fait des cycles saisonniers haut/moyens bas.
Axe 4 : climats polaire
-----------------------
- problème principal : températures des basses couches beaucoup trop froides.
Biais de 10K sur océans.
Coeff de mélange à tuner.
LW descendant trop faible en hiver.
- Biais autour de l'Antartique du forçage des nuages.
SW descendant trop fort malgré des nuages trop présents.
Axe 5 : Moyennes latitudes
--------------------------
- problème essentiel : position des jets et storms-track
Couplages avec les billans d'énergie
- Un moyen de séparer : faire de temps en temps des simulations guidées.
- Circulation méridienne.
- une énigme : la thermohaline s'effondre dans des simulations très
guidées en vent.
- Besoin de bilans d'eau par bassins et régions océaniques.
Axe 6 : Haute atmosphère
------------------------
- Marion travaille sur la connexion LMDZ/retrobus.
Veut ajouter des sorties spécifiques.
- Travail en prévision sur la lecture des forçages d'ozone.
Avant : lecture de fichiers de INCA+Reprobus
Maintenant : Ozone "communautaire".
Problème d'interpollation verticale/lissage.
Nouvelles interpolations horizontales.
Vertcialement : le but est de recaler les tropopauses
par rapport aux données. Nécessite un étirement de la grille
verticale.
François premier utilisateur.
- A faire : prise en charge des choix de .def et évaluation sur
les simulations longues. Notamment diagnostics QBO.
Nouveaux diagnostics
--------------------
* Diagnostics en régimes, y compris les taux de chauffage (3D)
* Potentiel de vitesse : existe dans ncl. A mettre en place dans
climaf. Demander à Jérome Servonnat ?
* dynzon : marche sans doute mais avec des sorties un peu bidon.
Francis doit faire le point.
* Précipitations par bassins.
* Rapatrier des fichiers jour à jour sur ciclad sur toute la période.
pluc, plul, topl, topl0 + u,v,T,q à 850hPa
Idées transverses à tester :
----------------------------
- Prise en compte des chutes de glace dans le transfert radiatif
- Prise en compte d'une fraction de nuages effective pour les angles
rasants.
- Contrôle des couches limites stables, pour les biais froids mais
également les plantages.
Plantages :
-----------
Laurent continue à tourner sans planter avec la version de l'été.
Abderrahmane a extrait un plantage qui selon Jean-Yves semble concerner
une couche limite avec une tr_s forte inversion.
svn upgrade (pour mémoire ...)
------------------------------
Besoin de faire un svn upgrade avant svn update avec des svn trop récent.
mettre un set -e avant svn upgrade
et set +e
Consolidation (Inchangé depuis les dernier poihl)
-------------
Séparer gcm.def en gcm.def et vert.def
Faire le point au prochain poihl sur ces différents aspects :
* Problème d'évaluation des SLP : Marie
* Problème de fichiers qui ne sortent pas normalement sur
ada (Ionela et Marie-Alice)
* Calendrier.
Mettre au carré. Faire le lien avec le forçage d'Olivier.
* Vérifier qu'il n'y a pas un problème avec les variables SAVE
(notamment l0) dans yamada4.
* A FAIRE : Test XIOS versus IOIPSL
* AFAIRE : Protection de la lecture des fichiers .data n'co.
Reprendre la routine getin_p des planetos.
Le infotrac n'est pas protégé mais il est en dehors du
domaine openMP.
FAIT : Ehouarn a mis la routine getin_p dans les sources.
* 1+1=2 avec les dernières versions de la physique (Laurent)
Tourner avec plusieurs physiques.
* Sensibilité au pas de temps de la fermeture stochastique
(récréation pour Jean-Yves).
* Traceurs et isotopes.
Pour le moment les isotopes ne marchent pas.
* Nettoyer le calcul de ALE et ALP dans physiq.F90 avec le problème
particulier de zw2.
* Faire que le adjust tourne en mpi-omp. Vérifier que ca tournerait avec 1 omp.
* trusting : monter en puissance sur la fréquence et le nombre des configuration
testées.
* Conservation de l'énergie avec l'eau.
* Choses bizares sur le filtrage de la vorticité.
* Question d'Olivier : rétro-compatibilité des start avec le nouvel albedo.
Laurent transmet la question à Olivier Boucher
* Problème des valeurs de la gravité.
* Sur les états initiaux :
- Généraliser l'utilisation des pondérations barycentriques
(plutôt que des interpollation de points de grille)
* grossismx : si on passe de 1. à 1.+epsilon, on décalle d'un quart
de point de grille.
Nouveaux développements
-----------------------
* Mélange convectif dans Emanuel (pour repartir à iflag_mix=1)
* Efficacité de précipitation.
* Splitting des poches
* Entrainement convectif
* Drags de surface sur océan (Frédéric, Traore)
* Pdfs verticales pour les nuages bas (à tester)
* Tube dentifrice (à tester)
* Richardson minimum en surface.
* Eau avec source méthane.
Optimisation
------------
En suspend