Index: /issm/branches/trunk-larour-SLPS2020/test/NightlyRun/test2005.m
===================================================================
--- /issm/branches/trunk-larour-SLPS2020/test/NightlyRun/test2005.m	(revision 25594)
+++ /issm/branches/trunk-larour-SLPS2020/test/NightlyRun/test2005.m	(revision 25595)
@@ -9,5 +9,5 @@
 md.solidearth.surfaceload.icethicknesschange=zeros(md.mesh.numberofelements,1);
 md.solidearth.sealevel=zeros(md.mesh.numberofvertices,1);
-md.dsl.global_average_thermosteric_sea_level_change=[0;0];
+md.dsl.global_average_thermosteric_sea_level_change=[1;0];
 md.dsl.sea_surface_height_change_above_geoid=zeros(md.mesh.numberofvertices+1,1);
 md.dsl.sea_water_pressure_change_at_sea_floor=zeros(md.mesh.numberofvertices+1,1);
@@ -54,5 +54,5 @@
 
 %Miscellaneous
-md.miscellaneous.name='test2002';
+md.miscellaneous.name='test2005';
 
 %Solution parameters
@@ -99,5 +99,5 @@
 S5=md.results.TransientSolution(5).Sealevel;
 S10=md.results.TransientSolution(10).Sealevel;
-Seus10=md.results.TransientSolution(10).SealevelRSLEustatic;
+Seus10=md.results.TransientSolution(10).Bslr;
 
 %Fields and tolerances to track changes
