由于数据是批量处理,后面是纯数组,不附带任何经纬度信息,需要利用一个函数来将陆地部分的值处理成空值(nan值)
一.首先定义一个二维的mask函数
,将纯数组的第一维
## 当纯数组的维度为(lat,lon)时,将纯数组的第一维赋值纬度信息,第二维赋值经度信息
注:此时这里的分辨率为1°×1°,所以间隔为1,纬度是90°N-90°S,当用CDO中的invertlat函数倒转纬度以后,所有的数据是从北纬到南纬
二.调用定义好的mask函数
MD1是要处理的目标二维数组,-90, 90, 0, 359这四个参数是纬度和经度范围参数,要对全球的陆地处理的话就是这个范围
三.目标数组是三维的情况
CMIP6 数据常常涉及到批量处理,需要对多个模式进行处理,所以可以定义一个三维函数来对多个模式的陆地进行处理
r1为目标数组,第一维是模式数量,然后是纬度、经度,将所有模式的陆地全部处理成空值