Index: /issm/trunk-jpl/src/m/classes/calving.js
===================================================================
--- /issm/trunk-jpl/src/m/classes/calving.js	(revision 19902)
+++ /issm/trunk-jpl/src/m/classes/calving.js	(revision 19903)
@@ -28,5 +28,5 @@
 		if (solution!=TransientSolutionEnum() | md.trans.iscalving==0) return;
 
-		checkfield(md,'fieldname','calving.spclevelset','timeseries',1);
+		checkfield(md,'fieldname','calving.spclevelset','Inf',1,'timeseries',1);
 		checkfield(md,'fieldname','calving.stabilization','values',[0,1,2]);
 		checkfield(md,'fieldname','calving.calvingrate(1:md.mesh.numberofvertices,:)','>=',0,'timeseries',1,'NaN',1,'Inf',1);
Index: /issm/trunk-jpl/src/m/classes/calving.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/calving.m	(revision 19902)
+++ /issm/trunk-jpl/src/m/classes/calving.m	(revision 19903)
@@ -45,5 +45,5 @@
 			if (solution~=TransientSolutionEnum() | md.transient.iscalving==0), return; end
 
-			md = checkfield(md,'fieldname','calving.spclevelset','timeseries',1);
+			md = checkfield(md,'fieldname','calving.spclevelset','Inf',1,'timeseries',1);
 			md = checkfield(md,'fieldname','calving.stabilization','values',[0 1 2]);
 			md = checkfield(md,'fieldname','calving.calvingrate(1:md.mesh.numberofvertices,:)','>=',0,'timeseries',1,'NaN',1,'Inf',1);
Index: /issm/trunk-jpl/src/m/classes/calving.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/calving.py	(revision 19902)
+++ /issm/trunk-jpl/src/m/classes/calving.py	(revision 19903)
@@ -53,5 +53,5 @@
 			return md
 
-		md = checkfield(md,'fieldname','calving.spclevelset','timeseries',1)
+		md = checkfield(md,'fieldname','calving.spclevelset','Inf',1,'timeseries',1)
 		md = checkfield(md,'fieldname','calving.stabilization','values',[0,1,2]);
 		md = checkfield(md,'fieldname','calving.calvingrate','>=',0,'timeseries',1,'NaN',1,'Inf',1);
Index: /issm/trunk-jpl/src/m/classes/calvingdev.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/calvingdev.m	(revision 19902)
+++ /issm/trunk-jpl/src/m/classes/calvingdev.m	(revision 19903)
@@ -47,5 +47,5 @@
 			if (solution~=TransientSolutionEnum() | md.transient.iscalving==0), return; end
 
-			md = checkfield(md,'fieldname','calving.spclevelset','timeseries',1);
+			md = checkfield(md,'fieldname','calving.spclevelset','Inf',1,'timeseries',1);
 			md = checkfield(md,'fieldname','calving.stabilization','values',[0 1 2]);
 			md = checkfield(md,'fieldname','calving.coeff','>',0,'size',[md.mesh.numberofvertices 1]);
Index: /issm/trunk-jpl/src/m/classes/calvinglevermann.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/calvinglevermann.m	(revision 19902)
+++ /issm/trunk-jpl/src/m/classes/calvinglevermann.m	(revision 19903)
@@ -47,5 +47,5 @@
 			if (solution~=TransientSolutionEnum() | md.transient.iscalving==0), return; end
 
-			md = checkfield(md,'fieldname','calving.spclevelset','timeseries',1);
+			md = checkfield(md,'fieldname','calving.spclevelset','Inf',1,'timeseries',1);
 			md = checkfield(md,'fieldname','calving.stabilization','values',[0 1 2]);
 			md = checkfield(md,'fieldname','calving.coeff','>',0,'size',[md.mesh.numberofvertices 1]);
Index: /issm/trunk-jpl/src/m/classes/calvinglevermann.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/calvinglevermann.py	(revision 19902)
+++ /issm/trunk-jpl/src/m/classes/calvinglevermann.py	(revision 19903)
@@ -53,5 +53,5 @@
 			return md
 
