Aller au contenu. | Aller à la navigation

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

9/10/2017

Compte rendu du poihl du 9/10/2017

Durée :  1h40


ATTENTION : stopper les commissions sur la branche LMDZ5


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

Plantage/1+1=2/reproductibulité : RAS
Optimisation pour la 6.0.13 (cangement de arch et compilo 2017, curie)
Laurent a créé la branche LMDZ6 vendredi.
 - Va créer une branche 6013 à partir de ce tronc.
 - Va corriger creation_modipsl.sh
 - On bascule le nom LMDZ5 en LMDZ dans l'arborescence.
Relancer les simulations guidées et climatiques quand les
forçages seront définitifs.

Arrivé au bout de 20 ans en HR, avec un plantage résolu en diminuant
une fois les constantes de dissipation, avec un pas de temps physique
de 15 min et pas de temps dynamique avec le même CFL qu'en LR.

Contrôle des sources :
----------------------
Remplacer tous les call getin par getin_p en dehors de conf_phys.
Frédéric s'occupe des thermcell.

Préparation des configs :
-------------------------
Ionela s'occupe de configurations
1/ avec CRE éteints.
2/ aquaplanète
 Rm :  IF (read_climoz==-1) zslat=0. ! Imposing hemispheric symetry

Data request (responsable Laurent)
============================
S'y sont remis : Laurent, Ionela et Abderrahmane.
Passé de 178 à 160 à trouver (sur 340).

Question sur les niveaux de pression :
   Arnaud a fait des tests avec XIOS ou on line.
dynzon fonctionne et sort bien les fonctions de courant.
Yann est parti sur l'écriture de disagnostics équivalents sous
   XIOS. Depuis la dynamique ou la physique ?
   Etendre la discussion à Aymeric et Sébastien pour la planéto.
   Laurent envoie un message à Yann pour lui demander.
Ionela lance une discussion sur le calcul des termes de transport.

Articulation convection profonde / convection peu profonde
=============================================
Faire le point sur ce qui est vraiment dans le code
* Splitting des poches : équations imprésentable
* Entrainement convectif / nouvel épluchage
* Tube dentifrice (à tester)
* Dyamique de population de poches
* gestion des précipitations

Evolution physique/versions futures :
============================
* Eau avec source méthane
* sursaturation
* Augmentation de la résolution.
* grille plus groissère en bas pour la convection

Dynamico :
==========
Attention, le rotationel n'est pas passé à la physique.
Laurent s'est installé sur curie et a fait tourner  une
   première aquaplanète. A fait 10 ans avec la physiq 6012
   avec les sources de mai dernier.
   sauf rrtm=0 et iflag_ener_conser=-1 (car pas prévu).
Traore continue à avoir des problèmes pour démarrer.
L'aquaplanète "tout LMDZ" déposée sur Heat se plante.
Mise en place de config modipsl, sans doute au bocal de jeudi.

Tuning automatique :
====================
Une réunion vendredi².
On a réussi à faire tourner la vague 2 (grâce à une remarque
   de Ionela).
En discussion avec Toulouse pour échanger des codes "r".
Programmation d'un exposé itune en réunion climat.

Convergence 1D std / Dephy :
============================
Marie-Pierre a mis tous les cas Ayotte au format commun.
Le cas 24F plante. Plante sans doute avec certains physiq.def
Fleur a envoyé des fichiers de forçage pour AMMA, Gabls 4 et Fire.

Publications :
===========
- Stratocumulus (Arnaud Jam)
- Thermiques / stratocumulus / bords Est (Frédéric)
- Splitting (Jean-Yve)
- Couplage / Flux de surface / rafales (Traore)
- Tuning (Ionela/Frédéric)
- Articulation convection peu profonde / profonde ?

