Changeset 1707

Timestamp:

Aug 11, 2005, 2:49:45 PM (19 years ago)

Author:

ole

Message:

 

Location:

inundation/ga/storm_surge/validation/LWRU2

Files:

• create_mesh.py (modified) (2 diffs)
• lwru2.py (modified) (3 diffs)

inundation/ga/storm_surge/validation/LWRU2/create_mesh.py

@@ -75 +75 @@
     m.addVertsSegs(dict)
 
+
+
+    #Localised refined area for gulleys
+    dict = {}
+    dict['points'] = [ [4.5, 1.8],
+                       [5.0, 1.8],
+                       [5.0, 2.5],
+                       [4.5, 2.5] ]                       
+                       
+    dict['segments'] = [[0,1], [1,2], [2,3], [3,0]]
+    dict['segment_tags'] = ['', '', '', '']
+    m.addVertsSegs(dict)    
+
     
     base_resolution = 1 
@@ -84 +97 @@
     mid.setMaxArea(0.001*base_resolution)
     
+
     inner = m.addRegionEN(point_mr2[0]+1, ybottom+1)
-    inner.setMaxArea(0.00001*base_resolution)
+    inner.setMaxArea(0.0001*base_resolution)
+
+    gulleys = m.addRegionEN(dict['points'][0]+0.1, dict['points'][1]+0.1)
+    gulleys.setMaxArea(0.00001*base_resolution)    
+    
     
     m.generateMesh('pzq28.0za1000000a')

inundation/ga/storm_surge/validation/LWRU2/lwru2.py

@@ -81 +81 @@
 
 try:
-    raise Exception
+    #raise Exception
     fid = open(picklefile)
     print 'Read pickled domain'
@@ -130 +130 @@
                                                   points,
                                                   attributes['elevation'],
-                                                  alpha = 0,
+                                                  alpha = 0.0001,
                                                   verbose=True)
 
@@ -159 +159 @@
 
 #Set boundary conditions
-domain.set_boundary({'left': Bts, 'right': Br, 'bottom': Br, 'top': Br})
+#domain.set_boundary({'left': Bts, 'right': Br, 'bottom': Br, 'top': Br})
+domain.set_boundary({'wave': Bts, 'wall': Br})
 
 #Evolve
 import time
 t0 = time.time()
-domain.visualise = True
+#domain.visualise = True
 
 for t in domain.evolve(yieldstep = 0.05, finaltime = 22.5):