Missing cell_method in the v2022 version of the SAFRAN forcing

[aducharne]

The SAFRAN forcing has been revised in summer 2022 see
​http://forge.ipsl.jussieu.fr/orchidee/ticket/854

But there is no cell_methods attribute (needed by the new driver) in the files of v2022, while they were present in the v2021 version.

The mean results with the new driver seem OK (consistent with the changes between v2021 and v2022) but we haven't checked the diurnal cycle.

[aducharne]

After some exchanges with the SAFRAN producers at METEOFRANCE, the good cell_methods are :
Tair: cell_methods = "time: instantaneous"
Idem pour Qair et PSurf et Wind (en approximant une moyenne sur 10 minutes par une valeur instantanée)
Rainf: cell_methods = "time: mean(end)"
Idem pour Snowf, SWdown, LWdown (en ne faisant pas d'interpolation linéaire pour LWdown afin de conserver l’énergie, il y a un keyword pour ça)

This is based on the info by Pierre Etchevers on Wed, 28 Apr 2021
1) ce sont des heures UTC
2) les variables scalaires sont instantannées (sauf le vent : moyenne entre H-10 mn et H).
3) pour les flux, c'est la moyenne horaire, avant time (donc de H-1 à H).

Here are the nco command prepared by Lucia to add the attributes, in case of one specific year (2010) :

ncatted -O -a cell_methods,LWdown,o,c,"time: mean(end)" -a cell_methods,SWdown,o,c,"time: mean(end)" ORCHIDEE_France-8km-1hour_fromSAFRAN_2010.nc

ncatted -O -a cell_methods,Rainf,a,c,"time: mean(end)" -a cell_methods,Snowf,a,c,"time: mean(end)" ORCHIDEE_France-8km-1hour_fromSAFRAN_2010.nc

ncatted -O -a cell_methods,Tair,a,c,"time: instantaneous" -a cell_methods,Qair,a,c,"time: instantaneous" -a cell_methods,PSurf,a,c,"time: instantaneous" -a cell_methods,Wind,a,c,"time: instantaneous" ORCHIDEE_France-8km-1hour_fromSAFRAN_2010.nc

[aducharne]

Il y a eu bcp de travail sur ce sujet début 2023 avec Vlad, Nicolas Viovy et Peng Huang, cf ​http://forge.ipsl.jussieu.fr/orchidee/attachment/ticket/889/Axe%20des%20temps%20SAFRAN.pdf. Comme illustré sur ce document, on a des instabilités numériques sur l'ET simulée si on utilise mean(end) pour les flux, alors que ce n'est pas le cas avec mean(center), le défaut du nouveau driver.

On a donc décidé de garder mean(center) bien que ce ne soit pas cohérent avec la nature du forçage selon Météo France (l'interprétation mean(center) ayant bien été reconfirmée par Météo France).