wiki:RetourSurLesPings

Version 145 (modified by jgipsl, 6 years ago) (diff)

--

Validation Workflow CMIP6

Cette page permet de :

  • centraliser les informations relatives à la validation des sorties produites par le Workflow CMIP6 IPSL.
  • renseigner les erreurs/problèmes trouvés
  • recenser les développements/fonctionalités non finalisées
  • assurer un suivi de ces erreurs jusqu'à leur correction

Chemins d'accès des simulations à analyser (modèle IPSLCM6.0.14-LR_CMIP6)

Version précédent IPSLCM6.0.13-LR_CMIP6 => obsolète :

  • 29/11/2017 : Simulation CM6013.cmip6-pi-08 de 10 ans (A NE PLUS UTILISER !!!)
  • 07/12/2017 : Simulation CM6013.cmip6-pi-08f de 1 an
    • DR 00.15, DR2XML 0.23, XIOS2 rev 1309
    • Les fichiers ping*xml et dr2xml*xml utilisés sont là :
      • Curie : /ccc/scratch/cont003/gencmip6/p86caub/IGCM_OUT/IPSLCM6/DEVT/piControl/CM6013.cmip6-pi-08f/CMIP6
      • Ciclad : https://vesg.ipsl.upmc.fr/thredds/catalog/work_thredds/p86caub/IPSLCM6/DEVT/piControl/CM6013.cmip6-pi-08f/CMIP6/catalog.html
    • Les sorties à analyser sont là :
      • Curie : /ccc/scratch/cont003/gencmip6/p86caub/IGCM_OUT/IPSLCM6/DEVT/piControl/CM6013.cmip6-pi-08f/CMIP6
      • Ciclad : https://vesg.ipsl.upmc.fr/thredds/catalog/work_thredds/p86caub/IPSLCM6/DEVT/piControl/CM6013.cmip6-pi-08f/CMIP6/catalog.html

Version actuelle :

  • 22/12/2017 : Simulation CM6014.cmip6-pi-01 de 1 an
    • IPSLCM6.0.14-LR, DR 00.15, DR2XML 0.23, XIOS2 rev 1387
    • Les sorties à analyser et les fichiers ping*xml et dr2xml*xml utilisés sont là :
      • Curie : /ccc/scratch/cont003/gencmip6/p86caub/IGCM_OUT/IPSLCM6/DEVT/piControl/CM6014.cmip6-pi-01/CMIP6
      • Ciclad : /ccc/work/cont003/thredds/p86denv/CMIP6/
      • vesg :  https://vesg.ipsl.upmc.fr/thredds/catalog/work_thredds/p86denv/CMIP6/catalog.html
    • Les sorties équivalentes en mode "standard" (fichiers hist*, grid_T, sechiba_history,...) sont là :
      • Curie : /ccc/store/cont003/gencmip6/p86caub/IGCM_OUT/IPSLCM6/DEVT/piControl/CM6014.standard-pi-01
      • Ciclad : /ccc/store/cont003/thredds/p86caub/IPSLCM6/DEVT/piControl/CM6014.standard-pi-01

Outils d'aide à la validation

LMDZ: Multi-atlas réalisé à partir de toutes les variables Amon disponibles dans la dernière simulation: https://vesgint-data.ipsl.upmc.fr/thredds/fileServer/IPSLFS/jservon/C-ESM-EP/DR_CMIP6_all_ping_variables_jservon/Atmosphere/atlas_Atmosphere_DR_CMIP6_all_ping_variables.html

ORCHIDEE: Multi-atlas réalisé à partir de toutes les variables Lmon disponibles dans la dernière simulation: https://vesgint-data.ipsl.upmc.fr/thredds/fileServer/IPSLFS/jservon/C-ESM-EP/DR_CMIP6_all_ping_variables_jservon/LandSurfaces/atlas_LandSurfaces_DR_CMIP6_all_ping_variables.html

