NCEP风场求散度和涡度

zane · 发表于 2012-11-9 09:11:22

ShaG 发表于 2012-10-28 10:39

登录/注册后可看大图

哈哈又见楼主，不过我觉气象上计算涡度和散度的时候是不可以用MATLAB自带的函数的，应该要把地球的曲率也考 ...

求出散度涡度后，如何把他们画出来，用什么函数？？急

zane · 发表于 2012-11-9 09:17:31

好像要那个程序啊~ 欲求！！ lzzoe1992@163.com

维欻? · 发表于 2012-11-14 12:42:07

谢谢楼主呀，学习中

三江魂 · 发表于 2012-11-15 10:42:12

厉害，厉害

neinei421 · 发表于 2012-12-17 00:29:37

虽然没有看到代码，不过也学到了不少东西，呵呵，等我画出来了也给大家看看

夕阳薇 · 发表于 2012-12-22 22:50:06

求教画用divergence如何画出散度场的图来，级别太低不能看到源程序，但是我自己是实在想不出办法了。。。@wlzhongouc @mofangbao

wlzhongouc · 发表于 2012-12-23 12:48:36

夕阳薇发表于 2012-12-22 22:50

登录/注册后可看大图

求教画用divergence如何画出散度场的图来，级别太低不能看到源程序，但是我自己是实在想不出办法了。。。@w ...

我在你的帖子里放上用到的程序了啊

夕阳薇 · 发表于 2012-12-26 10:21:53

wlzhongouc 发表于 2012-12-23 12:48

登录/注册后可看大图

我在你的帖子里放上用到的程序了啊

我弄好啦，我的那个程序问题在于速度场不是二维的，多了一个时间维……谢谢楼主先

wlzhongouc · 发表于 2012-12-26 12:59:48

夕阳薇发表于 2012-12-26 10:21

登录/注册后可看大图

我弄好啦，我的那个程序问题在于速度场不是二维的，多了一个时间维……谢谢楼主先

斥鷃 · 发表于 2013-1-31 17:42:32

本帖最后由斥鷃于 2013-1-31 17:53 编辑

只做过散度分析，没有做过涡度的，这里提出我的一点浅见：
1.divergence函数不能用于地理散度计算，原因自然是经距会随纬度变化这个问题了，divergence函数只适用于横纵坐标轴均匀的坐标系散度计算。具体解释很简单，看open divergence就行了，源代码写的很简单，个人觉得这个matlab函数写的有点水……重要部分我贴在下面（引自divergence函数源代码）：
if isempty(w)  % 2D

  [msg x y] = xyuvcheck(x,y,u,v);  error(msg)
  if isempty(x)
[px junk] = gradient(u); %#ok
[junk qy] = gradient(v); %#ok
  else
hx = x(1,:);
hy = y(:,1);
[px junk] = gradient(u, hx, hy); %#ok
[junk qy] = gradient(v, hx, hy); %#ok
  end
  div = px+qy;
2.49楼正解！顶一个，我在处理的时候直接用矩阵加减法做的中心差分求，做出来结果和grads绘出图形一致，不过代码不见得有49楼（他的方法去掉边界后也应该是中心差分）好，就不献丑了。但是另一方面49楼同学的做法在极地投影下也会出问题，这个的解决方法我按LZ和无径之林的方法试了试，比较笨~没成功~

以后什么时候有空再研究研究~

		自动登录	找回密码
密码			立即注册

[源程序] NCEP风场求散度和涡度