Index: /issm/trunk/src/c/parallel/GradJCompute.cpp
===================================================================
--- /issm/trunk/src/c/parallel/GradJCompute.cpp	(revision 1205)
+++ /issm/trunk/src/c/parallel/GradJCompute.cpp	(revision 1206)
@@ -49,4 +49,5 @@
 	femmodel->parameters->FindParam((void*)&sub_analysis_type,"sub_analysis_type");
 	femmodel->parameters->FindParam((void*)&numberofnodes,"numberofnodes");
+	femmodel->parameters->FindParam((void*)&numberofdofspernode,"numberofdofspernode");
 	femmodel->parameters->FindParam((void*)&solverstring,"solverstring");
 	femmodel->parameters->FindParam((void*)&control_type,"control_type");
@@ -57,5 +58,5 @@
 	diagnostic_core_nonlinear(&u_g,&K_ff0,&K_fs0,femmodel,inputs,analysis_type,sub_analysis_type);
 	VecToMPISerial(&u_g_double,u_g); VecFree(&u_g);
-	inputs->Add("velocity",u_g_double,2,numberofnodes);
+	inputs->Add("velocity",u_g_double,numberofdofspernode,numberofnodes);
 
 	_printf_("%s\n","      buid Du, difference between observed velocity and model velocity:");
@@ -76,5 +77,5 @@
 	_printf_("%s\n","      compute gradJ:");
 	VecToMPISerial(&lambda_g_double,lambda_g);VecFree(&lambda_g);
-	inputs->Add("adjoint",lambda_g_double,2,numberofnodes);
+	inputs->Add("adjoint",lambda_g_double,numberofdofspernode,numberofnodes);
 
 	Gradjx( &grad_g, numberofnodes,femmodel->elements,femmodel->nodes, femmodel->loads, femmodel->materials,