Index: /issm/trunk/test/Verification/test28_IceSheetIceFrontM2dDakotaPartition/DomainOutline.exp
===================================================================
--- /issm/trunk/test/Verification/test28_IceSheetIceFrontM2dDakotaPartition/DomainOutline.exp	(revision 4612)
+++ /issm/trunk/test/Verification/test28_IceSheetIceFrontM2dDakotaPartition/DomainOutline.exp	(revision 4612)
@@ -0,0 +1,10 @@
+## Name:domainoutline
+## Icon:0
+# Points Count  Value
+5 1.
+# X pos Y pos
+0 0
+1000000 0
+1000000 1000000
+0 1000000
+0 0
Index: /issm/trunk/test/Verification/test28_IceSheetIceFrontM2dDakotaPartition/Front.exp
===================================================================
--- /issm/trunk/test/Verification/test28_IceSheetIceFrontM2dDakotaPartition/Front.exp	(revision 4612)
+++ /issm/trunk/test/Verification/test28_IceSheetIceFrontM2dDakotaPartition/Front.exp	(revision 4612)
@@ -0,0 +1,10 @@
+## Name:icefront
+## Icon:0
+# Points Count  Value
+5 1.
+# X pos Y pos
+-1000 900000
+-1000 1100000
+1100000 1100000
+1100000 900000
+-1000 900000
Index: /issm/trunk/test/Verification/test28_IceSheetIceFrontM2dDakotaPartition/Lrel28.m
===================================================================
--- /issm/trunk/test/Verification/test28_IceSheetIceFrontM2dDakotaPartition/Lrel28.m	(revision 4612)
+++ /issm/trunk/test/Verification/test28_IceSheetIceFrontM2dDakotaPartition/Lrel28.m	(revision 4612)
@@ -0,0 +1,23 @@
+%Dakota options
+md.variables.nuv=normal_uncertain.empty();
+md.variables.nuv(end+1)=normal_uncertain('RhoIce',md.rho_ice,0.01);
+md.variables.nuv(end+1)=normal_uncertain('DragCoefficient',1,0.01);
+
+md.responses=struct();
+md.responses.rf =response_function.empty();
+md.responses.rf (end+1)=response_function('max_vel',[],[0.0001 0.001 0.01 0.25 0.5 0.75 0.99 0.999 0.9999]);
+
+md.qmu_method     =dakota_method('nond_l');
+
+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.tabular_graphics_data=true;
+
+%partitioning
+md.npart=20;
+md.part=partitioner(md,'package','chaco','npart',md.npart,'weighting','on')-1;
+md.qmu_analysis=1;
+
+md.eps_rel=10^-10; %tighten for qmu analysese
Index: /issm/trunk/test/Verification/test28_IceSheetIceFrontM2dDakotaPartition/Square.par
===================================================================
--- /issm/trunk/test/Verification/test28_IceSheetIceFrontM2dDakotaPartition/Square.par	(revision 4612)
+++ /issm/trunk/test/Verification/test28_IceSheetIceFrontM2dDakotaPartition/Square.par	(revision 4612)
@@ -0,0 +1,42 @@
+%Start defining model parameters here
+
+%dynamics
+md.verbose=0;
+md.dt=1; %1 year
+md.ndt=md.dt*10; 
+md.artificial_diffusivity=0;
+md.mem_debug=0;
+
+hmin=300;
+hmax=1000;
+ymin=min(md.y);
+ymax=max(md.y);
+md.thickness=hmax+(hmin-hmax)*(md.y-ymin)/(ymax-ymin);
+md.bed=-md.rho_ice/md.rho_water*md.thickness+100;
+md.surface=md.bed+md.thickness;
+
+md.drag_type=2; %0 none 1 plastic 2 viscous
+md.drag_coefficient=20*ones(md.numberofgrids,1); %q=1.
+%Take care of iceshelves: no basal drag
+pos=find(md.elementoniceshelf);
+md.drag_coefficient(md.elements(pos,:))=0;
+md.drag_p=ones(md.numberofelements,1);
+md.drag_q=ones(md.numberofelements,1);
+md.viscosity_overshoot=0.3;
+
+md.observed_temperature=(273-20)*ones(md.numberofgrids,1);
+
+md.rheology_B=paterson(md.observed_temperature);
+md.rheology_n=3*ones(md.numberofelements,1);
+
+%Deal with boundary conditions:
+md=SetMarineIceSheetBC(md,'Front.exp');
+
+%Parallel options
+md.np=8;
+md.time=50;
+md.waitonlock=30;
+
+
+%Dakota:
+Lrel28
Index: /issm/trunk/test/Verification/test28_IceSheetIceFrontM2dDakotaPartition/configuration.m
===================================================================
--- /issm/trunk/test/Verification/test28_IceSheetIceFrontM2dDakotaPartition/configuration.m	(revision 4612)
+++ /issm/trunk/test/Verification/test28_IceSheetIceFrontM2dDakotaPartition/configuration.m	(revision 4612)
@@ -0,0 +1,26 @@
+%name
+testname='test28_IceSheetIceFrontM2dDakotaPartition';
+
+%meshing
+resolution=100000;
+domainname='DomainOutline.exp';
+riftname='';
+rifts=0;
+
+%geography
+iceshelfname='';
+icesheetname='';
+
+%parameterization
+parametername='Square.par';
+			
+%elements type
+elementstype='''macayeal'',''all''';
+
+%extrusion
+numlayers=0;
+extrusionexponent=NaN;
+
+%sequences     analysis    sub_analysis  qmu   control  control_fit parallel
+sequences={  {DiagnosticSolutionEnum,  NoneAnalysisEnum,      1 ,     0,      'nan'        1    }
+	};
Index: /issm/trunk/test/Verification/test28_IceSheetIceFrontM2dDakotaPartition/testpostsolve.m
===================================================================
--- /issm/trunk/test/Verification/test28_IceSheetIceFrontM2dDakotaPartition/testpostsolve.m	(revision 4612)
+++ /issm/trunk/test/Verification/test28_IceSheetIceFrontM2dDakotaPartition/testpostsolve.m	(revision 4612)
@@ -0,0 +1,2 @@
+md=tres(md,'dakota');
+md.results.dakota.importancefactors=importancefactors(md,'drag','max_vel');
Index: /issm/trunk/test/Verification/test28_IceSheetIceFrontM2dDakotaPartition/testpresolve.m
===================================================================
--- /issm/trunk/test/Verification/test28_IceSheetIceFrontM2dDakotaPartition/testpresolve.m	(revision 4612)
+++ /issm/trunk/test/Verification/test28_IceSheetIceFrontM2dDakotaPartition/testpresolve.m	(revision 4612)
@@ -0,0 +1,9 @@
+%initial velocities
+if md.numberofgrids==340
+	load Velocities; md.vx=0.5*vx; md.vy=0.5*vy;
+end
+
+%rm existing results
+if exist('./qmu')
+	system('rm -rf qmu/');
+end