Index: /issm/trunk-jpl/src/m/classes/initialization.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/initialization.m	(revision 26006)
+++ /issm/trunk-jpl/src/m/classes/initialization.m	(revision 26007)
@@ -20,5 +20,5 @@
 		hydraulic_potential = NaN;
 		channelarea         = NaN;
-        sample              = NaN;
+        %sample              = NaN;
 	end
 	methods
@@ -112,9 +112,9 @@
                 end
             end
-            if ismember('SamplingAnalysis',analyses) & ~(strcmp(solution,'TransientSolution') & md.transient.issampling == 0),
-                if ~isnan(md.initialization.sample)
-                    md = checkfield(md,'fieldname','initialization.sample','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
-                end
-			end
+%             if ismember('SamplingAnalysis',analyses) & ~(strcmp(solution,'TransientSolution') & md.transient.issampling == 0),
+%                 if ~isnan(md.initialization.sample)
+%                     md = checkfield(md,'fieldname','initialization.sample','NaN',1,'Inf',1,'size',[md.mesh.numberofvertices 1]);
+%                 end
+% 			end
 		end % }}}
 		function disp(self) % {{{
@@ -135,5 +135,5 @@
 			fielddisplay(self,'hydraulic_potential','Hydraulic potential (for GlaDS) [Pa]');
 			fielddisplay(self,'channelarea','subglacial water channel area (for GlaDS) [m2]');
-            fielddisplay(self,'sample','Realization of a Gaussian random field');
+            %fielddisplay(self,'sample','Realization of a Gaussian random field');
 		end % }}}
 		function marshall(self,prefix,md,fid) % {{{
@@ -153,5 +153,5 @@
 			WriteData(fid,prefix,'object',self,'fieldname','channelarea','format','DoubleMat','mattype',1);
 			WriteData(fid,prefix,'object',self,'fieldname','hydraulic_potential','format','DoubleMat','mattype',1);
-            WriteData(fid,prefix,'object',self,'fieldname','sample','format','DoubleMat','mattype',1);
+            %WriteData(fid,prefix,'object',self,'fieldname','sample','format','DoubleMat','mattype',1);
             
 			if md.thermal.isenthalpy,
@@ -184,5 +184,5 @@
 			writejs1Darray(fid,[modelname '.initialization.hydraulic_potential'],self.hydraulic_potential);
 			writejs1Darray(fid,[modelname '.initialization.channel'],self.channelarea);
-            writejs1Darray(fid,[modelname '.initialization.sample'],self.sample);
+           % writejs1Darray(fid,[modelname '.initialization.sample'],self.sample);
             
 		end % }}}
Index: /issm/trunk-jpl/src/m/classes/transient.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/transient.m	(revision 26006)
+++ /issm/trunk-jpl/src/m/classes/transient.m	(revision 26007)
@@ -16,5 +16,5 @@
 		ismovingfront     = 0;
 		ishydrology       = 0;
-        issampling        = 0;
+        %issampling        = 0;
 		isslc             = 0;
 		iscoupler         = 0;
@@ -45,5 +45,5 @@
 			self.ismovingfront   =0;
 			self.ishydrology     = 0;
-            self.issampling      = 0;
+            %self.issampling      = 0;
 			self.isslc           = 0;
 			self.isoceancoupling = 0;
@@ -67,5 +67,5 @@
 			self.ismovingfront   = 0;
 			self.ishydrology     = 0;
-            self.issampling      = 0;
+            %self.issampling      = 0;
 			self.isslc           = 0;
 			self.isoceancoupling = 0;
@@ -102,5 +102,5 @@
 			md = checkfield(md,'fieldname','transient.isoceancoupling','numel',[1],'values',[0 1]);
 			md = checkfield(md,'fieldname','transient.iscoupler','numel',[1],'values',[0 1]);
-            md = checkfield(md,'fieldname','transient.issampling','numel',[1],'values',[0 1]);  
+            %md = checkfield(md,'fieldname','transient.issampling','numel',[1],'values',[0 1]);  
 			md = checkfield(md,'fieldname','transient.amr_frequency','numel',[1],'>=',0,'NaN',1,'Inf',1);
 
@@ -125,5 +125,5 @@
 			fielddisplay(self,'ismovingfront','indicates whether a moving front capability is used in the transient');
 			fielddisplay(self,'ishydrology','indicates whether an hydrology model is used');
-            fielddisplay(self,'issampling','indicates whether sampling is used in the transient')
+            %fielddisplay(self,'issampling','indicates whether sampling is used in the transient')
 			fielddisplay(self,'isslc','indicates whether a sea-level change solution is used in the transient');
 			fielddisplay(self,'isoceancoupling','indicates whether a coupling with an ocean model is used in the transient');
@@ -144,5 +144,5 @@
 			WriteData(fid,prefix,'object',self,'fieldname','ishydrology','format','Boolean');
 			WriteData(fid,prefix,'object',self,'fieldname','ismovingfront','format','Boolean');
-            WriteData(fid,prefix,'object',self,'fieldname','issampling','format','Boolean'); 
+            %WriteData(fid,prefix,'object',self,'fieldname','issampling','format','Boolean'); 
 			WriteData(fid,prefix,'object',self,'fieldname','isslc','format','Boolean');
 			WriteData(fid,prefix,'object',self,'fieldname','isoceancoupling','format','Boolean');
@@ -171,5 +171,5 @@
 			writejsdouble(fid,[modelname '.trans.ismovingfront'],self.ismovingfront);
 			writejsdouble(fid,[modelname '.trans.ishydrology'],self.ishydrology);
-            writejsdouble(fid,[modelname '.trans.issampling'],self.issampling); 
+            %writejsdouble(fid,[modelname '.trans.issampling'],self.issampling); 
 			writejsdouble(fid,[modelname '.trans.isslc'],self.isslc);
 			writejsdouble(fid,[modelname '.trans.isoceancoupling'],self.isoceancoupling);