cfl tolerance

BenMql · Jul 5, 2024 · fa29a59 · fa29a59
1 parent 5f12ee7
commit fa29a59
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/src/plane_layer/PL_IMEX_timestepping_marching.f90 b/src/plane_layer/PL_IMEX_timestepping_marching.f90
@@ -221,6 +221,7 @@ subroutine compute_cfl_based_timestep(self)
    v_over_dy_max = maxVal(dAbs(self%linear_variables(2)%phys))
    v_over_dy_max = v_over_dy_max * self%geometry%NYAA / self%geometry%Ly
    w_over_dz_max = v_over_dy_max
+   if (self%cargo%cflFactor_along_z > 1.d-10) then
    do ix = 1, self%geometry%phys%local_NX
    do iy = 1, self%geometry%phys%local_NY
    w_over_dz_max = dMax1(w_over_dz_max,&
@@ -229,6 +230,7 @@ subroutine compute_cfl_based_timestep(self)
                                        self%cargo%cflFactor_along_z   ) )
    end do
    end do
+   end if
 
 
    self%cargo%cfl_based_DT = 1._dp/ dMax1 ( u_over_dx_max,&