-		md = checkfield(md,'fieldname','calving.spclevelset','timeseries',1)
+		md = checkfield(md,'fieldname','calving.spclevelset','Inf',1,'timeseries',1)
 		md = checkfield(md,'fieldname','calving.stabilization','values',[0,1,2]);
 		md = checkfield(md,'fieldname','calving.coeff','size',[md.mesh.numberofvertices],'>',0)
Index: /issm/trunk-jpl/src/m/classes/calvingpi.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/calvingpi.m	(revision 19902)
+++ /issm/trunk-jpl/src/m/classes/calvingpi.m	(revision 19903)
@@ -47,5 +47,5 @@
 			if (solution~=TransientSolutionEnum() | md.transient.iscalving==0), return; end
 
-			md = checkfield(md,'fieldname','calving.spclevelset','timeseries',1);
+			md = checkfield(md,'fieldname','calving.spclevelset','Inf',1,'timeseries',1);
 			md = checkfield(md,'fieldname','calving.stabilization','values',[0 1 2]);
 			md = checkfield(md,'fieldname','calving.coeff','>',0,'size',[md.mesh.numberofvertices 1]);
Index: /issm/trunk-jpl/src/m/classes/damage.js
===================================================================
--- /issm/trunk-jpl/src/m/classes/damage.js	(revision 19902)
+++ /issm/trunk-jpl/src/m/classes/damage.js	(revision 19903)
@@ -71,5 +71,5 @@
 				checkfield(md,'fieldname','damage.law','numel',[1],'values',[0,1,2]);
 				checkfield(md,'fieldname','damage.D','>=',0,'<=',self.max_damage,'size',[md.mesh.numberofvertices ,1]);
-				checkfield(md,'fieldname','damage.spcdamage','timeseries',1);
+				checkfield(md,'fieldname','damage.spcdamage','Inf',1,'timeseries',1);
 				checkfield(md,'fieldname','damage.max_damage','<',1,'>=',0);
 				checkfield(md,'fieldname','damage.stabilization','numel',[1],'values',[0, 1, 2, 4]);
Index: /issm/trunk-jpl/src/m/classes/damage.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/damage.m	(revision 19902)
+++ /issm/trunk-jpl/src/m/classes/damage.m	(revision 19903)
@@ -140,5 +140,5 @@
 				md = checkfield(md,'fieldname','damage.law','numel',[1],'values',[0,1,2]);
 				md = checkfield(md,'fieldname','damage.D','>=',0,'<=',self.max_damage,'size',[md.mesh.numberofvertices 1]);
-				md = checkfield(md,'fieldname','damage.spcdamage','timeseries',1);
+				md = checkfield(md,'fieldname','damage.spcdamage','Inf',1,'timeseries',1);
 				md = checkfield(md,'fieldname','damage.max_damage','<',1,'>=',0);
 				md = checkfield(md,'fieldname','damage.stabilization','numel',[1],'values',[0 1 2 4]);
Index: /issm/trunk-jpl/src/m/classes/damage.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/damage.py	(revision 19902)
+++ /issm/trunk-jpl/src/m/classes/damage.py	(revision 19903)
@@ -125,5 +125,5 @@
 			md = checkfield(md,'fieldname','damage.max_damage','<',1,'>=',0)
 			md = checkfield(md,'fieldname','damage.law','numel',[1],'values',[0,1,2,3])
-			md = checkfield(md,'fieldname','damage.spcdamage','timeseries',1)
+			md = checkfield(md,'fieldname','damage.spcdamage','Inf',1,'timeseries',1)
 			md = checkfield(md,'fieldname','damage.stabilization','numel',[1],'values',[0,1,2,4])
 			md = checkfield(md,'fieldname','damage.maxiter','>=0',0)
Index: /issm/trunk-jpl/src/m/classes/frictionwaterlayer.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/frictionwaterlayer.m	(revision 19902)
+++ /issm/trunk-jpl/src/m/classes/frictionwaterlayer.m	(revision 19903)
@@ -35,5 +35,5 @@
 			md = checkfield(md,'fieldname','friction.q','NaN',1,'Inf',1,'size',[md.mesh.numberofelements 1]);
 			md = checkfield(md,'fieldname','friction.p','NaN',1,'Inf',1,'size',[md.mesh.numberofelements 1]);
