This diagnostic routine calculates the hydrostatic consistency in each T-point and each layer. is defined as:
For the numerical calculation it is used here that and can be cancelled out each. For , the grid box is hydrostatically consistent, else it is called hydrostatically inconsistent. In the latter case, numerical problems can be expected for terrain-following coordinates when stratification is strong.
is stored in the 3d netcdf output file. USES:
use domain, only: imin,imax,jmin,jmax,kmax,az,au,av,HU,HV use variables_3d, only: hn,hun,hvn,hcc IMPLICIT NONEREVISION HISTORY:
Original author(s): Karsten Bolding & Hans BurchardLOCAL VARIABLES:
integer :: i,j,k REALTYPE :: du1,du2,dv1,dv2 REALTYPE :: x,y