Index: /issm/trunk/src/m/solutions/jpl/ControlInitialization.m
===================================================================
--- /issm/trunk/src/m/solutions/jpl/ControlInitialization.m	(revision 3961)
+++ /issm/trunk/src/m/solutions/jpl/ControlInitialization.m	(revision 3962)
@@ -66,5 +66,5 @@
 m_ds.y_g=zeros(m_ds.nodesets.gsize,1);
 m_ds.y_g(dofsetgen([1,2,3],4,m_ds.nodesets.gsize))=u_g;
-[m_ds.ys m_ds.ys0]=Reducevectorgtos(m_ds.y_g,m_ds.nodesets);
+m_ds.ys=Reducevectorgtos(m_ds.y_g,m_ds.nodesets);
 
 %Compute Stokes velocities once to have a reasonably good ug in input
Index: /issm/trunk/src/m/solutions/jpl/CreateFemModel.m
===================================================================
--- /issm/trunk/src/m/solutions/jpl/CreateFemModel.m	(revision 3961)
+++ /issm/trunk/src/m/solutions/jpl/CreateFemModel.m	(revision 3962)
@@ -24,5 +24,5 @@
 
 	displaystring(md.verbose,'%s','   reducing single point constraints vector...');
-	[m.ys m.ys0]=Reducevectorgtos(m.yg.vector,m.nodesets);
+	m.ys=Reducevectorgtos(m.yg.vector,m.nodesets);
 	
 	displaystring(md.verbose,'%s','   normalizing rigid body constraints matrix...');
Index: /issm/trunk/src/m/solutions/jpl/diagnostic_core.m
===================================================================
--- /issm/trunk/src/m/solutions/jpl/diagnostic_core.m	(revision 3961)
+++ /issm/trunk/src/m/solutions/jpl/diagnostic_core.m	(revision 3962)
@@ -44,5 +44,5 @@
 	if ismacayealpattyn,
 		model.dh.y_g=u_g;
-		[model.dh.ys model.dh.ys0]=Reducevectorgtos(model.dh.y_g,model.dh.nodesets);
+		model.dh.ys=Reducevectorgtos(model.dh.y_g,model.dh.nodesets);
 	end
 
@@ -90,5 +90,5 @@
 		model.ds.y_g(dofsetgen([1,2],4,model.ds.nodesets.gsize))=u_g;
 		model.ds.y_g(dofsetgen([3],4,model.ds.nodesets.gsize))=u_g_vert;
-		[model.ds.ys model.ds.ys0]=Reducevectorgtos(model.ds.y_g,model.ds.nodesets);
+		model.ds.ys=Reducevectorgtos(model.ds.y_g,model.ds.nodesets);
 
 		displaystring(verbose,'\n%s',['computing stokes velocities and pressure ...']);
Index: /issm/trunk/src/m/solutions/jpl/gradjcompute_core.m
===================================================================
--- /issm/trunk/src/m/solutions/jpl/gradjcompute_core.m	(revision 3961)
+++ /issm/trunk/src/m/solutions/jpl/gradjcompute_core.m	(revision 3962)
@@ -24,5 +24,5 @@
 
 %Reduce adjoint load from g-set to f-set
-[Du_f] = Reduceloadfromgtof( Du_g, m.Gmn, K_fs0, m.ys0, m.nodesets);
+[Du_f] = Reduceloadfromgtof( Du_g, m.Gmn, K_fs0, m.ys, m.nodesets,1);% 1 because we want ys0
 
 %Solve for adjoint vector: 
@@ -31,5 +31,5 @@
 
 %Merge back to g set
-lambda_g= Mergesolutionfromftog( lambda_f, m.Gmn, m.ys0, m.nodesets ); 
+lambda_g= Mergesolutionfromftog( lambda_f, m.Gmn, m.ys, m.nodesets,1); % 1 because we want ys0
 inputs=add(inputs,'adjoint',lambda_g,'doublevec',m.parameters.NumberOfDofsPerNode,m.parameters.NumberOfNodes);