source: issm/oecreview/Archive/14312-15392/ISSM-14841-14842.diff

../trunk-jpl/test/NightlyRun/test2001.m

@@ -3 +3 @@
 md=setmask(md,'','');
 md=parameterize(md,'../Par/SquareSheetConstrained.par');
 
+%% indicate what you want to compute 
+md.gia.output_rates=1;           % just want "w" solution 
+md.gia.cross_section_shape=1;    % for square-edged x-section 
+
 %% define loading history 
 md.timestepping.start_time=2400000; %2,400 kyr
-md.timestepping.final_time=2500000; %2,400 kyr
+md.timestepping.final_time=2500000; %2,500 kyr
+% In order to run Ivins1999 benchmarks, for now, we do things manually as follows: 
+% 1. numtimes is fixed to be five 
+% 2. final_time given above is hard-coded in "GiaDeflectionCorex.cpp"
+% 3. evaluation time = start_time 
+% 4. define new load times, according to the chosen evaluation time (care that general load history is not altered). This is really important, because of #1 and also becasue of the fact that futre loading times are not recognized.
+% 5. Ice thickness assoiated with evaluation time is actually the same at final_time
 md.geometry.thickness=[...
-        zeros(md.mesh.numberofvertices+1,1),...
-        [md.geometry.thickness; 1.0e-4],...
+        [md.geometry.thickness*0.0; 0.0],...
+        [md.geometry.thickness/2.0; 0.1],...
+        [md.geometry.thickness; 0.2],...
+        [md.geometry.thickness; 1.0],...
         [md.geometry.thickness; md.timestepping.start_time],...
         ];