Changeset 1931

Timestamp:

Oct 14, 2005, 6:16:39 PM (20 years ago)

Author:

steve

Message:

 

File:

• inundation/analytical solutions/Analytical solution_oblique_shock_01.py (modified) (4 diffs)

inundation/analytical solutions/Analytical solution_oblique_shock_01.py

@@ -73 +73 @@
 
 def shock_hb(t):
-    return 2.5
+    return 1.5
 
 def shock_hh(t):
-    return 2.5
+    return 1.5
 
 def shock_pb(t):
-    return 10
+    return 5
 
 def shock_pp(t):
-    return 10
+    return 5
 
 
@@ -126 +126 @@
 domain.set_quantity('xmomentum',shock_p )
 
+class SetValueWhere:
 
+    def __init__(self,quantity,value=0,xrange=0.0):
+        self.quantity =  quantity
+        self.value = value
+        self.xrange = xrange
+
+    def __call__(self,x,y):
+        n = x.shape[0]
+        w = self.quantity.centroid_values
+        for i in range(n):
+            if x[i]<self.xrange:
+                w[i] = self.value
+        return w
+
+
+Stage = domain.quantities['stage']
+Xmom =  domain.quantities['xmomentum']
+Ymom =  domain.quantities['ymomentum']
 
 #for id, face in domain.boundary:
@@ -135 +153 @@
 import time
 t0 = time.time()
-for t in domain.evolve(yieldstep = 0.5, finaltime = 5.0):
+for t in domain.evolve(yieldstep = 0.1, finaltime = 6.0):
     domain.write_time()
     id = 3399
-    print domain.quantities['stage'].get_values(location='centroids',indices=[id])
-    print domain.quantities['xmomentum'].get_values(location='centroids',indices=[id])
-    print domain.quantities['ymomentum'].get_values(location='centroids',indices=[id])
+    print Stage.get_values(location='centroids',indices=[id])
+    print Xmom.get_values(location='centroids',indices=[id])
+    print Ymom.get_values(location='centroids',indices=[id])
+    vstage =  Stage.get_values(location='centroids',indices=[id])
+    vxmom  =  Xmom.get_values(location='centroids',indices=[id])
     id = 12719
-    print domain.quantities['stage'].get_values(location='centroids',indices=[id])
-    print domain.quantities['xmomentum'].get_values(location='centroids',indices=[id])
-    print domain.quantities['ymomentum'].get_values(location='centroids',indices=[id])
+    print Stage.get_values(location='centroids',indices=[id])
+    print Xmom.get_values(location='centroids',indices=[id])
+    print Ymom.get_values(location='centroids',indices=[id])
+    tclean = 2.0
+#    if t > tclean-0.09 and t < tclean+0.01:
+#        print 'Cleaning up Initial profile'
+#        setstage = SetValueWhere(quantity=Stage,value=vstage[0],xrange=10)
+#        setxmom  = SetValueWhere(quantity=Xmom,value=vxmom[0],xrange=10)
+#        Stage.set_values(setstage,location='centroids')
+#        Xmom.set_values(setxmom,location='centroids')
+#        Dnew = Dirichlet_boundary([vstage[0] , vxmom[0], 0.0])
+#        domain.set_boundary({'left': Dnew, 'right': T, 'top': R, 'bottom': R})
 
 print 'That took %.2f seconds' %(time.time()-t0)
@@ -155 +184 @@
 #im = ImageGrab.grab()
 #im.save("ccube.eps")
+