Changeset 3636

Timestamp:

Sep 20, 2006, 4:22:56 PM (17 years ago)

Author:

ole

Message:

Work on Okushiri validation (towards recreating 18 Aug 2005 results)

Location:

anuga_validation/okushiri_2005

Files:

• create_mesh.py (modified) (2 diffs)
• extract_timeseries.py (modified) (2 diffs)

anuga_validation/okushiri_2005/create_mesh.py

@@ -113 +113 @@
     
 #
-
-    #Fine resolution. Needs alpha >= 0.001    
+    
     base_resolution = 1
 
@@ -122 +121 @@
     mid = m.addRegionEN(point_mbottom[0]+.1, ybottom+.1)
     mid.setMaxArea(0.0005*base_resolution)
-    
+  
 
     inner = m.addRegionEN(island_3[0]+.1, island_3[1]+.1)
-    inner.setMaxArea(0.00007*base_resolution)
+    inner.setMaxArea(0.0003*base_resolution)
 
 
     gulleys = m.addRegionEN(p0[0]+0.1, p0[1]+0.1)
-    gulleys.setMaxArea(0.00002*base_resolution)
+    gulleys.setMaxArea(0.00005*base_resolution)
 
 
+    # From r 1709 11 August 2005
+    #base_resolution = 100
+    #
+    #ocean = m.addRegionEN(xleft+.1, ybottom+.1)
+    #ocean.setMaxArea(0.01*base_resolution)
+    #
+    #mid = m.addRegionEN(point_mbottom[0]+.1, ybottom+.1)
+    #mid.setMaxArea(0.001*base_resolution)
 
+    #inner = m.addRegionEN(island_3[0]+.1, island_3[1]+.1)
+    #inner.setMaxArea(0.0001*base_resolution)
+
+
+    #gulleys = m.addRegionEN(p0[0]+0.1, p0[1]+0.1)
+    #gulleys.setMaxArea(0.00001*base_resolution)        
+    
+    
     m.generateMesh('pzq28.0za1000000a')
 

anuga_validation/okushiri_2005/extract_timeseries.py

@@ -4 +4 @@
 import sys
 from os import sep
-from caching import cache
+from anuga.caching import cache
 
 sys.path.append('..'+sep+'..'+sep)
@@ -51 +51 @@
 
 
-from anuga.pyvolution.util import file_function
+from anuga.abstract_2d_finite_volumes.util import file_function
 from anuga.utilities.numerical_tools import ensure_numeric
 gauge_values = [ensure_numeric(input_stage),