%function []=vecplt2(cdf_names,start,stop);
start=[1990 10 23 15 0 0];
stop=[1991 1 14 0 0 0];
cdf_names=['3591-a1h.cdf'
           '3521-a1h.cdf'
           '3551-a1h.cdf'
           'u204ent2.cdf'
           'u304ent2.cdf'
           'u704ent2.cdf'];
hold off
%clear w
%[m,n]=size(cdf_names);
%load harry.ll
%for i=1:m
%  cdf=cdf_names(i,:);
%  lat0(i)=mcagt(cdf,'NC_GLOBAL','latitude');
%  lon0(i)=mcagt(cdf,'NC_GLOBAL','longitude');
%  w(:,i)=gvlp(cdf,start,stop);
%end
axis('square');
lat_mid=42.3;
lat_max=42.6;
lon_mid=-70.4;
dlon=1.2;
lat_fac=cos(lat_mid*pi/180);
dlat=dlon*lat_fac;
axis([lon_mid-.5*dlon lon_mid+.5*dlon lat_mid-.5*dlat lat_mid+.5*dlat]);
plot(harry(:,1),harry(:,2));
hold on
[n,m]=size(w);
i=1;
u=real(w(i,:));
v=imag(w(i,:));
%
% factor for arrows
%
xfac=.01;
yfac=xfac*lat_fac;
arrows(lon0,lat0,u,v,xfac,yfac,'g-');
uold=u;
vold=v;
for i=2:n;
  arrows(lon0,lat0,u,v,xfac,yfac,'i-');
  u=real(w(i,:));
  v=imag(w(i,:));
  arrows(lon0,lat0,u,v,xfac,yfac,'g-');
  uold=u;
  vold=v;
end

