У меня есть набор данных U, размер которого [ns, ny, nx] в сигма-уровне с единицами давления, и я хотел бы объединить его в другой стандартный уровень давления ([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
Однако я нашел, что это очень неэффективно при работе с сеткой данных по сетке.Я думаю, что более быстрый способ сделать это параллельно в х и у сетки.Но мне всегда не удавалось сделать это в модуле многопроцессорной обработки.
любое предложение о многопроцессорной обработке или какой-то другой способ, который более эффективен или эффективен.
Большое спасибо