求助，帮忙看一个程序！

clc; clear;
load sjy.txt
sjy=flipud(sjy);
%
ncols      = 177;
nrows      = 64;
xllcorner  =  89.401306152344;
yllcorner   = 31.62291656;
cellsize   =  0.07272727;
NODATA_value =-9999;
%
lon=[xllcorner:cellsize:xllcorner+(ncols-1)*cellsize];
lat=[yllcorner:cellsize:yllcorner+(nrows-1)*cellsize];
[sjy_lon, sjy_lat]=meshgrid(lon, lat); clear lon lat
%-----
% Pixelsize = 0.07272727;
glb_ncols      = 4950;
glb_nrows      = 2091;
UL=      -179.99568;
LR =       -62.849808;
lon=[UL:cellsize:UL+(glb_ncols-1)*cellsize];
lat=[LR:cellsize:LR+(glb_nrows-1)*cellsize];
[ndvi_lon, ndvi_lat]=meshgrid(lon, lat); clear lon lat
%
mm_list=['jan';'feb';'mar';'apr';'may';'jun';'jul';'aug';'sep';'oct';'nov';'dec'];
k=0;
for yy=1982:2006
    yr=num2str(yy)
    yr=yr(3:4);
    for mm=1:12
        for tt=1:2
            k=k+1;
            if tt==1
                fna=['D:\matfiles\' yr  mm_list(mm, :) '15a.mat'];
            else
                fna=['D:\matfiles\' yr  mm_list(mm, :) '15b.mat'];
            end
            ndvidata=load(fna); ndvidata=ndvidata.ndvidata;
            data=flipud(double(ndvidata));  clear ndvidata
            flg=find(data<0); data(flg)=nan; clear flg
            dist=(ndvi_lon-sjy_lon(1,1)).^2+(ndvi_lat-sjy_lat(1,1)).^2;
            [xx yy]=find(~(dist>min(min(dist))));
            sjy_ndvi=data(xx:xx+nrows-1, yy:yy+ncols-1);
            flg=find(~(sjy==1));
            sjy_ndvi(flg)=nan;
            if k==1
                ndvi=sjy_ndvi;
            else
                ndvi(:,:,k)=sjy_ndvi;
            end
        end
    end
end
sjy_ndvi=ndvi; clear ndvi
%
hp=pcolor(sjy_lon, sjy_lat, mean(sjy_ndvi, 3));
set(hp,'edgecolor','none'); colorbar