%                  sticker.m
%            gets individual cruise data from files gh#va.tid
%            and plots means, shears
%

choose=input('type 1 to load files ');
if choose==1
ncruise=input('Cruise number ');
ncruise=int2str(ncruise);

!ncd ghdata
eval(['load gh',ncruise,'va.tid']);
eval(['load gh',ncruise,'nav.bin']);
!ncd ghmat

lat=latbin;
lon=lonbin;

end

again=1;

while again==1

nlevel=input('velocity level (1-5), 0 for shear (1 minus 5) ');
if nlevel==0
    ucoef=r3d(uacoef,nrow,0,0,1)-r3d(uacoef,nrow,0,0,5);
    vcoef=r3d(vacoef,nrow,0,0,1)-r3d(vacoef,nrow,0,0,5);
else
    ucoef=r3d(uacoef,nrow,0,nlevel,0);
    vcoef=r3d(vacoef,nrow,0,nlevel,0);
end
stickgrd
again=input('type 1 for another plot ');

end
