Changeset 18267

Timestamp:

07/18/14 14:38:56 (11 years ago)

Author:

schlegel

Message:

NEW: add icebridge tutorial

Location:

issm/trunk-jpl/examples

Files:

• IceBridge (added)
• IceBridge/DomainOutline.exp (added)
• IceBridge/Greenland.par (added)
• IceBridge/Greenland_cheatsheet.par (added)
• IceBridge/Greenland_noOIB.par (added)
• IceBridge/Jak_grounded.exp (added)
• IceBridge/Jak_outline.exp (added)
• IceBridge/Models (added)
• IceBridge/data_gaps.exp (added)
• IceBridge/runme.m (added)
• Pig/runme.m (modified) (2 diffs)
• UncertaintyQuantification/runme.m (modified) (2 diffs)

issm/trunk-jpl/examples/Pig/runme.m

@@ -188 +188 @@
 
         % Load Model
-
-        % Disable inversion
+        md = loadmodel('./Models/PIG.Control_drag');
+        md.inversion.iscontrol=0;
+
+        disp('   Extruding mesh')
+        number_of_layers=3;
+        md=extrude(md,number_of_layers,0.9);
 
         % Extrude Mesh
 
-        % Set Flowequation
-
         % Solve
+        md=solve(md,StressbalanceSolutionEnum);
 
         % Save Model
+        save ./Models/PIG.ModelHO md;
 
 end
@@ -215 +219 @@
                                                 'ylim#all',[min(mdHO.mesh.y) max(mdHO.mesh.y)]/10^3,...
                                                 'FontSize#all',12,...
-                                                'data',mdHO.initialization.vel,'title','Observed velocity',...
+                                                'data',mdSSA.initialization.vel,'title','Observed velocity',...
                                                 'data',(mdHO.results.StressbalanceSolution.Vel(surf)-mdHO.initialization.vel(surf)),'title','(HO-observed) velocities',...
                                                 'data',mdSSA.results.StressbalanceSolution.Vel,'title','Modeled SSA Velocity',...
                                                 'data',(mdHO.results.StressbalanceSolution.Vel(surf)-mdSSA.results.StressbalanceSolution.Vel),'title','(HO-SSA) velocities',...
-                                                'data',mdHO.results.StressbalanceSolution.Vel,'title','Modeled HO surface Velocities',...
+                                                'data',project2d(mdHO,mdHO.results.StressbalanceSolution.Vel,1),'title','Modeled HO surface Velocities',...
                                                 'data',(mdHO.results.StressbalanceSolution.Vel(surf)-mdHO.results.StressbalanceSolution.Vel(basal)),'title','(HOsurf-HO base) velocities',...
                                                 'caxis#1',([1.5,4000]),'caxis#3',([1.5,4000]),'caxis#5',([1.5,4000]),...
                                                 'colorbar#all','on','view#all',2,...
                                                 'colorbartitle#all','[m/yr]',...
-                                                'layer#5',1, 'log#1', 10,'log#3', 10,'log#5', 10);
-end
-% }}}
+                                                'log#1', 10,'log#3', 10,'log#5', 10,'gridded#all',1);
+end
+% }}}

issm/trunk-jpl/examples/UncertaintyQuantification/runme.m

@@ -1 +1 @@
 %PIG Uncertainty Quantification Application
-steps=[1:3]; 
+steps=[7]; 
 
 if any(steps==1) 
@@ -162 +162 @@
 
         %variables
+        md.qmu.variables.Thickness=normal_uncertain('scaled_Thickness',1,0.05);
         md.qmu.variables.DragCoefficient=normal_uncertain(...
         'scaled_FrictionCoefficient',1,0.05);
         md.qmu.variables.rheology_B=normal_uncertain(...
         'scaled_MaterialsRheologyB',1,0.05);
-        md.qmu.variables.Thickness=normal_uncertain('scaled_Thickness',1,0.05);
 
         %responses