dir{1}='/cove/usr4/jhyatt/plume/orig/ecom3d.cdf';
cdf=dir{1};                                      
var='salt'
time=130
jindex=30
[u,x,z,u2,x2,z2] = ecom_jslice(cdf,var,time,jindex);

lev=[25 32];

figure(1)
pcolor(x2,z2,u2)
caxis(lev)
cdot(x,z,u,jet,20,0,lev)
line(x,z,'marker','.','linestyle','none','markersize',10,'color','black');
axis([1 40 -40 1])

 
%figure(2)
%pslice(x2,z2,u2,lev)
%shading faceted
%axis([1 40 -20 1])

