Index: /issm/trunk/test/NightlyRun/test235.m
===================================================================
--- /issm/trunk/test/NightlyRun/test235.m	(revision 9651)
+++ /issm/trunk/test/NightlyRun/test235.m	(revision 9652)
@@ -56,27 +56,27 @@
 
 %partitioning
-md.npart=md.numberofnodes;
-md=partitioner(md,'package','linear','npart',md.npart);
-md.part=md.part-1;
+md.qmu.numberofpartitions=md.numberofnodes;
+md.qmu.partitionitioner(md,'package','linear','npart',md.qmu.numberofpartitions);
+md.qmu.partition=md.qmu.partition-1;
 
 %Dakota options
 %variables
-md.variables.rheology_B=normal_uncertain('scaled_MaterialsRheologyB',1,.5);
+md.qmu.variables.rheology_B=normal_uncertain('scaled_MaterialsRheologyB',1,.5);
 
 %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.qmu.responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
 
 %method
-md.qmu_method     =dakota_method('nond_l');
+md.qmu.method     =dakota_method('nond_l');
 
 %parameters
-md.qmu_params.direct=true;
-md.qmu_params.analysis_driver='diagnostic';
-md.qmu_params.evaluation_concurrency=1;
-md.qmu_params.interval_type='forward';
+md.qmu.params.direct=true;
+md.qmu.params.analysis_driver='diagnostic';
+md.qmu.params.evaluation_concurrency=1;
+md.qmu.params.interval_type='forward';
 
 %imperative! 
 md.eps_rel=10^-10; %tighten for qmu analysese
-md.dakota_analysis=1;
+md.qmu.isdakota=1;
 
 %solve
Index: /issm/trunk/test/NightlyRun/test236.m
===================================================================
--- /issm/trunk/test/NightlyRun/test236.m	(revision 9651)
+++ /issm/trunk/test/NightlyRun/test236.m	(revision 9652)
@@ -56,27 +56,27 @@
 
 %partitioning
-md.npart=md.numberofnodes;
-md=partitioner(md,'package','linear','npart',md.npart);
-md.part=md.part-1;
+md.qmu.numberofpartitions=md.numberofnodes;
+md.qmu.partitionitioner(md,'package','linear','npart',md.qmu.numberofpartitions);
+md.qmu.partition=md.qmu.partition-1;
 
 %Dakota options
 %variables
-md.variables.rheology_B=normal_uncertain('scaled_MaterialsRheologyB',1,.05);
+md.qmu.variables.rheology_B=normal_uncertain('scaled_MaterialsRheologyB',1,.05);
 
 %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.qmu.responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
 
 %method
-md.qmu_method     =dakota_method('nond_l');
+md.qmu.method     =dakota_method('nond_l');
 
 %parameters
-md.qmu_params.direct=true;
-md.qmu_params.analysis_driver='diagnostic';
-md.qmu_params.evaluation_concurrency=1;
-md.qmu_params.interval_type='forward';
+md.qmu.params.direct=true;
+md.qmu.params.analysis_driver='diagnostic';
+md.qmu.params.evaluation_concurrency=1;
+md.qmu.params.interval_type='forward';
 
 %imperative! 
 md.eps_rel=10^-10; %tighten for qmu analysese
-md.dakota_analysis=1;
+md.qmu.isdakota=1;
 
 %solve
Index: /issm/trunk/test/NightlyRun/test423.m
===================================================================
--- /issm/trunk/test/NightlyRun/test423.m	(revision 9651)
+++ /issm/trunk/test/NightlyRun/test423.m	(revision 9652)
@@ -6,26 +6,26 @@
 
 %partitioning
-md.npart=md.numberofnodes;
-md=partitioner(md,'package','linear','npart',md.npart);
-md.part=md.part-1;
-md.dakota_analysis=1;
+md.qmu.numberofpartitions=md.numberofnodes;
+md.qmu.partitionitioner(md,'package','linear','npart',md.qmu.numberofpartitions);
+md.qmu.partition=md.qmu.partition-1;
+md.qmu.isdakota=1;
 
 %Dakota options
 %variables
