请教NCL绘制风场的辐散图

黃小言

最近在学习绘制风场的辐散图，我用的是ERA Interim的UV风场资料，打算绘制风场的辐散图，
我用了二种方法（如下表示）进行绘制，绘出来的图是一样的，可是我另外用ERA Interim的divengence资料绘制，出来的结果是完全不一样，请教各位高手，是什么原因造成如此大的差别？那一张才是正确的？？？

以UV风场资料绘制：
f    = addfile ("netcdf-web219-20141125045259-24015-72242.nc", "r")
u    =short2flt(f->u)
v    =short2flt(f->v)


uu=dim_avg_n_Wrap(u,0)
vv=dim_avg_n_Wrap(v,0)                 


;*************************************************


方法一：
  dv = uv2dvG_Wrap(uu,vv)                ; u , v ==> divergence
  uvd= dv2uvG_Wrap(dv)                   ; dv  ==> divergent  wind components
                                         ;     ==>uvd(0,nlat,mlon) -->u
                                         ;     ==>uvd(1,nlat,mlon) -->v


;*************************************************
方法二：
  dv    = new ( dimsizes(uu), typeof(uu),uu@_FillValue)
  ud    = new ( dimsizes(dv), typeof(dv),dv@_FillValue)
  vd    = new ( dimsizes(dv), typeof(dv),dv@_FillValue)


  uv2dvg(uu,vv,dv) 　　　　　　　; u,v  ==> divergence
  dv2uvg(dv,ud,vd) 　　　　　　　; dv  ==> divergent  wind components


;*************************************************
...以下略

以divengence资料绘制：
f    = addfile ("netcdf-web223-20141209083223-9912-2142.nc", "r")
d    =short2flt(f->d)


dv=dim_avg_n_Wrap(d,0)


;*************************************************
方法一：
  uvd= dv2uvG_Wrap(dv)                   ; dv  ==> divergent  wind components
                                         ;     ==>uvd(0,nlat,mlon) -->u
                                         ;     ==>uvd(1,nlat,mlon) -->v


;*************************************************
...以下略

0708011

input values must be in ascending latitude order
input array must be on a global grid
看看你的资料纬度是-90...90 这样的吗？printVarSummary看看。如果不是要u=u(:,::-1,:)

黃小言

终于都对上了～！真是受教了，谢谢～

songyan_711

谢谢～

oubababa

学习学习谢谢

sunnygirlbo

想问一下，楼主最后问题解决了没?吧风场分为辐散风是用哪一种方法呢？我用第二种方法算了一下，感觉和文章里面的不太对，麻烦楼主了

黃小言

sunnygirlbo 发表于 2015-12-27 15:21
想问一下，楼主最后问题解决了没?吧风场分为辐散风是用哪一种方法呢？我用第二种方法算了一下，感觉和文章 ...

是的,后來不管用哪一种方法结果都一样~

人生路

怎么感觉你那个大小不一致呢，方向一样

黃小言

人生路 发表于 2016-6-12 18:27
怎么感觉你那个大小不一致呢，方向一样

是的,所以後來修正後就一樣了~

yeah...

请问楼主最后是用哪种方法计算的？我还想请教下楼主那个divergence component wind指的是用势函数求出的散度风吗