Could anybody say to me how to get the vorticity in PHOE? Is there any pre-defined variable where it is stored? or do I need to add my own ground code? Thanks .

You need to add your own ground code....or export the data to matlab and calculate post processing....this is okay for steady state solutions...but for unsteardy ones better to calc as you go in ground..

Cheers Dan

Dan is correct about the use of GROUND, but note that the file c:\phoenics\d_earth\d_opt\d_advmph\gxintp.for contains subroutine GXVORT which computes the mean vorticity vector for cartesian and cylindrical meshes. This subroutine could be called from the user's GROUND. Some discussion is given in the comments at the top of this subroutine. The file could also readily be extended to handle BFC meshes.