| 2 | |
| 3 | ===Thursday 15th September=== |
| 4 | Réunion pirates 15 septembre 2022 |
| 5 | Présents: Yann, Claire, Myriam, Julie, Juliette, Olivier B, Arnaud, Brady, Christian, Claire, Laurent F, Martin, Olivier M, Patricia, Thibault, Guillaume G,Thomas D |
| 6 | |
| 7 | DYNAMICO |
| 8 | Rappel: https://forge.ipsl.jussieu.fr/igcmg/wiki/IPSLCM6/IPSL-CM6A-ATM-ICO#Validation |
| 9 | |
| 10 | comparaison de |
| 11 | |
| 12 | CM61-LR-pi-03 : piControl CMIP6, réalisé avec IPSLCM6.1 |
| 13 | |
| 14 | TGCC : (/ccc/store/cont003) gencmip6/p86maf/IGCM_OUT/IPSLCM6/PROD/piControl/CM61-LR-pi-03 |
| 15 | |
| 16 | VALID-CM622-LR.01 : 50 ans, piControl IPSLCM6.2, version à jour LMDZ (physique v6.2) |
| 17 | |
| 18 | TGCC : gencmip6/p86caub/IGCM_OUT/IPSLCM6/DEVT/piControl/VALID-CM622-LR.01 |
| 19 | |
| 20 | 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 |
| 21 | |
| 22 | TGCC : gencmip6/p86caub/IGCM_OUT/IPSLCM7/DEVT/piControl/TEST-CM72-SIMPLE-ROUTING |
| 23 | |
| 24 | TEST-CM72-SIMPLE-ROUTING.06 : valeur cld* corrigée 0,00080 pour recalage température globale |
| 25 | |
| 26 | TGCC : gencmip6/p86caub/IGCM_OUT/IPSLCM7/DEVT/piControl/TEST-CM72-SIMPLE-ROUTING.06 |
| 27 | |
| 28 | |
| 29 | Intermonitoring : https://thredds-su.ipsl.fr/thredds/fileServer/tgcc_thredds/work/p86caub/INTERMONITORING/intermonit_VALIDATION_DYNAMICO_COUPLE/index.html |
| 30 | |
| 31 | 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 |
| 32 | |
| 33 | ressemble bug masques Myriam VLR: Olivier Marti Arnaud et Yann verifient: ce probleme n'etait pas la dans les versions précédentes |
| 34 | |
| 35 | 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 euvre dans un 2e temps. |
| 36 | |
| 37 | |
| 38 | wfo (E-P-R sous la glace) signe différent par rapport à IPSLCM6ALR |
| 39 | |
| 40 | 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 ? |
| 41 | |
| 42 | Difficile de raisonner ainsi, friver pas indépendant de evap, pb ssh prioritaire |
| 43 | |
| 44 | 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. |
| 45 | |
| 46 | 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 |
| 47 | |
| 48 | 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. |
| 49 | |
| 50 | 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 |
| 51 | |
| 52 | 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 |
| 53 | |
| 54 | |
| 55 | |
| 56 | 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 |
| 57 | |
| 58 | PB Atlas ocean: cartes SST et SSS non exploitables |
| 59 | |
| 60 | [Arnaud]: ca marche quand la référence est OBS plutot que un autre run |
| 61 | |
| 62 | Cartes ATM montrent biais chauds au sud de 55°S et au nord de 60°N, en hiver. |
| 63 | |
| 64 | |
| 65 | [FH] Explications possibles: |
| 66 | |
| 67 | 1) un problème propre à Dynamico, mais je n'y crois pas trop. |
| 68 | |
| 69 | 2) une dé-compensation d'erreur structurelle par rapport à la version à résolution équivalente de LMDZlonlat (ou LMDZlola, c'est plus mignon). |
| 70 | |
| 71 | 3) entre les deux : les deux dynamiques ont leurs problèmes et le tuning de la physique est capable de compenser ça en partie. |
| 72 | |
| 73 | 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. |
| 74 | |
| 75 | |
| 76 | [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?? |
| 77 | |
| 78 | |
| 79 | [Myriam] Ces biais viennent ils de la dynamique plutot? Origine plutot equatoriale (hypothèse ) |
| 80 | |
| 81 | |
| 82 | [Yann] Autres différences fondamentales Dynamico vs LMDZlola: |
| 83 | |
| 84 | Dynamico a une moins bonne résolution aux poles pour nombre de points équivalents (mis meilleur résolution à l'equateur) |
| 85 | |
| 86 | [FH]: Le défaut biais chaud se corrige de façon très systématique pour LMDZlola avec l'augmentation de la résolution https://web.lmd.jussieu.fr/~hourdin/PUBLIS/LMDZ5.pdf pour une étude sytématique avec LMDZ5A Fig 9 de https://web.lmd.jussieu.fr/~hourdin/PUBLIS/2019MS001892.pdf pour une comparaison entre LMDZ6ALR et HR (50km de résolution). Pour rappel pour celles et ceux qui n'avaient pas suivi les épisodes précédents, la raison du choix d'avoir autant de latitudes que de longitudes dans LMDZlola, c'est à dire une grille très anisotrope à l'équateur avec une grille deux fois plus fine en latitude, vient du fait que l'accroissement en résolution était beaucoup plus important que l'acroissement en longitude pour mieux positioner les jets. |
| 87 | [AC] atlas de comparaison des simulations HighResMIP forcées faites avec dynamique LMDZ et DYNAMICO aux différentes résolutionsh https://thredds-su.ipsl.fr/thredds/fileServer/tgcc_thredds/work/p86caub/C-ESM-EP/force_HRMIP_ico_lonlat_p86caub/Atmosphere_Surface/atlas_Atmosphere_Surface_force_HRMIP_ico_lonlat.html |
| 88 | Attention, dans ces figures, comparaison directe aux observations pas possible. On peut comparer les runs les uns par rapport aux autres. |
| 89 | [YM] Sur HS: résolution semble aller dans le bon sens. Effet de la résolution par rapport à ce que disait Yann plus haut. |
| 90 | |
| 91 | Dissipation. : est on plus ou moins dissipé que dans LMDZlola? |
| 92 | |
| 93 | 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 ? |
| 94 | |
| 95 | [AC] Mêmes valeurs de dissipation |
| 96 | |
| 97 | [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 |
| 98 | |
| 99 | Action: Définir une série de test avec moins de dissipation? (Thomas, Yann, Arnaud, ?) |
| 100 | |
| 101 | Commencer par du forcé. |
| 102 | |
| 103 | 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. |
| 104 | |
| 105 | Orographie: cf anomalies de vent méridien dans le Pacifique qui semblent améliorer. |
| 106 | |
| 107 | |
| 108 | En forcé? |
| 109 | |
| 110 | atlas de différence entre une version LMDZico et LMDZlola |
| 111 | |
| 112 | Intermonitoring : https://thredds-su.ipsl.fr/thredds/fileServer/tgcc_thredds/work/p86caub/INTERMONITORING/intermonit_VALIDATION_DYNAMICO_FORCE/index.html |
| 113 | |
| 114 | Atlas CESMEP 10 ans: https://thredds-su.ipsl.fr/thredds/fileServer/tgcc_thredds/work/p86caub/C-ESM-EP/force_ico_lonlat_p86caub/C-ESM-EP_force_ico_lonlat.html |
| 115 | |
| 116 | avec |
| 117 | |
| 118 | ICOLOR-SROUT.02 : 10 ans expérience clim, DYNAMICO-LMDZ-ORCHIDEE |
| 119 | |
| 120 | TGCC : gencmip6/p86caub/IGCM_OUT/LMDZOR/DEVT/clim/LOR-SROUT.02 |
| 121 | |
| 122 | LOR-SROUT.03 : 10 ans expérience clim, LMDZ-ORCHIDEE |
| 123 | |
| 124 | TGCC : gencmip6/p86caub/IGCM_OUT/ICOLMDZOR/DEVT/test/ICOLOR-SROUT.03 |
| 125 | |
| 126 | |
| 127 | |
| 128 | 1) sur les biais en température (mais avec des SSTs imposées ce n'est pas évident) |
| 129 | |
| 130 | Oui clairement sur continent, même ordre de grandeur, max en hiver. |
| 131 | |
| 132 | 2) sur les flux à la surface (comme c'est le cas sur les problèmes de bord Est) |
| 133 | |
| 134 | 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. |
| 135 | |
| 136 | [FH] pense que c'est ca la source. |
| 137 | |
| 138 | [JJ] Raisonnement uniquement HS, quid HN? |
| 139 | |
| 140 | 3) sur la structure méridienne moyenne (latitude-pression) de la température et de la pression. |
| 141 | |
| 142 | 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 |
| 143 | |
| 144 | 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? |
| 145 | |
| 146 | |
| 147 | Installation des biais? |
| 148 | |
| 149 | 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"). |
| 150 | |
| 151 | Tuning automatique / history matching: quel effet sur ces configs, quelles efficacité? |
| 152 | |
| 153 | |
| 154 | Analyse: (ceci est un échange de mail OB FH remis ici pour mémoire) |
| 155 | |
| 156 | [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... |
| 157 | |
| 158 | |
| 159 | other topics : |
| 160 | |
| 161 | NEMO en climat froid |
| 162 | cf pres Olivier Marti |
| 163 | IPSLCM6 plante systématiquement en glaciaire |
| 164 | Tests NEMO4.2 IPSLCM6.5 |
| 165 | Tests froids en CO2 seul, sans les calottes, namelist NEMO4.2 fournie par Guillaume Gachon |
| 166 | -> ça tient. |
| 167 | épaisseur de glace plafonne à 20m dans Arctique. |
| 168 | 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. |
| 169 | |
| 170 | 4xCO2 tient aussi |
| 171 | |
| 172 | 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) |
| 173 | |
| 174 | |
| 175 | futur de CMIP7 - ou du moins ce qu'on en sait... |
| 176 | OB: ca commence à s'organiser. |
| 177 | |
| 178 | Thibault volontaire pour prendre en charge forçages. |
| 179 | |
| 180 | C'st à nous de nous orgaiser pour définir ce que l'on veut. |
| 181 | |
| 182 | A discuter à la retraite? |
| 183 | |
| 184 | OB sera au NCAR en novembre pour discuter avec JF Lamarque |
| 185 | |
| 186 | Discussions avant l'été au moment de workshop WCRP. Vers des ensemles moins opportunistes? |
| 187 | |
| 188 | A l'IPSL, calendrier CMIP7 risque de se téléscoper avec changement calculateur. |
| 189 | |
| 190 | Prochain calculateur TGCC attendu 2024-2025, va demander du temps d'adaptation. Cela va impliquer des choix par rapport à CMIP7. |
| 191 | |
| 192 | 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. |
| 193 | |
| 194 | Reflexion sur nb de MIP, nb de simulations etc à avoir. |
| 195 | |
| 196 | |
| 197 | infos sur retraite des 17-18 novembre |
| 198 | date 17/18 novembre maintenue - à confirmer dans qqs jours après vérifcations de présence de certains membres. domaine St Paul |
| 199 | |
| 200 | toutes suggestions d'atelier bienvenues ! |
| 201 | |
| 202 | |
| 203 | Objectifs: |
| 204 | |
| 205 | Remobiliser un peu plus large que le groupe pirates (ca 80 personnes) |
| 206 | |
| 207 | Faire une réunon utile pour nos développements |
| 208 | |
| 209 | |
| 210 | Il y a 2 Questions à traiter et à séparer: |
| 211 | |
| 212 | Quel modèle veut on? |
| 213 | |
| 214 | Que veut on faire pour CMIP7 |
| 215 | |
| 216 | 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? |
| 217 | |
| 218 | Sessions de présentations générales de l'avancée des modèles incontrounables pour remobiliser et informer |
| 219 | Sessions parrallèles par ex sur les groupes de travail des pirates. et autres. |
| 220 | |
| 221 | Présentations par MIP et leur configuration? Au moins pour les MIP traditionnels. |
| 222 | |
| 223 | |
| 224 | |
| 225 | but / organisation des reunions pirates |
| 226 | |
| 227 | 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? |
| 228 | |
| 229 | LF ne se sent pas référent scientifique |
| 230 | |
| 231 | |
| 232 | bilan du jour positif. |
| 233 | |
| 234 | Orchidee: présence à renforcer |
| 235 | |
| 236 | Fréquence: tous les 15 jours. |
| 237 | |
| 238 | Prochaine réunion: ESM? |
| 239 | |
| 240 | définition ESM? intégration des cycles? Bcp de différentes définitions. |
| 241 | |
| 242 | Patricia: biais très similaires entre modele physique et modele-cycle |
| 243 | |
| 244 | Myriam: un groupe de travail strato-REPROBUS démarre. Un run couplé REPROBUS existe. |
| 245 | |
| 246 | pb qu'il faudra discuter; Dynamico: quelle version d'orchidee? |
| 247 | |
| 248 | Quelle config d'orchidee en fonction des besoin? |
| 249 | |
| 250 | |
| 251 | - > Prochain ordre du jour NEMOv4 |
| 252 | |
| 253 | |
| 254 | |
| 255 | TRACCS ? |
| 256 | |
| 257 | Masa plutot en première ligne |
| 258 | |
| 259 | accepté avec budget réduit - 50 millions euros sur 8 ans échelle nationale |
| 260 | |
| 261 | 10 projets ciblés (PC) prédéfinis (70% budget) + 1 AAP géré par ANR (25%) |
| 262 | |
| 263 | Conflit d'interet entre coordinateur de PC et déposants de ce type de projets ANR? A priori non. |
| 264 | |
| 265 | PC: outilsCDD longs avec objectif (espoir?) pérennisation |
| 266 | |
| 267 | Discussions organisationnelles en cours: quelle tutelle prend quel PC? |
| 268 | |
| 269 | PC sur calcul intensif sera moin coupé que les autres car gros besoins. |
| 270 | |
| 271 | Petit étalement dans le temps à prévoir car surcharge admin |
| 272 | |
| 273 | |
| 274 | Mail avec liste PC + responsable à prévoir -> Olivier |
| 275 | |
| 276 | |
| 277 | |
| 278 | |
| 279 | |
| 280 | |
| 281 | |