Changeset 1742 for inundation/validation/LWRU2/create_mesh.py

Timestamp:

Aug 23, 2005, 2:56:21 PM (20 years ago)

Author:

ole

Message:

Recovered validation stuff from 18 August + wrapping up

File:

• inundation/validation/LWRU2/create_mesh.py (modified) (5 diffs)

inundation/validation/LWRU2/create_mesh.py

@@ -14 +14 @@
 
     #Basic geometry
-
+    
     xleft   = 0
     xright  = 5.448
@@ -20 +20 @@
     ytop    = 3.402
 
-    #Outline
+    #Outline 
     point_sw = [xleft, ybottom]
     point_se = [xright, ybottom]
-    point_nw = [xleft, ytop]
+    point_nw = [xleft, ytop]    
     point_ne = [xright, ytop]
 
     #Midway points (left)
-    point_ml1 = [xleft + (xright-xleft)/3+1, ytop]
-    point_ml2 = [xleft + (xright-xleft)/3+1, ybottom]
-
-    #Midway points (right)
-    point_mr1 = [xleft + 2*(xright-xleft)/3+0.5, ytop]
-    point_mr2 = [xleft + 2*(xright-xleft)/3+0.5, ybottom]
+    point_mtop = [xleft + (xright-xleft)/3+1, ytop]
+    point_mbottom = [xleft + (xright-xleft)/3+1, ybottom]
 
 
     geo = Geo_reference(xllcorner = xleft,
                         yllcorner = ybottom)
-
-
+    
+                        
     print "***********************"
     print "geo ref", geo
     print "***********************"
-
+    
     m = Mesh(geo_reference=geo)
 
@@ -49 +45 @@
     dict['points'] = [point_se,   #se
                       point_ne,
-                      point_mr1,
-                      point_ml1,
+                      point_mtop,
                       point_nw,
                       point_sw,
-                      point_ml2,
-                      point_mr2]
+                      point_mbottom]
 
-
+    
     dict['segments'] = [[0,1], [1,2], [2,3], [3,4],
-                        [4,5], [5,6], [6,7], [7,0], #The outer border
-                        [2,7], [3,6]]               #Separators
-
+                        [4,5], [5,0],  #The outer border
+                        [2,5]]         #Separator
+    
     dict['segment_tags'] = ['wall',
-                            'wall',
                             'wall',
                             'wall',
@@ -68 +61 @@
                             'wall',
                             'wall',
-                            'wall',
-                            '',           #Interior
                             '']           #Interior
 
+        
+    m.addVertsSegs(dict)
 
-    m.addVertsSegs(dict)
+
 
 
@@ -87 +80 @@
     p2 = [xr, yt]
     p3 = [xl, yt]
+    
+    dict['points'] = [p0, p1, p2, p3]
+    dict['segments'] = [[0,1], [1,2], [2,3], [3,0]]
+    dict['segment_tags'] = ['', '', '', '']
+    m.addVertsSegs(dict)    
 
-    dict['points'] = [p0, p1, p2, p3]
+    #Island area
+    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      
+    
+
+    dict['points'] = [island_0, island_1, island_2,
+                      island_3, island_4, island_5,
+                      #p0, p3]                      
+                      island_6, island_7]
+
+                      
+    dict['segments'] = [[0,1], [1,2], [2,3], [3,4],
+                        [4,5], [5,6], [6,7], [7,0]]
+                        
+    dict['segment_tags'] = ['', '', '', '', '', '', '', '']
 
 
-    dict['segments'] = [[0,1], [1,2], [2,3], [3,0]]
-    dict['segment_tags'] = ['', '', '', '']
-    m.addVertsSegs(dict)
+    m.addVertsSegs(dict)    
 
+    
+#
+    
+    base_resolution = 1
 
-    base_resolution = 100
+    ocean = m.addRegionEN(xleft+.1, ybottom+.1)
+    ocean.setMaxArea(0.1*base_resolution)
 
-    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.0005*base_resolution)
+    
 
-    mid = m.addRegionEN(point_ml2[0]+1, ybottom+1)
-    mid.setMaxArea(0.001*base_resolution)
-
-
-    inner = m.addRegionEN(point_mr2[0]+1, ybottom+1)
-    inner.setMaxArea(0.0001*base_resolution)
+    inner = m.addRegionEN(island_3[0]+.1, island_3[1]+.1)
+    inner.setMaxArea(0.00007*base_resolution)
 
 
     gulleys = m.addRegionEN(p0[0]+0.1, p0[1]+0.1)
-    gulleys.setMaxArea(0.00001*base_resolution)
-
-
+    gulleys.setMaxArea(0.00002*base_resolution)    
+    
+    
     m.generateMesh('pzq28.0za1000000a')
 
     import filenames
     m.export_mesh_file(filenames.mesh_filename)
+