p00=p0/max(abs(p0(:))+.001);
subplot(122);imagesc([1:nx]*dx,[1:nx]*dx,(p00+BVEL)); colorbar;
hold on;plot([0 nx*dx],[xs+recdepth xs+recdepth]*dx,'-w');
title(['Snapshot of Acoustic Waves: T = ',num2str(dt*it),' s']);
xlabel('X (km)')
ylabel('Z (km)');
plot(dx*zs,dx*xs,'*w');pause(.01);hold off
subplot(121)
a=REALDATA';a=a/max(abs(a(:)));clim=[-.672 .672];
imagesc([1:nx]*dx,[1:nt]*dt,a,clim);
title('Seismograms of Acoustic Waves'); xlabel('X (km)'); ylabel('Time (s)')
xlabel('X (km)')
ylabel('Time (s)'); pause(.1);

