你好,百度了一下ncl的一个相关问题,看到提到你的,想跟你请教一个问题。conform的应用,感觉没错,但是还是提示了“fatal:Number of dimensions on right hand side do not match number of dimension in left hand side”
begin
diri = "/mnt/hgfs/E/WUTIP/dat/"
fili1 = "ERA.temp.0920-1005.nc"
fili2 = "ERA.uv.0920-1005.nc"
fili3 = "ERA.rh.0920-1005.nc"
fili4 = "ERA.surf.0920-1005.nc"
diro = "/mnt/hgfs/E/WUTIP/dat/pic/"
filo = "tse_lat-time"
f1 = addfile(diri+fili1,"r")
f2 = addfile(diri+fili2,"r")
f3 = addfile(diri+fili3,"r")
f4 = addfile(diri+fili4,"r")
tk = short2flt(f1 ->t(:,{1000:200},{20.25},{110.25}))
u = short2flt(f2 ->u(:,{1000:200},{20.25},{110.25}))
v = short2flt(f2 ->v(:,{1000:200},{20.25},{110.25}))
rh = short2flt(f3 ->r(:,{1000:200},{20.25},{110.25}))
p = short2flt(f4 ->msl(:,{20.25},{110.25}))
u!0 ="time"
u!1 ="lev"
pressure_levels = (/1000.,975.,950.,925.,900.,875.,850.,825.,800.,775.,750.,700.,650.,600.,550.,500.,\
450.,400.,350.,300.,250.,225.,200./)
nlevels = dimsizes(pressure_levels) ; number of pressure levels
;计算假相当位温
p = conform(u(time|:,lev|:),pressure_levels,1)