INTERFACE:
subroutine cfl_check()DESCRIPTION:
This routine loops over all horizontal grid points and calculated the maximum time step according to the shallow water criterium by Beckers and Deleersnijder (1993):
![]() |
(95) |
with the local Courant number
![]() |
(96) |
where is the gravitational acceleration and
is the local
bathymetry value. In case that the chosen micro time step
is larger than
, the program will be aborted. In any
the CFL diagnostics will be written to standard output.
USES:
use parameters, only: g use domain, only: imin,imax,jmin,jmax,H,az #if defined(SPHERICAL) || defined(CURVILINEAR) use domain, only: dyc,dxc #else use domain, only: dy,dx #endif use m2d, only: dtm IMPLICIT NONEREVISION HISTORY:
Original author(s): Karsten Bolding & Hans BurchardLOCAL VARIABLES:
integer :: pos(2),max_pos(2),rc,i,j REALTYPE :: h_max=-99.,c,max_dt,dtt logical, dimension(:,:), allocatable :: lmask