Add geopotential height variable

The continuous equation is solving the advection of geopotential height
actually (not geopotential depth), so add 'gh' variable in the model.