Changeset 3915 for anuga_validation/okushiri_2005/create_okushiri.py

Timestamp:

Nov 3, 2006, 5:15:13 PM (18 years ago)

Author:

ole

Message:

Okushiri demo

File:

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

anuga_validation/okushiri_2005/create_okushiri.py

@@ -153 +153 @@
 
     # 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.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.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_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,
@@ -179 +179 @@
 
     # 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]
+    #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]
-
 
     
     # Interior regions and creation of mesh
     interior_regions = [[rhs_region, 0.0005],
-                        [island, 0.0003*base_resolution],
-                        [gulleys, 0.00003*base_resolution]]    
+                        [island, 0.0002*base_resolution],
+                        [gulleys, 0.00002*base_resolution]]    
 
     meshname = project.mesh_filename + '.msh'
@@ -206 +205 @@
                                  interior_regions=interior_regions,
                                  filename=project.mesh_filename,
-                                 use_cache=True,
                                  verbose=True)