[MITgcm-support] discretization of curl: vorticity

[Marco Reale]

Dear MITGCM users,

first I want to thank Dave , Jean-Michel  and Ryan for their help and useful suggestions in the previous post about momentum balance: averaging over many time steps gave me a good agreement between TOTUTEND/86400 and Um_Advec+Um_Diss+Um_dPHdx+dEta/dx+Um_ExternalForcings. Now I’m facing the question about the computation of vorticity balance. To compute vorticity I need to apply the curl to each member of  momentum equation.

My question is : Should I apply the curl in the form that is applied to U,V velocity fields as the model does to each member of Momentum budget (Um_Advec+Um_Diss+Um_dPHdx+dEta/dx+Um_ExternalForcings and Vm_Advec+Vm_Diss+Vm_dPHdx+dEta/dy+Vm_ExternalForcings) or there exists another to discretize the curl which might more suitable for this kind of computations (for example in the case of pressure terms)?I’ve already applied the model formula to my U,V output to compute Z3 and it works great. 

Looking forward to a response to my mail

Cheers

Marco





***********************************
Marco Reale
ph.D student "Environmental and Industrial Fluid Mechanics"-XXVII cycle
University of Trieste-Department of Mathematics and Informatics
National Institute of Geophysics and Oceanography (OGS)
International Center for Theorethical Physics (ICTP)
Italian Geophysical Association (AGI)
Italian Meteorological Society (SMI)
European Geophysical Union (EGU)
web-page: http://phdfluidmechanics.appspot.com/stud.html?Reale
skype: mcreal82
phone number: +39 3270888069

***********************************

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20131211/f48cdd64/attachment-0001.htm>