200高度场每月距平

xueqixiang

家园各位老师，改写了200高度场每月距平，时间1987—2016，可以运行，但始终觉得存在问题。请各位老师帮忙看看。谢谢
程序：
hgt1 = ncread('D:\matlab\MATLAB\wavelet\level\hgt.mon.mean.nc','hgt');
%hgt2 = hgt(8:85,17:46,10,97:456);
ran = 360;
hgt_mn=zeros(ran,2340);
for i=1:ran
    squ=squeeze(hgt1(8:85,17:46,10,i));                              
    hgt_mn(i,:)=reshape(squ,1,2340);   
end
hgt_nan=isnan(hgt_mn);      
k=0;


for j=1:2340
    if sum(hgt_nan(:,j)==0)  
        k=k+1;                 
        hgt_(:,k)=hgt_mn(:,j);  
    end
end
HGT=zeros(size(hgt_));
% 求每月的距平
for i=1:12
   
    HGT(i:12:ran-12+i,:)=hgt_(i:12:end,:)-repmat(mean(hgt_(i:12:end,:),1),size(hgt_(i:12:end,:),1),1);
end