docs: def io_surface(timestep, time, fid, fld):
      """Output for surface files"""
      fid.write("{} {}".format(timestep, time))
      fid.writelines(["%10.2e" % item for item in fld[:]])
      fid.writelines(["\n"])