-			md = checkfield(md,'fieldname','thermal.spctemperature','timeseries',1,'>=',0.);
+			md = checkfield(md,'fieldname','thermal.spctemperature','Inf',1,'timeseries',1,'>=',0.);
 
 		end % }}}
Index: /issm/trunk-jpl/src/m/classes/hydrologydc.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/hydrologydc.m	(revision 19902)
+++ /issm/trunk-jpl/src/m/classes/hydrologydc.m	(revision 19903)
@@ -116,5 +116,5 @@
 			md = checkfield(md,'fieldname','hydrology.basal_moulin_input','NaN',1,'Inf',1,'timeseries',1);
 
-			md = checkfield(md,'fieldname','hydrology.spcsediment_head','timeseries',1);
+			md = checkfield(md,'fieldname','hydrology.spcsediment_head','Inf',1,'timeseries',1);
 			md = checkfield(md,'fieldname','hydrology.sediment_compressibility','>',0,'numel',1);
 			md = checkfield(md,'fieldname','hydrology.sediment_porosity','>',0,'numel',1);
@@ -123,5 +123,5 @@
 
 			if self.isefficientlayer==1,
-				md = checkfield(md,'fieldname','hydrology.spcepl_head','timeseries',1);
+				md = checkfield(md,'fieldname','hydrology.spcepl_head','Inf',1,'timeseries',1);
 				md = checkfield(md,'fieldname','hydrology.mask_eplactive_node','size',[md.mesh.numberofvertices 1],'values',[0 1]);
 				md = checkfield(md,'fieldname','hydrology.epl_compressibility','>',0,'numel',1);
Index: /issm/trunk-jpl/src/m/classes/hydrologydc.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/hydrologydc.py	(revision 19902)
+++ /issm/trunk-jpl/src/m/classes/hydrologydc.py	(revision 19903)
@@ -164,5 +164,5 @@
 
 		md = checkfield(md,'fieldname','hydrology.basal_moulin_input','NaN',1,'Inf',1,'timeseries',1)
-		md = checkfield(md,'fieldname','hydrology.spcsediment_head','timeseries',1)
+		md = checkfield(md,'fieldname','hydrology.spcsediment_head','Inf',1,'timeseries',1)
 		md = checkfield(md,'fieldname','hydrology.sediment_compressibility','>',0.,'numel',[1])
 		md = checkfield(md,'fieldname','hydrology.sediment_porosity','>',0.,'numel',[1])
@@ -170,5 +170,5 @@
 		md = checkfield(md,'fieldname','hydrology.sediment_transmitivity','>=',0,'size',[md.mesh.numberofvertices,1])
 		if self.isefficientlayer==1:
-			md = checkfield(md,'fieldname','hydrology.spcepl_head','timeseries',1)
+			md = checkfield(md,'fieldname','hydrology.spcepl_head','Inf',1,'timeseries',1)
 			md = checkfield(md,'fieldname','hydrology.mask_eplactive_node','size',[md.mesh.numberofvertices,1],'values',[0,1])
 			md = checkfield(md,'fieldname','hydrology.epl_compressibility','>',0.,'numel',[1])
Index: /issm/trunk-jpl/src/m/classes/hydrologyshreve.js
===================================================================
--- /issm/trunk-jpl/src/m/classes/hydrologyshreve.js	(revision 19902)
+++ /issm/trunk-jpl/src/m/classes/hydrologyshreve.js	(revision 19903)
@@ -32,5 +32,5 @@
 			if(!ArrayAnyEqual(ArrayIsMember(HydrologyShreveAnalysisEnum(),analyses),1)) return;
 
-			checkfield(md,'fieldname','hydrology.spcwatercolumn','timeseries',1);
+			checkfield(md,'fieldname','hydrology.spcwatercolumn','Inf',1,'timeseries',1);
 			checkfield(md,'fieldname','hydrology.stabilization','>=',0);
 
Index: /issm/trunk-jpl/src/m/classes/hydrologyshreve.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/hydrologyshreve.m	(revision 19902)
+++ /issm/trunk-jpl/src/m/classes/hydrologyshreve.m	(revision 19903)
@@ -45,5 +45,5 @@
 			end
 
