function gomviz(cdf,var,tstep,depth,vax,ax,plotv)
% function gomviz(cdf,var,tstep,depth,[vax],[ax],[plotv])

if(exist('vax')~=1),
 vax=[nan nan];
end
load /usr/users/rps/gom/dino/wgom_half.mat;
[long_name,units]=base_att(cdf,var);
[u,x,y,jd]=zslice(cdf,var,tstep,-abs(depth));
pslice(x/1000,y/1000,u,vax,units);...
line(coast(:,1),coast(:,2),'color','white');...
xlabel('km');ylabel('km');
date_string=greg2str(gregorian(jd));
depstr=int2str(abs(depth))
title(['Slice of ' long_name ' at ' depstr ' m depth on ' date_string(1:10)]);
if(exist('ax')==1),
  if(isempty(find(isnan(ax)))),
     axis(ax)
  end
end
if(exist('plotv')==1)
  [w,x,y]=zsliceuv(cdf,tstep,-abs(depth));
  psliceuv(x/1000,y/1000,w,plotv(1),plotv(2),'white','none');
end 
