Aller au contenu. | Aller à la navigation

Outils personnels
Se connecter
Sections
Vous êtes ici : Accueil Développeurs CR de réunions POIHL 29/5/2017

29/5/2017

Compte rendu du poihl du 29/5/2017

Durée : 1h25

Plantages :
===========

* Vérifier que le save suprimé dans Mellor & Yamada permet de retrouver
    la reproductibilité.
    Retester avec la version corrigée les histoires de précipitation en
    bandes quand on active la diffusion de la TKE par la TKE


Consolidation :
===============

Reproductibilité :
------------------
Le save intempestif dans Mellor et Yamada expliquait bien la
non reproductibilité.
Attention : le modèle faisait sans doute un peu n'importe quoi
dans ce cas.
Mais ne concerne que l'activation de new_yamada4=y.

Conservation de l'énergie
-------------------------
Jean-Yves s'y remet : quelques corrections à commettre.
Travail côté orchidee.

Travail sur l'utilisation des termes de production de la TKE pour
la conservation de l'énergie.
Nouveaux flags iflag_ener_conser prochainement commis :
Soit ou prend juste la transformation depuis l'énérgie cinétique
de la colonne (u,v -> TKE) soit aussi depuis l'énergie statique (T -> TKE)
3 : Rien                      1.470
4 : vdf u,v                   0.08057
5 : vdf u,v,T                -0.02686
6 : 5 + the u,v              -0.1998
7 : 5 + the u,v,T             0.4193
8 : 5 + the&con u,v           -0.4111
9 : 5 + the&con&oro u,v       -0.6901
10: 5 + the&con&oro&lif u,v   -0.7470
(tests effectués sur 10 ans en 80x80x79)
L'option actuelle (=1) correspond a 7 (u,v,T pour vdf + the)
+ l'ajout de u,v pour la convection.
C'est a dire ici #7 + ( #8 - #6 )  = 0.4 + ( -0.4 - -0.2 ) = 0.2 W/m2
L'ajoute de T pour "the" ajoute un terme positif ce qui est abhérant
   pour de la convection (qui crée de la TKE normalement).
   Le problème est sans doute liée à la question des changement phase.
   Je propose donc d'abandonner ce terme qui a priori fonctionne à l'envers.
En n'activant pas ce terme, mais en activant toutes les tendances
de u,v, on arrive a -0.74 W/m2 (choix 10).
Donc on voit qu'on explique presque 1 W/m2 sur les 2,5 qu'on a dans
les pieds dans le couplé.
On voit aussi que l'effet total de conversion d'energie cinetique
grande échelle est de
#4 - #3 + #10 - #5 = -1.55 (vdf) - 0.55 (autres) = -2.1 W/m2
A ne pas négliger

1+1=2 :
-------

Contrôle des sources :
----------------------
Vérifier que les save sont bien en tread_private
Jean-Yves a écrit un script systmatique et envoyé les sorties.
   Le bug du save de yamada4 était identifié dedans.
   Il faut s'en saisir.

Ionela a commencé à vérifier aussi que les variables de physiq_mod
qui sont passées à des sous programmes sont bien en save au cas où
les sous-programmes ne soient pas appelés à tous les pas de temps.

Règle :
Faire attention sur les appels tous les N pas de temps à ce que les
variables de sorties soient en save.
Et impératif que toutes les paramétrisations soient appelées au
premier appel à la physique.

Reflechir à faire des allocate et pouvoir mettre des NaN en initialisation.


Optimisation :
==============
En cours (Laurent)


Réglages :
==========

Tuning :
--------
Jean-Baptiste regarde avec Jean les PDFs sous-maille sur la verticale.
    Vont sans doute faire une proposition un peu améliorée.

6010 : tuning de base du couplé V6010.2 LR. fait.
6010fallv : re réglée avec un fallv plus faible pour retrouver des nuages hauts.
    fait. Mais trop froide de 1K. Proposition de re-réglage (cld_lc 0.00012 - 0.0001)
    fallva,b,c... ?
