Index: /issm/trunk-jpl/src/c/classes/FemModel.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 21521)
+++ /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 21522)
@@ -3208,5 +3208,5 @@
 
 		/*Copy vector in matrix*/
-		xMemCpy<IssmDouble>(&P0inputsold[i*numelementsold],vector,numelementsold);
+		for(j=0;j<numelementsold;j++) P0inputsold[j*numelementsold+i] = vector[j]
 		xDelete<IssmDouble>(vector);
 	}