Index: /issm/trunk/src/m/solutions/adjointbalancedthickness_core.m
===================================================================
--- /issm/trunk/src/m/solutions/adjointbalancedthickness_core.m	(revision 5411)
+++ /issm/trunk/src/m/solutions/adjointbalancedthickness_core.m	(revision 5412)
@@ -15,8 +15,14 @@
 	femmodel=solver_linear(femmodel);
 
+	%Call SurfaceAreax, because some it might be needed by PVector
+	femmodel.elements=SurfaceArea(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters);
+
+	%compute adjoint
 	displaystring('\n%s',['      computing adjoint']);
 	femmodel=SetCurrentConfiguration(femmodel,BalancedthicknessAnalysisEnum,AdjointBalancedthicknessAnalysisEnum);
 	femmodel=solver_adjoint_linear(femmodel);
 
+	%Save results
+	femmodel.elements=InputToResult(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters,AdjointEnum);
 	if(solution_type==AdjointSolutionEnum & ~control_analysis)
 		displaystring(verbose,'\n%s',['      saving results...']);
Index: /issm/trunk/src/m/solutions/adjointdiagnostic_core.m
===================================================================
--- /issm/trunk/src/m/solutions/adjointdiagnostic_core.m	(revision 5411)
+++ /issm/trunk/src/m/solutions/adjointdiagnostic_core.m	(revision 5412)
@@ -19,4 +19,8 @@
 	femmodel=solver_diagnostic_nonlinear(femmodel,conserve_loads);
 
+	%Call SurfaceAreax, because some it might be needed by PVector
+	femmodel.elements=SurfaceArea(femmodel.elements,femmodel.nodes,femmodel.vertices,femmodel.loads,femmodel.materials,femmodel.parameters);
+
+	%compute adjoint
 	displaystring('\n%s',['      computing adjoint']);
 	if(isstokes), femmodel=SetCurrentConfiguration(femmodel,DiagnosticStokesAnalysisEnum,AdjointStokesAnalysisEnum);
@@ -24,4 +28,5 @@
 	femmodel=solver_adjoint_linear(femmodel);
 
+	%save results
 	if(solution_type==AdjointSolutionEnum & ~control_analysis)
 		displaystring(verbose,'\n%s',['      saving results...']);