Changeset 14123

Timestamp:

12/10/12 17:33:17 (12 years ago)

Author:

helsen

Message:

workshop update

File:

• issm/workshop/2012/Talks/13_IceBridge/runme.m (modified) (5 diffs)

issm/workshop/2012/Talks/13_IceBridge/runme.m ¶

@@ -1 +1 @@
 clear all;
-steps=[1:5];
+steps=[1];
 
 %Location of SeaRISE dataset
@@ -21 +21 @@
 
         %Mesh greenland without refinement in Jak basin
-        %md=bamg(md,'hmax',400000,'hmin',5000,'gradation',1.7,'field',vel,'err',8);
-        %plotmodel (md,'data','mesh');
-        %return
+        md=bamg(md,'hmax',400000,'hmin',5000,'gradation',1.7,'field',vel,'err',8);
+        return
         
         %Refine mesh in the region of Jakobshavn (resolution = 1000 m)
@@ -30 +29 @@
                 './Exp_Par/Jak_outline.exp','node',1);
         hmaxVertices(find(in))=1000;
-        
-        %Adapt the mesh to minimize error in velocity interpolation
         md=bamg(md,'hmax',400000,'hmin',5000,'gradation',1.7,'field',vel,...
                 'err',8,'hmaxVertices',hmaxVertices);
-        plotmodel (md,'data','mesh');
 
         %convert x,y coordinates (Polar stereo) to lat/lon
@@ -47 +43 @@
 
         md = setmask(md,'','');
-        md = parameterize(md,'./Exp_Par/Greenland_solution.par');
+        md = parameterize(md,'./Exp_Par/Greenland.par');
         md = setflowequation(md,'macayeal','all');
 
-        save ./Models/Greenland.Parameterization md; 
+        save ./Models/Greenland.Parameterization2 md; 
 end
 
@@ -134 +130 @@
         plotmodel(md,'data',md.results.TransientSolution(end).Vel,'caxis',[1e-1 6000],...
                 'log', 10, 'title', 'Velocity (m/y)','gridded',1, ...
-                'data', md.results.TransientSolution(1).SurfaceforcingsMassBalance, ...
+                'data', md.results.TransientSolution(end).SurfaceforcingsMassBalance, ...
                 'title', 'Surface mass balance (m/y)','gridded',1, ...
                 'data',md.results.TransientSolution(end).Thickness,...