
%Store seismograms

for i=1:nsp
   x= xsr(i);z=zsr(i);
   seismo(it,i,is) = p(x,z);
   UU(it,i,is) = u(x,z);
   WW(it,i,is) = w(x,z);
end
