NCEP风场求散度和涡度

蒙奇·D·斯

无法下载（自动回复：请不要使用迅雷等下载工具，使用IE右键另存为就可以）啊，桑心中，只能自己摸索了

蒙奇·D·斯

居然还有这种事情，下到了

meiguidaozhu

能不能发给我一份呢？315549202@qq.com

njzqxt

楼主，读不了NC数据啊。

小蓝不蓝

随缘 发表于 2012-4-23 20:57

                               
登录/注册后可看大图

原来是猛牛作品啊，不错啊

老刘很会卖萌啊。。。

wlzhongouc

njzqxt 发表于 2012-9-13 09:34

                               
登录/注册后可看大图

楼主，读不了NC数据啊。

看你用的Matlab版本了 不同的版本用的nc插件还不太一样

muri

好吧，级别不够，还得自力更生

matlab里面有计算地球球面距离的函数，
distance
貌似可以省去很多距离计算的麻烦

muri

无法下载（自动回复：请不要使用迅雷等下载工具，使用IE右键另存为就可以）（自动回复：请不要使用迅雷等下载工具，使用IE右键另存为就可以）你的程序，自己编了中国区域850的水平散度，贴出来大家看看没有没毛病
R0=6371004;  pi=3.1415926;
   lat0=55;  lon0=74;
    for lati=55:-1:15
      latm(55-lati+1)= distance(lat0,lon0,lati,lon0)/180*pi*R0;
    end
   
   for loni=74:1:137
      lonm(loni-73)= distance(lat0,lon0,lat0,loni)/180*pi*R0;
   end
   
   [xx,yy]=meshgrid(lonm,latm);
   div=divergence(xx,yy,u850,v850);  

ShaG

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

% calculate curvature vorticity
R=6.3781e6; % earth's radius
[dx, ~]=gradient(longitude);
dx=dx*pi./180;
[~, dy]=gradient(latitude);
dy=dy*pi./180;
dy(1, :)=NaN;  
dy(73, :)=NaN; % 由于使用差分的缘故，把dy的头尾元素设定为缺测
[dv, ~]=gradient(v);
[~, du]=gradient(u.*cos(latitude*pi./180));
curlz=(dv./dx-du./dy)./(R.*cos(latitude*pi./180));

zane

请问诸位大侠，小生计算到这里，如何把散度场绘制出来？？
x=-5:0.1:5;
>> y=x;
>> u=x.*x;
>> v=-y.*y;
>> [xx,yy]=meshgrid(x,y);
>> [uu,vv]=meshgrid(u,v);
>> div=divergence(xx,yy,uu,vv);

谢过~~