5/11/2018

Compte rendu du poihl du 5/11/2018

===================
(après épilation de deux poihls)

Durée: 1h30

Prochaine formation LMDZ :
==========================
17-18-19 décembre
20 incsrits (en comptant les deux high-tune en attente de réponse)
C. Rio pas là
Jean-Louis et Jean-Baptiste présents

Consolidation LMDZ6 :
=====================

Problème de réstartabilité résolu: c'était un problème
de configuration pour les simulations CLIM. Pour résoudre le problème,
il suffit de mettre VEGET_UPDATE_at_start=n dans orchidee.card

Vérifier si les problèmes de reproductibilité sur ada avec
install_lmdz.sh sont résolus. (F&F) Qui semblaient être là même
sans orchidee.

Inchangé en-dessous

Problèmes de reproductibilité sur curie (Abderrahmane).
   Résolu avec la branche svn3353 plutôt que la trunk 3360.
   Abderrahmane s'en occupe avec Laurent.

Inchangé ci dessous :
* Ionela a des problèmes équivalents quand elle met les calculs
       directement dans pbl ou dans un sous programme.
* Plantages bizares à la ligne de changement de date dans les aquaplanetes
de Max.


Problème de diagnostics pour T2m :
----------------------------------
Ionela :
Si on a T2m et T2m_max on peut identifier les points bizares
    à partir d'une simulation corrigée qu'on utilise pour
    déterminer M=max( T2m_max-T2m ) en chaque point de grille.
Sur une simulation buggée, on corrige les points où T2m_max-T2m>M.

On remplace la valeur abhérante par une interpolation
temporelle sur 2, 4 ou 6 jours et apres j'ai applique egalement une
interpolation spatiale sur 4 points (nord-sud+est-ouest) pour combler
les points abherants restants pour le premier et/ou le dernier jour
de simulation.
Privilégier la correction au même point plutôt que les moyennes spatiales.

J'ai teste sur le run Amip et un run historical r1i1p1f1. Ca a l'air de bien
marcher.
Pour Amip j'ai utilise comme reference le run Amip corrige on-line sur
les points terrestres (simulation PatchTER).
Pour le run historical j'ai pris comme reference une simulation couplee
corrigee on-line en global du membre 2 (r2i1p1f1).

Le vent max a 10m ne semble pas affecte par des valeurs "abherantes", mais
je n'ai pas verifie que les situations abherantes correspondent a des vents
tres faibles. A faire.

L'humidite specifique ne semble pas affectee mais elle doit l'etre
puisqu'en tracant l'humidite relative minimale (hursmin) j'ai trouve
beaucoup (je dirais beaucoup plus si l'on compare aux valeurs abherantes de
tasmax) de jours avec des valeurs negatives.
A cause de ces valeurs negatives qui risquent d'etre presentes plusieurs pas
de temps dans la journee (a verifier quand meme), je dirais qu'il est
impossible de corriger huss et hurs.
   Diagnostiquer la fraction de points avec des rh2m_min négatifs.
   Regarder si on a sorti systématiqument rh2m_min quand on sortait t2m_max.
   Regarder l'erreur sur rh2m, rh2m_min et q2m (par comparaison des simulations
    jumelles)

