CMIP6 模式数据处理,
从ECMWF网站上先注册和下载对应的数据
下载好了,要处理时间分辨率和空间分辨率,推荐用cdo,首先要pip install cdo
用双曲线性插值把空间分辨率统一化,先制作grid文件
grid 文件信息如下
根据grid信息双曲线性插值命令,这里假设原nc文件为pr_day_MRI-ESM2-0_ssp585_r1i1p1f1_gn_20780622-20991231.nc, 输出文件名为MRI-ESM2-0_ssp585_20780622-20991231.nc
运行
cdo remapbil,grid pr_day_MRI-ESM2-0_ssp585_r1i1p1f1_gn_20780622-20991231.nc MRI-ESM2-0_ssp585_20780622-20991231.nc
随后调整时间为标准时刻,这里的模式设置每年都是365天,在python读取很不方面。
cdo setcalendar,standard MRI-ESM2-0_ssp585_20780622-20991231.nc MRI-ESM2-0_ssp585_20780622-20991231_Standard.nc
这样调整好了,用python再读取每个数据的全球空间范围的降雨数值,空间分辨率1X1°,由于单位不是mm/天,把每个数值乘以86400,得到每一天的数值。
记录一下,处理数据是个体力活。