Thanks for the suggestion tiny tim.
I got my plotting software [mathematica] to automatically calculate the gradient of the z function so I don't have to worry about taking dz/dx...
This is the code i used: [mathematica]
(gradf[x_, y_, t_] = {D[z[x, y, t], x], D[z[x, y, t], y]...