if exist('index')~=1
      index=input('index array to plot  '); 
end
bhmap(42.325, 42.35)
arrows(lon(index),lat(index),uva(index),vva(index),'-g',.0001,.00007)
arrows(-70.98,42.345,50,0,'-g',.0001,.00007)
text(-70.98,42.343,'50 cm/s')
timestart=gregorian(jd(index(1)));
timeend=gregorian(jd(index(length(index))));
datestring=[int2str(timestart(1)),'/',int2str(timestart(2)),'/90  '];
timestring=[int2str2(timestart(3),4),'-',int2str2(timeend(3),4)];
title(['vertically averaged current, ',datestring,timestring,' GMT'])

