Index: /issm/trunk/test/NightlyRun/test235.m
===================================================================
--- /issm/trunk/test/NightlyRun/test235.m	(revision 5497)
+++ /issm/trunk/test/NightlyRun/test235.m	(revision 5498)
@@ -61,5 +61,5 @@
 %Dakota options
 %variables
-md.variables.rheology_B=normal_uncertain('RheologyBbar',1,.5);
+md.variables.rheology_B=normal_uncertain('scaled_RheologyBbar',1,.5);
 
 %responses
@@ -74,5 +74,4 @@
 md.qmu_params.evaluation_concurrency=1;
 md.qmu_params.interval_type='forward';
-md.qmu_params.tabular_graphics_data=true;
 
 %imperative! 
@@ -85,5 +84,5 @@
 %Fields and tolerances to track changes
 md=tres(md,'dakota');
-md.results.dakota.importancefactors=importancefactors(md,'RheologyBbar','MaxVel')';
+md.results.dakota.importancefactors=importancefactors(md,'scaled_RheologyBbar','MaxVel')';
 field_names     ={'importancefactors'};
 field_tolerances={1e-11};
Index: /issm/trunk/test/NightlyRun/test236.m
===================================================================
--- /issm/trunk/test/NightlyRun/test236.m	(revision 5497)
+++ /issm/trunk/test/NightlyRun/test236.m	(revision 5498)
@@ -61,5 +61,5 @@
 %Dakota options
 %variables
-md.variables.rheology_B=normal_uncertain('RheologyBbar',1,.5);
+md.variables.rheology_B=normal_uncertain('scaled_RheologyBbar',1,.05);
 
 %responses
@@ -74,5 +74,4 @@
 md.qmu_params.evaluation_concurrency=1;
 md.qmu_params.interval_type='forward';
-md.qmu_params.tabular_graphics_data=true;
 
 %imperative! 
@@ -85,5 +84,5 @@
 %Fields and tolerances to track changes
 md=tres(md,'dakota');
-md.results.dakota.importancefactors=importancefactors(md,'RheologyBbar','MaxVel')';
+md.results.dakota.importancefactors=importancefactors(md,'scaled_RheologyBbar','MaxVel')';
 field_names     ={'importancefactors'};
 field_tolerances={1e-13};
Index: /issm/trunk/test/NightlyRun/test423.m
===================================================================
--- /issm/trunk/test/NightlyRun/test423.m	(revision 5497)
+++ /issm/trunk/test/NightlyRun/test423.m	(revision 5498)
@@ -14,7 +14,6 @@
 %variables
 md.variables.rho_ice=normal_uncertain('RhoIce',md.rho_ice,0.01);
-
-%md.variables.drag_coefficient=normal_uncertain('DragCoefficient',1,.01);
-md.variables.drag_coefficient=normal_uncertain('DragCoefficient',1,.01*ones(md.npart,1));
+md.variables.drag_coefficient=normal_uncertain('scaled_DragCoefficient',1,.01);
+%md.variables.drag_coefficient=normal_uncertain('scaled_DragCoefficient',1,.01*ones(md.npart,1));
 
 %responses
@@ -29,5 +28,4 @@
 md.qmu_params.evaluation_concurrency=1;
 md.qmu_params.interval_type='forward';
-md.qmu_params.tabular_graphics_data=true;
 
 
@@ -40,5 +38,5 @@
 %Fields and tolerances to track changes
 md=tres(md,'dakota');
-md.results.dakota.importancefactors=importancefactors(md,'DragCoefficient','MaxVel')';
+md.results.dakota.importancefactors=importancefactors(md,'scaled_DragCoefficient','MaxVel')';
 field_names     ={'importancefactors'};
 field_tolerances={1e-11};
Index: /issm/trunk/test/NightlyRun/test424.m
===================================================================
--- /issm/trunk/test/NightlyRun/test424.m	(revision 5497)
+++ /issm/trunk/test/NightlyRun/test424.m	(revision 5498)
@@ -15,5 +15,5 @@
 %variables
 md.variables.rho_ice=normal_uncertain('RhoIce',md.rho_ice,0.01);
-md.variables.drag_coefficient=normal_uncertain('DragCoefficient',1,0.01);
+md.variables.drag_coefficient=normal_uncertain('scaled_DragCoefficient',1,0.01);
 
 %responses
