我有一组数据U，其尺寸为[ns，ny，nx]，单位为sigma级别，并且希望将其插入另一个标准压力水平（[nz，ny，nx]）这样我必须按照以下网格处理这个问题：
U2 = np.empty((nz, ny, nx))
for ix in range(NX):
for iy in range(NY):
tmp = np.interp(plev_new, plev_old[:,iy,ix] U[:,iy,ix], left=filled_value, right=filled_value)
U2[:,iy,ix] = tmp
但是，我发现按网格处理数据网格效率非常低。 我认为更快的方法是在x和y网格中并行。但我总是在模块多处理中没有这样做。
任何有关多处理的建议或其他更好或更有效的方法。
非常感谢