22/6/2015

Compte Rendu du Point Hebdomadaire LMDZ du 22/6/2015

Durée : 1h25 minutes


Discussion des priorités
========================


Prochain poihl :
----------------
Parler de l'interface physique/dynamique

Plantages
---------

Forcés :
  * RRTM : a fait 4 ans après correction du LWdown.
    Prolongée de 5 ans mes a fini par planté avec hgardfou dans vdf latitude 87.
    Tester avec un pas de temps plus court (passer de 720 à 960 pas de temps/jour)
    A tester : aérosols et SW6bandes.
  * La 17d longue plantait. Relancée avec avec le random_trig_max=0.95.
    A fait 3 ans de plus. Mais plante dans la convection.
  * La simulation avec grille corrigée tourne. Le problème de démarrage venait
     du fait que la création de l'état initial tournait en parallèle doonc
     les états initiaux étaient bidons.
     Lancé en 17d. A tourné 2 ans. Plante dans les lsc dans les tropiques.
     Regarder s'il s'agit d'un plantage himalayen.
     Relancer avec un pas de temps plus court.
     Sortir un shade des points de grille pour comparer les reliefs avec et
       sans décalage.
     Tirer au clair ces histoires
  * Pas mal de plantages convectifs dans les régions froides.
      Seuil introduit sur les probabilités de déclenchement pour
      éviter des déclenchement rares sur des points où le déclenchement
      est très improbable.
  * Jean-Yves a codé un ajustement sec avant l'appel de la convection
      sur le profile interne aux poches.
      Mais résultats bizares en comparant l'ajustement sec d'Emanuel
      et celui de LMDZ.
  * Plantages couplés : non convectif.
      Forte vitesse verticale. Sans doute un plantage dynamique.
      Passage day_step : 720 -> 960
      Deux simulations se sont prolongées 1 an ou 2 mais plantes
        sur la convection ( 17Ca et 17ctrue )
      Plante en juille 1877. Plantage convectif. tmax_fonte_cv=288.15
        Tester de changer la température du max de fussion de la convection.
      Une troisieme se plante au bout de 4 ans (21bReg02). Est repartie
        en augmentant la dissipation pendant 1 mois.
        Se plante au bout de 19 ans dans la convection.
  * sortir systématiquement dtcon[k=@max]*86400 et proba_notrig sur les
        fichiers 1D des plantages.

Réglage
-------
  * Une simulation a tourné avec les bons paramètres pour les thermiques.
      Le net au sommet est 4 W/m2 plus grand que dans la 17c avec
      une modification éventuelle du SW CS de 2 W/m2.
  * A lancé les simulations avec un seuil sur la probabilité de déclenchement
      à seuil.
      Mise à disposition. Evaluer les déclenchement effectifs de la convection.
  * A intégré les modifications d'Olivier pour pouvoir travailler avec le
      SW 2bandes du nouveau rayonnement et les aérosols.
  * RRTM améliore les bilan radiatifs ciel clair.
      Mais le forçage radiatif SW des nuages est plus faible.
      Regarder comment les tailles de gouttes et cristaux sont gérés dans RRTM.

4 choses à regarder :
  * Les nuages hauts, et le choix des rayons des cristaux, et bilan radiatifs
     des hautes latitudes (Jean-Baptiste, Ionela, Marine)
      Regarder la position des jets.
      Frédéric a regardé les W500 dans les Reg01/Reg02
          Il semble qu'il y ait bien une différence sur W500 mais pas sur W intégré
      Regarder les tests de sensibilité à fallv de Ionela
      Se donner des diag : corrélations cldh / LW CRE
                           les simulateurs d'observations de rayons de cristaux
  * Evaluation des nuages bas (Jean-Louis, Abderrahmane)
      Se donner les moyens d'évaluer.
      Regarder les dernières simulations d'Abderrahmane
  * Dsitributions de précipitations (Marine et al.)
  * Flux air/mer (Frédéric)

Diagnostics à activer :
  * Réactiver dynzon.
  * Repasser à des simulations de 10 ans.
       Abderrahmane teste encore certaines choses.
 -> Commencer à travailler sur la position des jets.

Test :

Réunion LMDZ6
-------------
   Le 6 juillet
   Ionela doit présenter des choses sur cldh/crelw

