Changeset 3913 for anuga_validation/okushiri_2005/create_okushiri.py

Timestamp:

Nov 3, 2006, 11:35:32 AM (17 years ago)

Author:

ole

Message:

Got similar number of triangles as Aug 2005 okushiri validation

File:

• anuga_validation/okushiri_2005/create_okushiri.py (modified) (3 diffs)

anuga_validation/okushiri_2005/create_okushiri.py

@@ -152 +152 @@
     
 
-    # Island area and drawdown region
+    # Island area and drawdown region (original)
     island_0 = [xleft + 2*(xright-xleft)/3+1.2, ytop-0.5]
     island_1 = [xleft + 2*(xright-xleft)/3+0.5, ybottom + 2*(ytop-ybottom)/3]
-    island_2 = [xleft + (xright-xleft)/2+0.4, ybottom + 2*(ytop-ybottom)/3-0.3]
-    island_3 = [xleft + (xright-xleft)/2+0.4, ybottom + (ytop-ybottom)/3+0.3]
-    island_4 = [xleft + 2*(xright-xleft)/3+0.4, ybottom + (ytop-ybottom)/3-0.3]
-    island_5 = [xleft + 2*(xright-xleft)/3+1.2, ybottom+0.8]
-    island_6 = [xl-.01, yb]  # Keep right edge just off the gulleys
-    island_7 = [xl-.01, yt]
+    island_2 = [xleft + (xright-xleft)/2+0.3, ybottom + 2*(ytop-ybottom)/3-0.3]
+    island_3 = [xleft + (xright-xleft)/2+0.3, ybottom + (ytop-ybottom)/3+0.3]
+    island_4 = [xleft + 2*(xright-xleft)/3+0.4, ybottom + (ytop-ybottom)/3-0.3] 
+    island_5 = [xleft + 2*(xright-xleft)/3+1.2, ybottom+0.2]
+    island_6 = [xl-.01, yb]  #OK
+    island_7 = [xl-.01, yt]  #OK      
+
+
+    # Island area and drawdown region
+    #island_0 = [xleft + 2*(xright-xleft)/3+1.2, ytop-0.5]
+    #island_1 = [xleft + 2*(xright-xleft)/3+0.5, ybottom + 2*(ytop-ybottom)/3]
+    #island_2 = [xleft + (xright-xleft)/2+0.4, ybottom + 2*(ytop-ybottom)/3-0.3]
+    #island_3 = [xleft + (xright-xleft)/2+0.4, ybottom + (ytop-ybottom)/3+0.3]
+    #island_4 = [xleft + 2*(xright-xleft)/3+0.4, ybottom + (ytop-ybottom)/3-0.3]
+    #island_5 = [xleft + 2*(xright-xleft)/3+1.2, ybottom+0.8]
+    #island_6 = [xl-.01, yb]  # Keep right edge just off the gulleys
+    #island_7 = [xl-.01, yt]
  
     island = [island_0, island_1, island_2,
@@ -167 +178 @@
 
 
+    # Region spanning half right hand side of domain just inside boundary (org)
+    rhs_nw = [xleft + (xright-xleft)/3+1, ytop-0.02]
+    rhs_sw = [xleft + (xright-xleft)/3+1, ybottom+0.02]
+    rhs_se = [xright-0.02, ybottom+0.02]
+    rhs_ne = [xright-0.02, ytop-0.02]
+
     # Region spanning half right hand side of domain just inside boundary
-    rhs_nw = [xleft + (xright-xleft)/3+1, ytop-1.4]
-    rhs_sw = [xleft + (xright-xleft)/3+1, ybottom+0.5]
-    rhs_se = [xright-0.1, ybottom+0.2]
-    rhs_ne = [xright-0.1, ytop-0.2]    
+    #rhs_nw = [xleft + (xright-xleft)/3+1, ytop-1.4]
+    #rhs_sw = [xleft + (xright-xleft)/3+1, ybottom+0.5]
+    #rhs_se = [xright-0.1, ybottom+0.2]
+    #rhs_ne = [xright-0.1, ytop-0.2]        
 
     rhs_region = [rhs_nw, rhs_ne, rhs_se, rhs_sw]
@@ -179 +196 @@
     # Interior regions and creation of mesh
     interior_regions = [[rhs_region, 0.0005],
-                        [gulleys, 0.00002*base_resolution],
-                        [island, 0.00007*base_resolution]] 
+                        [island, 0.0003*base_resolution],
+                        [gulleys, 0.00003*base_resolution]]    
 
     meshname = project.mesh_filename + '.msh'