Aller au contenu. | Aller à la navigation

Outils personnels

Navigation

24/8/2015

Compte Rendu du Point Hebdomadaire LMDZ du 24/8/2015

Durée : 1h25 minutes


Tour de table (divers)
======================
Marie-Alice :
   Fait le point sur les (très faibles) évolutions du dossier calcul/CMIP6
   Fait la pub pour le rapport de Nalanda (fille de Balaji) sur la
   configuration interactive des grilles zoomées de LMDZ.

Lionel :
   Fait état d'un calcul erroné aux pôles dans ce0l (toutes les mailles
     autour du pôle ayant la même masse (faux en version zoomée).

Ehouarn :
   a pas mal avancé sur le nettoyage de l'interface physique/dynamique
   Règles de programmation :
       - Les deux seules interfaces entre la dynamique et la physique sont
       calfis et iniphysiq. Dans la physique, on ne doit faire de
       référence (include, use, appels de procédures) qu'au contenu de
       phylmd et phy_common. phy_common contient tout ce qui est commun
       aux physiques (phydev,phylmd,phymar).
       - un point délicat qui reste à élucider/trancher: l'instant
       (calendaire) à transmettre à la physique lors de l'appel à la
       dynamique; important en
       particulier pour le transfert radiatif.
   Attention à la sous-consommation sur Ada (seulement 52% des 800000
       heures Ada allouées pour le projet LMDZT)

Ionela, Marie-Pierre, Abderrahmane, Frédéric :
   RAS

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

Plantages
---------
 L'EXCELLENTE NOUVELLE DE L'ETE, confirmé par Laurent après le poihl :
   Laurent a fait tourné deux couplés cet été (versions v5.17g et v5.17h)
   sans aucun plantage modèle !!!!
 Le moment de lancer un 4xCO2 pour voir si on tient aussi dans ces
   conditions et se faire une première idée de la sensibilité
   du nouveau modèle (n'est-ce pas Marie-Alice ?).
 Reste t'il des plantages avec RRTM ?

Réglage
-------
  * Les réglages de dernières minutes avant l'été semblent donner des
      résultats pas trop mauvais (les 17g et 17h).
  * Il faudra essayer de mettre sur le papier la stratégie suivie pour
      ces derniers réglages. Frédéric s'est beaucoup basé sur les
      diagnistics de surface, en travaillant sur les structures, et en
      n'ajoutant que des modifications qui ne détérioraient pas notamment
      les biais de flux en contraste tropiques/bord-Est.
  * Et la circulation thermoaline est bien maintenue (20 svdp pour la
      nadw)
  * Les biais de bords Est semblent attennués (c'était une cible du
      réglage).
  * Parmi les défauts notables : un biais chaud circum antartique lié à un
      biais radiatif SW déjà identifié.
      Ce biais est apparu avec la nouvelle nouvelle physique.
      Il n'était là ni dans les versions anciennes physique, ni même
         dans la v5.32 (nouvelle physique sur la grille L79).
         Ce serait bien de lancer une v5.32 avec les derniers sources pour
         confirmer.
  * Autre soucis à considéré : la non conservation de l'énergie.
      2W/m2 entre le sommet et la surface en global.
      Sans doute arrivé en grande partie avec la thermo de la glace.

A 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.
  * Evaluation des nuages bas (Jean-Louis, Abderrahmane)
      Se donner les moyens d'évaluer.
      Regarder les dernières simulations d'Abderrahmane
  * Position des jets
  * Dsitributions de précipitations (Marine et al.)
  * Flux air/mer (Frédéric)
  * Documenter la sensibilité au pas de temps radiatif et tourner avec
      16 appels / jour

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

Consolidation
-------------
Faire le point au prochain poihl sur ces différents aspects :
  * 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.
  * 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.
     Se donner les moyens de vérifier que les fichiers initiaux
     et de limites sont corrects.
  * 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)
  * 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
  * Evaluer les flux ciel clair une fois qu'on aura les nouveaux albedo
  * 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.

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

En septembre  ...
-----------------
Parler de la labélisation et de la formation.



Proposition de commission significative
=======================================
RAS

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)

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

Tempêtes en points de grille
============================
C'est résolu, non ?

LMDZ1D/MAR
==========
Etienne Vignon a repris le 1D / phymar

lmdz-info
=========
A imprimer.

 

Ajouter un commentaire

Vous pouvez ajouter un commentaire en complétant le formulaire ci-dessous. Le format doit être plain text. Les commentaires sont modérés.

Enter the word