求助下面程序用grads怎么写 非常感谢

waylee2008

matlab程序
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');

river

你这是要计算什么，vhurl涡度？说清楚了才能帮你嘛，你给个这样的程序别人比一定能看懂啊。就像我去matlab区发一个grads的脚本文件，问你用matlab怎么编程，你也不一定能看懂啊

waylee2008

olr.mon.mean.nc文件
求涡度场 公式如下

river

waylee2008 发表于 2012-11-21 15:57

                               
登录/注册后可看大图

olr.mon.mean.nc文件
求涡度场 公式如下

看来你是一点儿不懂grads，也不懂你用的这个资料啊。你给的程序里面只用到了olr这个资料，它不是向外长波辐射的资料吗？你给的公式里面怎么又都成涡度了。不懂！提问之前先把自己的思路理清楚，不明不白的这样得问好几次别人才能知道你要干什么，效率太低