function plot2(NX,NZ,NT,dx,dt,c,RES,BVEL,res,REALDATA,SYNDATA)
  
subplot(111);colormap('gray');imagesc([1:NX]*dx,[1:NT]*dt,RES');
   title('Residual Pressure Seismograms');
xlabel('Offset (ft)');ylabel('Time (s)');colorbar;pause(2);
  print -dps fig2a.ps;
   colormap('gray');subplot(111);imagesc([1:NX]*dx,[1:NZ]*dx,BVEL');
   title('Reconstructed Velocity (ft/s)');xlabel('Offset (ft)');
   ylabel('Depth (ft)');colorbar;pause(5);
  print -dps fig2b.ps;
figure(2);subplot(111);
plot([1:11],res(2:12),'*',[1:11],res(2:12),'-');axis([1 11 0 .14]);ylabel('Residual');
xlabel('Iteration Number');title('CSG Residual Norm vs Iteration No.');
print -dps fig3a.ps;
plot([150:NT]*dt,REALDATA(NX/2,150:NT)/max(REALDATA(NX/2,150:NT) ),'--');hold on;ylabel('Amplitude');
plot([150:NT]*dt,SYNDATA(NX/2,150:NT)/max(SYNDATA(NX/2,150:NT)),'-');ylabel('Amplitude'); xlabel('Time (s)')
;title('Observed (--) vs Predicted (-) Seismograms');  axis([150*dt NT*dt -1 1]);hold off;
print -dps fig3b.ps;
plot([1:NZ]*dx,BVEL(NX/2,:));hold on;
plot([1:NZ]*dx,c(NX/2,:),'--');
xlabel('Depth (ft)');title('Actual (--) vs Reconstructed (-) Velocity');
ylabel('Velocity (ft/s)'); axis([0 NZ*dx 4400 5300]);hold off;
print -dps fig3c.ps; 
