[原创]MATLAB矢量地理数据相关计算（涡度、散度、平流...）

sxq

楼主，我想请问一下，求curlz时公式后面的u.*tan()是什么意思

我不是谁wu

楼主，计算风应力旋度的程序，u,v,是用风速还是风应力呢？

kc121

最想学的是楼主的那个图用matlab如何画的？

lqe2004

谢谢分享！正好需要呢

muziyili

楼主赞啊！

muziyili

可惜我的金钱现在还不够下载的

muziyili

错误使用 gradient (line 49)
输出参数太多。

出错 dx_atmos (line 3)
    [dx, ~]=gradient(longitude);

出错 curlz_atmos (line 3)
    dx=dx_atmos(longitude, latitude);
这是怎么回事呢？

d1268y

谢谢分享!!!!

2daidai

楼主帮我看看程序吧

1. clc
   
2. clear
   
3. close all;
   
4. 
   
5. u=ncread('E:\study\monsoon\NCEP\uwnd\uwnd.2014.nc','uwnd');
   
6. v=ncread('E:\study\monsoon\NCEP\vwnd\vwnd.2014.nc','vwnd');
   
7. ncdisp E:\study\monsoon\NCEP\vwnd\vwnd.2014.nc
   
8. u_lon=ncread('E:\study\monsoon\NCEP\uwnd\uwnd.2014.nc','lon');
   
9. u_lat=ncread('E:\study\monsoon\NCEP\uwnd\uwnd.2014.nc','lat');
   
10. 
    
11. U=squeeze(u(43:51,30:34,3,:));
    
12. V=squeeze(v(43:51,30:34,3,:));
    
13. [X Y]=meshgrid(u_lon(43:51),u_lat(30:34));
    
14. for time=1:365;
    
15.     div1(:,:,time)=divh_atmos(X',Y',U(:,:,time),V(:,:,time));%850hap
    
16. end

复制代码

用的ncep的数据怎么div1全是NaN或者INF且U，V中没有缺省值。。不用divh_atmos,用divergence结果也一样。。非常悲伤

2daidai

muziyili 发表于 2017-4-18 10:35
错误使用 gradient (line 49)
输出参数太多。

好像用meshgrid重新把经纬度构建一下才可以[X Y]=meshgrid(u_lon(43:51),u_lat(30:34));div1(:,:,time)=divh_atmos(X',Y',U(:,:,time),V(:,:,time))我是这样改了下就没有那个报错了。但是算出来的结果又问题