Changeset 4709

Timestamp:

Sep 6, 2007, 1:29:36 PM (17 years ago)

Author:

sexton

Message:

updates

Location:

anuga_validation/convergence_study

Files:

• animatesww2d_alt.py (modified) (4 diffs)
• convergence.py (modified) (2 diffs)

anuga_validation/convergence_study/animatesww2d_alt.py

@@ -71 +71 @@
     min_stages = zeros((n,m), Float)
     min_stage = 100.
-
+    fid = open('stage.csv', 'w')
     for k, location in enumerate(x_points):
         for i, t in enumerate(f.get_time()):
@@ -86 +86 @@
 
             stage[i,k] = w
+            s = '%.2f,' %(w)
+            fid.write(s)
             elevation[i,k] = z
             xmom[i,k] = uh
@@ -96 +98 @@
             max_stages[i,k] = max_stage
             min_stages[i,k] = min_stage
+        s = '\n'
+        fid.write(s)
 
     min_stage = min(min(min_stages))
@@ -106 +110 @@
     #max_vec = zeros(100, Float)
 
-    for i in range(0,n,20):
-        #max_vec[i] = max(stage[i,:])
-        plot(x_points,stage[i,:]) #x_points,max_vec,'+') 
-        xlabel('x (m)')
-        ylabel('stage (m)')
-        axis(stage_axis)
-        name = 'time_%i' %i
-        savefig(name)
-        title(name)
-        figs.append(name)
+    if sys.platform == 'win32':
+        for i in range(0,n,20):
+            #max_vec[i] = max(stage[i,:])
+            plot(x_points,stage[i,:]) #x_points,max_vec,'+') 
+            xlabel('x (m)')
+            ylabel('stage (m)')
+            axis(stage_axis)
+            name = 'time_%i' %i
+            savefig(name)
+            title(name)
+            figs.append(name)
 
     close('all')
-    
-animatesww2d('myexample2.sww')
+
+swwfile = join(getenv('INUNDATIONHOME'),'data','validation','convergence_study','5000','myexample2.sww')
+animatesww2d(swwfile)
 
 # once png files have been created, they can then be dragged into

anuga_validation/convergence_study/convergence.py

@@ -44 +44 @@
                          boundary_tags={'left': [0], 'top': [1],
                                         'right': [2], 'bottom': [3]},
-                         maximum_triangle_area = 25000,
+                         maximum_triangle_area = 1000,
                          filename=meshname)
 
@@ -85 +85 @@
 #------------------------------------------------------------------------------
 
-for t in domain.evolve(yieldstep = 5.0, finaltime = 40*60.):
+for t in domain.evolve(yieldstep = 20.0, finaltime = 40*60.):
     domain.write_time()