Changeset 23762

Timestamp:

02/27/19 17:09:05 (6 years ago)

Author:

schlegel

Message:

CHG: add some complexity to timestepping for restart test

Location:

issm/trunk-jpl/test/NightlyRun

Files:

• test125.m (modified) (3 diffs)
• test125.py (modified) (3 diffs)

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

@@ -8 +8 @@
 
 md.verbose=verbose('solution',true);
-md.settings.recording_frequency=5;
+md.settings.recording_frequency=4;
 
 % time steps and resolution
-md.timestepping.final_time=8;
+md.timestepping.final_time=19;
+md.settings.output_frequency=2;
 
 md=solve(md,'Transient');
@@ -22 +23 @@
                                                 1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13,1e-13};
 field_values={...
-        (md.results.TransientSolution(6).Vx)-(md2.results.TransientSolution(6).Vx),...
-        (md.results.TransientSolution(6).Vy)-(md2.results.TransientSolution(6).Vy),...
-        (md.results.TransientSolution(6).Vel)-(md2.results.TransientSolution(6).Vel),...
-        (md.results.TransientSolution(6).TotalSmb)-(md2.results.TransientSolution(6).TotalSmb),...
-        (md.results.TransientSolution(6).Base)-(md2.results.TransientSolution(6).Base),...
-        (md.results.TransientSolution(6).Surface)-(md2.results.TransientSolution(6).Surface),...
-        (md.results.TransientSolution(6).Thickness)-(md2.results.TransientSolution(6).Thickness),...
-        (md.results.TransientSolution(6).IceVolume)-(md2.results.TransientSolution(6).IceVolume),...
         (md.results.TransientSolution(7).Vx)-(md2.results.TransientSolution(7).Vx),...
         (md.results.TransientSolution(7).Vy)-(md2.results.TransientSolution(7).Vy),...
@@ -46 +39 @@
         (md.results.TransientSolution(8).Thickness)-(md2.results.TransientSolution(8).Thickness),...
         (md.results.TransientSolution(8).IceVolume)-(md2.results.TransientSolution(8).IceVolume),...
+        (md.results.TransientSolution(9).Vx)-(md2.results.TransientSolution(9).Vx),...
+        (md.results.TransientSolution(9).Vy)-(md2.results.TransientSolution(9).Vy),...
+        (md.results.TransientSolution(9).Vel)-(md2.results.TransientSolution(9).Vel),...
+        (md.results.TransientSolution(9).TotalSmb)-(md2.results.TransientSolution(9).TotalSmb),...
+        (md.results.TransientSolution(9).Base)-(md2.results.TransientSolution(9).Base),...
+        (md.results.TransientSolution(9).Surface)-(md2.results.TransientSolution(9).Surface),...
+        (md.results.TransientSolution(9).Thickness)-(md2.results.TransientSolution(9).Thickness),...
+        (md.results.TransientSolution(9).IceVolume)-(md2.results.TransientSolution(9).IceVolume),...
         };

issm/trunk-jpl/test/NightlyRun/test125.py

@@ -20 +20 @@
 
 md.verbose=verbose('solution',1)
-md.settings.recording_frequency=5
+md.settings.recording_frequency=4
 
 # time steps and resolution
-md.timestepping.final_time=8
+md.timestepping.final_time=19
+md.settings.output_frequency=2
 
 md=solve(md,'Transient')
@@ -36 +37 @@
                 1e-13,1e-13,1e-13,1e-13,1e-13,1e-13]
 field_values=[\
-                md2.results.TransientSolution[5].Vx-md.results.TransientSolution[5].Vx,\
-                md2.results.TransientSolution[5].Vy-md.results.TransientSolution[5].Vy,\
-                md2.results.TransientSolution[5].Vel-md.results.TransientSolution[5].Vel,\
-                md2.results.TransientSolution[5].TotalSmb-md.results.TransientSolution[5].TotalSmb,\
-                md2.results.TransientSolution[5].Base-md.results.TransientSolution[5].Base,\
-                md2.results.TransientSolution[5].Surface-md.results.TransientSolution[5].Surface,\
-                md2.results.TransientSolution[5].Thickness-md.results.TransientSolution[5].Thickness,\
-                md2.results.TransientSolution[5].IceVolume-md.results.TransientSolution[5].IceVolume,\
                 md2.results.TransientSolution[6].Vx-md.results.TransientSolution[6].Vx,\
                 md2.results.TransientSolution[6].Vy-md.results.TransientSolution[6].Vy,\
@@ -59 +52 @@
                 md2.results.TransientSolution[7].Surface-md.results.TransientSolution[7].Surface,\
                 md2.results.TransientSolution[7].Thickness-md.results.TransientSolution[7].Thickness,\
-                md2.results.TransientSolution[7].IceVolume-md.results.TransientSolution[7].IceVolume\
+                md2.results.TransientSolution[7].IceVolume-md.results.TransientSolution[7].IceVolume,\
+                md2.results.TransientSolution[8].Vx-md.results.TransientSolution[8].Vx,\
+                md2.results.TransientSolution[8].Vy-md.results.TransientSolution[8].Vy,\
+                md2.results.TransientSolution[8].Vel-md.results.TransientSolution[8].Vel,\
+                md2.results.TransientSolution[8].TotalSmb-md.results.TransientSolution[8].TotalSmb,\
+                md2.results.TransientSolution[8].Base-md.results.TransientSolution[8].Base,\
+                md2.results.TransientSolution[8].Surface-md.results.TransientSolution[8].Surface,\
+                md2.results.TransientSolution[8].Thickness-md.results.TransientSolution[8].Thickness,\
+                md2.results.TransientSolution[8].IceVolume-md.results.TransientSolution[8].IceVolume\
                 ]