**INTERFACE:**

subroutine hcc_check()

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 NONE

Original author(s): Karsten Bolding & Hans Burchard

integer :: i,j,k REALTYPE :: du1,du2,dv1,dv2 REALTYPE :: x,y