-			md = checkfield(md,'fieldname','hydrology.spcwatercolumn','timeseries',1);
+			md = checkfield(md,'fieldname','hydrology.spcwatercolumn','Inf',1,'timeseries',1);
 			md = checkfield(md,'fieldname','hydrology.stabilization','>=',0);
 		end % }}}
Index: /issm/trunk-jpl/src/m/classes/hydrologyshreve.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/hydrologyshreve.py	(revision 19902)
+++ /issm/trunk-jpl/src/m/classes/hydrologyshreve.py	(revision 19903)
@@ -43,5 +43,5 @@
 			return md
 
-		md = checkfield(md,'fieldname','hydrology.spcwatercolumn','timeseries',1)
+		md = checkfield(md,'fieldname','hydrology.spcwatercolumn','Inf',1,'timeseries',1)
 		md = checkfield(md,'fieldname','hydrology.stabilization','>=',0)
 
Index: /issm/trunk-jpl/src/m/classes/masstransport.js
===================================================================
--- /issm/trunk-jpl/src/m/classes/masstransport.js	(revision 19902)
+++ /issm/trunk-jpl/src/m/classes/masstransport.js	(revision 19903)
@@ -46,5 +46,5 @@
 			if(!ArrayAnyEqual(ArrayIsMember(HydrologyShreveAnalysisEnum(),analyses),1) | (solution==TransientSolutionEnum() & md.trans.ismasstransport==0)) return; 
 
-			checkfield(md,'fieldname','masstransport.spcthickness','timeseries',1);
+			checkfield(md,'fieldname','masstransport.spcthickness','Inf',1,'timeseries',1);
 			checkfield(md,'fieldname','masstransport.isfreesurface','values',[0 ,1]);
 			checkfield(md,'fieldname','masstransport.hydrostatic_adjustment','values',['Absolute', 'Incremental']);
Index: /issm/trunk-jpl/src/m/classes/masstransport.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/masstransport.m	(revision 19902)
+++ /issm/trunk-jpl/src/m/classes/masstransport.m	(revision 19903)
@@ -89,5 +89,5 @@
 			if ~ismember(MasstransportAnalysisEnum(),analyses) |  (solution==TransientSolutionEnum() & md.transient.ismasstransport==0), return; end
 
-			md = checkfield(md,'fieldname','masstransport.spcthickness','timeseries',1);
+			md = checkfield(md,'fieldname','masstransport.spcthickness','Inf',1,'timeseries',1);
 			md = checkfield(md,'fieldname','masstransport.isfreesurface','values',[0 1]);
 			md = checkfield(md,'fieldname','masstransport.hydrostatic_adjustment','values',{'Absolute' 'Incremental'});
Index: /issm/trunk-jpl/src/m/classes/masstransport.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/masstransport.py	(revision 19902)
+++ /issm/trunk-jpl/src/m/classes/masstransport.py	(revision 19903)
@@ -72,5 +72,5 @@
 			return md
 
-		md = checkfield(md,'fieldname','masstransport.spcthickness','timeseries',1)
+		md = checkfield(md,'fieldname','masstransport.spcthickness','Inf',1,'timeseries',1)
 		md = checkfield(md,'fieldname','masstransport.isfreesurface','values',[0,1])
 		md = checkfield(md,'fieldname','masstransport.hydrostatic_adjustment','values',['Absolute','Incremental'])
Index: /issm/trunk-jpl/src/m/classes/stressbalance.js
===================================================================
--- /issm/trunk-jpl/src/m/classes/stressbalance.js	(revision 19902)
+++ /issm/trunk-jpl/src/m/classes/stressbalance.js	(revision 19903)
@@ -81,7 +81,7 @@
 			if(ArrayAnyEqual(ArrayIsMember(StressbalanceAnalysisEnum(),analyses),0))return;
 
