Changeset 1710

Timestamp:

Aug 11, 2005, 3:45:58 PM (18 years ago)

Author:

steve

Message:

Analytical solution_wave_runup.py can now run strang file

File:

• inundation/ga/storm_surge/analytical solutions/Analytical solution_wave_runup.py (modified) (9 diffs)

inundation/ga/storm_surge/analytical solutions/Analytical solution_wave_runup.py

@@ -5 +5 @@
    Christopher Zoppou, Stephen Roberts, Ole Nielsen, Duncan Gray
    Geoscience Australia
-   
+
 Specific methods pertaining to the 2D shallow water equation
 are imported from shallow_water
@@ -15 +15 @@
 
 ######################
-# Module imports 
+# Module imports
 import sys
 from os import sep
@@ -39 +39 @@
 #points, elements = strang_mesh('Run-up.pt')
 points, elements = strang_mesh('strang_7389.pt')
-domain = Domain(points, elements) 
+domain = Domain(points, elements)
 
 domain.default_order = 2
@@ -49 +49 @@
 for id, face in domain.boundary:
     domain.boundary[(id,face)] = 'external'
-    
-domain.boundary[(0,0)] = 'left'
+    if domain.get_vertex_coordinate(id,(face+1)%3)[0] < -200.0+1.0e-10:
+        domain.boundary[(id,face)] = 'left'
+
+
 
 
@@ -60 +62 @@
 print "Number of triangles = ", len(domain)
 
-#Reduction operation for get_vertex_values              
+#Reduction operation for get_vertex_values
 from util import mean
-domain.reduction = mean 
+domain.reduction = mean
 #domain.reduction = min  #Looks better near steep slopes
 
-#Define the boundary condition 
+#Define the boundary condition
 def stage_setup(x,t_star1):
     vh = 0
@@ -84 +86 @@
         sigma_prime = sigma0/a
         const = (1.-ii*lambda_prime)**2+sigma_prime**2
-          
+
         v1 = 8.*eta/a*imag(1./const**(3./2.) \
             -3./4.*(1.-ii*lambda_prime)/const**(5./2.))
-            
+
         x1 = -v1**2/2.-a**2*sigma_prime**2/16.+eta \
             *real(1.-2.*(5./4.-ii*lambda_prime) \
@@ -149 +151 @@
             h = 0
             w[i] = z[i]
-    return w   
+    return w
 
 domain.set_quantity('stage', stage)
@@ -161 +163 @@
 
 
+
+#domain.visualise = True
+
 ######################
 #Evolution
@@ -168 +173 @@
     domain.write_time()
     print boundary_stage(domain.time)
-    
+
 print 'That took %.2f seconds' %(time.time()-t0)
 
-