Aller au contenu. | Aller à la navigation

Outils personnels
Se connecter
Sections
Vous êtes ici : Accueil Développeurs CR de réunions POIHL Compte-rendus PoiHL 2015 29/6/2015

29/6/2015

Compte Rendu du Point Hebdomadaire LMDZ du 29/6/2015
Durée : 1h25 minutes


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


Prochain poihl :
----------------
Parler de l'interface physique/dynamique

Plantages
---------

Forcés :
  * Olivier a tourné avec la NPv5.17c et Orca 1 sans planter.
  * Runs couplés de Laurent
      Tout relancé avec random_notrig_max=0.95, tmax_fonte=288.15, les
      modifs d'Arnaud
      La v5.17true a fait 17 ans de plus et la 21breg02 a fait 21 ans de
      plus. Mais la 17true a planté dans la convection.
      Les fichiers sortis par Orchidee sont vérolés. Arrivé au moment où
      on a changé le pas de temps dynamique.
      Laurent va relancer un couplé propre avec la bonne grille, l'ozone,
      les aérosols.
  * Jean-Yves a commis les modifications sur l'ajustement sec dans la
      convection.
  * RRTM : a fait 4 ans après correction du LWdown.
       Prolongée de 5 ans mes a fini par planter avec hgardfou dans vdf
       latitude 87.
     Tester avec un pas de temps plus court (passer de 720 à 960 pas de
        temps/jour)
     Replante au même endroit.
     Relancée avec les mofications de l'ajustement de la convection :
     ok_adj_wk=y)
     A retourné 14 ans. (3 plantages de suite à latitude 88, longitude 73)
  * 17d longue : plante au même point avec les paramètres de la 17e
     Relancée avec les modifs de Jean-Yves. A fait 2 ans. Continue à
     tourner.
  * sortir systématiquement dtcon[k=@max]*86400 et proba_notrig sur les
        fichiers 1D des plantages.

Réglage
-------
  * Nouvelle série de simulations avec les modifications
       Ionela lance une série de tests de sensibilité.
       La 17e avait encore un problème de grille
      Une 17e01 arrive.
  * Reglage des nuages hauts
      Jean-Baptiste a regardé ce qui se faisait dans RRTM.
        Pour les nuages hauts : newmicro est appelé pour les rayons.
        Les fonctions qui relient contenu en eau et opacité sont dans RRTM
        Mais la taille des cristaux est pilotée par newmicro.
        Mais des commentaires laissent planer des doutes (notamment comme
          un facteur 2 sur rayon et et diametre).
      Ionela a regardé le rapport entre LWCRE et nuages hauts
          dans 21breg01/02 et les obs.
          La pente de la relation CRELW / nuages hauts Calipso
                 est trop faible dans LMDZ ancien rayonnement
                 augmente en diminuant la taille des cristaux
                 mais est déjà trop forte avec RRTM.
      Tirer au clair de savoir qui a changé la PDF de W500/Wmoyen
          Regarder sur les simulations de sensibilité tou
  * Pas de temps pour l'appel au rayonnement.
      Entre 16 et 32 appels (3 ou 1h30) on a le même temps de calcul.
  * dynzon ne fonctionne pas en parallele.
      Apparamment ca a marché sur la 4.12 en MIP/OpenMP à 59 niveaux.
      Passer à XIOS
  * 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
   Frédéric sur les flux air-mer

Consolidation
-------------
  * 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.
  * 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 : 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
  * 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.
  * Laurent regarde le problème de sensibilité au nombre de process openMP sur ADA.
       Quand on repasse à une physique ancienne ou 3.2, problème.

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

En septembre  ...
-----------------
Parler de la lablisation et de la formation.



Tour de table (divers)
======================
Infos calcul :
   2M heures sur ada. Nouveaux logins.
   10M heures sur curie.
Marie-Pierre a mis bomex dans le 1D. Répare Cindy
Refaire une nouvelle testing.
   Josefine a fait un test en mode debug avec orhcidee. Ca plante.
   Correction faite.
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
Ehouarn : 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.

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
=========
Relecture puis finalisation au retour de Ionela.

 

Ajouter un commentaire

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

Enter the word