A scalar unknown is added to the node_dof records. The scalar is the residue in the differential equations. Also derivatives should be initialized.
