Changeset 3636 for anuga_validation/okushiri_2005/create_mesh.py

Timestamp:

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

Author:

ole

Message:

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

File:

• anuga_validation/okushiri_2005/create_mesh.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')