figure(1);plot(zz,p2,xxx,yyy,'--');
axis([0 2000 -.25 1.2]);
title(['PRESSURE SNAPSHOT  T= ',num2str(dt*it),' sec']);xlabel('Depth (ft)'); pause(.01)