Changeset 20754

Timestamp:

06/16/16 12:36:19 (9 years ago)

Author:

seroussi

Message:

CHG: finished Pig runme

Location:

issm/trunk-jpl/examples/Pig

Files:

• Help.m (modified) (2 diffs)
• Pig.par (modified) (1 diff)
• runme.m (modified) (6 diffs)

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

@@ -1 +1 @@
 if step==6
-        md = loadmodel('./Models/PIG.Control_drag');
+        md = loadmodel('./Models/PIG_Control_drag');
         md.inversion.iscontrol=0;
 
@@ -12 +12 @@
         md=solve(md,StressbalanceSolutionEnum);
 
-        save ./Models/PIG.ModelHO md;
+        save ./Models/PIG_ModelHO md;
 end

issm/trunk-jpl/examples/Pig/Pig.par

@@ -44 +44 @@
 md.geometry.surface=md.geometry.thickness+md.geometry.base;
 md.geometry.bed=md.geometry.base;
-md.geometry.bed(pos)=md.geometry.base(pos)-100;
+md.geometry.bed(pos)=md.geometry.base(pos)-200;
 
 %Initialization parameters

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

@@ -1 @@
-step=[5];
+step=[7];
 
 if step==1   %Mesh Generation #1
@@ -53 +53 @@
 
         % Save model
-        save ./Models/PIG.Mesh_generation md;
+        save ./Models/PIG_Mesh_generation md;
 end
 
 if step==2  %Masks #2
 
-        md = loadmodel('./Models/PIG.Mesh_generation'); 
+        md = loadmodel('./Models/PIG_Mesh_generation'); 
 
         % Load SeaRISe dataset for Antarctica  
@@ -82 +82 @@
         
         % Save model
-        save ./Models/PIG.SetMask md;
+        save ./Models/PIG_SetMask md;
 end
 
 if step==3  %Parameterization #3
 
-        md = loadmodel('./Models/PIG.SetMask');
+        md = loadmodel('./Models/PIG_SetMask');
         md = parameterize(md,'./Pig.par');
 
@@ -94 +94 @@
         
         % Save model
-        save ./Models/PIG.Parameterization md;
+        save ./Models/PIG_Parameterization md;
 end
 
 if step==4  %Control Method #4
 
-        md = loadmodel('./Models/PIG.Parameterization');
+        md = loadmodel('./Models/PIG_Parameterization');
 
         % Control general
@@ -137 +137 @@
 
         % Save model
-        save ./Models/PIG.Control_drag md;
+        save ./Models/PIG_Control_drag md;
 end
 
 if step==5 %Plot #5
 
-        md = loadmodel('./Models/PIG.Control_drag');
+        md = loadmodel('./Models/PIG_Control_drag');
 
         plotmodel(md,'nlines',2,'ncols',2,'unit#all','km','axis#all','equal',...
@@ -175 +175 @@
 if step==7  %Plot #7
 
-        mdHO = loadmodel('./Models/PIG.ModelHO');
-        mdSSA = loadmodel('./Models/PIG.Control_drag');
+        mdHO = loadmodel('./Models/PIG_ModelHO');
+        mdSSA = loadmodel('./Models/PIG_Control_drag');
 
         basal=find(mdHO.mesh.vertexonbase);