07/04/2014
Compte Rendu du Point Hebdomadaire LMDZ du 07/04/2014
CR POIHL 2014/04/07 ------------------- Durée : 1h30 Tour de table, anticyclonique (dans l'hémisphère Nord) Lluis ----- - premiers résultats simulations WRF physique LMDZ - étude d'une simulation dans la méditerranée, exemple de figure montrant le vent à 10m dans l'ancienne et la nouvelle physique (NPv3.0) simulation 40km résolution avec 39 niveaux verticaux (grille WRF) et sponge layer à 10 hPa; simulation 5 jours non guidée (enfin sauf aux bords); initialisation ERA-Interim; version WRF 3.3 - pour l'instant, run séquentiel; parallélisation à venir - comparaison outputs WRF et physique LMDZ en cours
- résultatshttp://www.lmd.jussieu.fr/~lflmd/WRF_LMDZ_medic950116/medicane_19950115_WRF-LMDZ.html
Jean-Yves --------- - terminé la correction pour conserver l'eau; mêmes résultats qu'avant en 1D mais cette fois eau conservée à 10^(-12) près (attention: il faut mettre dqimpl à 0 par défaut dans les thermiques, i.e. thermcell_main.F90; dqimpl était par défaut à -1 et non conservation dans ce cas) - problème cependant dans la convection avec niveau de condensation dans la couche d'alimentation -> division par 0 lorsque condensation au sol (brouillard) -> dernières versions marchent pas -> correction en cours - vérifié également la conservation des traceurs avec le béryllium - mise au propre en cours de la convection: cv3p_mixing.F90 est fait; suppression des if(cvflag_grav) dans cv3_routines.F90 en cours (à la place, on met grav à 10m/s2 si cvflag_grav est .F. et à RG si cvflag_grav est .T.); suppression à venir du cas scalaire (subroutine wake_scal) dans les poches froides pour ne laisser que le code vectoriel. Laurent ------- - mise en place en cours de simulations couplées en 144x142x59 pour IPSL-CM6-RC0 avec 128 processeurs (27 MPI x 4 OMP, et 20 procs MPI sur océans) - soucis pour l'instant avec les interpolations du coupleur OASIS - il faut essuyer les plâtres; SUPPORTONS LAURENT (au sens de supporter usité au XVème siècle, i.e. soutenons Laurent!!) >>> DISCUSSION SUR LE PARALLÈLE: il faudrait progresser ensemble sur les méthodes et conventions adoptées pour le parallèle -> journée tutoriel/debugging, session code sprint (à la campagne!) pour remise au propre de certaines parties du code? ex. de thèmes abordés: comment gérer OMP master, les broadcasts (envoi info d'un proc à tous les autres), barrières OMP et MPI Marie-Pierre ------------ - nouvelle version install.sh 26/03 (basée sur SVN 2000) - nouvelle version instal1d.sh 28/03 avec correction du pb répertoire phy1d_orig - participation à l'intercomparaison CIRC (http://circ.gsfc.nasa.gov/, Continual Intercomparison of Radiation Codes) avec modèle 1D ancien rayonnement et RRTM >>> PARENTHÈSE RRTM: problèmes évoquées par Sunghye sur le forum résolus avec Lionel; modèle tourne avec résultats identiques en MPI et OMP; tests 1D pour flux ciel clair en cours Frédéric (arrivé en cours de réunion) -------- e-contribution: - A fait les commissions de Jean-Baptiste pour le contrôle de la gamme de température de la phase mixte liquide/glace pour l'eau nuageuse. svn 2006 Tests de réglage en cours en aquaplanète sur ada. Devrait proposer un premier jeu de paramètres mardi. - A fait les commissions de Jean-Yves pour la conservation des traceurs. Mais la version commise ne tourne pas en version nouvelle physique. Il s'agit de la svn 2007 La version modipsl.20140404.tar qui contenait ces deux modifications a donc été retirée. A voir si on décommet la 2007 en attendant de trouver le bug réponse lors de la réunion sur ce dernier point: on met un flag par défaut à 0 le temps que Jean-Yves résolve le problème également évoqué: LE DEBUGGEUR GDB DE GFORTRAN EST TRÈS PRATIQUE! Marie-Alice ----------- - préparation de la formation libiGCM dans une salle magnifique à l'IDRIS qui pourrait servir plus souvent (30 postes avec environnement créé au cas par cas par Patrick Corde) Ehouarn ------- - avait des doutes sur la reproductibilité (relancer le même exécutable n fois) des simulations sur ADA; tests plus poussés montrent que tout va bien - par contre, les résultats continuent de différer entre dyn3dmem et dyn3dpar; en résumé: Avec dyn3dpar : seq == MPI =/= OpenMP == Mixte MPI/OpenMP Par contre pour dyn3dmem: seq =/= MPI =/= OpenMP =/= mixte MPI/OpenMP - travaille sur XiOS et sorties station >>> PARENTHÈSE SUR DYN3DMEM: recommander dyn3dmem mais attention au guidage qui ne fonctionne pas avec dyn3dmem en général et avec OMP (aussi bien avec dyn3dmem que dyn3dpar). seul dyn3dpar/MPI marche. >>> DISCUSSION SUR XIOS: doit-on pousser le XML? à l'heure actuelle, sortie de 300 variables dans 5/6 fichiers -> passage à XML pourrait simplifier le code, mais il faut trouver un juste milieu Ionela ------ - runs à différentes résolutions horizontales (allant de 96x95 à 280x280) ancienne physique 39 niveaux -> écart avec obs de la précip globale (en mm/jour) augmente avec la résolution - va redégrader les résultats haute résolution pour tout comparer en 96x95 (augmentation précip vient sans doute en partie de là) Lionel ------ - travail sur la fonction d'Exner pour simplifier le code (apparaît en plusieurs endroits à l'heure actuelle) - travail sur le passage de la vorticité de la dynamique à la physique pour amélioration par François Lott de certaines paramétrisations >>> DISCUSSION INTERFACE DYN/PHYS: comment passer au mieux davantage de variables de la dynamique à la physique (par exemple, vorticité, champs en moyenne zonale, tendances vent du dv...)? comment bien interfacer dynamique et physique (par un module?); dépend de la dynamique utilisée; réflexion à poursuivre; Jean-Baptiste ------------- - a ajouté les modifications dans LMDZ pour changer les tmin, tmax et exposant de la fraction de glace dans les .def (flag "iflag_t_glace") -> voir log version SVN 2006 - a cherché à tester les changements avec ok_newmicro=n mais nuage.F90 semble ne plus bien fonctionner -> est-ce-qu'on préserve nuage.F90? - va regarder le bilan radiatif dans le cas TWPICE et comparer aux observations dans l'idée de trouver des paramètres optimaux pour les nuages hauts (tmin, tmax, et exposant pour fraction glace) Abderrahmane (e-contribution) ------------ - Nouvelle simulation NPv4.10 avec Orchidee 11 couches a tournée : Les sorties sont sur gaya : /workgpfs/rech/lmd/rlmd093/IGCM_OUT/LMDZOR/TEST/clim/NPv4.10/ATM/Output/ - Simulation APv4.0 : physique AR4 avec la résolution 144x142x59 Sorties sur gaya : /workgpfs/rech/lmd/rlmd093/IGCM_OUT/LMDZOR/TEST/clim/APv4.0/ATM/Output Atlas : http://dodsp.idris.fr/rlmd093/LMDZOR/DEVT/clim/APv4.0/ATLAS/SE_1981_1983/ATM/ATM.html - Tests avec 40 procs MPI, 8 et 16 taches OpenMP : RealCpuTime : 8 taches 1351.69000 (s) = 23mn 16 taches 1456.86000 = 24 mn Ave 47 procs : 8 taches 1087.14000 = 18 mn - En cours : adapter la partie où sont lus les fichiers *.data avec OpenMP activer Cosp Lidia (e-contribution) ----- * simulations AMIP basse résolution 96x95x39 nouvelle physique [1958-2009] - simulation de contrôle - simulation avec guidage 3h, 30mn L'exécution de ces simulations est terminée. Le post-traitement manuel de ces simulations est fait. * simulations AMIP moyenne résolution 144x142x39 nouvelle physique - simulation de contrôle - simulation avec guidage 3h, 30mn L'exécution de ces simulations est faite pour 10 ans [1979-1988], et le post-traitement correspondant est fait. * simulations tests AMIP basse résolution nouvelle physique, avec guidage 9h et 12h, et analyses des résultats avec Jean-Louis Dufresne. Lancement d'une simulation basse résolution nouvelle physique [1958-2009] avec guidage 12h, à partir de ces tests. L'exécution de cette simulation est terminée, et le post-traitement est fait. * la nouvelle physique La nouvelle physique définie dans les simulations AMIP basse résolution est la nouvelle physique CM5B (NPv3.1) de CMIP5. La nouvelle physique définie dans les simulations AMIP moyenne résolution est la nouvelle physique NPv3.2, c'est-à-dire la NPv3.1 avec la stabilité numérique du modèle. Réunion utilisateurs ==================== Idée de faire la réunion la semaine du 30 juin (voir calendrier universitaire http://www.upmc.fr/modules/resources/download/default/d_formation/calendrier_universitaire.pdf) Marie-Pierre prend ça en charge. LMDZ6 ===== Prochaine réunion LMDZ6 le 28 avril 14h30 Réglage des .def dans les lancements de la NPv4 ------------------------------------- 144x142xL59 disgmin=0.125 # dans gcm.def par exemple. day_step=960 ou 720 iphysiq=5 Configuration physique NPv4 : ------------------------------ + Fermeture stochastique * iflag_trig=1 * tau=8000 * Ajouter lels modifications de la fermeture. * Vérifier la taille critique des thermiques pour le déclenchement. + Transition liquide/glace * Sources à commettre. + Stratocumulus d'Arnaud * iflag_thermals_ed=8 * iflag_thermals=18 Forums : -------- http://forum.lmd.jussieu.fr/forums-lmdz/lmdz6 On se force à discuter la mise au point de la version beta de LMDZ6 sur le forum. On fait vivre les 2 fils de discussion pour commencer : - Configurations de référence - RTTM Profiling --------- A faire. Se souvenir de faire des tests de stabilité avec un abrupt 4xCO2 Basculer au maximum sur dyn3dmem. Faut-il rerégler les paramètres de l'orographie ? Qui prend en charge XIOS ? Rayonnement : commettre les sources rapidement. Réactiver les tests thermodynamique de la glace / déclenchement stochastique Regarder le changement de densité des poches. Jean-Yves fait remarquer qu'on pourrait mettre un klev_max. Avoir 1+1=2 avec la nouvelle physique. lmdz-info --------- On relance. lmdz-en-bref ------------ Nouvelle testing. Journées utilisateurs lmdz-info pleins de choses à raconter! Ehouarn s'en charge Conservation avec la nouvelle physique : ---------------------------------------- 2013/04/07 : conservation de l'eau résolue 2013/04/13 : nouvelle routine pour la conservation de l'énergie basée sur les termes source de l'équation pronostique de la TKE. 2012/03/05 : travail en cours A regarder. Lien avec les problèmes identifiés en 1D. Problème identifié dans les sorties, en 1D. 1+1=2 : ------- 2014/01/06 : le problème est un peu différent. C'est day_step qui n'est pas un multiple de idissip. Que fait-on ? 2013/09/02 : probleme identifié pour la dissipation : la période d'appel à la dissipation effective idissip n'était pas automatiquement un diviseur de day_step Résoudre la partie poches et convection de la nouvelle phys. Ehouarn reprend la partie dissipation. OK avec ancienne physique à partir de svn1626