tickets trac :
=========
37 en cours.
Une dizaine à fermer.
Jean-Baptiste: Concernant le ticket #51 "Vents niveaux standard dans le BENCH" ,
j'ai regardé, ce n'est pas un bug, mais qqchose qu'il est bon de rappeler;
Le ticket disait:
------------------------------------
le histmth.nc du BENCH du install.sh sort des vents aux niveaux standard (par ex. u850) absurdes, de l'ordre de 10000 m/s
ce n'est pas le cas dans les simulations type AMIP d'Abderrahmane par exemple;
c'est peut-être propre au BENCH, j'ai un peu creusé mais je n'ai rien vu pour le moment
------------------------------------
et Jean-Yves explique la raison de cela dans un mail de la même époque (fin 2015) que je reporte ci-dessous; Jean-Yves y soulève aussi une question intéressante qui vaut le coup sans doute d'être rediscutée:

    Les sorties sur les niveaux standards (telles u100, u500, ...) peuvent présenter des valeurs parfaitement fantaisistes. Si je comprends bien, ceci arrive lorsque la simulation est interrompue après une durée qui n'est pas un multiple de la période de sortie sur le fichier considéré. Ce n'est pas forcément un bug, simplement une caractéristique dont il faut se souvenir.
    Ce qui semble plus proche d'un bug, c'est que ceci n'est vrai que pour les fichiers de sortie numéro 1 et 2. Dans calcul_STDlev.h, la sortie de valeurs moyennées n'est programmée que pour les fichiers 1 et 2. Pour tous les autres, ce sont les valeurs instantannées qui sont sorties. Ceci est confirmé par les commentaites qui annoncent bien que les fichiers 1 et 2 comportent des sorties de type "day" ou "mth".
    Tout ça ne semble pas compatible avec la flexibilité qui existe en principe maintenant sur les natures des diverses fichiers.
   Alors, bug ou pas ??
Ionela, concernant la flexibilite des sorties sur les niveaux standards (telles u100, u500, ...).
Actuellement elle existe lorsque l'on utilise les sorties par XIOS (ok_all_xml=y dans config.def), mais pas si l'on utilise les sorties via IOIPSL (ok_all_xml=n).
Ceci date de l'epoque de CMIP5 ou ces sorties etaient faites par un mécanisme différent des autres fichiers hist*.nc.
J'avais code en dur les frequences de sortie et les types de sorties (moyenne ou instantannee) selon ce qui etait demande par CMIP5.
On a des valeurs moyenees pour les fichiers no. 1 (histmthNMC) et 2 (histdayNMC) et des valeurs instantannees pour le fichier 3 (histhfNMC). C'est toujours le cas si on met ok_all_xml=n.
Bien sur je peux changer cela si besoin.
Par contre si l'on met ok_all_xml=y alors on peut changer les frequences de sorties et les types de sorties (inst ou ave) via les parametres freq_outNMC__1, freq_outNMC__2, freq_outNMC__3 et phys_out_filetypes__7, phys_out_filetypes__8 et phys_out_filetypes__9.
freq_outNMC__1 -> histmthNMC
freq_outNMC__2 -> histdayNMC
freq_outNMC__3 -> histhfNMC
phys_out_filetypes__7 -> histmthNMC
phys_out_filetypes__8 -> histdayNMC
phys_out_filetypes__9 -> histhfNMC

Les numeros sont un peu tordus et c'est du au fait que maintenant l'ecriture des fichiers se fait par la meme routine phys_output_write_mod.F90
Les 6 premiers fichiers de sorties sont les hist* habituel et les 3 suivants (7, 8 et 9) sont les NMC, puis il y a le histstrataer.

    phys_out_filenames(1) = 'histmth'
    phys_out_filenames(2) = 'histday'
    phys_out_filenames(3) = 'histhf6h'
    phys_out_filenames(4) = 'histhf3h'
    phys_out_filenames(5) = 'histhf3hm'
    phys_out_filenames(6) = 'histstn'
    phys_out_filenames(7) = 'histmthNMC'
    phys_out_filenames(8) = 'histdayNMC'
    phys_out_filenames(9) = 'histhfNMC'
    phys_out_filenames(10)= 'histstrataer'


Formation LMDZ :
==============
16 inscriptions.
Marie-Pierre retrouve les infos échangées après la formation pour
les faire circuler.

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

Marie-Pierre :
--------------
install 3002 OK.

Jean-Yves :
-----------
Test de la ALP conditionnel en cours par Abderrahmane.

Ehouarn :
---------
Changement des restart avec une commission récente. Mais convergence
numérique OK.
A intégré les nouvelles routines slab de Francis.
  (Vérifier que c'est bien dans la LMDZ6)

Ionela :
--------
A pris la 3002 pour travailler les configurations.
Est tombée sur des problèmes de grilles dans les sorties xml de Cosp.
En attente de corrections d'Abderrahmane.
Vérification des guidages.

Laurent :
---------
RAS


Jean-Baptiste :
---------------
RAS

Frédéric :
----------
RAS

Ajouter un commentaire

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

Enter the word