在 NetCDF 文件 - 编程之家

如何解决在 NetCDF 文件

我有一个非常大的 NetCDF 文件，我想在域的某个部分设置 density、speed 的 sound 和 temperature 值，比方说从 x=-120 到 x=0。我没有以前的插值解决方案。你能帮助我吗？我现在正在尝试对解决方案进行插值，其中我希望仅对坐标 x

solutionfield = Field('data_for_interp/restart.pval')

dummygrid = Grid('mesh/dome_lox.grid')

injector_sol = Field("out/dome_lox.pval.1")

interpolated_field = solutionfield.interpolate(solutiongrid,dummygrid,\
                       nneighbours=10,method='inverse_distance',\
                       exponent=3.00)

interpolated_field.reset_iteration_counter()

injector_ids = np.where(dummygrid.pnt[0,:] < 0.0)[0]

for var in interpolated_field.vars.keys():
    interpolated_field.vars[var][injector_ids] = injector_sol.vars[var][injector_ids]

interpolated_field.write('restart_dome.pval')

...