Use an interface for output_write*d

- generic function output_write now handles arrays of
  all 3 dimensions.
Merge pull request #14 from has2k1/discovered-interface

Use an interface for output_write*d