output cleanup + aux vars handling in hdf5
trying to fix reindexing in courant()