%
OPATH=['../fcast_archives/',runid,'/OBS/']
UVOF=[OPATH,'UV_INTEG.m3d']
%
[UVO,gridname,year,ncol,gmt]=read_obs(UVOF);


gridname='../MESH/bank150';
gridstruct=loadgrid(gridname);

figure;

%

plotbnd(gridstruct);
% axis([150000 250000 -250000 -150000]);
axis('square');set(gca,'XTick',[]);set(gca,'YTick',[]);
hold on
V=UVO(ii,:);
hv=vecplot2(V(:,2),V(:,3),V(:,4),V(:,5),0.5,'m/s',230000.,-160000.);
set(hv,'color','blue');
title(verbatim(['Observations: RMS = ',num2str(rms(V(:,4)+sqrt(-1.0)*V(:,5))),' m/s']))
