- 积分
- 992
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2019-10-2
- 最后登录
- 1970-1-1
|
发表于 2022-4-13 22:35:00
|
显示全部楼层
clc,clear,clear off
ncdisp('F:\Data\hgt.mon.mean.nc')
hgt=ncread('F:\Data\hgt.mon.mean.nc','hgt');
time=ncread('F:\Data\hgt.mon.mean.nc','time');
lat=ncread('F:\Data\hgt.mon.mean.nc','lat');
lon=ncread('F:\Data\hgt.mon.mean.nc','lon');
format = 'mm dd, yyyy HH:MM:SS.FFF AM';
time = double(time);
dstr = datestr((datenum('1800-01-01') + time./24),format);
TM = datevec(dstr);
hgt=squeeze(hgt(:,:,3,385:884));%数据从1980年1月开始截取;截取至2020年12月 ;
hgt(:,:,501:504)=nan;
%%
hgt_year=reshape(hgt,[size(hgt,1),size(hgt,2),12,floor(size(hgt,3)/12)]);
hgt_mean=nanmean(hgt_year,4);
hgt_ano=hgt_year-hgt_mean;
% 截取数据
lon_need1=find(lon>=115&lon<=150);
lat_need1=find(lat>=15&lat<=25);
index=hgt_ano(lon_need1,lat_need1,6:8,:);
%%
%求出指数
WPSH_Osindex=squeeze(nanmean(nanmean(nanmean(index,1),2),3)); |
|