function repeated_wave_propagation_multi(iters,T::DataType=Float64;Nx::Int=128,λ=1.55e-6,w0=1e-2,dx1=1e-3,dxn=dx1,dz=5e1,Nz::Int=21,ng=3,kwargs...)
function repeated_wave_propagation_multi(T::DataType=Float64;iters=100,Nx::Int=128,λ=1.55e-6,w0=1e-2,dx1=1e-3,dxn=dx1,dz=5e1,Nz::Int=21,ng=3,kwargs...)