-			checkfield(md,'fieldname','stressbalance.spcvx','timeseries',1);
-			checkfield(md,'fieldname','stressbalance.spcvy','timeseries',1);
-			checkfield(md,'fieldname','stressbalance.spcvz','timeseries',1);
+			checkfield(md,'fieldname','stressbalance.spcvx','Inf',1,'timeseries',1);
+			checkfield(md,'fieldname','stressbalance.spcvy','Inf',1,'timeseries',1);
+			checkfield(md,'fieldname','stressbalance.spcvz','Inf',1,'timeseries',1);
 			checkfield(md,'fieldname','stressbalance.restol','size',[1, 1],'>',0,'NaN',1,'Inf',1);
 			checkfield(md,'fieldname','stressbalance.reltol','size',[1, 1]);
Index: /issm/trunk-jpl/src/m/classes/stressbalance.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/stressbalance.m	(revision 19902)
+++ /issm/trunk-jpl/src/m/classes/stressbalance.m	(revision 19903)
@@ -134,7 +134,7 @@
 			if ~ismember(StressbalanceAnalysisEnum(),analyses), return; end
 
-			md = checkfield(md,'fieldname','stressbalance.spcvx','timeseries',1);
-			md = checkfield(md,'fieldname','stressbalance.spcvy','timeseries',1);
-			md = checkfield(md,'fieldname','stressbalance.spcvz','timeseries',1);
+			md = checkfield(md,'fieldname','stressbalance.spcvx','Inf',1,'timeseries',1);
+			md = checkfield(md,'fieldname','stressbalance.spcvy','Inf',1,'timeseries',1);
+			md = checkfield(md,'fieldname','stressbalance.spcvz','Inf',1,'timeseries',1);
 			md = checkfield(md,'fieldname','stressbalance.restol','size',[1 1],'>',0,'NaN',1,'Inf',1);
 			md = checkfield(md,'fieldname','stressbalance.reltol','size',[1 1]);
Index: /issm/trunk-jpl/src/m/classes/stressbalance.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/stressbalance.py	(revision 19902)
+++ /issm/trunk-jpl/src/m/classes/stressbalance.py	(revision 19903)
@@ -133,8 +133,8 @@
 			return md
 
-		md = checkfield(md,'fieldname','stressbalance.spcvx','timeseries',1)
-		md = checkfield(md,'fieldname','stressbalance.spcvy','timeseries',1)
+		md = checkfield(md,'fieldname','stressbalance.spcvx','Inf',1,'timeseries',1)
+		md = checkfield(md,'fieldname','stressbalance.spcvy','Inf',1,'timeseries',1)
 		if m.strcmp(md.mesh.domaintype(),'3D'):
-			md = checkfield(md,'fieldname','stressbalance.spcvz','timeseries',1)
+			md = checkfield(md,'fieldname','stressbalance.spcvz','Inf',1,'timeseries',1)
 		md = checkfield(md,'fieldname','stressbalance.restol','size',[1],'>',0)
 		md = checkfield(md,'fieldname','stressbalance.reltol','size',[1])
Index: /issm/trunk-jpl/src/m/classes/thermal.js
===================================================================
--- /issm/trunk-jpl/src/m/classes/thermal.js	(revision 19902)
+++ /issm/trunk-jpl/src/m/classes/thermal.js	(revision 19903)
@@ -58,5 +58,5 @@
 
 		checkfield(md,'fieldname','thermal.stabilization','numel',[1],'values',[0 ,1, 2]);