-md.variables.rho_ice=normal_uncertain('RhoIce',md.materials.rho_ice,0.01);
-md.variables.drag_coefficient=normal_uncertain('scaled_FrictionCoefficient',1,.01);
-%md.variables.drag_coefficient=normal_uncertain('scaled_FrictionCoefficient',1,.01*ones(md.npart,1));
+md.qmu.variables.rho_ice=normal_uncertain('RhoIce',md.materials.rho_ice,0.01);
+md.qmu.variables.drag_coefficient=normal_uncertain('scaled_FrictionCoefficient',1,.01);
+%md.qmu.variables.drag_coefficient=normal_uncertain('scaled_FrictionCoefficient',1,.01*ones(md.qmu.numberofpartitions,1));
 
 %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.qmu.responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
 
 %method
-md.qmu_method     =dakota_method('nond_l');
+md.qmu.method     =dakota_method('nond_l');
 
 %parameters
-md.qmu_params.direct=true;
-md.qmu_params.analysis_driver='diagnostic';
-md.qmu_params.evaluation_concurrency=1;
-md.qmu_params.interval_type='forward';
+md.qmu.params.direct=true;
+md.qmu.params.analysis_driver='diagnostic';
+md.qmu.params.evaluation_concurrency=1;
+md.qmu.params.interval_type='forward';
 
 
Index: /issm/trunk/test/NightlyRun/test424.m
===================================================================
--- /issm/trunk/test/NightlyRun/test424.m	(revision 9651)
+++ /issm/trunk/test/NightlyRun/test424.m	(revision 9652)
@@ -6,26 +6,26 @@
 
 %partitioning
-md.npart=md.numberofnodes;
-md=partitioner(md,'package','linear','npart',md.npart);
-md.part=md.part-1;
-md.dakota_analysis=1;
+md.qmu.numberofpartitions=md.numberofnodes;
+md.qmu.partitionitioner(md,'package','linear','npart',md.qmu.numberofpartitions);
+md.qmu.partition=md.qmu.partition-1;
+md.qmu.isdakota=1;
 
 
 %Dakota options
 %variables
