Accumulated stress decreases when no interactions

Caltech-OCTO · Jun 27, 2024 · 2f6b16b · 2f6b16b
1 parent 4df98b9
commit 2f6b16b
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 4 deletions.
diff --git a/src/physical_processes/types.jl b/src/physical_processes/types.jl
@@ -33,7 +33,7 @@ Fields:
  that stress_accum = stress_accum(1-λ) + stress_instant(λ)
  α <AbstractFloat> Adjusts ellipse in stress space by raising area ratio to the α
 Note:
- τ is used in calc_stress!, whereas α is used in determine_fractures.
+ λ is used in calc_stress!, whereas α is used in determine_fractures.
 """
 @kwdef struct DecayAreaScaledCalculator{FT<:AbstractFloat} <: AbstractStressCalculator{FT}
  λ::FT = 0.2

diff --git a/src/physical_processes/update_floe.jl b/src/physical_processes/update_floe.jl
@@ -500,9 +500,8 @@ function timestep_floe_properties!(
  cforce = floes.collision_force[i]
  ctrq = floes.collision_trq[i]
  # Update stress
- if floes.num_inters[i] > 0
- calc_stress!(get_floe(floes, i), floe_settings, Δt)
- end
+ calc_stress!(get_floe(floes, i), floe_settings, Δt)
+
  # Ensure no extreem values due to model instability
  if floes.height[i] > floe_settings.max_floe_height
  @warn "Reducing height to $(floe_settings.max_floe_height) m"