Changeset 1708

Timestamp:

Aug 11, 2005, 3:15:54 PM (18 years ago)

Author:

ole

Message:

Added refined region around gulleys to mesh

Location:

inundation/ga/storm_surge/validation/LWRU2

Files:

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

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

@@ -79 +79 @@
     #Localised refined area for gulleys
     dict = {}
-    dict['points'] = [ [4.5, 1.8],
-                       [5.0, 1.8],
-                       [5.0, 2.5],
-                       [4.5, 2.5] ]                       
-                       
+    xl = 4.8
+    xr = 5.3
+    yb = 1.6
+    yt = 2.3
+    p0 = [xl, yb]
+    p1 = [xr, yb]
+    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'] = ['', '', '', '']
@@ -89 +96 @@
 
     
-    base_resolution = 1 
+    base_resolution = 100 
 
     ocean = m.addRegionEN(xleft+1, ybottom+1)
@@ -101 +108 @@
     inner.setMaxArea(0.0001*base_resolution)
 
-    gulleys = m.addRegionEN(dict['points'][0]+0.1, dict['points'][1]+0.1)
+
+    gulleys = m.addRegionEN(p0[0]+0.1, p0[1]+0.1)
     gulleys.setMaxArea(0.00001*base_resolution)    
     

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

@@ -81 +81 @@
 
 try:
-    #raise Exception
+    raise Exception
     fid = open(picklefile)
     print 'Read pickled domain'
@@ -169 +169 @@
 for t in domain.evolve(yieldstep = 0.05, finaltime = 22.5):
     domain.write_time()
-    Stage_cv = domain.quantities['stage'].centroid_values
-    print Stage_cv[0]
-    print 'Function value',function(t)
+    #Stage_cv = domain.quantities['stage'].centroid_values
+    #print Stage_cv[0]
+    #print 'Function value',function(t)
 
 print 'That took %.2f seconds' %(time.time()-t0)