In theory, thicknesses (ice and snow), should be advected with the div_Uℵh operator, which scales the thickness by sea ice then renormalizes by the current thickness. However, at the moment we use horizontal_div_Uc because the div_Uℵh is very unstable and sensitive to regions of low concentration.
I opened this issue as a reminder to verify and improve the robustness of the div_Uℵh operator such that we can use it in practical simulations
In theory, thicknesses (ice and snow), should be advected with the
div_Uℵhoperator, which scales the thickness by sea ice then renormalizes by the current thickness. However, at the moment we usehorizontal_div_Ucbecause thediv_Uℵhis very unstable and sensitive to regions of low concentration.I opened this issue as a reminder to verify and improve the robustness of the
div_Uℵhoperator such that we can use it in practical simulations