pism/pism — Issue #26: inward scheme for calculation of principal strain rates not rotationally invariant in stable version
Author: talbrecht
Created: 2012-02-03 10:21:44+00:00
State: open
Labels: bug

The inward-scheme for the calculation of strain rates in stressbalance/SSA.cc in SSA::compute_principal_strain_rates() still asks for the ice thickness of the neighbors, but not for all four directions. In the dev version I committed a mask-based solution. Since eigencalving uses principla components of the strain rates in the penultimate grid cell upstream from the front (offset=2), this should not impact the results so much.