Changeset 25217

Timestamp:

07/06/20 19:26:14 (5 years ago)

Author:

Eric.Larour

Message:

CHG: dak/slr related tests.

Location:

issm/trunk-jpl/test/SandBox

Files:

• Makefile (modified) (1 diff)
• bp2002.m (added)
• corr2002.m (added)
• dak2002.m (added)
• geoid2002.m (added)
• gia2002.m (added)
• hier2002.m (added)
• interpBedmap2.m (added)
• mmedak.m (added)
• res2002.m (added)
• testbayes.m (modified) (4 diffs)
• tran2002.m (added)
• val2002.m (added)

issm/trunk-jpl/test/SandBox/Makefile

@@ -1 +1 @@
 clean: 
-        rm -rf *.tar.gz *.errlog *.outlog *.bin *.qmu.in *.qmu.out *qmu.err *.queue *.toolkits
+        rm -rf *.tar.gz *.errlog *.outlog *.bin *.qmu.in *.qmu.out *qmu.err *.queue *.toolkits *.out

issm/trunk-jpl/test/SandBox/testbayes.m

@@ -5 +5 @@
 md=setflowequation(md,'SSA','all');
 md.cluster=generic('name',oshostname(),'np',3);
-md.materials.rho_ice=10^7; %involved in the mass flux, make it easy
-md.geometry.thickness(:)=1; %make it easy
-md.geometry.surface=md.geometry.base+md.geometry.thickness;
 
 %constrain all velocities to 1 m/yr, in the y-direction
-md.stressbalance.spcvx(:)=0;
-md.stressbalance.spcvy(:)=1;
+pos=find(md.mesh.y==0); 
+md.stressbalance.spcvx(pos)=0;
+md.stressbalance.spcvy(pos)=0;
+%md.stressbalance.spcvy(:)=1;
 md.stressbalance.spcvz(:)=0;
 
@@ -24 +23 @@
 
 %input/output:
-md.qmu.variables.drag_coefficient=normal_uncertain('scaled_FrictionCoefficient',1,0.01);
+md.qmu.variables.drag_coefficient=normal_uncertain('scaled_FrictionCoefficient',1,0.10);
 md.qmu.responses.misfit=calibration_function('MaxVel');
 %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]);
@@ -32 +31 @@
 md.qmu.method(end)=dmeth_params_set(md.qmu.method(end),...
 'seed',1234,...
-'samples',20,...
+'samples',2000,...
 'queso',true,...
 'metropolis_hastings',true,...
@@ -53 +52 @@
 %Fields and tolerances to track changes
 md.qmu.results=md.results.dakota;
-
-error;
-
-%ok, mass flux of 3 profiles should be -3 Gt/yr -3 Gt/yr and the sum, which is -6 Gt/yr
-%we recover those mass fluxes through the mean of the response.
-%also, we recover the max velo, which should be 1m/yr. 
-%we put all that data in the montecarlo field, which we will use to test for success.
-%also, check that the stddev are 0.
-md.results.dakota.montecarlo=[];
-for i=1:8,
-        md.results.dakota.montecarlo=[md.results.dakota.montecarlo md.results.dakota.dresp_out(i).mean];
-end
-for i=1:8,
-        md.results.dakota.montecarlo=[md.results.dakota.montecarlo md.results.dakota.dresp_out(i).stddev];
-end
-field_names     ={'montecarlo'};
-field_tolerances={1e-11};
-field_values={...
-         md.results.dakota.montecarlo,...
-        };