Index: /issm/trunk-jpl/src/m/classes/qmu/normal_uncertain.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/qmu/normal_uncertain.m	(revision 25024)
+++ /issm/trunk-jpl/src/m/classes/qmu/normal_uncertain.m	(revision 25025)
@@ -18,5 +18,5 @@
 		stddev		= NaN;
 		partition	= [];
-		nsteps        = 0;
+		nsteps		= 0;
 	end
 	methods
@@ -24,5 +24,5 @@
 
 			%recover options:
-			options = pairoptions(varargin{:});
+			options=pairoptions(varargin{:});
 
 			%initialize fields:
@@ -106,5 +106,5 @@
 		%virtual functions needed by qmu processing algorithms:
 		%implemented:
-		function [desc]  =prop_desc(nuv,dstr) % {{{
+		function [desc]=prop_desc(nuv,dstr) % {{{
 			desc=cell(1,numel(nuv));
 			for i=1:numel(nuv)
@@ -121,5 +121,5 @@
 			desc=allempty(desc);
 		end %}}}
-		function [mean]  =prop_mean(nuv) % {{{
+		function [mean]=prop_mean(nuv) % {{{
 			mean=zeros(1,numel(nuv));
 			for i=1:numel(nuv)
@@ -134,11 +134,11 @@
 		end % }}}
 		%default
-		function [abscissas] =prop_abscissas(hbu) % {{{
+		function [abscissas]=prop_abscissas(hbu) % {{{
 			abscissas=[];
 		end % }}}
-		function [counts] =prop_counts(hbu) % {{{
+		function [counts]=prop_counts(hbu) % {{{
 			counts=[];
 		end % }}}
-		function [pairs_per_variable] =prop_pairs_per_variable(hbu) % {{{
+		function [pairs_per_variable]=prop_pairs_per_variable(hbu) % {{{
 			pairs_per_variable=[];
 		end % }}}
@@ -146,8 +146,8 @@
 			initpt=[];
 		end % }}}
-		function [lower]  =prop_lower(nuv) % {{{
+		function [lower]=prop_lower(nuv) % {{{
 			lower=[];
 		end % }}}
-		function [upper]  =prop_upper(nuv) % {{{
+		function [upper]=prop_upper(nuv) % {{{
 			upper=[];
 		end % }}}
@@ -155,12 +155,12 @@
 			initst=[];
 		end % }}}
-		function [stype] =prop_stype(nuv) % {{{
+		function [stype]=prop_stype(nuv) % {{{
 			stype={};
 		end % }}}
-		function [scale] =prop_scale(nuv) % {{{
+		function [scale]=prop_scale(nuv) % {{{
 			scale=[];
 		end % }}}
 		%new methods:
-		function scaled =isscaled(self) % {{{
+		function scaled=isscaled(self) % {{{
 			if strncmp(self.descriptor,'scaled_',7),
 				scaled=1;
Index: /issm/trunk-jpl/src/m/classes/qmu/uniform_uncertain.m
===================================================================
--- /issm/trunk-jpl/src/m/classes/qmu/uniform_uncertain.m	(revision 25024)
+++ /issm/trunk-jpl/src/m/classes/qmu/uniform_uncertain.m	(revision 25025)
@@ -18,4 +18,5 @@
 		upper		= Inf;
 		partition	= [];
+		nsteps		= 0;
 	end
 	methods
@@ -23,5 +24,5 @@
 
 			%recover options: 
-			options = pairoptions(varargin{:});
+			options=pairoptions(varargin{:});
 
 			%initialize fields: 
@@ -34,11 +35,18 @@
 			%partitions as the lower and upper vectors:
 			if self.isscaled(),
-				self.partition=getfieldvalue(options,'partition'); 
-				npart=partition_npart(self.partition);
-				if npart~=length(self.upper),
-					error(['uniform_uncertain constructor: for the scaled variable' self.descriptor ' the upper field is not currently a vector of values for all the partitions described in the partition vector']);
+				self.partition=getfieldvalue(options,'partition');
+				self.nsteps=getfieldvalue(options,'nsteps',1);
+				npart=qmupart2npart(self.partition);
+				if npart~=length(self.upper,1),
+					error(['uniform_uncertain constructor: for the scaled variable' self.descriptor ' the row size of the upper field should be identical to the number of partitions']);
 				end
-				if npart~=length(self.lower),
-					error(['uniform_uncertain constructor: for the scaled variable' self.descriptor ' the lower field is not currently a vector of values for all the partitions described in the partition vector']);
+				if npart~=length(self.lower,1),
+					error(['uniform_uncertain constructor: for the scaled variable' self.descriptor ' the row size of the lower field should be identical to the number of partitions']);
+				end
+				if self.nsteps~=size(self.upper,2),
+					error(['uniform_uncertain constructor: for the scaled variable ' self.descriptor ' the col size of the upper field should be identical to the number of time steps']);
+				end
+				if self.nsteps~=size(self.lower,2),
+					error(['uniform_uncertain constructor: for the scaled variable ' self.descriptor ' the col size of the lower field should be identical to the number of time steps']);
 				end
 			end
@@ -54,4 +62,5 @@
 				fielddisplay(self,'partition','partition vector defining where sampling will occur');
 			end
+			fielddisplay(self,'nsteps','number of time steps');
 		end 
 		%}}}
@@ -65,9 +74,16 @@
 				end
 				%better have a partition vector that has as many partitions as upper and lower's size: 
-				if length(self.upper)~=partition_npart(self.partition),
-					error('uniform_uncertain error message: upper and partition should be vectors of identical size');
+				if length(self.upper,1)~=partition_npart(self.partition),
+					error('uniform_uncertain error message: row size of upper and partition size should be identical');
 				end
-				if length(self.lower)~=partition_npart(self.partition),
-					error('uniform_uncertain error message: lower and partition should be vectors of identical size');
+				if length(self.lower,1)~=partition_npart(self.partition),
+					error('uniform_uncertain error message: row size of lower and partition size should be identical');
+				end
+				%we need as steps in stddev and mean as there are time steps: 
+				if size(self.upper,2)~=self.nsteps,
+					error('uniform_uncertain error message: col size of upper and number of time steps should be identical');
+				end
+				if size(self.lower,2)~=self.nsteps,
+					error('uniform_uncertain error message: col size of lower and number of time steps should be identical');
 				end
 
@@ -90,5 +106,5 @@
 		%virtual functions needed by qmu processing algorithms:
 		%implemented:
-		function [desc] = prop_desc(uuv,dstr) % {{{
+		function [desc]=prop_desc(uuv,dstr) % {{{
 			desc=cell(1,numel(uuv));
 			for i=1:numel(uuv)
@@ -105,5 +121,5 @@
 			desc=allempty(desc);
 		end %}}}
-		function [lower] = prop_lower(uuv) % {{{
+		function [lower]=prop_lower(uuv) % {{{
 			lower=zeros(1,numel(uuv));
 			for i=1:numel(uuv)
@@ -112,5 +128,5 @@
 			lower=allequal(lower,-Inf);
 		end %}}}
-		function [upper] = prop_upper(uuv) % {{{
+		function [upper]=prop_upper(uuv) % {{{
 			upper=zeros(1,numel(uuv));
 			for i=1:numel(uuv)
@@ -123,5 +139,5 @@
 			stddev=[];
 		end % }}}
-		function [mean]  =prop_mean(nuv) % {{{
+		function [mean]=prop_mean(nuv) % {{{
 			mean=[];
 		end % }}}
@@ -132,21 +148,21 @@
 			initst=[];
 		end %}}}
-		function [stype] =prop_stype(uuv) %{{{
+		function [stype]=prop_stype(uuv) %{{{
 			stype={};
 		end %}}}
-		function [scale] =prop_scale(uuv) %{{{
+		function [scale]=prop_scale(uuv) %{{{
 			scale=[];
 		end %}}}
-		function [abscissas] =prop_abscissas(hbu) % {{{
+		function [abscissas]=prop_abscissas(hbu) % {{{
 			abscissas=[]; 
 		end % }}}
-		function [counts] =prop_counts(hbu) % {{{
+		function [counts]=prop_counts(hbu) % {{{
 			counts=[]; 
 		end % }}}
-		function [pairs_per_variable] =prop_pairs_per_variable(hbu) % {{{
+		function [pairs_per_variable]=prop_pairs_per_variable(hbu) % {{{
 			pairs_per_variable=[];
 		end % }}}
 		%new methods: 
-		function scaled =isscaled(self) % {{{
+		function scaled=isscaled(self) % {{{
 			if strncmp(self.descriptor,'scaled_',7),
 				scaled=1;
Index: /issm/trunk-jpl/src/m/qmu/setupdesign/QmuSetupVariables.m
===================================================================
--- /issm/trunk-jpl/src/m/qmu/setupdesign/QmuSetupVariables.m	(revision 25024)
+++ /issm/trunk-jpl/src/m/qmu/setupdesign/QmuSetupVariables.m	(revision 25025)
@@ -17,7 +17,13 @@
 
 	if isa(variables,'uniform_uncertain'),
-		nlower=length(variables.lower); nupper=length(variables.upper);
+		nlower=length(variables.lower,1);
+		nupper=length(variables.upper,1);
 		if (nlower ~= npart || nupper ~=npart),
-			error('QmuSetupVariables error message: upper and lower fields should be same size as the number of partitions');
+			error('QmuSetupVariables error message: upper and lower fields should have the same number of rows as the number of partitions');
+		end
+		nlower=length(variables.lower,2);
+		nupper=length(variables.upper,2);
+		if (nstddev ~= nt || nmean ~=nt),
+			error('QmuSetupVariables error message: upper and lower fields should have the same number of cols as the number of time steps');
 		end
 	elseif isa(variables,'normal_uncertain'),