From a version of this code that was floating around, added read_2d_restart
to Science.[ch]pp; this reads a 2D array from disk that was previously
written by a SPINS job and extends it along the y-dimension to fill the
given array; this is useful for extending 2D runs to 3D.