Changeset 1632 for inundation/ga/storm_surge/validation/LWRU1/lwru1.py

Timestamp:

Jul 22, 2005, 5:59:19 PM (20 years ago)

Author:

ole

Message:

Interpolation work and some refactoring

File:

• inundation/ga/storm_surge/validation/LWRU1/lwru1.py (modified) (5 diffs)

inundation/ga/storm_surge/validation/LWRU1/lwru1.py

@@ -32 +32 @@
 #(only about 200m is expected, though)
 
-#points, vertices, boundary = rectangular(400, 40,
+points, vertices, boundary = rectangular(400, 40,
+                                         len1=55000, len2=5000,
+                                         origin = (-5000, 0.0))
+
+#points, vertices, boundary = rectangular(100, 10,
 #                                         len1=55000, len2=5000,
 #                                         origin = (-5000, 0.0))
-
-points, vertices, boundary = rectangular(100, 10,
-                                         len1=55000, len2=5000,
-                                         origin = (-5000, 0.0))
 
 
@@ -90 +90 @@
         #plot(self.x, self.w)
         #show()
-
+        #import sys; sys.exit() 
 
     def __call__(self, x, y):
@@ -154 +154 @@
 
 #Set boundary conditions
-domain.set_boundary({'left': Bd, 'right': Bd, 'bottom': Bd, 'top': Bd})
+domain.set_boundary({'left': Br, 'right': Bt, 'bottom': Bt, 'top': Bt})
 
 
@@ -167 +167 @@
 xes = []
 y = 2500 
-x0 = -5000
+x0 = -100
 step = 50
 for i in range(1000):
@@ -176 +176 @@
 from pylab import *
 from pyvolution.least_squares import Interpolation
-I = Interpolation(domain.coordinates,
-                  domain.triangles,
+
+ 
+V = domain.get_vertex_coordinates(obj=True)
+T = domain.get_triangles(obj=True)
+
+I = Interpolation(V,
+                  T,
                   point_coordinates = pt,
                   verbose = True)
 
 
+f = domain.quantities['elevation'].vertex_values.flat
+z = I.interpolate( f )
+
 print 'xxxxx'
 
+ion()
 hold(False)
-for t in domain.evolve(yieldstep = 1, finaltime = 220.0):
-
-    f = domain.quantities['stage'].get_vertex_values(smooth = True)
+f = domain.quantities['stage'].vertex_values.flat
+y = I.interpolate( f )
+plot(xes, y, '-b', xes, z, '-k')
+set( gca(), Ylim=(-10,10) )
+
+raw_input('go')
+for t in domain.evolve(yieldstep = 1, finaltime = 300.0):
+
+    f = domain.quantities['stage'].vertex_values.flat
     y = I.interpolate( f )
-
+    ioff()
+    plot(xes, y, '-b', xes, z, '-k')
     ion()
-    plot(xes, y)    
+    set( gca(), Ylim=(-10,10) )
+    
+    
+    #print y[:], y.shape
+
+
     domain.write_time()