How to implement vorticity calculations
 

Hello,

I made a code in fortran for solving LBM equations (the enemies of CFD :)).

Anyway, I have the three velocity components for node and I would like to plot the vorticity.

I know the matrix which define the vorticity, but how to implement in your opinion this to a fortran file. I mean the derivatives like (dw/dy-dv/dz)=w_x (w_x vorticity along x).

How can I implement it in a code:
- taking the velocity in z-direction at node i and i-1 divided by the distance between the nodes (which is constant) and calculate the vorticity?

Thank you ver much

Kind regards

Antonio