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