matlab关于画风场的问题

夕阳薇

{:soso_e154:}我有一个850hPa的全球风速场（Uwnd850，Vwnd850）的数据。如何绘制出大气上下层的风速场和散度场？初学matlab，自己研究了好久都调试不出来。请大家帮帮忙，谢谢!

夕阳薇

问题补充，我的程序是这样的：
%% prepare:
clear all;close all;clc;
%% open nc files:
fuwnd=netcdf('uwnd.mon.mean.nc','nowrite');
fvwnd=netcdf('vwnd.mon.mean.nc','nowrite');
%% data:
%850hPa&200hPa U_wind&V_wind
Uwnd850=squeeze(fuwnd{'uwnd'}(:,3,:,:)*fuwnd{'uwnd'}.scale_factor+fuwnd{'uwnd'}.add_offset);
Vwnd850=squeeze(fvwnd{'vwnd'}(:,3,:,:)*fvwnd{'vwnd'}.scale_factor+fvwnd{'vwnd'}.add_offset);
%longitude&latitude
LON=fuwnd{'lon'}(:);
LAT=fuwnd{'lat'}(:);
%%
[x,y]=meshgrid(LON,LAT);
quiver(x,y,Uwnd850,Vwnd850);
div = divergence(x,y,Uwnd850,Vwnd850);
由于我的数据是nc的，我不太会弄。。。求改错，求优化。。。

夕阳薇

求各路神仙指导啊。。{:soso_e154:}

夕阳薇

@wlzhongouc看了你的程序，觉得你很牛。请你帮我看看我的程序哪里出问题了好么？

jacky

非常感谢分享。

wlzhongouc

夕阳薇 发表于 2012-12-21 22:27

                               
登录/注册后可看大图

@wlzhongouc看了你的程序，觉得你很牛。请你帮我看看我的程序哪里出问题了好么？

没有我也在学习  计算的时候应该用实际的网格距离
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
        
[xx,yy]=meshgrid(dy,dx);

div=divergence(xx,yy,Uwnd850,Vwnd850);

这样还是会有一些问题  有的地方跟我想的不一样 我还没有找出问题

夕阳薇

wlzhongouc 发表于 2012-12-22 09:39

                               
登录/注册后可看大图

没有我也在学习  计算的时候应该用实际的网格距离
R0=6371004;
pi=3.1415926;

我修改之后只能画出风场，他的散度是一堆数，没有图。。。

wlzhongouc

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

                               
登录/注册后可看大图

我修改之后只能画出风场，他的散度是一堆数，没有图。。。

看明白你的问题是什么意思

313722504

matlab出血中

tsl08

菜鸟进来学习学习~