Point sur la DR CMIP6 pour LMDZ

Point sur les sorties LMDZ6 / CMIP6 atmosphère

Résumé discussion du 12/01/2018
avec Laurent F, Ionela M, Abderrahmane I, Jean-Yves G, Sandrine B

Variables renseignées dans le ping mais non sorties (https://forge.ipsl.jussieu.fr/igcmg/wiki/RetourSurLesPings) :

areacelli (vérif avec dernière version de la DR)
→ should be areacella (cf http://clipc-mfl970103/services.ceda.ac.uk/dreq/index/CMORvar.html)

cfadDbze94 (COSP)
clcalipso2 (COSP)
climodis (COSP)
clmisr (MISR SIMULATOR, COSP?)
cltmodis (COSP)
clwmodis (COSP)
jpdftaureicemodis (COSP)
jpdftaureliqmodis (COSP)
→ variables COSP qui n’étaient pas sorties car simulateurs CloudSat, MODIS et MISR n’étaient pas activés → seront activés (AI)

Le fichier namelist cosp_output_nl.txt, permet de définir les clés logiques pour activer ou non à la fois les simulateurs et la sortie des variables diagnostiques correspondants. Dans ce fichier, actuellement, ne sont activés par défaut que 2 simulateurs et la sortie de champs de fractions nuageuses correspondant.
- Ce fichier est mis à jour de façon à activer les 5 simulateurs (Calipso, Cloudsat, Isccp, Misr et Modis) et la sortie de l'ensemble des variables demandées dans CMIP6 (version svn 3155).
- Dans les sources COSP, l'activation des champs à sortir en fonction de ce qui est demandé dans les fichiers .xml est codé. Les tests sont en cours. L'équivalent quand on utilise IOIPSL n'est pas encore fait (à y réfléchir).

parasolRefl (COSP)
→ vérifier pourquoi parasolRefl n’est pas déjà sorti vu que le simulator Parasol et en principe appelé en même temps que le simulateur Calipso (AI)

Dans le fichier ping_LMDZ.xml, le nom de la variable CMIP6 est parasolRefl_sea (non parasolRefl) correspondant pour LMDZ à parasol_refl*fract_oce

cldicemxrat (cldicemxrat27)
cldwatmxrat (cldwatmxrat27)
→ variables demandées par CMIP6 pour certaines fréquences ; pas de besoin particulier exprimé par F. Chéruy (cf discussion) mais peuvent être laissées puisque variables CMIP6 ; ces variables avaient été codées pour le projet CAUSES
Après vérification, ces variables sont bien sorties

cldncl => désactivée car pb
→ réglé ?

co2mass => désactivée car pb
→ réglé ? Oui, à contrôler au prochain test

mcd => désactivée car pb
→ défini « positive down » dans CMIP6
réglé, à contrôler au prochain test

hfdsn (Downward heat flux into snow where land over land)
→ à retirer a priori

lialb (vérif avec dernière version de la DR)
lialbIs (vérif avec dernière version de la DR)
→ kesako ? (il semble que ce ne soit pas des variables CMIP6)

mrroLi ( Land Ice Runoff Flux [mon: Temporal mean, Antarctic Polar Stereographic Grid] )
→ sortie Orchidée
Variable codée dans LMDZ mais non sortie dans le test du workflow, à contrôler au prochain test

rlut4co2 (vu avec Ionela)
rlutcs4co2 (vu avec Ionela)
rsut4co2 (vu avec Ionela)
rsutcs4co2 (vu avec Ionela)
→ Calculs radiatifs 4xCO2 en cours de traitement (IM)

utendnogw => désactivée car pb
vtendnogw => désactivée car pb
tdps (2m dewpoint temperature [day: Temporal mean, Global field (single level))
→ réglé

tntogw (Zonal Mean pas encore implanté)
zmtnt => désactivée car pb
→ variables manquantes pour l’instant car moyennes zonales pas encore sorties par XIOS

Anomalies ou points à regarder dans la liste des sorties du workflow CMIP6 envoyée par Laurent le 12/01/18 :

ci : Fraction of time convection occurs
→ pourquoi le max des Amon et CF3hr est-il à 0.5 ? rapport avec l’appel 1 pas de temps sur 2 de la convection ? A vérifier [JYG, IM]
→ pb effectivement dans la sauvegarde des champs de clouds_gno qd la convection n’est pas appelée (voir Ionela et Jean-Yves pour les détails).

clc, clic, clis, clwc, clws : convective/stratiform cloud variables :
→ variables à retirer des outputs car n’ont pas de sens pour LMDZ (le rayt ne voit qu’une seule variable nuageuse qui est la combinaison conv/strat) [LF]

cldicemxrat27, cldwatmxrat27 :
→ variables demandées par CMIP6 pour certaines fréquences ; pas de besoin particulier exprimé par F. Chéruy (cf discussion) mais peuvent être laissées puisque variables CMIP6 ; ces variables avaient été codées pour le projet CAUSES

clisccp (E3hrPt) :
→ à vérifier [AI]

clisccp est un champs 5D (temps, tau, pression, lat, lon)
Mise en place d'un script python, permettant de calculer la moyenne, le min et le max. La somme est faite sur l'axe tau (épaisseurs optiques), et la moyenne, le min ou le max pour les autres axes. Le script marche en interactif pour les fichiers Cfday et Cfmon et donne des résultats raisonnables mais plante pour E3hrPt vue la taille du fichier.
Pour Cfday :
('Moyenne clisccp = ', 9.2559375856955253)
('Min clisccp = ', 0.0)
('Max clisccp = ', 62.832307272810276)
Pour Cfmon :
('Moyenne clisccp = ', 9.3181915701942941)
('Min clisccp = ', 0.0)
('Max clisccp = ', 70.282876014709473)
Pour E3hrPt : le script est lancé en bash :
('Moyenne clisccp = ', 9.1206810673905565)
('Min clisccp = ', 0.0)
('Max clisccp = ', 61.010966415352975

hurs et huss (CF3hr et Cfsubhr) : near-surface humidity
→ valeurs minimales négatives. À vérifier et corriger. [IM]

mcu, mcd : convective mass fluxes
→ vérifier que les outputs ne sont pas biaisés par l’appel 1 fois sur 2 de la convection [JYG]

prcprof, prsn, prrc, prsn, etc : Convective/stratiform rainfall/snow fluxes
→ vérifier les variables demandées et fournies [JYG & SB]

sci : Fraction of time shallow convection occurs :
→ à diagnostiquer ; prendre la fréquence d’appel des thermiques ? (devrait duffir car peu de chance que convect soit appelé seul sans thermiques dans les cas de convection peu profonde) [JYG, IM, SB]
→ envoyer mail à Steve Klein / Mark Webb pour connaître leur définition des outputs de shallow convection [SB]

Autres :

sorties COSP au dernier pas de temps manquantes :
→ à vérifier [AI]

Corrigé : rajout de l'attribut XIOS freq_offset="0ts" pour toutes les sorties Cosp.

→ plot des principales variables de tous les simulateurs [AI]

http://lmdz.lmd.jussieu.fr/Members/aidelkadi/cosp/#Validation

wind10m, u10m, v10m : diagnostique du vent à 10m
→ calcul à revoir [IM]

Flux radiatifs 4xCO2 :
→ appel à traiter [IM]

Flux radiatifs ciel-clair ou nuageux demandés que dans certaines tables:
→ à dupliquer éventuellement [SB] → voir ci-après :
→ Dupliquer les champs suivants (pour faciliter l’analyse):
rsds_day en rsds_CFday
rsus_day en rsus_CFday
rsdt_day en rsdt_CFday
rsut_day en rsut_CFday
rlus_day en rlus_CFday
rlds_day en rlds_CFday
rlut_day en rlut_CFday

Missing variables listées par Laurent (non fournies à CMIP6) (liste ici):

<field id="CMIP6_tntscp" field_ref="dummy_XYA" />  </field_group> → n’est-ce pas la variable LMDZ dtlsc+lteva aussi connue sous le nom dthlsc ?