Index: /issm/branches/trunk-larour-SLPS2020/test/NightlyRun/test2006.m
===================================================================
--- /issm/branches/trunk-larour-SLPS2020/test/NightlyRun/test2006.m	(revision 25595)
+++ /issm/branches/trunk-larour-SLPS2020/test/NightlyRun/test2006.m	(revision 25595)
@@ -0,0 +1,156 @@
+%Test Name: EarthSlr Dakota Sampling glaciers.
+
+%mesh earth:
+md=model;
+md.cluster=generic('name',oshostname(),'np',5);
+md.mesh=gmshplanet('radius',6.371012*10^3,'resolution',700.); %700 km resolution mesh
+
+%parameterize solidearth solution:
+%solidearth loading:  {{{
+md.solidearth.surfaceload.icethicknesschange=zeros(md.mesh.numberofelements,1);
+md.solidearth.sealevel=zeros(md.mesh.numberofvertices,1);
+md.dsl.global_average_thermosteric_sea_level_change=[1;0];
+md.dsl.sea_surface_height_change_above_geoid=zeros(md.mesh.numberofvertices+1,1);
+md.dsl.sea_water_pressure_change_at_sea_floor=zeros(md.mesh.numberofvertices+1,1);
+%antarctica
+late=sum(md.mesh.lat(md.mesh.elements),2)/3;
+longe=sum(md.mesh.long(md.mesh.elements),2)/3;
+pos=find(late < -80);
+md.solidearth.surfaceload.icethicknesschange(pos)=-100;
+%greenland
+pos=find(late>70 & late<80 & longe>-60 & longe<-30);
+md.solidearth.surfaceload.icethicknesschange(pos)=-100;
+
+%elastic loading from love numbers:
+md.solidearth.lovenumbers=lovenumbers('maxdeg',100);
+
+%}}}
+%mask:  {{{
+mask=gmtmask(md.mesh.lat,md.mesh.long);
+icemask=ones(md.mesh.numberofvertices,1);
+pos=find(mask==0);
+icemask(pos)=-1;
+pos=find(sum(mask(md.mesh.elements),2)<3);
+icemask(md.mesh.elements(pos,:))=-1;
+md.mask.ice_levelset=icemask;
+md.mask.ocean_levelset=-icemask;
+
+%make sure that the elements that have loads are fully grounded:
+pos=find(md.solidearth.surfaceload.icethicknesschange);
+md.mask.ocean_levelset(md.mesh.elements(pos,:))=1;
+
+%make sure wherever there is an ice load, that the mask is set to ice:
+pos=find(md.solidearth.surfaceload.icethicknesschange);
+md.mask.ice_levelset(md.mesh.elements(pos,:))=-1;
+% }}}
+
+md.solidearth.settings.ocean_area_scaling=0;
+
+%Geometry for the bed, arbitrary: 
+md.geometry.bed=-ones(md.mesh.numberofvertices,1);
+
+%Materials: 
+md.materials=materials('hydro');
+
+
+%Miscellaneous
+md.miscellaneous.name='test2005';
+
+%Solution parameters
+md.solidearth.settings.reltol=NaN;
+md.solidearth.settings.abstol=1e-3;
+md.solidearth.settings.computesealevelchange=1;
+
+% max number of iteration reverted back to 10 (i.e., the original default value)
+md.solidearth.settings.maxiter=10;
+
+%eustatic + rigid + elastic + rotation run:
+md.solidearth.settings.rigid=1;
+md.solidearth.settings.elastic=1;
+md.solidearth.settings.rotation=1;
+
+%transient settings: 
+md.timestepping.start_time=0;
+md.timestepping.final_time=10;
+md.timestepping.time_step=1;
+md.transient.isslr=1;
+md.transient.issmb=0;
+md.transient.isgia=1;
+md.transient.ismasstransport=0;
+md.transient.isstressbalance=0;
+md.transient.isthermal=0;
+dh=md.solidearth.surfaceload.icethicknesschange;
+deltathickness=zeros(md.mesh.numberofelements+1,10);
+for i=1:10,
+	deltathickness(1:end-1,i)=dh*i;
+end
+deltathickness(end,:)=0:1:9;
+md.solidearth.surfaceload.icethicknesschange=deltathickness;
+
+%hack: 
+md.geometry.surface=zeros(md.mesh.numberofvertices,1);
+md.geometry.thickness=ones(md.mesh.numberofvertices,1);
+md.geometry.base=-ones(md.mesh.numberofvertices,1);
+md.geometry.bed=md.geometry.base;
+
+
+%Uncertainty Quantification
+%ice sheets {{{
+npart=1; nt=1;
+partition=-ones(md.mesh.numberofelements,1);
+pos=find(late < -80); partition(pos)=0;
+pos=find(late>70 & late<80 & longe>-60 & longe<-30); partition(pos)=0;
+
+%variables: 
+qmuvar.surfaceload=normal_uncertain('descriptor','scaled_SurfaceloadIceThicknessChange',...
+	'mean',1*ones(npart,nt),...
+	'stddev',.1*ones(npart,nt),... %10% standard deviation
+	'partition',partition,...
+	'transient','on',...
+	'nsteps',nt);
+%}}}
+
+%correlation:
+md.qmu.correlation_matrix=[];
+
+%variables final declaration:
+md.qmu.variables=struct();
+md.qmu.variables.surfaceload=qmuvar.surfaceload;
+
+%responses  % {{{
+md.qmu.responses.sealevel1=response_function('descriptor','Outputdefinition1');
+%output definitions: 
+md.outputdefinition.definitions={nodalvalue('name','SNode','definitionstring','Outputdefinition1', ...
+	'model_string','Sealevel','node',1)};  %Pick vertex index at 1.
+% }}}
+	%algorithm:  % {{{
+	md.qmu.method     =dakota_method('nond_samp');
+	md.qmu.method(end)=dmeth_params_set(md.qmu.method(end),...
+	'seed',1234,...
+	'samples',10,...
+	'sample_type','random');
+	md.qmu.output=1; 
+	%}}}
+	%parameters % {{{
+	md.qmu.params.direct=true;
+	md.qmu.params.interval_type='forward';
+	md.qmu.params.analysis_driver='matlab';
+	md.qmu.params.evaluation_scheduling='master';
+	md.qmu.params.processors_per_evaluation=2;
+	md.qmu.params.tabular_graphics_data=true; 
+	md.qmu.isdakota=1;
+	md.verbose=verbose(0); md.verbose.qmu=1;
+	% }}}
+
+%run transient dakota solution: 
+md=solve(md,'Transient');
+
+%sealevel mean and standard deviation:
+Smean=md.results.dakota.dresp_out.mean;
+Sstd=md.results.dakota.dresp_out.stddev;
+
+%Fields and tolerances to track changes
+field_names={'Sealevel Mean','Sealevel Standard Deviation'};
+field_tolerances={1e-13,1e-13};
+field_values={Smean,Sstd};
+