Changeset 25290 for issm/trunk-jpl/test/NightlyRun/test2051.m

Timestamp:

07/18/20 23:19:58 (5 years ago)

Author:

adhikari

Message:

CHG: GIA ivins nightly merged, 2051,52,53 => 2051; 2071,72,73=> 2052; 2081,2082,2083 => 2053

File:

• issm/trunk-jpl/test/NightlyRun/test2051.m (modified) (3 diffs)

issm/trunk-jpl/test/NightlyRun/test2051.m

@@ -1 @@
-%Test Name: GiaIvinsBenchmarksAB2dA1
+%Test Name: GiaIvinsBenchmarksAB2dA
 % Benchmark experiments (Figure A2a Ivins and James, 1999, Geophys. J. Int.) 
 md=triangle(model(),'../Exp/RoundFrontEISMINT.exp',200000.);
@@ -8 +8 @@
 md.gia.cross_section_shape=1;    % for square-edged x-section 
 
-%% define loading history 
+% evaluation time (termed start_time) 
 md.timestepping.start_time=2002100; % after 2 kyr of deglaciation 
 md.timestepping.final_time=2500000; % 2,500 kyr
+
+%% define loading history {{{ 
 md.geometry.thickness=[...
         [md.geometry.thickness*0.0; 0.0],...
@@ -18 +20 @@
         [md.geometry.thickness*0.0; md.timestepping.start_time],...
         ];
+% }}} 
+
+% find out elements that have zero loads throughout the loading history.
+pos = find(sum(abs(md.geometry.thickness(1:end-1,:)),2)==0);
+md.mask.ice_levelset(pos)=1; % no-ice.
+
+md.cluster=generic('name',oshostname(),'np',3);
+md.verbose=verbose('1111111');
 
 %% solve for GIA deflection 
-md.cluster=generic('name',oshostname(),'np',3);
-md.verbose=verbose('1111111');
 md=solve(md,'Gia');
 
+%Test Name: GiaIvinsBenchmarksAB2dA1
+U_AB2dA1 = md.results.GiaSolution.UGia; 
+URate_AB2dA1 = md.results.GiaSolution.UGiaRate; 
+
+%Test Name: GiaIvinsBenchmarksAB2dA2
+%% different evaluation time. {{{ 
+md.timestepping.start_time=2005100; % after 5 kyr of deglaciation 
+md.geometry.thickness(end,end) = md.timestepping.start_time; 
+
+md=solve(md,'Gia');
+
+U_AB2dA2 = md.results.GiaSolution.UGia; 
+URate_AB2dA2 = md.results.GiaSolution.UGiaRate; 
+% }}} 
+
+%Test Name: GiaIvinsBenchmarksAB2dA3
+%% different evaluation time. {{{ 
+md.timestepping.start_time=2010100; % after 10 kyr of deglaciation 
+md.geometry.thickness(end,end) = md.timestepping.start_time; 
+
+md=solve(md,'Gia');
+
+U_AB2dA3 = md.results.GiaSolution.UGia; 
+URate_AB2dA3 = md.results.GiaSolution.UGiaRate; 
+% }}} 
+
 %Fields and tolerances to track changes
-field_names     ={'UGia','UGiaRate'};
-field_tolerances={1e-13,1e-13};
-field_values={...
-        (md.results.GiaSolution.UGia),...
-        (md.results.GiaSolution.UGiaRate),...
-        };
+field_names     ={'U_AB2dA1','URate_AB2dA1','U_AB2dA2','URate_AB2dA2','U_AB2dA3','URate_AB2dA3'};
+field_tolerances={1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
+field_values={U_AB2dA1,URate_AB2dA1,U_AB2dA2,URate_AB2dA2,U_AB2dA3,URate_AB2dA3}; 
+