Stress balance free drift model #72
Conversation
Codecov ReportAttention: Patch coverage is
Additional details and impacted files@@ Coverage Diff @@
## main #72 +/- ##
==========================================
- Coverage 70.44% 69.24% -1.20%
==========================================
Files 27 28 +1
Lines 829 865 +36
==========================================
+ Hits 584 599 +15
- Misses 245 266 +21 ☔ View full report in Codecov by Sentry. 🚀 New features to boost your workflow:
|
| and an explicit part ``τaˣ = τaˣₑ + u * τaˣᵢ`` and ``τoˣ = τoˣₑ + u * τoˣᵢ`` (and similarly for the y-component) such that | ||
| ``` | ||
| uᶠ = (τoˣₑ - τaˣₑ) / (τoˣᵢ - τaˣᵢ) | ||
| vᶠ = (τoʸₑ - τaʸₑ) / (τoʸᵢ - τaʸᵢ) |
|
|
||
| A free drift parameterization that computes the free drift velocities as a balance between | ||
| top and bottom stresses ``τa ≈ τo`` where we split the stresses into a linear implicit part | ||
| and an explicit part ``τaˣ = τaˣₑ + u * τaˣᵢ`` and ``τoˣ = τoˣₑ + u * τoˣᵢ`` (and similarly for the y-component) such that |
There was a problem hiding this comment.
τ has units m^2 / s^2, u is m/s, and so if you want something of the form τ = something * u, then something has units m/s --- its a velocity.
Is that what you want? Not a quadratic term?
There was a problem hiding this comment.
right, I will change the docs and the names in the code to make sure it is more understandable.
I am not sure this is necessarily what we want, this is how we treat the stress in the SeaIceMomentumEquations (semi-implicitly), but maybe here we can treat it quadratically. I ll need to add a couple more methods.
There was a problem hiding this comment.
We probably need to solve a nonlinear system of equations if we want the whole implicit quadratic drag
There was a problem hiding this comment.
It might be simpler if one of the two stresses does not depend on the ice velocity, so that we can write an explicit formula. I ll add a method for when the top stress is independent on the ice velocity and when the bottom stress is independent of the ice velocity. We can leave the case in which both stresses are implicit for which we need a non-linear solve for the future.
|
An example of the dynamics = StressBalanceFreeDrift(top_stress, bottom_stress)sea_ice_advected_by_anticyclone.mp4 |
This PR adds a stress balance free drift calculation that computes the velocity as a balance between the atmospheric stress at the top and the ocean stress at the bottom.
This addition warrants the change of the module name from
SeaIceMomentumEquationstoSeaIceDynamicssinceSeaIceMomentumEquationbecomes only one of the dynamics model that we can use