Changeset 25938

Timestamp:

01/15/21 11:52:26 (4 years ago)

Author:

adhikari

Message:

BUG: transient runs fixed

File:

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

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

@@ -2 +2 @@
 addpath('../Data','../Functions');
 
-steps=[6];
+steps=[7];
 
 if any(steps==1) %{{{
@@ -199 +199 @@
         md.transient.isstressbalance=0;
         md.transient.isthermal=0;
+        md.transient.isgia=1; 
         md.transient.isslr=1;
 
@@ -210 +211 @@
         md.verbose=verbose('111111111');
 
+        md.solidearth.requested_outputs = {'Sealevel','SealevelRSL'};
+        
         md=solve(md,'tr');
 
@@ -219 +222 @@
         md = loadmodel('./Models/SlrGRACE_Transient');
 
-        time = md.slr.deltathickness(end,:);
+        time = md.solidearth.surfaceload.icethicknesschange(end,:);
 
         for tt=1:length(time)
-                gmsl(tt) = md.results.TransientSolution(tt).SealevelEustatic*1000;      % GMSL rate mm/yr
-                sol1(:,tt) = md.slr.deltathickness(1:end-1,tt)*100;                                     % ice equivalent height [cm/yr]
-                sol2(:,tt) = md.results.TransientSolution(tt+1).Sealevel*1000;          % mm/yr
+                gmsl(tt) = md.results.TransientSolution(tt).Bslr*1000; % GMSL rate mm/yr
+                sol1(:,tt) = md.solidearth.surfaceload.icethicknesschange(1:end-1,tt)*100;      % ice equivalent height [cm/yr]
+                sol2(:,tt) = md.results.TransientSolution(tt+1).SealevelRSL*1000;       % mm/yr
         end
         sol_name = {'Change in water equivalent height [cm]', 'Relative sea-level [mm]'};