6010blcv : avec les changements pour avoir plus de variabilité des pluies
    et les changements pour la couche limite stable.
    (en gros la 60103b des tests documentés plus hauts)
    En cours de réchauffement.
6010gust : juste les rafales sur le blcv.
6010blcvfallv : pour gagner un peu de nuages hauts.

Ionela a lancé aussi des tests de sensibilité aux tailles des cristaux.
    Tests croisants des changements des tailles min et max.

Travail sur les diagnostics de tailles de cristaux et de gouttes.
    Vérifier que ce qui sort de newmicro et qui est utilisé dans le simulateur,
    est vraiment utilisé dans LMDZ.
    En tous cas, les tests de Ionela montre qu'on contrôle le radiatif avec les
      rayons min et max des cristaux dans les .def.

Faire un test f_cdrag_oce=0.9 et 1 sur la 6010fallv.

Lancer les tests de Jean-Yves sur la même simulation.

Regarder la pluie et P-E sur la côte Est des Etats-Unis.
   Un biais positif pourrait contribuer à envoyer beaucoup d'eau
   douce dans l'Atlantique Nord et défaboriser la THC
   et favoriser la banquise.

Tests de Ionela en cours sur fallv plus faible et ratqs plus faible.
Aucun test ne change l'OLR ciel clair.

Tempéartures continentales :
----------------------------
RAS

Couches limites stables :
-------------------------
RAS

Convection / précipitations tropicales :
----------------------------------------
Travail en cours sur les analyses RCE avec les LES de Caroline.
Ionela est sur un re-réglage avec les valeurs de elcrit/tlcrit de la NPv5.5
Frédéric a constaté qu'en compensant une augmentation des nuages convectifs
    (liée au changement de elcrit/tlcrit) par un changement de epmax,
    on revenait à une variabilité faible.
Les rafales seront activés dedans.
En attente du 3eme réglage de Ionela pour lancer les simulations de
    sensibilité à iflag_adj_wk.


Diagnostics :
=============
Nouveaux diagnostics Gabls 4

Abderrahmane a pas mal travaillé sur les multi-atlas
   Les LMA est diags de Marine ne semblent plus tournés.
   Les métriques PCMDI ne sortent toujours pas automatiquement.
   C'est urgent ?
   Les LMA remarchent.

Multi-atlas différentiels à finaliser (Ionela).

Automatisation des diagnostics en régime devrait marcher.
   Va mettre par défaut sans les régimes sur les nuages.

A ajouté les nuages de lmdz dans les comparaisons à Calipso (quand
   le simulateur n'est pas activé).



Tour de table :
===============

Adriana :
---------
passe

Marie-Pierre :
--------------
A inclus des diagnostics de Jean-Baptiste dans le 1D.
+ des diags de Frédérique sur dice.

Jean-Yves :
-----------
Se bat avec la conservation de l'énergie dans la convection.
A compris que dans le schéma d'Emanuel, l'eau condensée est
   à la température de l'environnement (ce qui permet de calculer
   une fois pour toutes les Lv qui dépendent de la température).
A enlevé les poches, la glace et l'entrainement.
Mais ca ne conserve toujours pas.
Erreur trouvé sur l'ascendance adiabatique.
S'était trompé sur le mélange sur les couches d'alimentation en bas.
Modifications à commettre.

Ionela :
--------
En plus : a tracé les rayons des gouttes d'eau.
Modis voit des rayons plus gros dans les tropiques (15-20 microns)
au contraire de LMDZ (5 microns).
Remonter vers Olivier et Jean-Louis sur ces questions.

Ehouarn :
---------
RAS

Abderrahmane :
--------------
A commencé à vérifier les sorties Cosp demandées dans CMIP6.
   Simulations d'un an et contrôle des sorties.
   Problème d'unité en fraction ou % des sorties.

Frédéric :
----------
RDP

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.

Ajouter un commentaire

You can add a comment by filling out the form below. Plain text formatting. Comments are moderated.

Enter the word