Index: /issm/trunk-jpl/src/m/classes/initialization.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/initialization.m	(revision 16430)
+++ /issm/trunk-jpl/src/m/classes/initialization.m	(revision 16431)
@@ -51,5 +51,7 @@
 				md = checkfield(md,'initialization.vx','NaN',1,'size',[md.mesh.numberofvertices 1]);
 				md = checkfield(md,'initialization.vy','NaN',1,'size',[md.mesh.numberofvertices 1]);
-				md = checkfield(md,'initialization.vz','NaN',1,'size',[md.mesh.numberofvertices 1]);
+				if strcmp(meshtype(md.mesh),'3D')
+					md = checkfield(md,'initialization.vz','NaN',1,'size',[md.mesh.numberofvertices 1]);
+				end
 				md = checkfield(md,'initialization.pressure','NaN',1,'size',[md.mesh.numberofvertices 1]);
 			end
Index: /issm/trunk-jpl/src/m/classes/initialization.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/initialization.py	(revision 16430)
+++ /issm/trunk-jpl/src/m/classes/initialization.py	(revision 16431)
@@ -65,5 +65,6 @@
 			md = checkfield(md,'initialization.vx','NaN',1,'size',[md.mesh.numberofvertices])
 			md = checkfield(md,'initialization.vy','NaN',1,'size',[md.mesh.numberofvertices])
-			md = checkfield(md,'initialization.vz','NaN',1,'size',[md.mesh.numberofvertices])
+			if strcmp(md.mesh.meshtype(),'3D'):
+				md = checkfield(md,'initialization.vz','NaN',1,'size',[md.mesh.numberofvertices])
 			md = checkfield(md,'initialization.pressure','NaN',1,'size',[md.mesh.numberofvertices])
 		if (EnthalpyAnalysisEnum() in analyses and md.thermal.isenthalpy) or solution==EnthalpySolutionEnum():