MATLAB程序用grads怎么写 恳请各位大侠帮帮忙！

waylee2008

clc,clear,close

ncload olr.mon.mean.nc;

for it=1:456;
   for iy=1:73;
       for ix=1:144;
               if(ix==1);%west
                   x1=144;
               else
                   x1=ix-1;
               end
               if(ix==144);%east
                   x2=1;
               else
                   x2=ix+1;
               end
               if(iy==1);%south
                   y1=73;
               else
                   y1=iy-1;
               end
               if(iy==73);%north
                   y2=1;
               else
                   y2=iy+1;
               end

           vhurl(it,iy,ix)=(4*olr(it,iy,ix)-olr(it,iy,x1)-olr(it,iy,x2)-olr(it,y1,ix)-olr(it,y2,ix))/4;
       end
   end
  % disp(it)
end
save Vhurl.mat vhurl
save lonlat.mat lon lat


clc,clear,close
load Vhurl.mat
load lonlat.mat
t=1
% vhurl(1,34,35)
lat=lat(5:end-4);vhurl=vhurl(:,5:end-4,:);
vhurl1=squeeze(vhurl(t,:,:));

[Lon,Lat]=meshgrid(lon,lat);
m_proj('miller','lat',[lat(end) lat(1)],'lon',[lon(1) lon(end)]);
  m_pcolor(Lon,Lat,vhurl1);shading interp
%   hold on
%   m_contour(Lon,Lat,vhurl1,[-400:100:400],'linewidth',2,'linecolor','k');
    bar=colorbar;
          set(bar,'fontname','Times New Roman','fontweight','bold');
         m_coast('patch',[0.6 0.6 0.6]);
          m_grid('linewi',1.7,'linest','none','tickdir','in','fontsize',11,...
                'fontname','Times New Roman','fontweight','bold');

Aires

楼主不能奢望别人会给你写程序吧