fields: iDof internalE linearKE momentum q rhoM rhoN because the boundary values are required to record surface data.