而对于不规则的范围,数据的匹配精度有一定要求,譬如,需要严格按照某个特定区域的shapefile文件来截取数据。...虽然,NCL官网提供了可行的解决方案,但是 shapefile_mask_data(包含在shapefile_utils.ncl中,官网有提供)也仅仅是较好地适用于2维的Lat-Lon数据,对于3维或者更高维度的数据...)
34 f_eccc = addfile(filp1+"ECCC_TP_hr.nc","r")
35 f_ecmwf = addfile(filp1+"ECMWF_TP_hr.nc","r"...)
36 f_hmcr = addfile(filp1+"HMCR_TP_hr.nc","r")
37 f_ukmo = addfile(filp1+"UKMO_TP_hr.nc","r")...0,0,:,:),shp_filename,True) ; 构建母版——最基础的2D mask范围
2 mask_io = where(ismissing(mask_in), 0, 1) ;将所需范围内外数据的分离开