- added scalar time series 'timestamp' to merged pism files
Well-known reference for PETSc, out of date.
White-space changes.
Tweaking PISMNC.

This makes it a little easier to write a field with a custom set of
attributes.
Fix basal B.C. for enthalpy (port to 0.5).

Now PISM will always use the insulating Neumann B.C. for enthalpy at the
base *if* the base is temperate. We no longer follow the decision chart
in [Aschwanden et al].

Previously PISM set B.C. in the temperate base case as follows:

- insulating Neumann B.C. if a non-zero-thickness layer of basal
  temperate ice is present, and
- Dirichlet B.C. otherwise
update routing test script to include thermo diagnostics