Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

a bug in module_cu_gf_wrfdrv.F ?

This post was from a previous version of the WRF&MPAS-A Support Forum. New replies have been disabled and if you have follow up questions related to this post, then please start a new thread from the forum home page.


New member
In ./phys/module_cu_gf_wrfdrv.F, the vertical velocities omega (pa/s) was calculated by omeg(I,K)= -g*rho(i,k,j)*w(i,k,j).
The w was obtained from grid%w_2, which was not at the mass grid point.
Therefore, is it more reasonable to calculate Omega after interplating w to the mass grid point?
You are right that W (hence omega) is vertically staggered, meaning it is located on the interface, not in the middle of the layer. However, omega is only used to compute the flux of Qv and mconv is summed in the vertical, so I think that is ok. Let me know if you think otherwise.