function []=vecplt(cdf1,cdf2);
clg
lat0(1)=mcagt(cdf1,'NC_GLOBAL','latitude');
lon0(1)=mcagt(cdf1,'NC_GLOBAL','longitude');
lat0(2)=mcagt(cdf2,'NC_GLOBAL','latitude');
lon0(2)=mcagt(cdf2,'NC_GLOBAL','longitude');
load harry.ll
[r,t,w(:,1),w(:,2)]=bcomcor(cdf1,'east','north',cdf2,'east','north');
axis('square');
axis([-71.1 -70.1 41.7 42.4]);
plot(harry(:,1),harry(:,2));
[n,m]=size(w);
i=1;
u=real(w(i,:));
v=imag(w(i,:));
arrows(lon0,lat0,u,v,.01,.007,'g-');
hold on
uold=u;
vold=v;
for i=2:n;
  arrows(lon0,lat0,u,v,.01,.007,'i-');
  u=real(w(i,:));
  v=imag(w(i,:));
  arrows(lon0,lat0,u,v,.01,.007,'g-');
  uold=u;
  vold=v;
end

