INTERFACE:
subroutine cfl_check()DESCRIPTION:
This routine loops over all horizontal grid points and calculated the maximum time step according to the CFL criterium by Beckers and Deleersnijder (1993):
(95) |
with the local shallow water wave speed
(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 case 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