Hi, I try to use NDSolve to do this partial differential equations, but the result is wrong... The source term is a Gaussian pulse as time and an exponential decay with the depth z. The result should be going to some small value at a large time, but it keeps "plat" when it goes to a large time. I attach my notebook in the attachment, thank you!