Multi-atlas C-ESM-EP simplifié (atlas d'évaluation classique): https://vesg.ipsl.upmc.fr/thredds/fileServer/IPSLFS/jservon/C-ESM-EP/DR_CMIP6_check_quick_atlas_jservon/C-ESM-EP_DR_CMIP6_check_quick_atlas.html

Résultats d'un petit script utilisant CDO pour sortir des stats simples sur tous les fichiers d'un répertoire:

LMDZ

AC : Variables sorties vs variables renseignées dans le ping ("désactivée car pb" signifie que vous avez été averti du pb et que la variable a été désactivée pour cette simulation)

  • Variables renseignées dans le ping : 232
  • Variables sorties : 204
  • Variables renseignées dans le ping mais non sorties
    areacelli (vérif avec dernière version de la DR)
    cfadDbze94 (COSP)
    clcalipso2  (COSP)
    cldicemxrat  (cldicemxrat27)
    cldncl  => désactivée car pb
    cldwatmxrat  (cldwatmxrat27)
    climodis (COSP)
    clmisr  (MISR SIMULATOR, COSP?)
    cltmodis (COSP)
    clwmodis (COSP)
    co2mass => désactivée car pb (a priori réglé avec LMDZ r3145)
    hfdsn (Downward heat flux into snow where land over land)
    jpdftaureicemodis (COSP)
    jpdftaureliqmodis (COSP)
    lialb (vérif avec dernière version de la DR)
    lialbIs (vérif avec dernière version de la DR)
    mcd => désactivée car pb
    mrroLi ( Land Ice Runoff Flux [mon: Temporal mean, Antarctic Polar Stereographic Grid] )
    parasolRefl (COSP)
    rlut4co2 (vu avec Ionela)
    rlutcs4co2 (vu avec Ionela)
    rsut4co2 (vu avec Ionela)
    rsutcs4co2 (vu avec Ionela)
    tdps (2m dewpoint temperature [day: Temporal mean, Global field (single level))
    tntogw (Zonal Mean pas encore implanté)
    utendnogw => désactivée car pb (a priori réglé avec LMDZ r3139)
    vtendnogw => désactivée car pb (a priori réglé avec LMDZ r3139)
    zmtnt => désactivée car pb
    
Rapporteur Variable Fréquence Simulation Problème/Tâche Outil responsable Statut Simulation avec correction
OB cdnc AERmon CM6013.cmip6-pi-08 Valeurs trop petites ping file changer pour scdnc/lcc3d où on moyenne les deux champs avant d'en faire le rapport. A FAIRE
OB ch4globalClim, n2oglobalClim Amon CM6013.cmip6-pi-08 redondant avec ch4global et n2oglobal ping file devrait disparaître avec DRv1.00.18. Toujours là dans 6.0.14 ?'''
OB cl, clc Amon CM6013.cmip6-pi-08 0-1 au lieu de 0-100% Corrigé ping file
OB cldnvi Eday, Emon CM6013.cmip6-pi-08 Valeurs qui me semblent bizarres. Diag LMDz à discuter. LMDz Je pense que cldnvi doit être divisé par le diag lcc après moyennage (=<cldnvi>/<lcc>). Implémenté dans 6.0.14 mais trop de valeurs indéfinies ==> pb
AC hfls, hfss AMon CM6013.cmip6-pi-08 Erreur de signe ping file Corrigé (commit LMDZ rev 3096) CM6013.cmip6-pi-08f
OB hus, huss, hus7h CM6013.cmip6-pi-08 L'unité du fichier ping (1) semble erronée, kg kg-1 plus approprié CDR Github issue ouvert. Martin Juckes l'a fermé. Sortie LMDz corrigé. On laisse comme cela pour CMIP6 même si pas terrible.
OB sftlf CM6013.cmip6-pi-08 pourc_ter ou pourc_ter+pourc_lic ?? Masa me confirme que c'est plutôt pourc_ter+pourc_lic ping Corrigé ping file Validé dans 6.0.14
AC ua,ta,va,... AMon CM6013.cmip6-pi-08 Pb interpolation verticale niveaux de pression DR2XML + XIOS Corrigé (correction DR2XML) CM6013.cmip6-pi-08f
OB ztp AERmon CM6013.cmip6-pi-08 valeurs trop faibles LMDz Corrigé LMDz rev 3119 puis dans rev 3123
SB clt 3hr CM6013.cmip6-pi-08 Cette valeur paraît tres faible pour une total cloud fraction. Les valeurs de la clt pour les Amon et day paraissent plus raisonnables. Corrige pb unité
SB hur Amon CM6013.cmip6-pi-08 Le max suggere que l'humide relative varie ici entre 0 et 1 et non entre 0 et 100%. Idem pour tous les champs hur suivants (mais pas les hurs) Corrigé ping unité
SB rlu CFmon CM6013.cmip6-pi-08 ce champ (rlu) semble identique au champ ciel clair (rlucs) ??
SB rlucs CFmon CM6013.cmip6-pi-08 ce champ (rlucs) semble identique au champ all-sky (rlu)/ Les min/max ressemblent aussi a ceux du champ rlus, mais c'est peut-etre normal si les min/max se trouvent a la surface. ??
SB rlutcs CFday CM6013.cmip6-pi-08 la valeur moyenne de l'OLR ciel clair ne devrait pas etre identique a la valeur moyenne de l'OLR all-sky (champ suivant). ??
SB rsdt Amon CM6013.cmip6-pi-08 Ici j'imagine que l'on donne toujours l'insolation par une valeur positive, meme si c'est un flux downwelling? ??
SB tas 3hr CM6013.cmip6-pi-08 minimum: C'est froid pour une near-surface air temperature. C'est en Antarctique? ??
SB ts Amon CM6013.cmip6-pi-08 minimum froid ??
SB wap Amon CM6013.cmip6-pi-08 max: Ca fait 648 hPa/day...c'est fort pour une subsidence a 500 hPa! Est-ce que c'est dû à une interpolation sur un relief? ??
JS ta, hus AMon CM6013.cmip6-pi-08f Pas de donnée ?? Corrigé dans CM6014.cmip6-pi-01
JS hur AMon CM6013.cmip6-pi-08f unité en fraction (0,1) mais demandé en % dans le ping Corrigé ping file
AC zg AMon CM6013.cmip6-pi-08f Pb vu avec multi-atlas cmip6 vs standard ping file en cours
OB rld4co2, rlucs4co2, ... CM6014 0 et bandes Ne peut pas fonctionner si co2_ppm_per n est pas défini ==> pb de config ?
JS sci AMon CM6014.cmip6-pi-1 Pas de donnée ?? En cours

Variables regardées:

  • JS: toutes les variables Amon disponibles dans la dernière simulation CM6014.cmip6-pi-01:

https://vesgint-data.ipsl.upmc.fr/thredds/fileServer/IPSLFS/jservon/C-ESM-EP/DR_CMIP6_all_ping_variables_jservon/Atmosphere/atlas_Atmosphere_DR_CMIP6_all_ping_variables.html

  • OB: solbnd, ztp, tatp, ptp, sza, sftlf, sconcss, sconcso4, sconcdust, hur, ec550aer, ch4global, ch4globalclim, n2oglobal, n2oglobalclim, cfc11global, cfc12global, abs550aer, aod550volso4, cldnvi (pb), conc*, load*, od*, lwtoacsaer, sconc* disponibles dans le run 6.0.14

Manque tous les flux 4xco2, il faudrait modifier le code pour activer le diagnostic LMDz automatiquement si le diag est demandé mais à ne pas demander systématiquement car coûteux.

Remarque de Sandrine: Il y a un rlutcs_CFday et un rlut_day mais pas de rlut_CFday ni de rlutcs_day. C'est sans doute que pour les champs TOA, les champs CFday et day sont identiques. Mais pour l'analyse, ca va pas mal compliquer les choses. Serait-il possible de rajouter (ou dupliquer) les champs pour que l'on ait le ciel-clair et le all-sky pour chaque table?

ORCHIDEE

  • AC : Variables sorties vs variables renseignées dans le ping ("désactivée car pb" signifie que vous avez été averti du pb et que la variable a été désactivée pour cette simulation)
    • Variables renseignées dans le ping : 153
    • Variables sorties : 150
    • Variables renseignées dans le ping mais non sorties
      areacella (sorti par ATM)
      nudgincsm (Nudging Increment of Water in Soil Mositure)   JG : c'est normal. La variable est désactivé car l'option ok_nudge_mc=false.
      nudgincswe (Nudging Increment of Water in Snow)           JG : c'est normal. La variable est désactivé car l'option ok_nudge_snow=false.
      
Rapporteur Variable Fréquence Simulation Problème/Tâche Outil responsable Statut Simulation avec correction
JS cProduct, fFire Lmon CM6013.cmip6-pi-08f Champ vide PC : ok (run de contrôle donc pas de land use change et feux pas activés) A faire (PC : à verifier dans un run transitoire) A faire (PC : à verifier dans un run transitoire)
JS fFire Lmon CM6013.cmip6-pi-08f Champ vide PC : ok (run de contrôle donc pas de land use change et feux pas activés) A faire (PC : à verifier dans un run transitoire) JG : fFire n'est plus donné car les feux sont pas activé (dummy_not_provided)
PC cLitter Lmon et Eyr CM6014.cmip6-pi-01 valeurs cLitter_Lmon (l=1:12@ave) ≠ cLitter_Eyr  
PC cProduct Lmon et Eyr CM6014.cmip6-pi-01 valeurs/definition cette variable doit être nulle dans une simulation de contrôle - pas de land use change (et donc pas de wood harvest non plus), ni de land cover change) vérifier les flags dans la config
JS mrtws Emon CM6014.cmip6-pi-01 warning: facteur 5 avec la version précédente A confirmer ?? A confirmer --
JS hflsLut, hfssLut, fFire (PC : feux pas activés dans CMIP6) Emon CM6014.cmip6-pi-01 pas de fichier sorti A faire A faire --
PC cSoil Emon et Eyr CM6014.cmip6-pi-01 valeurs cSoil_Emon (l=1:12@ave) ≠ cSoil_Eyr  
PC cVeg Lmon et Eyr CM6014.cmip6-pi-01 valeurs cVeg_Lmon (l=1:12@ave) ≠ cVeg_Eyr  
PC ra E3hr CM6014.cmip6-pi-01 frequences de sortie sorties journalières au lieu des sorties demandées toutes les 3 heures  
AD sandfrac Efx CM6014.cmip6-pi-01 mauvaises valeurs : 100 partout alors que ça doit valoir entre 0 et 1 orchidee (slowproc.f90) Pb identifié, lié au restart (AD) Commit fait rev 4871. Il faut ajouter soils_param.nc si on démarre avec un vieux restart
PC rhSoil Emon CM6014.cmip6-pi-01 mauvaises valeurs  
PC rhLitter Emon CM6014.cmip6-pi-01 mauvaises valeurs  

Variables regardées:

  • JS (mise à jour le 09-01-2018): toutes les variables *mon disponibles dans la dernière simulation CM6014.cmip6-pi-01:

https://vesgint-data.ipsl.upmc.fr/thredds/fileServer/IPSLFS/jservon/C-ESM-EP/DR_CMIP6_all_ping_variables_jservon/LandSurfaces/atlas_LandSurfaces_DR_CMIP6_all_ping_variables.html Note: cProduct semble ok, il faut confirmer si c'est normal que l'on ait pas le fichier fFire en sortie

  • PC (CM6014.cmip6-pi-01) : baresoilFrac (Lmon, Eyr), cLitterAbove, cLitterBelow, cLitter (Lmon, Eyr), cLitterGrass, cLitterTree, cLitterSubSurf, cLitterSurf, cProduct (Lmon, Eyr), cSoil (Emon, Eyr), cSoilGrass, cSoilTree, cSoilFast, cSoilMedium, cSoilSlow, cVeg (Lmon, Eyr), cVegGrass, cVegTree, gpp (E3hr, Lmon), gppGrass, gppTree, ra (E3hr, Lmon), raGrass, raTree, rh (E3hr, Lmon), rhGrass, rhLitter, rhSoil, rhTree

remarque : créer éventuellement une variable cLitterAgricultural pour vérifier que le bilan est bien fermé. Validation visuelle effectuée avec le masque crop (cropFracC3+cropFracC4).

NEMO

AC : Variables sorties vs variables renseignées dans le ping ("désactivée car pb" signifie que vous avez été averti du pb et que la variable a été désactivée pour cette simulation)

  • Variables renseignées dans le ping : 253
  • Variables sorties : 232
  • Variables renseignées dans le ping mais non sorties
    bigthetaoga => désactivée car pb
    cfc11 (Moles Per Unit Mass of CFC-11 in sea water)
    cfc12 (Moles Per Unit Mass of CFC-12 in sea water)
    deptho => désactivée car pb
    difvho (OClim présent, manque le Oyr)
    difvmo => désactivée car pb
    difvmto (Oclim manquant)
    difvtrto (Oclim manquant)
    dissicos (vérif avec dernière version de la DR)
    fgcfc11 (Surface Downward CFC11 flux [mon: Temporal mean, Global field (single level))
    fgcfc12 (Surface Downward CFC12 flux [mon: Temporal mean, Global field (single level))
    fgsf6 ( Surface Downward SF6 flux)
    fsitherm (Water Flux into Sea Water due to Sea Ice Thermodynamics)
    hfx => désactivée car pb
    hfy => désactivée car pb
    intpcalcite ( Calcite Production )
    poc (vérif avec dernière version de la DR. Devenu?intpoc {Omon}: Particulate Organic Carbon Content)
    pocos (vérif avec dernière version de la DR)
    rsdo (Downwelling Shortwave Radiation in Sea Water)
    sfdsi (attention désigne deux "choses". sfdsi {SImon}: Salt flux from sea ice [mon: Temporal mean, Global field (single level) [XY-na] [amnsia-twm]] (1)
    sfdsi {Omon}: Downward Sea Ice Basal Salt Flux [mon: Temporal mean, Global field (single level) [XY-na] [amse-tmn]] (8))
    tomint (integral wrt depth of product of sea water density and prognostic temperature)
    
Rapporteur Variable Fréquence Simulation Problème/Tâche Outil responsable Statut Simulation avec correction
AC Variables groupe SBC AMon CM6013.cmip6-pi-08 11 records par an (au lieu de 12) DR2XML + XIOS Corrigé (correction DR2XML) CM6013.cmip6-pi-08f okJM
AC Variables PISCES manquantes ?? CM6013.cmip6-pi-08 Pas de sorties PISCES DR2XML Corrigé (correction setting DR2XML) CM6013.cmip6-pi-08f
AC msftbarot OMon CM6013.cmip6-pi-08f Pb avec somme cumulée suivant axe j XIOS + NEMO En cours
AC ?? ?? CM6013.cmip6-pi-08f Transects XIOS + NEMO En cours
JM+JD *ga Omon CM6013.cmip6-pi-08 valeurs a 0 DR2XML corrigé (correction DR2XML) CM6014.cmip6-pi-01 okJM
JM+JD hcont300 Omon CM6013.cmip6-pi-08 unités NEMO corrigé (code): CM6013.cmip6-pi-08f ok en UNITES mais code NEMO diahth.F90 pas pris en compte: a recompiler. Devrait etre ok dans CM6014 CM6014.cmip6-pi-01 okJM
JM+JD salinity fields Omon CM6013.cmip6-pi-08 valeurs * 35 XML corrigé (XML) ok CM6013.cmip6-pi-08f
JM thetaot300 thetato700 thetaot2000 Omon CM6013.cmip6-pi-08 no valid data XIOS Corrigé CM6014.cmip6-pi-01 okJM
JM msftyyz + hfbasin etc Omon CM6013.cmip6-pi-08 ok but only 3 basins needed rather than 5 for the moment in the output (1:glb, 2:atl, 5:indoPac) + information on basin numbers to be added somewhere in metadata NEMO en cours Julie
JM hfbasin Omon CM6013.cmip6-pi-08f values = 0. Par contre hfbasinpmadv ok. sltbasin a des valeurs NEMO? en cours Julie...
JM sltovgyre et sltbasins hfbasinpmadv (m>1) Omon CM6013.cmip6-pi-08f suspicieux: mauvais masque de bassin. htovgyre probablement faux aussi mais ne se voit pas a cause de hfbasin. NEMO en cours Julie: diaptr.F90: btmsk30 a rajouter dans calcul des composantes totales. update 21 Dec: ok diaptr.F90 chez JD
JM sltnortha CM6013.cmip6-pi-08f lon x lat x time devrait etre lat x time seulement? XIOS corrigé. Devrait être pris en compte dans CM6014 correct dimensions in CM6014.cmip6-pi-01 but = 0
JM tos Omon Oday 3hr CM6013.cmip6-pi-08 given in °C but indicated unit is K. After emails exchanges with Martin Jukes, he indicated that the DR will eventually be changes with tos to be given in °C. So we leave the field as it is given but units will have to be changed in agreement with the evolution of the DR
JM masso,volo Omon CM6013.cmip6-pi-08f lon x lat x time, should depend on time only. XIOS corrigé. Devrait être pris en compte dans CM6014 CM6014.cmip6-pi-01 okJM
JM *_Oclim Oclim CM6014.cmip6-pi-01 Oclim = 365 pas de temps. ok? to be checked by XIOS
JM dispkevfo Oclim CM6014.cmip6-pi-01 repartition bizarre de valeurs selon latitude et longitude NEMO MPI issue ? to be checked by JM&JD
JM sftof Ofx CM6014.cmip6-pi-01 suspicious values 0 in ocean, 100% under sea ice Ping corrigé rev 200 de shaconemo: 100-this -> Arnaud: verifier syntaxe stp?

Variables regardées:

  • JS: variables océan Omon disponibles dans l'atlas C-ESM-EP
  • JM: areacello, tauuo tauvo thetao, variables SBC, uo, vo, umo, vmo, wo, wmo, obvfsq, pbo, mlotst_Eday, mlostmax, mlostmin, omldamax, rsntds, hfevapds, hfds, hfgeou, hfrainds, masscello, hfx,hfy, sltbasin, htovgyre, htvovort, rsdo, rsntds, so, sos, somint, sob, sossq, t20d, thkcello, thetao, tob, tos, tossq, thetaot,
  • JM CM6014.cmip6-pi-01 : basin, diftrblo2d (ordre de grandeur?), diftreblo2d (=diftrblo2d), diftrxylo2d (constant dans le temps et constante scalée par la geometrie), difvho (noisy because evd activation at some places and depths), difvso, dispkevfo, dispkexyfo (ordre de grandeur?), ficeberg, flandice, friver, hfibthermds2d, hfrunoffds2d, hfsnthermds2d, masso, mlost,
  • JM CM6014.cmip6-pi-01 : ocontempdif (ordre de grandeur?), ocontemppadvect, ocontemppmdiff, ocontemprmadvect, ocontemptend (remarque ordre de grandeur: 12 W/m2 = 0.7 K/mois -> ok?),
  • JM CM6014.cmip6-pi-01 : osaltdiff (strong values in Hudson bay), osaltpadvect, osaltpmdiff, osaltrmadvect, osalttend (4e-3 kg / m2 / s = 1 psu / mois -> ok?)