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