Frédérique :
* Dans les cas qui plantent, la longueur de Monin Obukov et Richardson ne sont
pas du même signe. Peut-être lié à l'humidité.
   Si on prend tout sec on n'a plus de problème.
   Ou si on prend l'humidité de l'air patout (et pas beta qsat d'Orhcidee)
      on n'a plus de problème non plus.
   Le beta qsurf est très fort et beaucoup plus fort que l'humidité de l'air.
   Est-ce le cas ?
* Vérifier que l'utilisation de l'humidité à la surface beta qsat est bien
    cohérente et applicable à toutes les sous surface (Frédérique)
* Vérifier qu'on ne viole pas la 2ème principe (Jean-Yves)


Haute résolution + dynamico :
-----------------------------
Note sur la comparaison dynamico/lon-lat avec la physque 6A.

Contrôle Qualité / trunk / scripts
==================================
(inchangé)

install_lmdz.sh marche en local et sur ciclad avec ifort (ancien ~ 2013).
      Il faudrait tourner avec 2018 pour pouvoir initialiser à NaN
      Mais toute la chaine pour lancer est différente.
      En plan pour le moment parce qu'il faudra faire en même temps la
bascule du lancement parallèle par libigcm de poe vers mpirun.

Passer au nouvel orchidee par défaut.
      Fait en option.

La compilation avec XIOS n'est pas à jour pour Orchidee.

Version 6B
==========
Projet PRACE déposé pour porter entre autre un modèle
MR025 : grille MR dans l'atmosphère et résolution 1/4° océan.
Début de l'allocation avril 2019.
Définir d'ici là, et tuner, une version nouvelle de la nouvelle physique.
* Correction T2m
A compléter ...

Articulation convection profonde / convection peu profonde
==========================================================
(inchangé)
* Jean-Yves : toutes les équations sont simples et claires en cours de publi.
* Modèle pronostic de densité de poches : en test
* Splitting des poches : en test
* Entrainement convectif / nouvel épluchage (à tester)
* Tube dentifrice (à tester)

Evolution physique/versions futures :
=====================================
(inchangés ci-dessous)

Frédéric a fait des tests sur les stratocumulus, sur la convection
profonde (densité de poches interactive et épluchage modifié de Catherine)

* Jean-Yves, Jean-Philippe et Catherine pensent qu'il faut éliminer une
partie des précipitations dans les ascendances adiabatiques.
* Eau avec source méthane
* sursaturation
* grille plus groissère en bas pour la convection

Tuning automatique :
====================
(inchangés ci-dessous)
Marie-Pierre a un script pour interpoller LES et 1D sur la même grille
verticale. Important pour les EOFs. A mettre sur svn (Marie-Pierre et Ionela)

Convergence 1D std / Dephy :
============================


Publications :
==============
(inchangé)
Frédérique confirme qu'elle va écrire une publi sur le forçage à la surface
Biais de Bords Est (Frédéric, en attente de correction pat Catherine)
Stratocumulus (entre les mains d'Arnaud)
Gusts/z0 (Traore, figures en cours).
Nuages (demander à Jean-Baptiste où il en est)
Splitting (Nicolas et Jean-Yves)
Convection/poches (Catherine)
Reste description générale, tuning forcé et couplé, stratégie des
lancement ajustement des historiques.

Prospective LMDZ :
==================
(inchangé)
Organiser une prospective LMDZ en janvier/février.
Une grosse réunion ?
Des réunions régulières avec visio ?


tickets trac :
==============
#9     Modifications pour une dynamique commune (terre – planete)
#11     Menage sur les prints
#21     "day_ini" dans "create_etat0_limit" et "gcm"
#28     Initialisation problems in cloudth ?
#30     Pb valeurs nulles/non définies sur configs linux
#32     Concervation de la masse des traceurs dans la dynamique
#34     Calving faux en zoome
#44     getin_dump et name
#52     readchlorophyll en OpenMP?
#54     Erreur sur la vorticité en mode zoomé
#55     Constantes non uniformes dans le modèle
#56     Vérifier le calcul des nombres de jours de pluie
#61     Calendrier
#66     Sensibilité au pas de temps de la fermeture stochastique
#67     Traceurs et isotopes.
#70     Conservation de l'énergie avec l'eau
#71     Filtrage de la vorticité
#75     ppm3d non fonctionnel
#76     Sortie epmax
#82     Nouveau calcul de wind10m
#84     Use getin_p instead of getin
#85     Rayonnement dans le 1D
#86     Get rid of "stop" comands in phylmd
#88     routines cldrag
#89     Problèmes grille zoomée :
#90     Nbre d'appel à cdrag buggé?
#91     iflag_con et restartabilité du modèle
#92     Incohérence paramètres d'overlap des nuages
#93     Plantage à la ligne de changement de date
#94     Problème de reproductibilité
#95     Problème de reproductibilité

le #94 doit pouvoir être fermé


Tour de table :
===============
JYG
---
A regardé travailler Marie-Pierre et Ehouarn.
Le 1D fonctionne en 365 jours.
Problème de notations pour le papier splting.
    Coefs A, B et K pour le couplage avec la surface.
    Phi = K ( h1 - hs )
    K est arbitaire et redondant avec A et B.
Jean-Yves fait circuler un texte extrait de son texte splitting.

Abderrahmane :
--------------
RAS

Adriana :
---------
RAS

Ionela :
--------
Runs avec la branche IPSL-CM5A mais .def nouvelle physique.
On a aussi les soucis de T2m. Y compris avec Orchidee2

Marie-Pierre :
--------------
Voudrait commettre les modifications pour avoir le calendrier 365 en 1D.
A préparé un document scientifique pour le MTS (ex MEDDE).

Frédéric :
----------
A travaillé sur une revue de la paramétrisation de la convection avec Catherine.
A écrit deux rapports convergence : 1 sur dynamico
http://www.lmd.jussieu.fr/~lmdz/DOCUMENTS/TRAVAIL/D2.1b.pdf
 et 1 sur les simulations à 50km.
http://www.lmd.jussieu.fr/~lmdz/DOCUMENTS/TRAVAIL/D2.5a.pdf
A contribué à l'écriture du projet PRACE soumis.
http://www.lmd.jussieu.fr/~lmdz/DOCUMENTS/TRAVAIL/Demande_Prace.pdf