Changeset 5741 for anuga_work/development/demos

Timestamp:

Sep 5, 2008, 9:16:45 PM (17 years ago)

Author:

steve

Message:

Updated timestepping

File:

• anuga_work/development/demos/netherlands.py (modified) (5 diffs)

anuga_work/development/demos/netherlands.py

@@ -27 +27 @@
 N = 130 # size = 33800
 N = 600 # Size = 720000
-N = 100
+N = 50
 
 points, elements, boundary = rectangular_cross(N, N)
@@ -34 +34 @@
 domain.check_integrity()
 
+import sys
+base = os.path.basename(sys.argv[0])
+domain.simulation_name = 'netherlands'
 domain.set_name(os.path.splitext(__file__)[0])
-domain.set_timestepping_method('rk3')
+domain.set_timestepping_method('rk2')
 domain.set_default_order(2)
 domain.set_store_vertices_uniquely(True) # Store as internally represented
@@ -130 +133 @@
 domain.set_boundary({'left': Bd, 'right': Br, 'bottom': Br, 'top': Br})
 
+#-------------------------------------------------------------------------------
+# Copy scripts to time stamped output directory and capture screen
+# output to file
+#-------------------------------------------------------------------------------
+time = strftime('%Y%m%d_%H%M%S',localtime())
+
+output_dir = 'dam_break_'+time
+output_file = 'dam_break'
+
+copy_code_files(output_dir,__file__)
+#start_screen_catcher(output_dir+'_')
+
 
 #------------------------------------------------------------------------------
@@ -139 +154 @@
     # Initialise real-time visualiser 
 
-    pass
-    #vis = RealtimeVisualiser(domain)
-    #vis.render_quantity_height("elevation", dynamic=False)
-    #vis.render_quantity_height("stage", dynamic=True)
-    #vis.colour_height_quantity('stage', (0.0, 0.0, 0.8))
-    #vis.start()
+    #pass
+    visualise = True
+    if visualise:
+        from anuga.visualiser import RealtimeVisualiser
+        vis = RealtimeVisualiser(domain)
+        vis.render_quantity_height("elevation", offset=0.01, dynamic=False)
+        vis.render_quantity_height("stage", dynamic=True)
+        vis.colour_height_quantity('stage', (0.2, 0.2, 0.8))
+        vis.start()
+        import time
+        time.sleep(2.0)
+    
+    
 
 
@@ -151 +173 @@
 t0 = time.time()
 
-for t in domain.evolve(yieldstep = 0.005, finaltime = None):
+for t in domain.evolve(yieldstep = 0.05, finaltime = 5.0):
     print domain.timestepping_statistics()
-    print domain.quantities['stage'].get_values(location='centroids',
-                                                indices=[0])
-                                                
-    #vis.update()                                               
-    #time.sleep(0.1)
-    #raw_input('pause>')
-    #V.update_quantity('stage')
-    #rpdb.set_active()
-    #integral_label.text='Integral=%10.5e'%domain.quantities['stage'].get_integral()
-#vis.evolveFinished()
+    #print domain.quantities['stage'].get_values(location='centroids',
+                                                #indices=[0])
+    if visualise:
+        vis.update()                                            
+
+        
+
+if visualise: vis.evolveFinished()
 
 print 'That took %.2f seconds' %(time.time()-t0)