wiki:Pirates2022

Version 67 (modified by cetlod, 19 months ago) (diff)

--

Informal working group on the development of IPSL-CM7

Thursday 29th September

Présents: Olivier B, Brady, Julie, Laurent F, Christian, Juliette, Myriam, Thomas D, Patricia C, Arnaud C, Guillaume, G, Claire, Ionela, Masa, Nicolas V, Philippe P

point info VLR == Myriam

  • Point de depart : IPSL-CM6.2-LR puis 6.5-LR degradées en resolution (ORCA2, 96x95 dans l'atmosphere), debuggage NEMOv4 ORCA2 (mers fermées), Correction biais vents moyen (position jets) - approche F, Cheruy.
  • Par ailleurs modifs namelist NEMO (a la casimir - Namelist2 cf ci dessous - nn_etau passé à 0) -> sans modifs vents, MOC se maintient (mean -2Sv par rapport à ORCA1), mais avec variabilité différente (spikes).
  • D'où viennent ces pics?
    • Variabilité temporelle couche de mélange hivernale dans VLR correcte, variabilité interannuelle.
    • Variabilité spatiale: perte de la convection Labrador et en grand partie mers Nordiques. Convection essentiellement dans Irminger et sud Islande, manque de filtre BF.
    • Banquise correcte dans HN (mais quid exactement dans Labrador? Probablement trop de glace? il faudra regarder cartes spatiales) mais absente dans HS.
  • NB: il y a aussi pas mal de changement dans namelist glace (Clément, notamment amax) pour coller à ORCA2 + tunings cdlc aussi.
  • Performances:
    • sur Irene ROME-SKL : 50-60 ans par jour
    • Fréquence de couplage passée à 9600. Permet de gagner en performances. Pas de temps de l'océan baissé à 4800. Glace couplée tous les 3 pas de temps de l'océan.
  • Actions:
    • faire tourner et faire circuler CESMEP pour poursuivre analyse
    • integrer une VLR dans atlas croisé des tests en cours.
  • Remarque: Différences potentielles avec la version Namelist2 (voir ci dessous) (en plus de pas de temps de couplage etc)
    • coefficients de diffusion ajoutés en dur par Casimir
    • cartes de diffusion de tmx différents

point sur simulations (couplé) en cours avec NEMO v4

  • intermonitoring des simus tuning :

​https://thredds-su.ipsl.fr/thredds/fileServer/ipsl_thredds/cetlod/interMonitoring/CM65v42work_IDRIS_InterMonit/index.html

  • CM62-LR-pi-JZ-01 = proche de 6ALR, ie NEMO 36 - 250 yr
  • CM65v420-LR-SKL-pi-05 = ref Christian avec NEMO v4 Namelist par défaut - 500 yr
  • Namelist 1 = new ref NEMO v4 - 200 yr
  • Namelist 2 = test etau = 1NEMO v4 - 100 yr
  • CM65v420-LR-CdL-pi-01 = modifs-casimir mixing - 400 yr
  • tableau namelists NEMO :

​https://docs.google.com/spreadsheets/d/1HBqlJ7Wfm9Co06qATAL9gEZ9M9RF1BgIJ7Z-Io_PXII/edit?usp=sharing

  • AMOC:
    • Namelist 2 : AMOC s'effondre. très cohérent avec tests de Myriam
    • Namelist 1 et CdL: acceptable, mais ne fait pas beaucoup mieux que 6.5 ie NEMOv4 par défaut. En AMOC. Quid autres diagnostiques? Blue Spot Atl N, biais chauds SO, : A priori, à la fois Namelist2 et 6.5 sont très proches de 6ALR.
  • ATTENTION: Gros changement NEMOv4 = passage à SI3! Il faut donc regarder et tuner en détail!

quelles prochaines simus ?

  • En machine actuellement: CdL (modifs Casimir). AMOC Se maintient voir augmente.
  • CM65v420-LR-SKL-pi-0: 500 ans, pas très itneressante par rapport à CdL.

tuning glace de mer ?

  • quels params à ajuster à la main ? (Martin ? Clement ?)
    • Voir les modifs à la main de Myriam?
  • est-ce qu'on lance le tuning automatique en couplé des params de glace du modèle couplé ? (Brady ? Laurent ? autre?)
    • Article soumis Martin et al tuning SI3 donne un premier range de paramètres (pour ORCA2 forcé) : ​https://gmd.copernicus.org/preprints/gmd-2022-170/
    • Quelles métriques? (à coder en XIOS)
    • L'idée serait HM en couplé avec des param de glace, en séparant les échelles de temps et en supposant qu'on peut pour le moment ignorer la BF

OB: on n'a pas de garantie que le tuning obtenu sera optimal en full couplé
JD et MK: oui en effet, mais cldlc doit etre tuné avec la glace de mer car fort couplage entre les 2.

protocole pdControl ?

  • est-ce qu'on modifie le protocole pour intégrer carte d'anomalies d'albedo plutot que anomalie globale (pour reproduire le heat uptake oceanique) ?
    • travail d'analyse preliminaire pour quantifier gains vs challenges ? Juliette et Brady
    • est-ce qu'on lance maintenant des simus pdControl (pour tuning glace de mer...) ? -> Oui il faudra etre en pd pour ca glace de mer

tuning ORCHIDEE ?

  • reflexion en cours dans le groupe orchidee sur la version du modèle à utiliser -> reunion le 14/10
  • travail actuel avec ORCHIDEE2. Ca n'est pas ce que le groupe orchidee préconise pour CMIP7, il faudra passer à ORCHIDEE4, meme si cette version n'est pas completement prete pour le couplé
  • ORCHIDEE2-4: gros saut technique, qui va a priori impliquer des changements en terme de condition au limite pour l'atmopshère, mais à voir.
  • Rappel contexte OPtimESM, projet européen début janvier qui peut faire jonction CMIP6-CMIP7 avec version modele intermediaire

-> inclure orchideev3 en cours de developpement? Ca permettrait d'experimenter et de valider. Mais ca ne doit pas décourager et freiner les dvts sur la v4.

  • quel(le) referent(e) pour les pirates ?

ouverture des cavités ? JULIE avait prévu de montrer qq slides de ce qui a été fait en forcé - manque de temps

retours sur reunion groupe plateforme: utilisation cylc vs libigcm ?

  • Laurent, Nicolas, Romain et Lola: explorent Cylc (ordonnanceur de tâches developpé par services météo NZ et Hadley Center) notamment sur SPIRIT
  • en attente d'une reponse des centres de calcul pour verifier si Cylc est compatible avec leurs contraintes de securité (reponse attendue en novembre pour TGCC). En effet, Cylc implique un démon python (assez classique pour un gestionnaire de workflow)-> négo en cours avec centre de calculs pour savoir si Cylc pourrait tourner chez eux. Si la réponse est non: fin des tests Cylc!
  • reste beaucoup de travail pour ordonnancer les taches de libIGCM par Cylc (travail dans Cylc + libIGCM) (NB pour les nuls: l'idée n'est pas de remplacer LibIGCM par Cylc mais d'utoliser Cylc pour mieux ordonnancer les taches de libIGCM)
    • qui prendrait en charge le travail ? embauche CDD à prévoir ?
    • Avantages attendues de Cylc: performances bof, plus propre oui, tuning & HM probablement

Organisation des reunions pirates

Prochaine reunion Couple DYNAMICO VENDREDI 14 OCTOBRE!

Autres ?

  • Thomas a des mauvaises nouvelles du NCAR.
    • Démission de J.F. Lamarque
    • Discussion autour de la ultra-haute résolution.

Thursday 15th September

Réunion pirates 15 septembre 2022

Présents: Yann, Claire, Myriam, Julie, Juliette, Olivier B, Arnaud, Brady, Christian, Claire, Laurent F, Martin, Olivier M, Patricia, Thibault, Guillaume G,Thomas D

DYNAMICO

Rappel: ​https://forge.ipsl.jussieu.fr/igcmg/wiki/IPSLCM6/IPSL-CM6A-ATM-ICO#Validation

  • Comparaison de
    • CM61-LR-pi-03 : piControl CMIP6, réalisé avec IPSLCM6.1
      • TGCC : (/ccc/store/cont003) gencmip6/p86maf/IGCM_OUT/IPSLCM6/PROD/piControl/CM61-LR-pi-03
    • VALID-CM622-LR.01 : 50 ans, piControl IPSLCM6.2, version à jour LMDZ (physique v6.2)
      • TGCC : gencmip6/p86caub/IGCM_OUT/IPSLCM6/DEVT/piControl/VALID-CM622-LR.01
    • TEST-CM72-SIMPLE-ROUTING : 50 ans, piControl DYNAMICO-LMDZ-ORCHIDEE-NEMO 3.6, version à jour LMDZ (physique v6.2), routage simplifié, valeur cld* par défaut 0,00065
      • TGCC : gencmip6/p86caub/IGCM_OUT/IPSLCM7/DEVT/piControl/TEST-CM72-SIMPLE-ROUTING
    • TEST-CM72-SIMPLE-ROUTING.06 : valeur cld* corrigée 0,00080 pour recalage température globale
      • TGCC : gencmip6/p86caub/IGCM_OUT/IPSLCM7/DEVT/piControl/TEST-CM72-SIMPLE-ROUTING.06
  • Intermonitoring : ​​https://thredds-su.ipsl.fr/thredds/fileServer/tgcc_thredds/work/p86caub/INTERMONITORING/intermonit_VALIDATION_DYNAMICO_COUPLE/index.html
    • SSH diminue lineairement, -30 cm en 50 ans -> pb bilan flux d'eau douce dans l'ocean à regler ? pas coherent avec excedent de friver... et pas de signal compensatoire en evap
      • Ressemble bug masques Myriam VLR: Olivier Marti Arnaud et Yann verifient: ce probleme n'etait pas la dans les versions précédentes
      • Autre hypothèse: on stocke de l'eau sur les calottes. Script Olivier Boucher, Christian sait faire, en refaisant tourner 20 ans pour ajouter tous les diagnostiques. Action à mettre en oeuvre dans un 2e temps.
    • wfo (E-P-R sous la glace) signe différent par rapport à IPSLCM6ALR
      • Cohérent ? avec friver bcp + fort dans DYNAMICO. biais diminué au moment du retuning, qui a diminué biais precip_land -> est-ce que c'est "grave" ce biais friver pour NEMO ?
      • Difficile de raisonner ainsi, friver pas indépendant de evap, pb ssh prioritaire
    • SSS acceptable dans les tropiques et eq. Moyennes latitudes: trop salé pdt 35 ans, meme si stabilisation puis diminution sur la fin. AMOC semble suivre (ou créer ?) ces 2 échelles de temps.
    • t2m_global : retuning recale bien. Mais sst 50S-50N devient trop froide. Biais chauds hautes latitudes fortement corrigé dans HN mais très peu dans HS
    • glace de mer HN et HS été et hiver insuffisante dans les 2 DYNAMICO, mieux dans le plus froid (TEST-CM72-SIMPLE-ROUTING.06) au bout de 50 ans.
    • AMOC s'effondre dans les 2 DYNAMICO (un peu moins dans la froide???), et formation NADW avec -> anomalies en profondeur T,S differentes de 6ALR
      • AMOC a besoin de plus de temps? Mais la disparistion de la formation d'eaxu dense est de toute facon mauvais signe. Lié aux problemes de stress de vent? A voir quand premiers pb corrigés
  • Atlas CESMEP 50 ans (reférence piControl CMIP6) : ​​https://thredds-su.ipsl.fr/thredds/fileServer/tgcc_thredds/work/p86caub/C-ESM-EP/couple_dynamico_ref_piControl_p86caub/C-ESM-EP_couple_dynamico_ref_piControl.html
    • PB Atlas ocean: cartes SST et SSS non exploitables
      • [Arnaud]: ca marche quand la référence est OBS plutot que un autre run
    • Cartes ATM montrent biais chauds au sud de 55°S et au nord de 60°N, en hiver.
      • [FH] Explications possibles:
        • 1) un problème propre à Dynamico, mais je n'y crois pas trop.
        • 2) une dé-compensation d'erreur structurelle par rapport à la version à résolution équivalente de LMDZlonlat (ou LMDZlola, c'est plus mignon).
        • 3) entre les deux : les deux dynamiques ont leurs problèmes et le tuning de la physique est capable de compenser ça en partie. C'était le cas par exemple entre les version 5A et 5B qui montraient des biais très différents avec le même coeur dynamique.
      • [Martin]: tout au début de IPSLCM6A-LR: vortex polaire était trop mou, associé à une différence tropique - pole trop faible. il y avait eu des ajustements au dessus de la calotte. Il y avait eu des ajustements d'albedo (Masa? JB Madeleine?). Est ce que ces ajustements ont été pris en compte ici? A priori oui, car même physique, et mêmes .def. -> Action: vérifier l'historique de ces modifications. Qui??
      • [Myriam] Ces biais viennent ils de la dynamique plutot? Origine plutot equatoriale (hypothèse )
      • [Yann] Autres différences fondamentales Dynamico vs LMDZlola:
        • Dynamico a une moins bonne résolution aux poles pour nombre de points équivalents (mis meilleur résolution à l'equateur)
        • Dissipation. : est on plus ou moins dissipé que dans LMDZlola?
          • Remarque [FH]: C'est curieux que dynamico amplifie le biais de position des jets par rapport à la dynamique longitude-latitude (mais on l'avait sans doute déjà vu). Je pensais qu'on dissipait moins le vent. Or cette dissipation tend généralement à ramener l'atmosphère vers une rotation solide, donc à ramener de la quantité de mouvement vers l'équateur. Est-ce que ca peut être une histoire de dissipation numérique ?
          • [AC] Mêmes valeurs de dissipation
          • [YM] Sébastien Fromang avait regardé les effets de ce paramètre. Retours? Tests en aquaplanète seulement? -> cf slides ajoutés en attaché sur cette page

