NCL处理wrfout散度

Bozhero

我在用NCL处理wrfout时，想画散度图，但是ncl散度函数不能用，自己想做中央差分计算，但是最后的图有问题，想请教一下大家有没有计算wrfout散度更好的方法
我的代码：
    u  = wrf_user_getvar(a,"ua",it)        
    v  = wrf_user_getvar(a,"va",it)   

     u_plane  = wrf_user_intrp3d( u,p,"h",pressure,0.,False)
     v_plane  = wrf_user_intrp3d( v,p,"h",pressure,0.,False)

      spd     = (u_plane*u_plane + v_plane*v_plane)^(0.5) ; m/sec
      spd@description = "Wind Speed"
      spd@units = "m/s"
      u_plane = u_plane*1.94386     ; kts
      v_plane = v_plane*1.94386     ; kts
      u_plane@units = "kts"
      v_plane@units = "kts"
          
;计算散度

          part1 = center_finite_diff_n(u_plane,lat,False,0,0)
          part2 = center_finite_diff_n(v_plane,lon,False,0,0)
          div = part1 + part2
          
          plot = gsn_csm_contour_map(wks,div,mpres)

放逐之刃

latm2d=a->XLAT(0,:,:)
lonm2d=a->XLONG(0,:,:)
fstart=min(latm2d)
fend=max(latm2d)
dims=dimsizes(latm2d)
lat = fspan(fstart, fend, dims(0))
fstart=min(lonm2d)
fend=max(lonm2d)
lon = fspan(fstart, fend, dims(1))

u_grid = rcm2rgrid(latm2d,lonm2d,u1_intrp,lat,lon,0)
u_grid!0="lev"
u_grid!1="lat"
u_grid!2="lon"
u_grid&lat=lat
u_grid&lon=lon
v_grid = rcm2rgrid(latm2d,lonm2d,v1_intrp,lat,lon,0)
v_grid!0="lev"
v_grid!1="lat"
v_grid!2="lon"
v_grid&lat=lat
v_grid&lon=lon
div_grid = uv2dv_cfd(u_grid, v_grid, lat, lon, 1)
dv = rgrid2rcm (lat, lon, div_grid, latm2d, lonm2d, 0)

Bozhero

自己顶一下

气象草鸟

请问您解决了吗，可以教我一下吗

Bozhero

最后使用了NCL里计算散度的函数解决的

aycc90

可以贴出来看一下吗

六间小可爱

您好～请问用哪个函数可以计算wrfout的散度呢？

好像一条dog

放逐之刃 发表于 2020-1-10 22:22
latm2d=a->XLAT(0,:,:)
lonm2d=a->XLONG(0,:,:)
fstart=min(latm2d)

请问u1_intrp是什么？

放逐之刃

好像一条dog 发表于 2020-12-17 14:25
请问u1_intrp是什么？

wrfout中的风速

KJK76

谢谢，帮了大忙