Consolidation
-------------
  * Calendrier.
      Mettre au carré. Faire le lien avec le forçage d'Olivier.
  * Initialisation :
     Pouvoir initialiser en parallele.
     Commencer par le rendre faussement parallèle en ne lançant que sur
     le maître. Ehouarn demande à David si il peut regarder ça.
     EN COURS : David a commis une version nettoyée de ce0l, avec séparation
       dynamique physique (etat0_dyn, etat0_phys)
       Plantage : Tsoil et Snow disparus ou à 0
     A FAIRE le point sur le filtrage des orographies
     Se donner les moyens de vérifier que les fichiers initiaux et de limites
       sont corrects.
  * EN COURS : Jean-Yves sur le front de fusion.
  * 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.
  * oa2.nc : Olivier a commencé à faire du Orca1 en LMDZ-LR-NPv5.17c
  * 1+1=2 avec les dernières versions de la physique (Laurent)
  * 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 : CA MARCHE sur ada. LMDZORV6. LMDZ96x95x39. tourne 1 mois.
      Le mettre sur la page web service.
      Mettre en place sur la v5.17e
  * 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
  * Evaluer les flux ciel clair une fois qu'on aura les nouveaux albedo
  * Faire en sorte de pouvoir sortir les dynzon.
  * Problème des valeurs de la gravité.
  * Sur les états initiaux :
     - Orographie non filtrée pour le calcul du géopotentiel.
      Filtrés pour les paramétrisations des ondes de gravité.
      Bizarrement, le point du pôle part du champ filtré.
     - Généraliser l'utilisation des pondérations barycentriques
      (plutôt que des interpollation de points de grille)
     - Est-ce que oldice=.true. sert à quelque chose
     - questions rugoro et frac_liq
  Priorités : commissions des corrections de l'initialisation.
  * grossismx : si on passe de 1. à 1.+epsilon, on décalle d'un quart
       de point de grille.

Evaluation
----------
* Récupérer les diagnostics des rayons effectifs des cristaux de glace
     Marine. Sans doute sur les histday.
     Voir si les sorties mensuelles suffisent.
* Distribution des precips.
     Marine a récupéré des données TRMM à 3h et 0.25°.0
* QBO
* Problème de diagnostic de la pression au niveau de la mer.
   La solution a priori : il faut exterapoler la temperature sous
   la surface avec un gradient constant de 6K/km

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 cours.

Réalisation CMIP6
-----------------
Sorties, ...

CU IDRIS
--------
hybride (mpi)x(mpi-omp) : marche sur curie mais pas sur ada.
bugs du nouveau compilateur.

En septembre  ...
-----------------
Parler de la lablisation et de la formation.



Tour de table (divers)
======================
Nouvelle testing juste avant le travail sur l'initialisation. svn2291
Nouveau jeu de commissions en attente sur l'interface physique/dynamique
   * passage propre de tous les paramètres entre les deux.
   * disparition des dimensions.h dans la physique.
François Lott a besoin de la vorticité pour la physique.
   Ehouarn est OK pour qu'on passe la vorticité dans les arguments
   de physiq.F90
Ehoaurn : Le modèle plante toujours systématiquement avec les
   dernières versions du compilateur untel, en debug mpi/openmp.
   Essayer sans physique parce que plante avec dyn3dmem mais pas dyn3dpar.
Lionel : Sur ada, dyn3dmem debug plante avec un message qui dit
   qu'il ne trouve pas gdb.
Ionela : question de Frédérique Cheruy sur les .def.
   Regarder omp_chunck.
Ionela : n'a pas réussi à faire tourner l'ozone traceur.


Proposition de commission significative
=======================================
(nouvelle rubrique pour décrire à l'avance les grosses évolutions
 à discuter au poihl préalablement aux commisssions)
D'abord les corrections de l'initialisation.
En attente des modifications de l'interface physique/dynamique.
   dimensions.h va disparaitre.


A FAIRE
=======
- Labélisation LMDZ (Relancer Cyrile Crevoisier).

- Créer un module warning_print_mod.F90
    qui réunirait les impression des différentes routines tous
    les N pas de temps, ou N jours.

- Envoyer lmdzinfo


Gestion des versions
====================

Trusting :
----------
    Lionel veut bien travailler sur un script automatique qui teste plein de choses.
    Laurent en a déjà un.
    Intégrer les tests du 1D de Marie-Pierre dans le cecking local (avec Frédéric)
    Ce serait bien aussi de prenser à des tests Inca / guidage

Utilisation des branches svn :
------------------------------
Peut-être utiliser davantage des branches pour des développements importants
    (pour éviter d'être bloqué ou sous pression comme ça a été le cas récemment avec
    les modifications de Jean-Yves sur le splitting)





Labélisation LMDZ
=================
Dossier de labélisation pour l'été.
Relire les documents existants à quelques uns (Frédéric, Marie-Pierre, Laurent ...)
puis réfléchir ensemblel au bout par le quel prendre les choses.
Préparer un envoie du dossier.
Se souvenir que le responsble Insu a dit à la réunion des entrants que LMDZ était
labélisée (sous la responsabilité)
On relit à quelques uns pour se rendre compte de ce qu'il faudrait faire.
Est-ce que l'enquête a été faite après ?
Et la labélisation ?

Standardisation des noms de variables
=====================================
? Comment en discute-t-on ?

Tempêtes en points de grille
============================
Comment l'aborde-t-on ?

LMDZ1D/MAR
==========
Etienne Vignon va relancer le 1D / phymar

lmdz-info
=========
Reléctures puis finalisation au retour de Ionela.