Changeset 18264

Timestamp:

07/18/14 12:33:21 (11 years ago)

Author:

schlegel

Message:

CHG: add plot execise

Location:

issm/trunk-jpl/examples/Greenland

Files:

• CheatSheet.m (modified) (4 diffs)
• runme.m (modified) (4 diffs)

issm/trunk-jpl/examples/Greenland/CheatSheet.m

@@ -9 +9 @@
 
         %Save surface mass balance, by looping through 200 years, or 1000 steps
-        surfmb=[]; for i=1:100; surfmb=[surfmb ...
+        surfmb=[]; for i=1:1000; surfmb=[surfmb ...
                 md.results.TransientSolution(i).SurfaceforcingsMassBalance]; end
 
@@ -19 +19 @@
 
         %Save velocity by looping through 1000 steps
-        vel=[]; for i=1:100; vel=[vel md.results.TransientSolution(i).Vel]; end
+        vel=[]; for i=1:1000; vel=[vel md.results.TransientSolution(i).Vel]; end
 
         %Plot velocity time series in second subplot
@@ -28 +28 @@
 
         %Save Ice Volume by looping through 1000 steps
-        volume=[]; for i=1:100; volume=[volume md.results.TransientSolution(i).IceVolume]; end
+        volume=[]; for i=1:1000; volume=[volume md.results.TransientSolution(i).IceVolume]; end
 
         %Plot volume time series in third subplot
@@ -35 +35 @@
         %Title this plot Mean Velocity and add an x label of years
         title('Ice Volume'); xlabel('years');
-
 end

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

@@ -146 +146 @@
 if any(steps==6)
         disp('   Step 6: Extract Box SMB');
+        md = loadmodel('./Models/Greenland.Transient');
 
         ncbox='../Data/Box_Greenland_SMB_monthly_1840-2012_5km_ver20140421.nc';
@@ -164 +165 @@
         %Area of grid for 5km box
         area_of_grid=5000*5000;
-        totalsmb=reshape(sum(sum(smbbox,1),2),length(t),1);
+        totalsmb=reshape(sum(sum(smbbox/1000,1),2),length(t),1)*area_of_grid;
 
         %save surface mass balance mat dataset
@@ -171 +172 @@
 
         %plot a time series of total SMB
-        plot(t,totalsmb); title('Total Surface mass balance, m^3');
-
-        clear smbbox smbmean totalsmb x1 y1 t
+        plot(t,totalsmb/1e9); title('Total Surface mass balance, Gt'); xlabel('year'); ylabel('Gt/yr');
+
+        clear smbbox
 
 end
@@ -211 +212 @@
 
         %load historic transient model
-        md = loadmodel('./Models/Greenland.HistoricTransient');
 
         %Create Line Plots of relaxation run.  Create a figure.
-        figure
+
 
         %Save surface mass balance, by looping through 200 years, or 1000 steps
-        surfmb=[]; for i=1:100; surfmb=[surfmb ...
-                md.results.TransientSolution(i).SurfaceforcingsMassBalance]; end
+
 
         %Plot surface mass balance time series in first subplot
-        subplot(3,1,1); plot([0.2:0.2:200],mean(surfmb)); 
+
         
         %Title this plot Mean surface mass balance
-        title('Mean Surface mass balance');
+
 
         %Save velocity by looping through 1000 steps
-        vel=[]; for i=1:100; vel=[vel md.results.TransientSolution(i).Vel]; end
+
 
         %Plot velocity time series in second subplot
-        subplot(3,1,2); plot([0.2:0.2:200],mean(vel)); 
+
         
         %Title this plot Mean Velocity
-        title('Mean Velocity');
+
 
         %Save Ice Volume by looping through 1000 steps
-        volume=[]; for i=1:100; volume=[volume md.results.TransientSolution(i).IceVolume]; end
+
 
         %Plot volume time series in third subplot
-        subplot(3,1,3); plot([0.2:0.2:200],volume); 
+
         
         %Title this plot Mean Velocity and add an x label of years
-        title('Ice Volume'); xlabel('years');
+
 end