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