Fixed a bug in the routine to fill the matrix for an elliptic system; added a procedure to project data defined at the elements to data defined at the nodes