Index: /issm/trunk/src/m/solutions/ice/diagnostic2d.m
===================================================================
--- /issm/trunk/src/m/solutions/ice/diagnostic2d.m	(revision 637)
+++ /issm/trunk/src/m/solutions/ice/diagnostic2d.m	(revision 638)
@@ -8,4 +8,7 @@
 %
 %   See also: ICEDIAGNOSTIC, ICEDIAGNOSTIC3D, ICEDIAGNOSTIC_CORE_NONLINEAR
+
+%global variables
+iceglobal
 
 %recover inputs
@@ -36,4 +39,9 @@
 	u_g=diagnostic_core_linear(m_dhu,'diagnostic_hutter',inputs);
 
+	if fem.ismacayealpattyn,
+		gridset=fem.m_dh.gridset;
+		fem.m_dh.ys=Reducevector_g(u_g);
+	end
+
 	%Computation of depth averaged pressure with MacAyeal's assumptions (P_bar=1/2*rho_ice*g*H in Pa)
 	u_g(4:6:m_dhu.gridset.gsize)=1/2*md.rho_ice*md.g*md.thickness;