Index: /issm/trunk-jpl/src/m/classes/frictionjosh.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/frictionjosh.m	(revision 24561)
+++ /issm/trunk-jpl/src/m/classes/frictionjosh.m	(revision 24562)
@@ -41,5 +41,5 @@
 
 			%Check that temperature is provided
-			md = checkfield(md,'fieldname','initialization.temperature','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
+			md = checkfield(md,'fieldname','initialization.temperature','NaN',1,'Inf',1,'size','universal');
 		end % }}}
 		function disp(self) % {{{
Index: /issm/trunk-jpl/src/m/classes/initialization.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/initialization.m	(revision 24561)
+++ /issm/trunk-jpl/src/m/classes/initialization.m	(revision 24562)
@@ -75,9 +75,9 @@
 				end
 				md = checkfield(md,'fieldname','initialization.pressure','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
-				md = checkfield(md,'fieldname','initialization.temperature','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
+				md = checkfield(md,'fieldname','initialization.temperature','NaN',1,'Inf',1,'size','universal');
 			end
 			if (ismember('EnthalpyAnalysis',analyses) & md.thermal.isenthalpy)
-				md = checkfield(md,'fieldname','initialization.waterfraction','>=',0,'size',[md.mesh.numberofvertices 1]);
-				md = checkfield(md,'fieldname','initialization.watercolumn'  ,'>=',0,'size',[md.mesh.numberofvertices 1]);
+				md = checkfield(md,'fieldname','initialization.waterfraction','>=',0,'size','universal');
+				md = checkfield(md,'fieldname','initialization.watercolumn'  ,'>=',0,'size','universal');
 				pos=find(md.initialization.waterfraction>0.);
 				if(~isempty(pos)),
Index: /issm/trunk-jpl/src/m/classes/thermal.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/thermal.m	(revision 24561)
+++ /issm/trunk-jpl/src/m/classes/thermal.m	(revision 24562)
@@ -93,5 +93,5 @@
 				md = checkfield(md,'fieldname','thermal.watercolumn_upperlimit','>=',0);
 
-				%Make sure the spc are less than melting point
+				%Make sure the spc are less than melting point (Josh commented out the next 2 lines)
 				TEMP=md.thermal.spctemperature(1:md.mesh.numberofvertices,:);
 				replicate=repmat(md.geometry.surface-md.mesh.z,1,size(md.thermal.spctemperature,2));