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

but I don't know how. I'm sure it's very simple, but I'm not used to doing this in python.
