格点数据算空间上（各格点）气候倾向率，for循环

奚若

格点数据，想算空间上（各格点）气候倾向率，其中，pre_116数据为51*71*116。求教各位；以下是脚本
clear all;
clc;
fdir='F:\毕业论文\MATLAB\Shanxi_CRU_ts3.25_1901-2016_pre_mon_interp-0.1_nc_mask.nc';
load('F:\毕业论文\MATLAB\结果\数据\time1901-2016');
%F:\毕业论文\MATLAB\Shanxi_CRU_ts3.25_1901-2016_pre_mon_interp-0.1_nc_mask.nc将字符串定义为fdir，方便操作，减少语句长度。
ncid=netcdf.open(fdir,'NC_NOWRITE');
ncdisp(fdir)
pre=ncread(fdir, 'pre');
pre1 = pre(:,:,1:1392);%将1392个月的格点的降水量赋值给pre1矩阵
[nlon, nlat,ntim]=size(pre1);%将1的属性返到前面
for i=1:nlon
    for j=1:nlat;%从格点（1,1）开始
        y=1;
        for t=1:12:ntim%第一年1月，每年开始月份间隔为12           
          pre_116(i,j,y)=sum(pre1(i,j,t));%将格点（i,j）的一年的季节平均降水赋值给前面矩阵
          y=y+1;
        end
    end
end
%这是求每一年的每个格点的季降水
for i=1:nlon
    for j=1:nlat;%从格点（1,1）开始
        o=1;
        for t=1:116%第一年1月，每年开始月份间隔为12           
          pre_116_t(i,j,t)=pre_116(i,j,t);
          o=o+1;
          pre_116_t_z=pre_116_t(i*j,t)%目的是让pre_116_t转为116*1
        end
        LI_cctr(time,pre_116_t)%LI_ccctr（）为调用函数
    end
end
%求各点气候倾向率

飞天猫

  pre_116(i,j,y)=sum(pre1(i,j,t)) 求出来的是 所有1月份的降水累计值