Temperature is initialized by assuming all surfaces are more or less
flat and sp01 defines the base elevation. Temperature is then calculated
as a steady state solution with T=0 at bottom and T=1 at ztemp.