TestC is a test on a square ice sheet.
The bed is fla in this test, but the basal drag has sinusoidal bumps in x and y.
To fit the model, we have to divide the drag coefficient by the effective pressure, 
and to do so, the effective must be stricly positive everywhere so we moved up the surface.
Periodic boundary conditions are applied on the boundaries.

