Added writeout routine, viscoelastic variables,
... writing viscoelastic fluxes (not finished)