-		checkfield(md,'fieldname','thermal.spctemperature','timeseries',1);
+		checkfield(md,'fieldname','thermal.spctemperature','Inf',1,'timeseries',1);
 		if(ArrayAnyEqual(ArrayIsMember(EnthalpyAnalysisEnum(),analyses),1) & md.thermal.isenthalpy & md.mesh.dimension() == 3){
 			
Index: /issm/trunk-jpl/src/m/classes/thermal.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/thermal.m	(revision 19902)
+++ /issm/trunk-jpl/src/m/classes/thermal.m	(revision 19903)
@@ -99,5 +99,5 @@
 
 			md = checkfield(md,'fieldname','thermal.stabilization','numel',[1],'values',[0 1 2]);
-			md = checkfield(md,'fieldname','thermal.spctemperature','timeseries',1);
+			md = checkfield(md,'fieldname','thermal.spctemperature','Inf',1,'timeseries',1);
 			if (ismember(EnthalpyAnalysisEnum(),analyses) & md.thermal.isenthalpy & dimension(md.mesh)==3),
 				pos=find(md.thermal.spctemperature(1:md.mesh.numberofvertices,:)~=NaN);
Index: /issm/trunk-jpl/src/m/classes/thermal.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/thermal.py	(revision 19902)
+++ /issm/trunk-jpl/src/m/classes/thermal.py	(revision 19903)
@@ -95,5 +95,5 @@
 
 		md = checkfield(md,'fieldname','thermal.stabilization','numel',[1],'values',[0,1,2])
-		md = checkfield(md,'fieldname','thermal.spctemperature','timeseries',1)
+		md = checkfield(md,'fieldname','thermal.spctemperature','Inf',1,'timeseries',1)
 		if EnthalpyAnalysisEnum() in analyses and md.thermal.isenthalpy and md.mesh.dimension()==3:
 			pos=numpy.nonzero(numpy.logical_not(numpy.isnan(md.thermal.spctemperature[0:md.mesh.numberofvertices])))
Index: /issm/trunk-jpl/src/py3/classes/calving.py
===================================================================
--- /issm/trunk-jpl/src/py3/classes/calving.py	(revision 19902)
+++ /issm/trunk-jpl/src/py3/classes/calving.py	(revision 19903)
@@ -53,5 +53,5 @@
 			return md
 
-		md = checkfield(md,'fieldname','calving.spclevelset','timeseries',1)
+		md = checkfield(md,'fieldname','calving.spclevelset','Inf',1,'timeseries',1)
 		md = checkfield(md,'fieldname','calving.stabilization','values',[0,1,2]);
 		md = checkfield(md,'fieldname','calving.calvingrate','>=',0,'timeseries',1,'NaN',1,'Inf',1);
Index: /issm/trunk-jpl/src/py3/classes/calvinglevermann.py
===================================================================
--- /issm/trunk-jpl/src/py3/classes/calvinglevermann.py	(revision 19902)
+++ /issm/trunk-jpl/src/py3/classes/calvinglevermann.py	(revision 19903)
@@ -53,5 +53,5 @@
 			return md
 
-		md = checkfield(md,'fieldname','calving.spclevelset','timeseries',1)
+		md = checkfield(md,'fieldname','calving.spclevelset','Inf',1,'timeseries',1)
 		md = checkfield(md,'fieldname','calving.stabilization','values',[0,1,2]);
 		md = checkfield(md,'fieldname','calving.coeff','size',[md.mesh.numberofvertices],'>',0)
Index: /issm/trunk-jpl/src/py3/classes/damage.py
===================================================================
--- /issm/trunk-jpl/src/py3/classes/damage.py	(revision 19902)
+++ /issm/trunk-jpl/src/py3/classes/damage.py	(revision 19903)
@@ -124,5 +124,5 @@
 			md = checkfield(md,'fieldname','damage.max_damage','<',1,'>=',0)
 			md = checkfield(md,'fieldname','damage.law','numel',[1],'values',[0,1,2,3])
-			md = checkfield(md,'fieldname','damage.spcdamage','timeseries',1)
+			md = checkfield(md,'fieldname','damage.spcdamage','Inf',1,'timeseries',1)
 			md = checkfield(md,'fieldname','damage.stabilization','numel',[1],'values',[0,1,2,4])
 			md = checkfield(md,'fieldname','damage.maxiter','>=0',0)
Index: /issm/trunk-jpl/src/py3/classes/hydrologydc.py
===================================================================
--- /issm/trunk-jpl/src/py3/classes/hydrologydc.py	(revision 19902)
+++ /issm/trunk-jpl/src/py3/classes/hydrologydc.py	(revision 19903)
@@ -164,5 +164,5 @@
 
 		md = checkfield(md,'fieldname','hydrology.basal_moulin_input','NaN',1,'Inf',1,'timeseries',1)
-		md = checkfield(md,'fieldname','hydrology.spcsediment_head','timeseries',1)
+		md = checkfield(md,'fieldname','hydrology.spcsediment_head','Inf',1,'timeseries',1)
 		md = checkfield(md,'fieldname','hydrology.sediment_compressibility','>',0.,'numel',[1])
 		md = checkfield(md,'fieldname','hydrology.sediment_porosity','>',0.,'numel',[1])
@@ -170,5 +170,5 @@
 		md = checkfield(md,'fieldname','hydrology.sediment_transmitivity','>=',0,'size',[md.mesh.numberofvertices,1])
 		if self.isefficientlayer==1:
-			md = checkfield(md,'fieldname','hydrology.spcepl_head','timeseries',1)
+			md = checkfield(md,'fieldname','hydrology.spcepl_head','Inf',1,'timeseries',1)
 			md = checkfield(md,'fieldname','hydrology.mask_eplactive_node','size',[md.mesh.numberofvertices,1],'values',[0,1])
 			md = checkfield(md,'fieldname','hydrology.epl_compressibility','>',0.,'numel',[1])
Index: /issm/trunk-jpl/src/py3/classes/hydrologyshreve.py
===================================================================
--- /issm/trunk-jpl/src/py3/classes/hydrologyshreve.py	(revision 19902)
+++ /issm/trunk-jpl/src/py3/classes/hydrologyshreve.py	(revision 19903)
@@ -43,5 +43,5 @@
 			return md
 
-		md = checkfield(md,'fieldname','hydrology.spcwatercolumn','timeseries',1)
+		md = checkfield(md,'fieldname','hydrology.spcwatercolumn','Inf',1,'timeseries',1)
 		md = checkfield(md,'fieldname','hydrology.stabilization','>=',0)
 
Index: /issm/trunk-jpl/src/py3/classes/masstransport.py
===================================================================
--- /issm/trunk-jpl/src/py3/classes/masstransport.py	(revision 19902)
+++ /issm/trunk-jpl/src/py3/classes/masstransport.py	(revision 19903)
@@ -72,5 +72,5 @@
 			return md
 
-		md = checkfield(md,'fieldname','masstransport.spcthickness','timeseries',1)
+		md = checkfield(md,'fieldname','masstransport.spcthickness','Inf',1,'timeseries',1)
 		md = checkfield(md,'fieldname','masstransport.isfreesurface','values',[0,1])
 		md = checkfield(md,'fieldname','masstransport.hydrostatic_adjustment','values',['Absolute','Incremental'])
Index: /issm/trunk-jpl/src/py3/classes/stressbalance.py
===================================================================
--- /issm/trunk-jpl/src/py3/classes/stressbalance.py	(revision 19902)
+++ /issm/trunk-jpl/src/py3/classes/stressbalance.py	(revision 19903)
@@ -133,8 +133,8 @@
 			return md
 
-		md = checkfield(md,'fieldname','stressbalance.spcvx','timeseries',1)
-		md = checkfield(md,'fieldname','stressbalance.spcvy','timeseries',1)
+		md = checkfield(md,'fieldname','stressbalance.spcvx','Inf',1,'timeseries',1)
+		md = checkfield(md,'fieldname','stressbalance.spcvy','Inf',1,'timeseries',1)
 		if m.strcmp(md.mesh.domaintype(),'3D'):
-			md = checkfield(md,'fieldname','stressbalance.spcvz','timeseries',1)
+			md = checkfield(md,'fieldname','stressbalance.spcvz','Inf',1,'timeseries',1)
 		md = checkfield(md,'fieldname','stressbalance.restol','size',[1],'>',0)
 		md = checkfield(md,'fieldname','stressbalance.reltol','size',[1])
Index: /issm/trunk-jpl/src/py3/classes/thermal.py
===================================================================
--- /issm/trunk-jpl/src/py3/classes/thermal.py	(revision 19902)
+++ /issm/trunk-jpl/src/py3/classes/thermal.py	(revision 19903)
@@ -95,5 +95,5 @@
 
 		md = checkfield(md,'fieldname','thermal.stabilization','numel',[1],'values',[0,1,2])
-		md = checkfield(md,'fieldname','thermal.spctemperature','timeseries',1)
+		md = checkfield(md,'fieldname','thermal.spctemperature','Inf',1,'timeseries',1)
 		if EnthalpyAnalysisEnum() in analyses and md.thermal.isenthalpy and md.mesh.dimension()==3:
 			pos=numpy.nonzero(numpy.logical_not(numpy.isnan(md.thermal.spctemperature[0:md.mesh.numberofvertices])))