@@ -28,5 +28,4 @@
 md.qmu_params.evaluation_concurrency=1;
 md.qmu_params.interval_type='forward';
-md.qmu_params.tabular_graphics_data=true;
 
 %imperative! 
@@ -38,5 +37,5 @@
 %Fields and tolerances to track changes
 md=tres(md,'dakota');
-md.results.dakota.importancefactors=importancefactors(md,'DragCoefficient','MaxVel')';
+md.results.dakota.importancefactors=importancefactors(md,'scaled_DragCoefficient','MaxVel')';
 field_names     ={'importancefactors'};
 field_tolerances={1e-10};
Index: /issm/trunk/test/NightlyRun/test425.m
===================================================================
--- /issm/trunk/test/NightlyRun/test425.m	(revision 5497)
+++ /issm/trunk/test/NightlyRun/test425.m	(revision 5498)
@@ -14,5 +14,5 @@
 %variables
 md.variables.rho_ice=normal_uncertain('RhoIce',md.rho_ice,0.01);
-md.variables.drag_coefficient=normal_uncertain('DragCoefficient',1,0.01);
+md.variables.drag_coefficient=normal_uncertain('scaled_DragCoefficient',1,0.01);
 
 %responses
@@ -27,5 +27,4 @@
 md.qmu_params.evaluation_concurrency=1;
 md.qmu_params.interval_type='forward';
-md.qmu_params.tabular_graphics_data=true;
 
 
@@ -39,5 +38,5 @@
 %Fields and tolerances to track changes
 md=tres(md,'dakota');
-md.results.dakota.importancefactors=importancefactors(md,'DragCoefficient','MaxVel')';
+md.results.dakota.importancefactors=importancefactors(md,'scaled_DragCoefficient','MaxVel')';
 field_names     ={'importancefactors'};
 field_tolerances={1e-11};
Index: /issm/trunk/test/NightlyRun/test426.m
===================================================================
--- /issm/trunk/test/NightlyRun/test426.m	(revision 5497)
+++ /issm/trunk/test/NightlyRun/test426.m	(revision 5498)
@@ -14,5 +14,5 @@
 %variables
 md.variables.rho_ice=normal_uncertain('RhoIce',md.rho_ice,0.01);
-md.variables.drag_coefficient=normal_uncertain('DragCoefficient',1,0.01);
+md.variables.drag_coefficient=normal_uncertain('scaled_DragCoefficient',1,0.01);
 
 %responses
@@ -27,5 +27,4 @@
 md.qmu_params.evaluation_concurrency=1;
 md.qmu_params.interval_type='forward';
-md.qmu_params.tabular_graphics_data=true;
 
 
@@ -39,5 +38,5 @@
 %Fields and tolerances to track changes
 md=tres(md,'dakota');
-md.results.dakota.importancefactors=importancefactors(md,'DragCoefficient','MaxVel')';
+md.results.dakota.importancefactors=importancefactors(md,'scaled_DragCoefficient','MaxVel')';
 field_names     ={'importancefactors'};
 field_tolerances={1e-10};
Index: /issm/trunk/test/NightlyRun/test427.m
===================================================================
--- /issm/trunk/test/NightlyRun/test427.m	(revision 5497)
+++ /issm/trunk/test/NightlyRun/test427.m	(revision 5498)
@@ -21,14 +21,14 @@
 
 %variables
-md.variables.drag_coefficient=normal_uncertain('DragCoefficient',1,0.01);
+md.variables.drag_coefficient=normal_uncertain('scaled_DragCoefficient',1,0.01);
 
 %responses
 md.responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
