matlab画2012年11月NCEP风场和散度场分享（基础版）

夕阳薇

第一次用matlab画出来的图，比较基础，没有地图。只是想和初学者分享一下。地图怎么弄还在研究，争取弄好也跟大家分享。还有要谢谢@wlzhongouc 给我的指导。 z5.m (1.18 KB, 下载次数: 233)

2012-12-26 10:40 上传

点击文件名下载附件
阅读权限: 20 至于数据怎么下载，网上有不少攻略。我是初学者，欢迎大家交流心得~

wlzhongouc

看了一下你的程序 你的x,y好像没有用实际距离，计算出来的散度量级上可能不对

Aires

楼主继续努力，希望分享更好的成果~

kongfeng0824

要加入地图数据吗？还是？可以利用load coast或者worldmap命令或者加载一个M_map命令。

夕阳薇

wlzhongouc 发表于 2012-12-26 13:03

                               
登录/注册后可看大图

看了一下你的程序 你的x,y好像没有用实际距离，计算出来的散度量级上可能不对

电脑一下坏了好些天，一直没来看看。实际距离是怎么操作呢呢？初学只是明白点儿皮毛，师兄之前给借鉴看看的程序，我还有没看明白的地方……

wlzhongouc

夕阳薇 发表于 2013-1-12 14:21

                               
登录/注册后可看大图

电脑一下坏了好些天，一直没来看看。实际距离是怎么操作呢呢？初学只是明白点儿皮毛，师兄之前给借鉴看看 ...

哪里没看明白？

夕阳薇

wlzhongouc 发表于 2013-1-12 17:44

                               
登录/注册后可看大图

哪里没看明白？

就是师兄提到的用实际距离的这一部分，我不明白这一段程序的原理或者说公式是什么样的，为什么是这样的。可能我理论知识学得不好。
R0=6371004;
pi=3.1415926;
[m n]=size(Uwnd850);
deltx(1:m)=2*pi*R0*cosd(lati)/360*mean(diff(longi));          % x lon (the actual distance of the resolution of data in longitude)
delty(1:n)=2*pi*R0/360*abs(mean(diff(lati)));             % y lat (the actual distance of the resolution of data in latitude)
for li=1:m
     dx(li)=sum(deltx(1:li),2);
end
for li=1:n
      dy(li)=sum(delty(1:li),2);
end

wlzhongouc

夕阳薇 发表于 2013-1-12 18:58

                               
登录/注册后可看大图

就是师兄提到的用实际距离的这一部分，我不明白这一段程序的原理或者说公式是什么样的，为什么是这样的。 ...

要计算实际距离你首先需要计算经纬度的网格实际间距吧，deltx(1:m)=2*pi*R0*cosd(lati)/360*mean(diff(longi));    和   delty(1:n)=2*pi*R0/360*abs(mean(diff(lati)));  就是计算这个的，然后你每隔一个网格间距再累加实际距离   dx(li)=sum(deltx(1:li),2);   和  dy(li)=sum(delty(1:li),2);

long_2013

初学  看看

liujianbao

我也是初学者，学习学习