[求助]关于速度势和辐散风计算

萤火飞侠

利用官网上的例子来计算多年速度势和辐散风，然后合成，发现速度势异常为正对应辐散而速度势异常为负对应辐合，和常理相反，请问各位大神有出现过这种情况吗？或者是我计算过程有什么处理没做吗？下面是我的程序，最后利用的是chi0和ud,vd三个变量
  f1    = addfile ("D:/necp/uwnd.mon.mean.nc", "r")
  f2    = addfile ("D:/necp/vwnd.mon.mean.nc", "r")
  ; 读取时间维度
  time = f1->time  ; time 为月数据，单位可以是 "months since 1961-01-01"
  ; 获取对应的年份和月份
  ytime  = cd_calendar(time, 0)
  itime  = ind((ytime(:,1).eq.3).and.(ytime(:,0).ge.1960).and.(ytime(:,0).le.2021))
  itime1  = ind((ytime(:,1).eq.4).and.(ytime(:,0).ge.1960).and.(ytime(:,0).le.2021))
  u1      = f1->uwnd(itime,:,:,:)
  v1      = f2->vwnd(itime,:,:,:)
  u2      = f1->uwnd(itime1,:,:,:)
  v2      = f2->vwnd(itime1,:,:,:)
  u       = (u1+u2)/2.
  v       = (v1+v2)/2.
  copy_VarMeta(u1, u)
  copy_VarMeta(v1, v)
  printVarSummary(u)
;*************************************************
; Calculate divergence and vortivity via spherical harmonics
; Use Wrap version to maintain meta data
;*************************************************

  div = uv2dvF_Wrap(u,v)           ; u,v ==> divergence 散度      在固定网格上给定 u 和 v 风分量（保留元数据），使用球谐计算发散度
  vrt = uv2vrF_Wrap(u,v)           ; u,v ==> vorticity 涡度       在固定网格上给定 u 和 v 风分量（保留元数据），通过球谐计算涡度
  printVarSummary(div)

;*************************************************
; Calculate velocity potential and stream function  
; 计算速度势和流函数
; Use Wrap version to maintain meta data
; 使用Wrap版本维护元数据
;*************************************************
  ;chi = ilapsG_Wrap ( div , 0)  ;高斯网格
  chi0 = ilapsF_Wrap ( div , 0)
  chi0@long_name = "velocity potential"
  chi0@units     = "m/s"
  printVarSummary(chi0)
;*************************************************
; Use spherical harmonics to derive divergent and rotational wind components
; 利用球谐函数导出辐散风分量和旋转风分量
; Use Wrap version to maintain meta data
;*************************************************   
  ;高斯网格
;  uvd = dv2uvG_Wrap(div) ; div  ==> divergent  wind components

  uvd = dv2uvF_Wrap(div) ; div  ==> divergent  wind components
  printVarSummary(uvd)
;*************************************************
; For clarity, explicitly extract the wind components
; 为了清晰，显式地提取风分量
; Add explicit lomg_name for each variable
; 为每个变量添加显式的lomg _ name
;*************************************************   
  ud  = uvd(0,:,:,:,:)     ;
  vd  = uvd(1,:,:,:,:)

  ud@long_name  = "Zonal Divergent Wind" ; 纬向辐散风
  vd@long_name  = "Meridional Divergent Wind" ; 经向辐散风

dieam

涡度和散度计算，要求纬度是递增的

千与千叶

涡度和散度计算，要求纬度是递增的

萤火飞侠

dieam 发表于 2024-12-10 23:20
涡度和散度计算，要求纬度是递增的

哎？不是要求从北到南吗

dieam

萤火飞侠 发表于 2024-12-11 10:42
哎？不是要求从北到南吗

不是哦，要求递增的哈

Arguments
u
v
wind components (input, arrays with two or more dimensions, rightmost two dimensions must be nlat x nlon)

input values must be in ascending latitude order
input array must be on a global grid