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.