用Grads求青藏高原区域平均温度报错，求助

bzz

http://bbs.06climate.com/forum.php?mod=attachment&aid=MTA4Njk4fDhlMmY2N2ZkNjVjZDBmZmZhMzFhMDgwNTg3NzMzMzdhfDE3MzIzMzYwMTg%3D&request=yes&_f=.nc我想求出青藏高原两米温度的面积平均值
当我将青藏高原温度作为第二个文件打开时，求区域平均报错
但是单独打开计算时却没有问题
请问各位大神，这个该如何解决？？或者说有没有其他办法？？
谢谢

附青藏高原温度nc文件和图，这个nc使用python重新写的
代码如下下：






import numpy as np
import xarray as xr

data1=xr.open_dataset(r'D:\t2m_1940_2022_year.nc')
tas01=data1.t2m
tas01=tas01[:,::-1,:]

#提取区域数据
a=np.where(tas01.longitude==66)[0][0]
b=np.where(tas01.longitude==106)[0][0]
tas01=tas01[:,:,a:b+1]
a=np.where(tas01.latitude==25)[0][0]
b=np.where(tas01.latitude==42)[0][0]
tas01=tas01[:,a:b+1,:]

#提取青藏高原数据
import geopandas as gpd
shp = r"D:\ship\TP\TPBoundary_new(2021)\TPBoundary_new(2021).shp"
tp_shp=gpd.read_file(shp)
tas01_tp=tas01.salem.roi(shape=tp_shp)
tas01_tp.to_netcdf(r'D:\t2m_TP_1940_2022_year.nc')