TestD is a test on a square ice sheet.
The bed is flat as in TestC, but the basal drag has sinusoidal variations in x.
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.
