Implemented in pism::stressbalance::BlatterMod, pism::stressbalance::SIAFD, and pism::stressbalance::ConstantInColumn.