-md.variables.rho_ice=normal_uncertain('RhoIce',md.materials.rho_ice,0.01);
-md.variables.drag_coefficient=normal_uncertain('scaled_FrictionCoefficient',1,0.01);
+md.qmu.variables.rho_ice=normal_uncertain('RhoIce',md.materials.rho_ice,0.01);
+md.qmu.variables.drag_coefficient=normal_uncertain('scaled_FrictionCoefficient',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.qmu.responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
 
 %method
-md.qmu_method     =dakota_method('nond_l');
+md.qmu.method     =dakota_method('nond_l');
 
 %parameters
-md.qmu_params.direct=true;
-md.qmu_params.analysis_driver='diagnostic';
-md.qmu_params.evaluation_concurrency=1;
-md.qmu_params.interval_type='forward';
+md.qmu.params.direct=true;
+md.qmu.params.analysis_driver='diagnostic';
+md.qmu.params.evaluation_concurrency=1;
+md.qmu.params.interval_type='forward';
 
 %imperative! 
Index: /issm/trunk/test/NightlyRun/test425.m
===================================================================
--- /issm/trunk/test/NightlyRun/test425.m	(revision 9651)
+++ /issm/trunk/test/NightlyRun/test425.m	(revision 9652)
@@ -8,28 +8,28 @@
 
 %partitioning
-md.npart=20;
-md=partitioner(md,'package','chaco','npart',md.npart,'weighting','on');
-md.part=md.part-1;
+md.qmu.numberofpartitions=20;
+md.qmu.partitionitioner(md,'package','chaco','npart',md.qmu.numberofpartitions,'weighting','on');
+md.qmu.partition=md.qmu.partition-1;
 
 %variables
-md.variables.rho_ice=normal_uncertain('RhoIce',md.materials.rho_ice,0.01);
-md.variables.drag_coefficient=normal_uncertain('scaled_FrictionCoefficient',1,0.01);
+md.qmu.variables.rho_ice=normal_uncertain('RhoIce',md.materials.rho_ice,0.01);
+md.qmu.variables.drag_coefficient=normal_uncertain('scaled_FrictionCoefficient',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.qmu.responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
 
 %method
-md.qmu_method     =dakota_method('nond_l');
+md.qmu.method     =dakota_method('nond_l');
 
 %parameters
-md.qmu_params.direct=true;
-md.qmu_params.analysis_driver='diagnostic';
-md.qmu_params.evaluation_concurrency=1;
-md.qmu_params.interval_type='forward';
+md.qmu.params.direct=true;
+md.qmu.params.analysis_driver='diagnostic';
+md.qmu.params.evaluation_concurrency=1;
+md.qmu.params.interval_type='forward';
 
 
 %imperative!
 md.eps_rel=10^-5; %tighten for qmu analyses
-md.dakota_analysis=1;
+md.qmu.isdakota=1;
 
 %solve
Index: /issm/trunk/test/NightlyRun/test426.m
===================================================================
--- /issm/trunk/test/NightlyRun/test426.m	(revision 9651)
+++ /issm/trunk/test/NightlyRun/test426.m	(revision 9652)
@@ -8,28 +8,28 @@
 
 %partitioning
-md.npart=20;
-md=partitioner(md,'package','chaco','npart',md.npart,'weighting','on');
-md.part=md.part-1;
+md.qmu.numberofpartitions=20;
+md.qmu.partitionitioner(md,'package','chaco','npart',md.qmu.numberofpartitions,'weighting','on');
+md.qmu.partition=md.qmu.partition-1;
 
 %variables
-md.variables.rho_ice=normal_uncertain('RhoIce',md.materials.rho_ice,0.01);
-md.variables.drag_coefficient=normal_uncertain('scaled_FrictionCoefficient',1,0.01);
+md.qmu.variables.rho_ice=normal_uncertain('RhoIce',md.materials.rho_ice,0.01);
+md.qmu.variables.drag_coefficient=normal_uncertain('scaled_FrictionCoefficient',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.qmu.responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
 
 %method
-md.qmu_method     =dakota_method('nond_l');
+md.qmu.method     =dakota_method('nond_l');
 
 %parameters
-md.qmu_params.direct=true;
-md.qmu_params.analysis_driver='diagnostic';
-md.qmu_params.evaluation_concurrency=1;
-md.qmu_params.interval_type='forward';
+md.qmu.params.direct=true;
+md.qmu.params.analysis_driver='diagnostic';
+md.qmu.params.evaluation_concurrency=1;
+md.qmu.params.interval_type='forward';
 
 
 %imperative!
 md.eps_rel=10^-5; %tighten for qmu analyses
-md.dakota_analysis=1;
+md.qmu.isdakota=1;
 
 %solve
Index: /issm/trunk/test/NightlyRun/test427.m
===================================================================
--- /issm/trunk/test/NightlyRun/test427.m	(revision 9651)
+++ /issm/trunk/test/NightlyRun/test427.m	(revision 9652)
@@ -15,33 +15,33 @@
 
 %partitioning
-md.npart=20;
-md=partitioner(md,'package','chaco','npart',md.npart,'weighting','on');
-md.part=md.part-1;
+md.qmu.numberofpartitions=20;
+md.qmu.partitionitioner(md,'package','chaco','npart',md.qmu.numberofpartitions,'weighting','on');
+md.qmu.partition=md.qmu.partition-1;
 
 %variables
-md.variables.drag_coefficient=normal_uncertain('scaled_FrictionCoefficient',1,0.01);
+md.qmu.variables.drag_coefficient=normal_uncertain('scaled_FrictionCoefficient',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('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]);
+md.qmu.responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
+md.qmu.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.qmu.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.qmu.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.qmu.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.qmu.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.qmu.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;
+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;
 
 %method
-md.qmu_method     =dakota_method('nond_l');
+md.qmu.method     =dakota_method('nond_l');
 
 %parameters
-md.qmu_params.direct=true;
-md.qmu_params.analysis_driver='diagnostic';
-md.qmu_params.evaluation_concurrency=1;
-md.qmu_params.interval_type='forward';
-md.dakota_analysis=1;
+md.qmu.params.direct=true;
+md.qmu.params.analysis_driver='diagnostic';
+md.qmu.params.evaluation_concurrency=1;
+md.qmu.params.interval_type='forward';
+md.qmu.isdakota=1;
 md.eps_rel=10^-5; %tighten for qmu analyses
 
Index: /issm/trunk/test/NightlyRun/test428.m
===================================================================
--- /issm/trunk/test/NightlyRun/test428.m	(revision 9651)
+++ /issm/trunk/test/NightlyRun/test428.m	(revision 9652)
@@ -15,33 +15,33 @@
 
 %partitioning
-md.npart=20;
-md=partitioner(md,'package','chaco','npart',md.npart,'weighting','on');
-md.part=md.part-1;
+md.qmu.numberofpartitions=20;
+md.qmu.partitionitioner(md,'package','chaco','npart',md.qmu.numberofpartitions,'weighting','on');
+md.qmu.partition=md.qmu.partition-1;
 
 %variables
-md.variables.drag_coefficient=normal_uncertain('scaled_FrictionCoefficient',1,0.01);
+md.qmu.variables.drag_coefficient=normal_uncertain('scaled_FrictionCoefficient',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('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]);
+md.qmu.responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
+md.qmu.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.qmu.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.qmu.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.qmu.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.qmu.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.qmu.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;
+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;
 
 %method
-md.qmu_method     =dakota_method('nond_l');
+md.qmu.method     =dakota_method('nond_l');
 
 %parameters
-md.qmu_params.direct=true;
-md.qmu_params.analysis_driver='diagnostic';
-md.qmu_params.evaluation_concurrency=1;
-md.qmu_params.interval_type='forward';
-md.dakota_analysis=1;
+md.qmu.params.direct=true;
+md.qmu.params.analysis_driver='diagnostic';
+md.qmu.params.evaluation_concurrency=1;
+md.qmu.params.interval_type='forward';
+md.qmu.isdakota=1;
 md.eps_rel=10^-5; %tighten for qmu analyses
 
Index: /issm/trunk/test/NightlyRun/test434.m
===================================================================
--- /issm/trunk/test/NightlyRun/test434.m	(revision 9651)
+++ /issm/trunk/test/NightlyRun/test434.m	(revision 9652)
@@ -13,23 +13,23 @@
 
 %Dakota options
-md.variables.drag_coefficient=normal_uncertain('scaled_FrictionCoefficient',1,0.01);
+md.qmu.variables.drag_coefficient=normal_uncertain('scaled_FrictionCoefficient',1,0.01);
 
-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]);
+md.qmu.responses.MaxVel=response_function('MaxVel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
+md.qmu.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.qmu.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.qmu.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.qmu.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.qmu.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.qmu.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;
+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;
 
 
 %%  nond_sampling study
 
-md.qmu_method     =dakota_method('nond_samp');
-md.qmu_method(end)=dmeth_params_set(md.qmu_method(end),...
+md.qmu.method     =dakota_method('nond_samp');
+md.qmu.method(end)=dmeth_params_set(md.qmu.method(end),...
 'seed',1234,...
 'samples',20,...
@@ -37,13 +37,13 @@
 
 %%  a variety of parameters
-md.qmu_params.evaluation_concurrency=1;
-md.qmu_params.analysis_driver='';
-md.qmu_params.analysis_components='';
+md.qmu.params.evaluation_concurrency=1;
+md.qmu.params.analysis_driver='';
+md.qmu.params.analysis_components='';
 
 %partitioning
-md.npart=20;
-md=partitioner(md,'package','chaco','npart',md.npart,'weighting','on');
-md.part=md.part-1;
-md.dakota_analysis=1;
+md.qmu.numberofpartitions=20;
+md.qmu.partitionitioner(md,'package','chaco','npart',md.qmu.numberofpartitions,'weighting','on');
+md.qmu.partition=md.qmu.partition-1;
+md.qmu.isdakota=1;
 
 md.eps_rel=10^-5; %tighten for qmu analyses
Index: /issm/trunk/test/NightlyRun/test435.m
===================================================================
--- /issm/trunk/test/NightlyRun/test435.m	(revision 9651)
+++ /issm/trunk/test/NightlyRun/test435.m	(revision 9652)
@@ -7,11 +7,11 @@
 
 %partitioning
-md.npart=100;
-md=partitioner(md,'package','chaco','npart',md.npart);
-md.part=md.part-1;
+md.qmu.numberofpartitions=100;
+md.qmu.partitionitioner(md,'package','chaco','npart',md.qmu.numberofpartitions);
+md.qmu.partition=md.qmu.partition-1;
 
 vector=(1:1:md.numberofnodes)';
 vector_on_partition=AreaAverageOntoPartition(md,vector);
-vector_on_nodes=vector_on_partition(md.part+1);
+vector_on_nodes=vector_on_partition(md.qmu.partition+1);
 
 field_names     ={'vector_on_nodes'};
Index: /issm/trunk/test/NightlyRun/test439.m
===================================================================
--- /issm/trunk/test/NightlyRun/test439.m	(revision 9651)
+++ /issm/trunk/test/NightlyRun/test439.m	(revision 9652)
@@ -6,24 +6,24 @@
 
 %partitioning
-md.npart=10;
-md=partitioner(md,'package','chaco','npart',md.npart);
-md.part=md.part-1;
-md.dakota_analysis=1;
+md.qmu.numberofpartitions=10;
+md.qmu.partitionitioner(md,'package','chaco','npart',md.qmu.numberofpartitions);
+md.qmu.partition=md.qmu.partition-1;
+md.qmu.isdakota=1;
 
 %Dakota options
 %variables
-md.variables.rho_ice=normal_uncertain('RhoIce',md.materials.rho_ice,0.01);
+md.qmu.variables.rho_ice=normal_uncertain('RhoIce',md.materials.rho_ice,0.01);
 
 %responses
-md.responses.MaxVel=response_function('scaled_Thickness',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
+md.qmu.responses.MaxVel=response_function('scaled_Thickness',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
 
 %method
-md.qmu_method     =dakota_method('nond_l');
+md.qmu.method     =dakota_method('nond_l');
 
 %parameters
-md.qmu_params.direct=true;
-md.qmu_params.analysis_driver='diagnostic';
-md.qmu_params.evaluation_concurrency=1;
-md.qmu_params.interval_type='forward';
+md.qmu.params.direct=true;
+md.qmu.params.analysis_driver='diagnostic';
+md.qmu.params.evaluation_concurrency=1;
+md.qmu.params.interval_type='forward';
 
 
@@ -36,11 +36,11 @@
 
 %test on thickness
-h=zeros(md.npart,1);
-for i=1:md.npart,
+h=zeros(md.qmu.numberofpartitions,1);
+for i=1:md.qmu.numberofpartitions,
 	h(i)=md.dakotaresults.dresp_out(i).mean;
 end
 
 %project onto grid
-thickness=h(md.part+1);
+thickness=h(md.qmu.partition+1);
 
 %Fields and tolerances to track changes
Index: /issm/trunk/test/NightlyRun/test440.m
===================================================================
--- /issm/trunk/test/NightlyRun/test440.m	(revision 9651)
+++ /issm/trunk/test/NightlyRun/test440.m	(revision 9652)
@@ -6,24 +6,24 @@
 
 %partitioning
-md.npart=10;
-md=partitioner(md,'package','chaco','npart',md.npart);
-md.part=md.part-1;
-md.dakota_analysis=1;
+md.qmu.numberofpartitions=10;
+md.qmu.partitionitioner(md,'package','chaco','npart',md.qmu.numberofpartitions);
+md.qmu.partition=md.qmu.partition-1;
+md.qmu.isdakota=1;
 
 %Dakota options
 %variables
-md.variables.rho_ice=normal_uncertain('RhoIce',md.materials.rho_ice,0.01);
+md.qmu.variables.rho_ice=normal_uncertain('RhoIce',md.materials.rho_ice,0.01);
 
 %responses
-md.responses.MaxVel=response_function('scaled_Thickness',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
+md.qmu.responses.MaxVel=response_function('scaled_Thickness',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
 
 %method
-md.qmu_method     =dakota_method('nond_l');
+md.qmu.method     =dakota_method('nond_l');
 
 %parameters
-md.qmu_params.direct=true;
-md.qmu_params.analysis_driver='diagnostic';
-md.qmu_params.evaluation_concurrency=1;
-md.qmu_params.interval_type='forward';
+md.qmu.params.direct=true;
+md.qmu.params.analysis_driver='diagnostic';
+md.qmu.params.evaluation_concurrency=1;
+md.qmu.params.interval_type='forward';
 
 
@@ -36,11 +36,11 @@
 
 %test on thickness
-h=zeros(md.npart,1);
-for i=1:md.npart,
+h=zeros(md.qmu.numberofpartitions,1);
+for i=1:md.qmu.numberofpartitions,
 	h(i)=md.dakotaresults.dresp_out(i).mean;
 end
 
 %project onto grid
-thickness=h(md.part+1);
+thickness=h(md.qmu.partition+1);
 
 %Fields and tolerances to track changes