Modularize some initial and boundary condition code.
Deal with hdf/asymptote issues.