p00=p2/max(abs(p2(:))+.001);

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
subplot(122);
imagesc([1:nx]*dx,[1:nz]*dz,(p00+BVEL)); 

% imagesc([1:nx]*dx, [1:nz]*dz,p2+BVEL/50, [-0.2 0.2]);
colorbar;
hold on;

%plot([1:nx]*dx,[1:nx]*0.+zr*dz,'-b');
title(['Snapshot of Acoustic Waves: T = ',num2str(dt*it),' s']);
xlabel('X (km)');
ylabel('Z (km)');
plot(dx*xs,dz*zs,'*y');
hold off;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
subplot(121);
a=REALDATA';
a=a/max(abs(a(:)));
clim=[-0.02 0.02];

imagesc([1:nx]*dx, [1:nt]*dt, a, clim);
colormap(gray);

title('Seismograms of Acoustic Waves'); 
xlabel('X (km)');
ylabel('Time (s)'); 
%colorbar;