file << "#"

  file << "# 1:iteration"

       << "time"

       << "2:time"

       << "level"

       << "3:level"

       << "component"

       << "4:component"

       << "i"

       << "5:i"

       << "j"

       << "6:j"

       << "k"

       << "7:k"

       << "x"

       << "8:x"

       << "y"

       << "9:y"

       << "z";

       << "10:z";
  int col = 11;

    file << "\t" << varname;

    file << "\t" << col++ << ":" << varname;