Changeset 27029

Timestamp:

05/30/22 10:22:28 (3 years ago)

Author:

adhikari

Message:

CHG: some fine tuning of load models

File:

• issm/trunk-jpl/examples/SlrGRACE/runme.m (modified) (8 diffs)

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

@@ -1 @@
-clear all;
+%clear all;
 addpath('../Data','../Functions');
 
-steps=[1:7];
+%steps=[1:7];
 
 if any(steps==1) % {{{ 
         disp('   Step 1: Global mesh creation');
 
-        numrefine=5;
+        numrefine=1;
         resolution=150*1e3;                     % inital resolution [m]
         radius = 6.371012*10^6;         % mean radius of Earth, m
@@ -50 +50 @@
         save ./Models/SlrGRACE_Mesh md;
 
-        plotmodel (md,'data',md.mask.ocean_levelset,'edgecolor','k');
+        plotmodel (md,'data',md.mask.ocean_levelset,'edgecolor','k','view',[45 45]);
+
 end % }}} 
 if any(steps==2) % {{{ 
@@ -70 +71 @@
         md.geometry.surface=md.geometry.bed+md.geometry.thickness;
 
-        md.masstransport.spcthickness = [md.geometry.thickness + ice_load; 0];
+        md.masstransport.spcthickness = repmat(md.geometry.thickness,1,2); 
+        md.masstransport.spcthickness(:,end) = md.masstransport.spcthickness(:,end) + ice_load; 
+        md.masstransport.spcthickness(end+1,:) = [0 1]; % dummy time. 
 
         md.smb.mass_balance=zeros(md.mesh.numberofvertices,1); 
@@ -140 +143 @@
         md = loadmodel('./Models/SlrGRACE_Solution');
 
-        sol1 = (md.masstransport.spcthickness(1:end-1)-md.geometry.thickness)*100; % [cm]
+        sol1 = (md.masstransport.spcthickness(1:end-1,2)-md.masstransport.spcthickness(1:end-1,1))*100; % [cm]
         sol2 = (md.results.TransientSolution.Sealevel-md.results.TransientSolution.Bed)*1000;   % [mm]
 
@@ -208 +211 @@
         % masstransport evalulates diff between the successive times, so we should cumsum. 
         num_time = size(ice_load,2); 
-        md.masstransport.spcthickness = [repmat(md.geometry.thickness,1,size(ice_load,2)) + ice_load];
-        md.masstransport.spcthickness =[ md.masstransport.spcthickness; 0:num_time-1]; 
+        md.masstransport.spcthickness = repmat(md.geometry.thickness,1,num_time+1); 
+        md.masstransport.spcthickness(:,2:end) = md.masstransport.spcthickness(:,2:end) + ice_load;
+        md.masstransport.spcthickness(end+1,:) = 0:num_time; % dummy time. 
 
         %Physics 
@@ -241 +245 @@
         md = loadmodel('./Models/SlrGRACE_Transient');
 
-        time = md.masstransport.spcthickness(end,:);
+        time = md.masstransport.spcthickness(end,2:end);
 
         for tt=1:length(time)
                 gmsl(tt) = md.results.TransientSolution(tt).CumBslc*1000; % [mm]
-                sol1(:,tt) = (md.masstransport.spcthickness(1:end-1,tt)-md.geometry.thickness)*100; % [cm]
+                sol1(:,tt) = (md.masstransport.spcthickness(1:end-1,tt+1)-md.geometry.thickness)*100; % [cm]
                 sol2(:,tt) = (md.results.TransientSolution(tt).Sealevel-md.results.TransientSolution(tt).Bed)*1000;     % [mm]
         end
         sol_name = {'Change in water equivalent height [cm]', 'Relative sea level [mm]'};
-        movie_name = {'Movie_dH.mp4','Movie_slr.mp4'};
+        movie_name = {'Movie_dH','Movie_slr'};
         
         res = 1.0;
@@ -273 +277 @@
                 vidObj = VideoWriter(movie_name{kk});
                 vidObj.FrameRate=2; % frames per second
-                vifObj.FileFormat='mp4';
                 open(vidObj);
 
@@ -311 +314 @@
                 close(vidObj);
         end
-        !open *.mp4;
 
         % plot GMSL time series