I have a function for changing time and depth, simplified to:
def f(z,t): return np.exp(-z)*np.sin(tz) z = np.linspace(0,3000,num=3001) t = np.arange(0,40000,4000)
I want to build the result for all z at each time step in t, which will lead to something like this:

but I donβt know how to do it. I am sure this is very simple, but I'm not used to doing this in python.
source share