I am running a general circulation model (EPIC, which uses isentropic vertical coordinate system, potential temperature) and need to add vertical diffusion to the zonal velocity U. The vertical coordinate is theta (Θ, potential temperature) and there is no vertical height coordinate in z but I...