Hi Sam, I looked over the ValidateSpatial briefly, but I'm not sure I understand it. However, > I could see us using a fix that does something like: > if len(data.pf.dimensions) > 2: > ds = div_fac * data['dz'].flat[0] > f += data["z-velocity"][1:-1,1:-1,sl_right]/ds > f -= data["z-velocity"][1:-1,1:-1,sl_left ]/ds seems pretty reasonable to me. j