-md.responses.MassFlux1=response_function('MassFlux1',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
-md.responses.MassFlux2=response_function('MassFlux2',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
-md.responses.MassFlux3=response_function('MassFlux3',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
-md.responses.MassFlux4=response_function('MassFlux4',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
-md.responses.MassFlux5=response_function('MassFlux5',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
-md.responses.massFlux6=response_function('MassFlux6',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
+md.responses.MassFlux1=response_function('indexed_MassFlux_1',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
+md.responses.MassFlux2=response_function('indexed_MassFlux_2',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
+md.responses.MassFlux3=response_function('indexed_MassFlux_3',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
+md.responses.MassFlux4=response_function('indexed_MassFlux_4',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
+md.responses.MassFlux5=response_function('indexed_MassFlux_5',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
+md.responses.massFlux6=response_function('indexed_MassFlux_6',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
 
 %mass flux profiles
@@ -44,5 +44,4 @@
 md.qmu_params.evaluation_concurrency=1;
 md.qmu_params.interval_type='forward';
-md.qmu_params.tabular_graphics_data=true;
 md.qmu_analysis=1;
 md.eps_rel=10^-5; %tighten for qmu analyses
Index: /issm/trunk/test/NightlyRun/test428.m
===================================================================
--- /issm/trunk/test/NightlyRun/test428.m	(revision 5497)
+++ /issm/trunk/test/NightlyRun/test428.m	(revision 5498)
@@ -21,14 +21,14 @@
 
 %variables
-md.variables.drag_coefficient=normal_uncertain('DragCoefficient',1,0.01);
+md.variables.drag_coefficient=normal_uncertain('scaled_DragCoefficient',1,0.01);
 
 %responses
 md.responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
-md.responses.MassFlux1=response_function('MassFlux1',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
-md.responses.MassFlux2=response_function('MassFlux2',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
-md.responses.MassFlux3=response_function('MassFlux3',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
-md.responses.MassFlux4=response_function('MassFlux4',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
-md.responses.MassFlux5=response_function('MassFlux5',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
-md.responses.massFlux6=response_function('MassFlux6',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
+md.responses.MassFlux1=response_function('indexed_MassFlux_1',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
+md.responses.MassFlux2=response_function('indexed_MassFlux_2',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
+md.responses.MassFlux3=response_function('indexed_MassFlux_3',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
+md.responses.MassFlux4=response_function('indexed_MassFlux_4',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
+md.responses.MassFlux5=response_function('indexed_MassFlux_5',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
+md.responses.MassFlux6=response_function('indexed_MassFlux_6',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
 
 %mass flux profiles
@@ -44,5 +44,4 @@
 md.qmu_params.evaluation_concurrency=1;
 md.qmu_params.interval_type='forward';
-md.qmu_params.tabular_graphics_data=true;
 md.qmu_analysis=1;
 md.eps_rel=10^-5; %tighten for qmu analyses
Index: /issm/trunk/test/NightlyRun/test434.m
===================================================================
--- /issm/trunk/test/NightlyRun/test434.m	(revision 5497)
+++ /issm/trunk/test/NightlyRun/test434.m	(revision 5498)
@@ -14,21 +14,17 @@
 
 %Dakota options
-md.variables.nuv=normal_uncertain.empty();
-md.variables.nuv(end+1)=normal_uncertain('DragCoefficient',1,0.01);
+md.variables.drag_coefficient=normal_uncertain('scaled_DragCoefficient',1,0.01);
 
-md.responses=struct();
-md.responses.rf =response_function.empty();
-md.responses.rf (end+1)=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
-md.responses.rf (end+1)=response_function('MassFlux1',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
-md.responses.rf (end+1)=response_function('MassFlux2',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
-md.responses.rf (end+1)=response_function('MassFlux3',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
-md.responses.rf (end+1)=response_function('MassFlux4',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
-md.responses.rf (end+1)=response_function('MassFlux5',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
-md.responses.rf (end+1)=response_function('MassFlux6',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
+md.responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
+md.responses.MassFlux1=response_function('indexed_MassFlux_1',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
+md.responses.MassFlux2=response_function('indexed_MassFlux_2',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
+md.responses.MassFlux3=response_function('indexed_MassFlux_3',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
+md.responses.MassFlux4=response_function('indexed_MassFlux_4',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
+md.responses.MassFlux5=response_function('indexed_MassFlux_5',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
+md.responses.MassFlux6=response_function('indexed_MassFlux_6',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
 
 %mass flux profiles
 md.qmu_mass_flux_profiles={'../Exp/MassFlux1.exp','../Exp/MassFlux2.exp','../Exp/MassFlux3.exp','../Exp/MassFlux4.exp','../Exp/MassFlux5.exp','../Exp/MassFlux6.exp'};
 md.qmu_mass_flux_profile_directory=pwd;
-