使用ncl脚本时出现警告:day_of_year: illegal calendar = 'proleptic_gregorian'
其原因是读取的降水nc文件是我手动合并生成,所以时间的calendar不是很标准,数据信息如下所示,可以发现Calendar
是proleptic_gregorian
,这个日历不被ncl脚本所识别。
在我检查cmip6模式资料时,发现他们的数据日历类型如下所示,calendar为365_day
本来是打算将calendar统一为模式资料的,但是感觉有点麻烦,而且他们的reftime都不一样。后来,我又看了一下网上下载的olr.daily资料,发现他们的格式是标准格式的,calendar=standard,这个看起来不错。
接着又对于我合成的降水资料的原始数据进行了检查,发现原本降水资料的calendar的类型就是gregorian,就是有点奇怪哈,所以合并之后才有点问题
为了解决上面的问题,我决定将合并的文件统一成olr中的标准格式。在这里使用cdo命令可以很方面的实现以上需求
cdo setcalendar,standard GPCP_data_interp_2.5x2.5.nc gpcp2.nc
运行之后,生成的文件信息如下所示,可以发现已经修改过来了。再次运行ncl脚本也不会有警告啦
参考
https://yidongwonyi.wordpress.com/linux-data-handling-netcdf-nc/cdo-change-calendar-type/
https://code.mpimet.mpg.de/boards/1/topics/3613
https://stackoverflow.com/questions/75771997/creating-an-ensemble-with-different-calendar-types
https://earth.bsc.es/gitlab/ces/cdo/-/commit/78b29b4ea1a8651863b20baddbb5bea85efd2a05