[Yt-dev] DivV in < 3D

j s oishi jsoishi at gmail.com
Tue Jul 12 17:05:07 PDT 2011


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



More information about the yt-dev mailing list