matlab编写南海夏季风指数程序

陈伯英

想问求南海夏季风指数的matlab程序怎么修改，参考了 李崇银南海夏季风指数matlab程序http://bbs.06climate.com/forum.p ... &fromuid=128407
(出处: 气象家园) 这个帖子，不是很看得明白

1. clc
   
2. clear
   
3. close all;
   
4. 
   
5. div1=zeros(1,365);
   
6. div2=zeros(1,365);
   
7. D=zeros(1,365);
   
8. Id=zeros(1,365);
   
9. u=ncread('I:\NCEPdata\pressure\uwnd.1986.nc','uwnd');
   
10. v=ncread('I:\NCEPdata\pressure\vwnd.1986.nc','vwnd');
    
11. ncdisp I:\NCEPdata\pressure\uwnd.1986.nc
    
12. ncdisp I:\NCEPdata\pressure\vwnd.1986.nc
    
13. u_lon=ncread('I:\NCEPdata\pressure\uwnd.1986.nc','lon');
    
14. u_lat=ncread('I:\NCEPdata\pressure\uwnd.1986.nc','lat');
    
15. 
    
16. U=squeeze(u(42:50,3:7,3,:));
    
17. V=squeeze(v(42:50,3:7,3,:));
    
18. [X,Y]=meshgrid(u_lon(42:50),u_lat(3:7));
    
19. for time=1:365
    
20.     div1(:,:,time)=divh_atmos(Y',X',U(:,:,time),V(:,:,time));%850hap
    
21. end
    
22. for time=1:365
    
23.     div2(:,:,time)=divh_atmos(Y',X',U1(:,:,time),V1(:,:,time));%200hap
    
24. end
    
25. for i=1:365
    
26.    D(i)=nanmean(nanmean((div1(:,:,i)-div2(:,:,i)),1),2);
    
27. end
    
28. 
    
29. for i=2:365
    
30.     Id(i)=D(i)/sqrt(sum(D(1:i).^2)/i);
    
31. end
    
32. Id(1)=0;

复制代码

wjy_ecnu

这个文件损坏了  重新下载试试

陈伯英

wjy_ecnu 发表于 2023-3-9 18:23
这个文件损坏了  重新下载试试

感谢，可以运行了