diff --git a/src/physical_processes/types.jl b/src/physical_processes/types.jl index a8e10ec..afaf82f 100644 --- 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(λ) α 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 index 474d26d..f7ba643 100644 --- 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"