From 4d0fefd231e62bb3d9d1c92009bd0ec1a2b06614 Mon Sep 17 00:00:00 2001 From: Christoph Kords Date: Wed, 6 Apr 2011 11:47:24 +0000 Subject: [PATCH] security factor of 1.2 for CFL condition in constitutive_nonlocal_dotState --- code/constitutive_nonlocal.f90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/code/constitutive_nonlocal.f90 b/code/constitutive_nonlocal.f90 index df501bef3..f4eb861da 100644 --- a/code/constitutive_nonlocal.f90 +++ b/code/constitutive_nonlocal.f90 @@ -1614,8 +1614,8 @@ forall (s = 1:ns, t = 1:4, rhoSgl(s,t+4) * constitutive_nonlocal_v(s,t,g,ip,el) !**************************************************************************** !*** check CFL condition for flux -if (any(minval(mesh_ipVolume(ip,el) / max(tiny(1.0_pReal),mesh_ipArea(:,ip,el))) & - < constitutive_nonlocal_v(1:ns,1:4,g,ip,el) * timestep)) then +if (any(1.2_pReal * constitutive_nonlocal_v(1:ns,1:4,g,ip,el) * timestep & ! security factor 1.2 + > mesh_ipVolume(ip,el) / maxval(mesh_ipArea(:,ip,el)))) then dotState%p = NaN return endif