* Action: Définir une série de test avec moins de dissipation? (Thomas, Yann, Arnaud, ?) Commencer par du forcé.

  • Gestion des ondes orographiques sous mailles -> Consulter Guillame Gastineau? Différence = algorithme de séparation résolu vs sous maille. Cette séparation a du être changée pour raisons techniques et en particulier faites en offline-préprocessing. Cette lecture pré-processing a été implémentée dans LMDZlola (pas dans version standard - Laurent F). David Cugnet avait vérifé que les 2 méthodes donnaient les mêmes résultats.
    • Orographie: cf anomalies de vent méridien dans le Pacifique qui semblent améliorer.

* En forcé?

avec

  • ICOLOR-SROUT.02 : 10 ans expérience clim, DYNAMICO-LMDZ-ORCHIDEE
    • TGCC : gencmip6/p86caub/IGCM_OUT/LMDZOR/DEVT/clim/LOR-SROUT.02
  • LOR-SROUT.03 : 10 ans expérience clim, LMDZ-ORCHIDEE
    • TGCC : gencmip6/p86caub/IGCM_OUT/ICOLMDZOR/DEVT/test/ICOLOR-SROUT.03
  • 1) sur les biais en température (mais avec des SSTs imposées ce n'est pas évident)
    • Oui clairement sur continent, même ordre de grandeur, max en hiver.
  • 2) sur les flux à la surface (comme c'est le cas sur les problèmes de bord Est)
    • On voit aussi une diminution du flux latent dans le circumpolaire, compté positivement pour une évaporation positive. -> réchauffera en couplé. Ce flux, légèrement positif sur l'océan circum antartique, devient encore plus faible (différence négative) avec Dynamico.
    • [FH] pense que c'est ca la source.
    • [JJ] Raisonnement uniquement HS, quid HN?
  • 3) sur la structure méridienne moyenne (latitude-pression) de la température et de la pression.
    • Déplacement du jet vers l'équateur en surface dans l'HS, visible sur les atlas de différence des tensions de vent zonal. diminution de la tension de vent zonal (d'ouest) dans ces latitudes peut expliquer moindre évaporation
    • Tension de vent méridienne vers le nord intensifiée dans les bassins nord, jet décalé vers le nord, ce qui devrait plutot améliorer l'océan (ventilation etc). Grosse anomalie + de SLP dans Pac N hiver. Mais ce changement corrige un fort biais de CM6ALR, qui avait basse pression trop à l'ouest et HP trop à l'est, le tout associé à fort biais chaud. Cebiais dynamique semble corrigé, peut etre réampliifié par autre effet?
  • Installation des biais?
    • est-ce qu'on voit ce biais se mettre en place sur une année ou quelques unes dans le couplé en partant de Levitus (ou autre "observation").
    • Tuning automatique / history matching: quel effet sur ces configs, quelles efficacité?
  • Analyse: (ceci est un échange de mail OB FH remis ici pour mémoire) [OB] : à AG de l'European Meteorological Society il y avait une session sur les vagues de chaleur: approche lagrangienne sur les masses d'air qui arrivent en un point donné, puis ils séparent l'intégrale du DT/Dt en plusieurs composantes dues à l'advection, aux processus adiabatiques et aux processus diabatiques. Ca ne donne pas la réalité, mais je me demandais si un truc similaire ne pourrait pas peut-être expliquer pourquoi deux modèles développent des biais différents à un endroit donné et pas ailleurs. Par contre ce n'est pas forcément simple à mettre en œuvre -> [FH] propose de coder ca en rétro, discussion en cours...

Other topics :

* NEMO en climat froid - cf pres Olivier Marti

IPSLCM6 plante systématiquement en glaciaire Tests NEMO4.2 IPSLCM6.5 Tests froids en CO2 seul, sans les calottes, namelist NEMO4.2 fournie par Guillaume Gachon -> ça tient. épaisseur de glace plafonne à 20m dans Arctique. Martin et Clément auraient une petite idée des modifs faites dans NEMO4.2 qui pourraient expliquer cette limitation qui évite les cathédrales de glace. 4xCO2 tient aussi Tests à faire en Dynamico régulièrement , quand on sera passé à NEMO4. Intégration de NEMO4 apres les premières validations/verifications (voir plus haut)

* futur de CMIP7 - ou du moins ce qu'on en sait... OB: ca commence à s'organiser.

  • Thibault volontaire pour prendre en charge forçages.
  • C'est à nous de nous orgaiser pour définir ce que l'on veut.
  • A discuter à la retraite?
  • OB sera au NCAR en novembre pour discuter avec JF Lamarque
  • Discussions avant l'été au moment de workshop WCRP. Vers des ensemles moins opportunistes?
  • A l'IPSL, calendrier CMIP7 risque de se téléscoper avec changement calculateur.
  • Prochain calculateur TGCC attendu 2024-2025, va demander du temps d'adaptation. Cela va impliquer des choix par rapport à CMIP7.
  • Mandat des WG IPCC actuels se termine au printemps 2023. Renouvellement du bureau à prévoir. Peut on attendre de réels changements? A priori pas trop, inertie.
  • Reflexion sur nb de MIP, nb de simulations etc à avoir.

* infos sur retraite des 17-18 novembre

  • date 17/18 novembre maintenue - à confirmer dans qqs jours après vérifcations de présence de certains membres. domaine St Paul
  • toutes suggestions d'atelier bienvenues !
  • Objectifs:
    • Remobiliser un peu plus large que le groupe pirates (ca 80 personnes)
    • Faire une réunon utile pour nos développements
  • Il y a 2 Questions à traiter et à séparer:
    • Quel modèle veut on?
    • Que veut on faire pour CMIP7
  • Sondage à mettre en place en amont? Avec ces questions + objectifs scientifiques des uns des autres -> volontaires? Inclure appel à volontaires dans mail de confirmation de la date?
  • Organisation:
    • Sessions de présentations générales de l'avancée des modèles incontournables pour remobiliser et informer
    • Sessions parrallèles par ex sur les groupes de travail des pirates. et autres.
    • Présentations par MIP et leur configuration? Au moins pour les MIP traditionnels.

* but / organisation des reunions pirates

  • on est toujours à la recherche d'un referent LMD(Z). Cela a été discuté en poihl. FH trop pris, Etienne potentiellement intéressé mais a de l'enseignement le jeudi ce trimestre. -> Envisager de changer de jour??? Julie et Juliette lancent un sondage?
  • LF ne se sent pas référent scientifique
  • bilan du jour positif.
  • Orchidee: présence à renforcer
  • Fréquence: tous les 15 jours.
  • Prochaine réunion: ESM?
    • Définition ESM? intégration des cycles? Bcp de différentes définitions.
    • Patricia: biais très similaires entre modele physique et modele-cycle
    • Myriam: un groupe de travail strato-REPROBUS démarre. Un run couplé REPROBUS existe.
    • pb qu'il faudra discuter; Dynamico: quelle version d'orchidee? Quelle config d'orchidee en fonction des besoin?

* - > Prochain ordre du jour NEMOv4

* TRACCS ?

  • Masa plutot en première ligne
  • accepté avec budget réduit - 50 millions euros sur 8 ans échelle nationale
  • 10 projets ciblés (PC) prédéfinis (70% budget) + 1 AAP géré par ANR (25%)
  • Conflit d'interet entre coordinateur de PC et déposants de ce type de projets ANR? A priori non.
  • PC: outilsCDD longs avec objectif (espoir?) pérennisation
  • Discussions organisationnelles en cours: quelle tutelle prend quel PC?
  • PC sur calcul intensif sera moin coupé que les autres car gros besoins.
  • Petit étalement dans le temps à prévoir car surcharge admin
  • Mail avec liste PC + responsable à prévoir -> Olivier

Thursday Septembre 1st 2022

Reunion pirates 01/09/2022

présents : Christiant Ethé, Laurent Fairhead, Olivier Boucher, Thibault Lurton, Juliette Mignot, Arnaud Caubel, Jérôme Servonnat, Etienne Vignon, Yann Merdeusoif, Thomas Dubos, Rémy Bonnet, Gurvan Madec, Julie Deshayes, Didier Swingedouw

  1. positionnement et organisation des réunions pirates

Identification par le CS restreint du besoin de relancer animation scientifique régulière.
-> 1ère réunion en juillet (Arnaud, Juliette, Christian, Julie, Olivier M, Yann) a permis de proposer un premier format, cf proposition en bas de cette page
-> Réunions ciblées, courtes (1h30), régulières (toutes les 2 semaines)
3 thématiques d'ores et déjà identifiées (proposition!)

[a] évolution de IPSL-CM “historique" (ie avec LMDZ) (+ tuning automatique?)
[b] mise en place des configurations “ICO” (avec Dynamico donc)
[c] évolution de l’ESM (ou tout sujet “cycles” au travers des différentes composantes).

Porteurs par labo : Olivier Marti pour LSCE, Julie Deshayes pour LOCEAN, XXX pour LMD -> XXX à identifier d'ici la prochaine réunion ([LF] responsable de cette action)

[OB]: coté CMC: réunions CS restreints régulières. + Retraite prévue fin nov-début décembre à confirmer.
[AC]: exp dvt CM6 a montré que ce genre de réunions fixes et ciblées fonctionne.
[TL]: Quelle articulation avec ESM2025, CMIP7, etc. y compris autres projets
[JD]: à construire. Il va y avoir TRACCS aussi. Mandat des "pirates"?
[OB]: ce groupe est complémentaire du CS restreint et des projets. Groupe de développement du modèle au sens large.
Un groupe séparé sur l'ESM a essayé de se lancer, mais peu de réunions (même si le dvt se fait). Peut etre bonne idée de l'intégrer ici. S'assurer que ESM ne diverge pas trop de la partie physique.
[JD]: même fonctionnement que CMIP6 sans la pression CMIP6. Réunions de travail effectif. Ne pas dupliquer d'autres groupes de travail qui existent par ailleurs. Groupe basé sur volontariat et envie/besoin de travailler ensemble et échanger.
[CE]: ce groupe va bénéficier des émulations liées aux projets divers.
[JD]: 30 min de chaque réunion dédiée à échange info sur projets et actions en cours.
[AC]: ne pas trop figer le plan de fonctionnement, on s'adaptera aux besoins.
+ Attention, toutes les 6 semaines du DYNAMICO (par ex) c'est trop peu.
[JD]: Indeed, quelle fréquence?
[TD]: Activité d'évaluation sera très utile pour tous ceux qui travaillent sur Dynamico. Important pour Dynamico de se confronter aux utilisateurs et autres développeurs du modèle.
Fréquence? Peut être coup de collier sur Dynamico au début?
Dynamico a mis en place réunions un peu informelles le vdi après midi. Peut être complémentaire de ce qui se passe dans le groupe ici.
[JD]: intermonitoring vs CESMEP vs codes maison?? Pour évaluation scientifique mais aussi identification éventuelle des bugs techniques.
[LF]: Souplesse dans le choix des thématiques chaque 2 semaines, selon les sujets brulants et dynamiques.
-> on garde toutes les 2 semaines jeudi matin 10h a priori.

  1. partage d’informations sur ce qui s’est passé pendant l’été

cf ​https://forge.ipsl.jussieu.fr/igcmg/wiki/Pirates2022[[BR]]

Partage d'info autour de l'implémentation de NEMO4, premiers tests namelists
NEMOv4 par défaut: état moyen pas mal, variabilité centennale.
namelist1: très proche du défaut.
namelist2: nnetau = 1 ie pénétration de tke sous la glace -> AMOC s'écroule. Pas de convection, ventilation le long des cotes écossaise, englacement. (remarque [GM] nnetau = 4: pénétration coupée dans l'HN et pas dans l'HS. )
-> on enterre namelist2. Pourtant on esperait regler pb de stratification trop forte et AMOC trop shallow avec cette pénétration.
Simu CdL: Dernières modifications de Casimir (intensite GM, intensite isoneutral): en cours chez Christian. A prolonger. puis a entrer dans une CESMEP avec des spread CMIP6 et séries temporelles. et la dynamico (voir plus bas) -> [Action JD]

Decide t on de faire tous nos tests en pd control?
Nécessite des abaques. Les a t on pour Dynamico? Pas de raison de ne pas les avoir.
pictrl permet de se comparer à CMIP6. Pb qu'on a pour l'instant = ordre 1 par rapport à des différences pd/pi...

Il y a du travail à faire sur glace de mer. amax est pour l'isntant le meme au nord et sud
[OB]: littérature: bilan énergetique TOA suggère bilan similaire dans les 2 hémisphères. Quid dans notre modèles? Rq: si on veut comparer aux obs, pd serait mieux.

Orchidee: [CE] a comparé à ce qu'on faisait dans CMIP6. Pas de surprise

CESMEP sur Dynamico. : [AC] a déjà fait tourné. Lien à partager. Les infos sont disponibles là : ​https://forge.ipsl.jussieu.fr/igcmg/wiki/IPSLCM6/IPSL-CM6A-ATM-ICO#Validation @tous: A REGARDER d'ici la prochaine réunion!!'''

[JS]: Stéphae Sénési a presque fini automatisé lancement CESMEP par LibIGCM
Dynamico: on a actuellement 50-100 ans en couplé avec du NEMO3.6.
Tests en climat froid par OM: cf wiki 5 aout. Take home message: ca pete pas.

  1. modalités de partage d’informations

mailing list cmip6_dev ?
Action pour [OB]:dupliquer cmip6_dev + nouveau nom: icmc_dev: dupliquer + corriger les noms (avec [JD]) + eventuellement appel général.
​https://forge.ipsl.jussieu.fr/igcmg/wiki/Pirates2022
slack ? Jugé pas indispensable pour le moment

Thursday August 5th 2022 : test froids

Olivier M. a fait deux tests froids IPSLCM 6.5 pour voir si la glace de mer se stabilise avec NEMO4 et SI3. On crée le froid en mettant dans LMDz le CO2 à 180 ppm (simulation CM65-OM01) et 100 ppm (simulation CM65-OM02).

Les deux simulations ont plus de 200 ans. Elles sont stabilisées en surface, glace de mer comprise. Les SST globales sont descendues à 16°C et 14,5 °C. Les t2m (température de l'air à 2m) sont à 11,5°C et 9,5°C. La simulation à 100 ppm a un contenu de chaleur océanique qui continue à baisser.

Les maxima d'épaisseurs de glace de mer atteignent 20m, mais très localement. A 100 ppm, la glace est arrivée en Bretagne nord. Elle pourrait gêner les manœuvres militaires chinoises et états-uniennes dans le détroit de Formose ...

A voir :

Les simulations sont toujours en machine, pour atteindre 400 ans.

Tuesday July 26th 2022

Olivier B

  • Mise à jour des forçages pour le nouveau rayonnement ECRAD: aérosols tropos (Lionel et Abderrahmane), aérosols stratos (Thibaut), solaire (Thibaut et Abderrahmane), en cours

Monday July 25th 2022

Julie, Olivier M, Guillaume G, Juliette

  • Tests namelist juillet 2022: colonne H tableau ci dessus ( a priori plus "saine" et à jour - dixit Seb Masson, Gurvan et Casimir en partie). Avec ca, 2 types de test: CM6.5_work
    • 1 simu pi lancées par Guillaume Gachon
    • 2 simus froides lancées par Olivier: 180 et 100ppm
      • elles ne crashent pas, apres 50-70 ans. Il faut attendre encore un peu (100 ans plutot). Volume de glace en Arctique augmente bien. Les simus restent en machine a priori 300 ans (mais Olivier doit les stopper pour mettre + de proc)
      • Remarque: la simu Namelist-1 de Guillaume Gachon part de scratch (car experiences utilisées aussi pour travail sur tuning), alors que la version svn part d'une ref SKL.
  • Prochains axes de tests:
    • jouer avec nn_etau? Actuellement, nn_etau = 0 ie on n'injecte rien sous la couche de melange. Or on est trop statifié, comme dans 6ALR -> Tester nn_etau = 1?
    • implémenter derniers développements Casimir -> voir avec Christian car changements de Casimir impliquent modification de code non distribués.
    • Routage flux d'eau douce AA et cavités (moyen terme)
    • Tuning Glace de mer SI3 pas encore abordé -> voir avec Martin à la rentree.
  • Diagnostiques?
    • Sur climat froid, enjeu = surtout stabilité du modèle, pas encore d'enjeu de diagnostique
    • Sur climat pi et chaud: se mettre à jour sur developpements CESMEP et autres (esmvaltool).
      • Remarque cesmep: un peu lourd (desactiver des choses a la demande), lien avec climaf pas evident.
      • Travail cesmep ete 2022 (Guillaume Gachon) : cartes de validation de base + monitoring propre AMOC, ACC, glace de mer, etc
  • Prochaine réunion:
    • 1er septembre - 10h-11h30, lancement des pirates 1 jeudi / 2
    • Olivier met un message sur slack au moment de son depart en vacances (15 aout) pour informer de la stabilite de ses simus

Tuesday July 19th 2022

e-mail de cadrage envoyé par Julie Deshayes 

Bonjour, 
hier (lundi 18/07) nous (Olivier Marti, Yann, Arnaud, Christian, Juliette et moi) avons discuté de ces réunions pirates, et imaginé la proposition ci-dessous d’animation dans ICMC, pour rendre effective la dimension “continue” de l’evolution des modèles de climat à l’IPSL :
 

[1] un créneau fixe dans la semaine (jeudi matin ? vendredi matin ? lundi matin ? à decider selon nos agendas — ie creneaux seminaires et autres —  et pour rendre le plus efficace le travail collaboratif), point de rdv de la construction continue des modèles de climat à l’IPSL (un peu comme le POIHL du lundi au LMD)

[2] une durée fixe aussi, pour ne pas chambouler le planning de chacun, fixée à 1h30 pour laisser 30 min de discussion aux affaires courantes et 1h à une thématique dédiée (voir point 4 ci-dessous) - au passage, l’information sur ces points de discussion sera partagée très largement, pour que tous les potentiels intéressés s’identifient et y participent.

[3] une fréquence fixe, de 1 reunion toutes les 2 semaines ***, tout au long de l'année : pour conserver le momentum dans la durée, avec prise de notes en live et partage sur une page de https://forge.ipsl.jussieu.fr/igcmg/wiki/Pirates2022 

[4] une thématique tournante pour couvrir les 3 sujets principaux qui occupent ICMC, par exemple en ce moment : 
	[a] evolution de IPSL-CM “historique" (i.e. avec LMDZ) 
	[b] mise en place des configurations “ICO” (avec Dynamico donc) 
	[c] évolution de l’ESM (ou tout sujet “cycles” au travers des differentes composantes).
Il faudra faire evoluer ces sujets au gré des actions en cours, par exemple [a] pourra inclure (à moyen ou long terme, selon les avancées) la validation et/ou le tuning automatique, tandis que les questions relatives aux differentes resolutions spatiales seront plutot traitées dans [b], et les histoires de couplage avec la calotte dans [c]. 

[5] trois pirates pour porter cette nouvelle animation transverse dans ICMC, issus des 3 “gros” labos de l’IPSL à savoir LSCE (Olivier Marti accepte d’être désigné volontaire), LOCEAN (je suis moi-même volontaire, cf ce message) et LMD (à determiner… après le 15 août, date à laquelle Etienne rentre de vacances ?  ).

qu’en pensez vous ? 
Julie

*** avec 2 réunions par semaine, on “parlera” des configs ICO toutes les 6 semaines. si on juge que  c’est insuffisant, il faudra passer à 1 réunion par semaine...

Monday July 18th 2022

Présents : Julie, Juliette, Christian, Yann, Olivier, Arnaud

Objectif : partage d'information autour du travail réalisé et à venir pour intégrer NEMO4 dans le couplé IPSL.

  • De IPSLCM6 à IPSLCM7
    • voir slides là : Reunion_NEMO4.pdf​
    • le changement de version d'une composante (tout comme l'activation d'un processus) doit rentrer dans un processus d'intégration continue. Besoin d'outils pour cela. Ce n'est pas comme cela qu'on fonctionne actuellement : alternance entre des phases d'intégration intenses et des phases plus calmes.
    • besoin de diffuser l'information de la feuille de route vers IPSLCM7. Plus généralement, il y aurait peut etre une utilité à avoir des réunions pleinières/élargies plus souvent.
  • NEMO 4
    • voir slides là : SimusCM65.pdf​
    • plusieurs versions de NEMO ont été testées en couplé océan-atmosphère. Un travail est fait aussi en forcé (simulations tests, atelier lecture namelist NEMO 4.2,...)
    • la dernière version avec NEMO 4.2 est moins englacée mais relativement stable. Les bilans d'eau douces sont similaires à ceux de IPSL-CM6A-LR.
    • même si c'est une version intermédiaire (étape vers IPSLCM7), il y a une utilité à pousser la validation : le travail incrémental sur les versions de NEMO a un interêt scientifique et il est à faire maintenant.
  • Test de NEMO4 en climat froid
    • le modèle IPSL-CM6A-LR est fragile en climat froid (~180 ppm de CO2) : formation de hauteur de glace aberrante en Arctique. Runs holocène : plantages réguliers en Arctique, instabilités de NEMO3.6.
    • choses intéressantes à tester dans NEMO4.2. Action : Olivier teste avec la namelist que lui envoie Julie.
    • besoin d'intégrer les climats froids dans les tests classiques du modèle.
  • Vers la plus haute résolution
    • voir slides là : Reunion_NEMO4.pdf​
    • objectif final du projet ESIWACE2 : couplé océan-atmosphere autour de 10km. Intérêt technique mais peu d'interêt scientifique pour cette résolution : on reste sur de la paramétrisation et le réglage du modèle peut couter tres cher.
    • la résolution intermédiaire à 25 km océan-atmosphère présente un interêt scientifique. Une version "technique" ICO-LMDZ-ORCHIDEE-NEMO4 tourne déjà à environ 6 SYPD sur 12 000 coeurs.
    • la complexité (activation de PISCES par ex) à 50km est plus intéressante qu'une version moins "complexe" à 10 km.
    • 20 SYPD minimum pour faire du climat ? En tout cas pour un exercice de type CMIP6 oui. Il doit être faisable d'améliorer les perfs de calcul actuelles pour parvenir à cette cible.
    • besoin de se mettre d'accord sur les sorties à regarder.
    • sacalabilité de NEMO4 meilleure que NEMO3.
  • Organisation pour avancer
    • besoin de réunions fréquentes. Proposition : 1 réunion par semaine (créneau fixé : jeudi matin ?) avec des thèmes tournants (ESM, DYNAMICO dans le couplé, CM6, NEMO4,...). Premier point fixé lundi 25 juillet à 10h pour les personnes présentes. Proposition à soumettre à Masa et Olivier.
    • Actions
      • Julie répond au mail de Masa pour lui proposer ce format. Un représentant par labo serait utile : Olivier pour le LSCE, Julie pour le LOCEAN et Etienne (?) pour le LMD?
      • Juliette crée une page wiki
      • Arnaud fait le CR de cette réunion.

Attachments (3)