How can I display the vorticity from a velocity field?
First load your vector field with v = loadvec(‘B00001.vec’); Then display it with the background option ‘rot’ showvec(v,’rot’); Another way to do this is to create a scalar field curl, curl = vec2scal(v,’rot’); and then to display it: showscal(curl); Before computing the vorticity field, you will probably want to filter your vector field: showvec(filterf(v,1),’rot’); Many other vector-to-scalar conversions (divergence, kinetic energy, velocity derivatives…) are available from vec2scal.
