Changes between Version 74 and Version 75 of Pirates2022


Ignore:
Timestamp:
10/14/22 12:56:38 (19 months ago)
Author:
jdlod
Comment:

--

Legend:

Unmodified
Added
Removed
Modified

  • Pirates2022

    v74 v75  
    22 
    33=== Friday 14th October === 
    4  
    5  * DYNAMICO 
    6   * Biais chaud Antarctique vu en couplé  
    7    * Intermonitoring : https://thredds-su.ipsl.fr/thredds/fileServer/tgcc_thredds/work/p86caub/INTERMONITORING/intermonit_VALIDATION_DYNAMICO_COUPLE_3/index.html 
    8    * Atlas : ​https://thredds-su.ipsl.fr/thredds/fileServer/tgcc_thredds/work/p86caub/C-ESM-EP/couple_dynamico_ref_piControl_2_p86caub/C-ESM-EP_couple_dynamico_ref_piControl_2.html 
    9    * Tests avec sponge=0 : https://thredds-su.ipsl.fr/thredds/fileServer/tgcc_thredds/work/p86caub/C-ESM-EP/force_ico_lonlat_sponge_p86caub/C-ESM-EP_force_ico_lonlat_sponge.html 
    10    * Tests dissipation et résolutions : ​https://thredds-su.ipsl.fr/thredds/fileServer/tgcc_thredds/work/p86caub/C-ESM-EP/force_ico_lonlat_dissip_3D_p86caub/C-ESM-EP_force_ico_lonlat_dissip_3D.html 
    11   * On perd de l'eau dans l'océan (cf intermonitoring ci dessus) 
    12    * Résultat script waterbudget Olivier B : https://thredds-su.ipsl.fr/thredds/fileServer/tgcc_thredds/work/p86caub/WATER_BUDGET_THREDDS/TEST-CM72-SIMPLE-ROUTING.10/waterbudget_TEST-CM72-SIMPLE-ROUTING.10_1860_1869.out 
    13     * le modèle est en train d'accumuler un peu de neige sur les continents à hauteur de 0.00267 Sv 
    14     * le modèle est en train d'accumuler de la neige sur les calottes continentales, à hauteur de 0.0559 Sv 
    15     * => accumulation de neige  ? état initial ?  
    16  
    174 
    185Présents:  
    … …  
    207 
    218 
    22 DYNAMICO 
    23 2 soucis ont été mis en lumière la dernière fois: biais chaud AA & perte eau, dérive ssh 
    24  
    25 * 1. Biais chaud AA: forcé, test en dissipation.  
    26 [Sebastien Fromang] Probleme identifie: Sponge:  ce qui se passe dans couche eponge au top de l'atmosphere. necessaire pour éviter spurious réflexions d'ondes dues à la raréfaction de l'air au sommet de l'atmosphère.  
    27 LMDZ: relaxation à la moyenne zonale vent 
    28 DYNAMICO: pas facile de faire relaxation moyenne zonale. Ce qu'il y avait par defaut = relaxation à 0. Pb aux poles: tue le vortex polaire et donc biais chaud. n'apparait qu'en hiver. 
    29 Sans couche éponge: biais chauds corrigés en hiver. Reste biais de vent zonal et biais froid.  
    30  
    31     Fred : Ne pas négliger ce biais froid résiduel notamment au nord. voyait aussi des anomalies symetriques dans DYNAMICO - les memes qu'identifiees ci-dessus ou d'autres qui s'ajoutent ? 
    32  
    33 Atlas couplés sans couche éponge: ​​https://thredds-su.ipsl.fr/thredds/fileServer/tgcc_thredds/work/p86caub/C-ESM-EP/couple_dynamico_ref_piControl_2_p86caub/C-ESM-EP_couple_dynamico_ref_piControl_2.html 
    34  
    35     Sponge = 0 derniere colonne. ATTENTION: 12 ans seulement! 
    36  
    37     -> correction du biais chaud en hiver au sud. Pas de correction en été (DJF) car pas de vortex, pas d'effet de sponge 
    38  
    39     + biais arctique DJF semble renforcé (attention comparaison à CM6A-LR pictrl. Peut etre qu'on corrige. A verifier...) . Bcp de changement DJF Subpolaire et polaire nord. Prolonger la simulation pour verifier impact AMOC.  
    40  
    41 -> comment corriger biais vent zonal DJF dans hemisphere sud ?  
    42  
    43  
    44 Simulations forcées à plus haute résolution: https://thredds-su.ipsl.fr/thredds/fileServer/tgcc_thredds/work/p86caub/C-ESM-EP/force_ico_lonlat_dissip_3D_p86caub/C-ESM-EP_force_ico_lonlat_dissip_3D.html 
    45 (NBP80 = 100km, vs 200km pour les autres resolutions. 3 dissipations differentes à chaque fois) 
    46 Amélioration en DJF avec la résolution, possiblement grâce à une amélioration de la topographie et/ou trait de cote AA, alors que en cette saison, le jet est décalé vers le sud et donc plus sensible à la topo. -> Action: refaire NBP80 and Sponge=0 
    47 Remarque surcout: NBP80 (80x80x10 soit 4x+ de points de grille que NBP40) = x8 (en forcé... à pondérer par cout relatif de l'atm dans le couplé pour les simus couplées) 
    48  
    49     [FH] avait l'impression que les biais persistent dans Runs HIghResMIP (50 km). [SF] verifie.  
    50  
    51     [FH] souligne que choix pour grille LMDZ favorise representation des ondes dans hautes latitudes hemisphere sud 
    52  
    53  
    54 Discussion: remettre une couche éponge avec relaxation à moyenne zonale? Cela necessite un petit developpement de code dans DYNAMICO. [FH]: à faire!  
    55 [SF]: faire un run forcé LMDZ SANS sponge, pour verifier de l'interet de la sponge.  
    56 [FH]: insiste sur le besoin de faire ce developpement dans DYNAMICO quand meme, car necessité acceptée largement dans communauté modelisation atm 
    57 [SF]: HighResMIP a été publié sans sponge.  
    58 [YM]: est ce que si on ajoutait des niveaux, on reglerait le pb?  
    59  
    60     [SF]: a priori non on dissipe moins les ondes, donc + de reflexion? 
    61  
    62  
    63 Actions:   
    64     Prolonger le couplé NBP40 qui a tourné 12 ans sans SPONGE (juste pour voir... cf besoin de tuning +/- automatique, ci-dessous) 
    65     Forcé LMDZ sans sponge 
    66     Forcé DYNAMICO NBP80 avec sponge  
    67     Tester aussi NBP60. Avantage d'etre resolution equivalente LMDZ aux moyennes latitudes.  
    68     Developper Sponge moyenne zonale dans DYNAMICO 
    69  
    70 [FH] Utiliser outils d'exploration paramétrique systématique avec HM sur DYNAMICO car physiq.def tres ajustée pour LMDZ, avec compensation de biais qui ne correspondent surement pas à DYNAMICO 
    71 [JD] est-il raisonnable d'étudier la sensibilite à la resolution sans faire retuning à chaque fois? 
    72 [FH] biais vents zonaux s'ameliorent systematiquement avec augmentation de la resolution, et de maniere plus flagrante que ce qu'on s'attend a corriger avec compensation d'erreurs 
    73  
    74 Remarque resolution:  
    75     NBP40 = 40x40x10 = 16 000 cellules 
    76     LMDZ LR = 20 000 points. eq NBP45.  
    77  
    78 [FH] rappelle qu'il a développé une param à la Gent & McWilliams pour corriger le transport zonal vers les poles.  
     9 * DYNAMICO : Biais chaud Antarctique vu en couplé   
     10 
     11   Intermonitoring : https://thredds-su.ipsl.fr/thredds/fileServer/tgcc_thredds/work/p86caub/INTERMONITORING/intermonit_VALIDATION_DYNAMICO_COUPLE_3/index.html 
     12     
     13   Atlas : ​https://thredds-su.ipsl.fr/thredds/fileServer/tgcc_thredds/work/p86caub/C-ESM-EP/couple_dynamico_ref_piControl_2_p86caub/C-ESM-EP_couple_dynamico_ref_piControl_2.html 
     14 
     15   Tests avec sponge=0 : https://thredds-su.ipsl.fr/thredds/fileServer/tgcc_thredds/work/p86caub/C-ESM-EP/force_ico_lonlat_sponge_p86caub/C-ESM-EP_force_ico_lonlat_sponge.html 
     16     
     17   Tests dissipation et résolutions : ​https://thredds-su.ipsl.fr/thredds/fileServer/tgcc_thredds/work/p86caub/C-ESM-EP/force_ico_lonlat_dissip_3D_p86caub/C-ESM-EP_force_ico_lonlat_dissip_3D.html 
     18 
     19[SF] Probleme identifie dans Sponge au TOA (ce qui se passe dans couche eponge au top de l'atmosphere, necessaire pour éviter spurious réflexions d'ondes dues à la raréfaction de l'air au sommet de l'atmosphère)\\ 
     20  LMDZ: relaxation à la moyenne zonale vent. \\ 
     21  DYNAMICO: pas facile de faire relaxation moyenne zonale, donc par defaut : relaxation à 0. Pb aux poles: tue le vortex polaire et crée donc biais chaud qui n'apparait qu'en hiver.\\ 
     22  Test sans couche éponge (derniere colonne dans atlas): biais chauds corrigés en hiver au sud, mais reste biais de vent zonal et biais froid. pas de correction en été au sud car pas de vortex donc pas d'effet de sponge. + biais arctique DJF semble renforcé (attention comparaison à CM6A-LR pictrl. Peut etre qu'on corrige. A verifier...) . attention : tests sur 12 ans seulement\\ 
     23  [FH] Ne pas négliger ce biais froid résiduel notamment au nord. voyait aussi des anomalies symetriques dans DYNAMICO - les memes qu'identifiees ci-dessus ou d'autres qui s'ajoutent ?\\ 
     24  -> comment corriger biais vent zonal DJF dans hemisphere sud ?  
     25 
     26Simulations forcées à plus haute résolution  
     27(NBP80 = 100km, vs 200km pour les autres resolutions,  3 dissipations differentes à chaque fois) :  
     28Amélioration en DJF avec la résolution, possiblement grâce à une amélioration de la topographie et/ou trait de cote AA, alors que en cette saison, le jet est décalé vers le sud et donc plus sensible à la topo.  
     29 
     30  -> Action: refaire NBP80 and Sponge=0. Remarque surcout: NBP80 (80x80x10 soit 4x+ de points de grille que NBP40) = x8 (en forcé... à pondérer par cout relatif de l'atm dans le couplé pour les simus couplées) 
     31 
     32  [FH] avait l'impression que les biais persistent dans Runs HIghResMIP (50 km). [SF] verifie.  
     33 
     34  [FH] souligne que choix pour grille LMDZ favorise representation des ondes dans hautes latitudes hemisphere sud 
     35 
     36Discussion: remettre une couche éponge avec relaxation à moyenne zonale? Cela necessite un petit developpement de code dans DYNAMICO. \\ 
     37  [FH]: à faire! \\ 
     38  [SF]: faire un run forcé LMDZ SANS sponge, pour verifier de l'interet de la sponge \\ 
     39  [FH]: insiste sur le besoin de faire ce developpement dans DYNAMICO quand meme, car necessité acceptée largement dans communauté modelisation atm\\ 
     40  [SF]: HighResMIP a été publié sans sponge \\ 
     41  [YM]: est ce que si on ajoutait des niveaux, on reglerait le pb? \\ 
     42  [SF]: a priori non on dissipe moins les ondes, donc + de reflexion? 
     43 
     44 
     45-> Actions:  \\ 
     46  Prolonger le couplé NBP40 qui a tourné 12 ans sans SPONGE (juste pour voir... cf besoin de tuning +/- automatique, ci-dessous)\\ 
     47  Forcé LMDZ sans sponge\\ 
     48  Forcé DYNAMICO NBP80 avec sponge \\ 
     49  Tester aussi NBP60 : Avantage d'etre resolution equivalente LMDZ aux moyennes latitudes. \\ 
     50  Developper Sponge moyenne zonale dans DYNAMICO 
     51 
     52[FH] insiste sur la pertinence d'utiliser outils d'exploration paramétrique systématique avec HM sur DYNAMICO car physiq.def tres ajustée pour LMDZ, avec compensation de biais qui ne correspondent surement pas à DYNAMICO\\ 
     53  [JD] est-il raisonnable d'étudier la sensibilite à la resolution sans faire retuning à chaque fois?\\ 
     54  [FH] biais vents zonaux s'ameliorent systematiquement avec augmentation de la resolution, et de maniere plus flagrante que ce qu'on s'attend a corriger avec compensation d'erreurs 
     55 
     56Remarque resolution: \\ 
     57  NBP40 = 40x40x10 = 16 000 cellules\\ 
     58  LMDZ LR = 20 000 points. eq NBP45.  
     59 
     60[FH] rappelle qu'il a développé une param à la GM (diffusion isopycnale dans l'ocean par terme advectif supplementaire) pour corriger le transport zonal vers les poles.  
    7961 
    8062 
    8163* 2. Perte en eau 
    82 test avec nouveau routage, avec poids mieux pris en compte a priori : produit la meme  
     64 
     65Dans DYNAMICO on perd de (toujours) de l'eau dans l'océan ... 
     66 
     67Résultat script waterbudget Olivier B : https://thredds-su.ipsl.fr/thredds/fileServer/tgcc_thredds/work/p86caub/WATER_BUDGET_THREDDS/TEST-CM72-SIMPLE-ROUTING.10/waterbudget_TEST-CM72-SIMPLE-ROUTING.10_1860_1869.out 
     68 
     69  le modèle est en train d'accumuler un peu de neige sur les continents à hauteur de 0.00267 Sv 
     70 
     71  le modèle est en train d'accumuler de la neige sur les calottes continentales, à hauteur de 0.0559 Sv 
     72 
     73  -> est-on en train de mettre à l'équilibre les différents reservoirs d'eau (hors ocean) ?  
     74 
     75  [OM] accumulation d'eau en neige dans orchidee, et en glace dans land ice, limité à 3 tonnes d'eau/m2 -> on met le surplus dans l'ocean. Peut etre que la perte se stabiliserait apres equilibrage orchidee et LIC. 
     76 
     77  [YM] Masa dit qu'elle a vu ce genre de pb dans simus froides et résolues par equilibrage+ correction neige.  
     78 
     79  [OM] diagnostiques à affiner pour voir où l'eau s'accumule 
     80 
     81  + clarifier le residual dans diags eau douce, de l'ordre de 0.03 Sv ce qui depasse le seuil critique (0.003) qu'on s'etait autorisé -> verifier que les diagnostiques sont bien adaptés pour DYNAMICO  
     82 
     83  [OM] verifier les bilans d'eau douce en calculant à partir des restarts et non des time-series 
     84 
     85test avec nouveau routage, avec poids mieux pris en compte a priori : produit le meme resultat que precedemment 
    8386https://thredds-su.ipsl.fr/thredds/fileServer/tgcc_thredds/work/p86caub/INTERMONITORING/intermonit_VALIDATION_DYNAMICO_COUPLE/images/OCE_sossheig_global_ave.jpg 
    8487 
    85 Script OB pour analyses flux d'eau 
    86 -> accumulation de neige dans orchidee, land ice et sea ice probablement car pas de restart orchidee et LIC dans ces simulations tests. Ocean (et Sechiba?): limite à 3tonnes d'eau / m2. ->  on met le surplus dans l'ocean. Peut etre que la perte se stabiliserait apres equilibrage orchidee et LIC.  
    87 Masa semble dire à Yann qu'elle a vu ce genre de pb dans simus froides et résolues par equilibrage+ correction neige.  
    88 Diagnostiques à affiner pour voir où eau s'accumule  
    89 + clarifier le residual dans diags eau douce, de l'ordre de 0.03 Sv ce qui depasse le seuil critique (0.003) qu'on s'etait autorisé -> verifier que les diagnostiques sont bien adaptés pour DYNAMICO  
    90 [OM] : verifier les bilans d'eau douce en calculant à partir des restarts et non des time-series 
    91  
    9288Actions:  
    9389 
    94     laisser tourner simu avec routage v12 pour vérifier que ça se stabilise (memo JD stabilisation à -60cm dans 6ALR) 
    95  
    96     Demander validation du diagnostique à OB car schéma de routage ici est un peu différent.  
    97  
    98     [AC] Test LMDZ avec restart non equilibrés pour verifier qu'on a la même dérive ? 
    99  
    100     [JD] ecrit à OB pour attirer son attention sur les questions autour du bilan d'eau douce automatique 
    101  
    102     [JD] lance aussi action pour verifier et si beosin adapter bilans d'eau douce automatiques aux sorties NEMO v4 (+ verif mers fermées en 4.2) 
    103  
     90  laisser tourner simu avec routage v12 pour vérifier que ça se stabilise (memo JD stabilisation à -60cm dans 6ALR) 
     91 
     92  Demander validation du diagnostique à OB car schéma de routage ici est un peu différent.  
     93 
     94  [AC] Test LMDZ avec restart non equilibrés pour verifier qu'on a la même dérive ? 
     95 
     96  [JD] ecrit à OB pour attirer son attention sur les questions autour du bilan d'eau douce automatique 
     97 
     98  [JD] lance aussi action pour verifier et si beosin adapter bilans d'eau douce automatiques aux sorties NEMO v4 (NB : leur SSH derive aussi, mais se stabilise autour de -30cm au bout de 50 ans...) 
     99 
     100Prochaine reunion vendredi 28 octobre 10h30 sur NEMOv4 (VLR et LR) 
    104101 
    105102