这与早期的 thread 密切相关，但我在用3d数组再现解决方案时遇到了一些困难。我的3d阵列有一个时间坐标和两个空间坐标。每个空间坐标本身就是沿空间维的指数的二维函数，就像当空间网格是曲线时的(经度、纬度)值一样。

目标:每次 ，我想找出坐标对应于空间中的最大值。以下是一些代码：

# Create cftime object to mimic my real data.
yr1=250; yr2=250
times = xr.cftime_range(start=format(yr1,'04d') + '-01',
                                        end  =format(yr2,'04d') + '-02',
                                        freq='1MS',calendar='noleap')
x=np.arange(6).reshape(2,3); y = x;
nt=len(times); ny=x.shape[0]; nx=x.shape[1]