shapefile_mask_data出错

synt

mask作图出错，shapefile_mask_data: Error: not a valid rectilinear, curvilinear, or unstructured grid，有大佬帮我看下哪里有错误吗
数据处理部分
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "d:/dataPRE/shapefile_utils.ncl"
begin


;------------------------读取数据------------------------
f11 = addfile("D:/dataPRE/PRE_MON_1961-2020.nc","r")
pre1 = f11->PRE_MON         ;(time=720, lat=72, lon=128)

;从数据上看这是1961-2020的月平均数据
;所以2018年的6-8月的数据坐标位置为：
;201806:(2018-1961)*12+6 = 690
;201808:(2018-1961)*12+8 = 692


time = ispan(201806,201808,1)

;pre_jja11 = dim_avg_n(pre1(690:692,:,:),0) ;夏季降水
pre_avg = new((/1/),"float")


          pre_mask   = shapefile_mask_data(pre1(690:692,:,:),"d:/dataPRE/shp/shp/tibetan.shp",True)
          pre_avg = dim_avg_n(pre_mask(690:692,:,:),0);区域平均

pre_jja11 = dim_rmvmean_Wrap(pre_avg)

AXD

这个函数似乎不能有时间维，把时间维循环起来就好了，不要直接放进去

@me

请问解决了吗？求个完整程序呢shapefile_utils.ncl文件是官网下载后直接用的么？