= 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 == Historique des simulations à analyser (modèle IPSLCM6.0.13-LR_CMIP6) == * 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}}} * '''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 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 à analyser 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 }}} * Variables sorties vs variables renseignées dans le ping ("désactivée car pb" signifie que vous êtes avez été averti du pb et que la variable a été désactivée pour cette simulation) * LMDZ * 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 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 vtendnogw => désactivée car pb zmtnt => désactivée car pb }}} * NEMO * 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 cfc12 deptho => désactivée car pb difvho difvmo => désactivée car pb difvmto difvtrto dissicos fgcfc11 fgcfc12 fgsf6 fsitherm hfx => désactivée car pb hfy => désactivée car pb intpcalcite poc pocos rsdo sfdsi tomint }}} * ORCHIDEE * Variables renseignées dans le ping : 153 * Variables sorties : 150 * Variables renseignées dans le ping mais non sorties {{{ areacella nudgincsm nudgincswe }}} == 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: * pour la simulation CM6013.cmip6-pi-08f * fichiers ATM: http://www.lmd.jussieu.fr/~fairhead/Labo/CMIP6_Validation/CMIP6_validation_CM6013.cmip6-pi-08f.ATM.html * fichiers OCE: http://www.lmd.jussieu.fr/~fairhead/Labo/CMIP6_Validation/CMIP6_validation_CM6013.cmip6-pi-08f.OCE.html * fichiers SRF: http://www.lmd.jussieu.fr/~fairhead/Labo/CMIP6_Validation/CMIP6_validation_CM6013.cmip6-pi-08f.SRF.html == LMDZ == ||'''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 (=/). 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 == ||'''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) || A faire (PC : à verifier dans un run transitoire)|| || PC || cLitter || Lmon/Eyr || CM6014.cmip6-pi-01 || valeurs || cLitter_Lmon (l=1:12@ave) ≠ cLitter_Eyr || ||  || || PC || cProduct || Lmon/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 || Emon || CM6014.cmip6-pi-01 || pas de fichier sorti || A faire || A faire || -- || 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 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 == ||'''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?)