Pass thv as argument to bu function
Subtract horizontal average pressure field