Changeset 647

Timestamp:

Dec 1, 2004, 2:25:43 PM (20 years ago)

Author:

chris

Message:

 

Location:

inundation/ga/storm_surge/analytical solutions

Files:

• Analytical_solution_MacDonald_import_mesh.py (modified) (3 diffs)
• Sydney.py (added)
• sydney_ubd.tsh (added)

inundation/ga/storm_surge/analytical solutions/Analytical_solution_MacDonald_import_mesh.py

@@ -8 +8 @@
    Copyright 2004
    Christopher Zoppou, Stephen Roberts, Ole Nielsen, Duncan Gray
-   Geoscience Australia
+   Geoscience Australia, ANU
    
 Specific methods pertaining to the 2D shallow water equation
@@ -83 +83 @@
 #--------------------
 # Boundary Conditions
-upstream_depth = 5.035-4.393
+upstream_depth = 5.035 - 4.393
+downstream_depth = 1.5 - 0
 discharge = 20
 tags = {}
-tags['Upstream'] = Dirichlet_boundary([5.035, 2, 0.0])
+tags['Upstream'] = Dirichlet_boundary([upstream_depth, 2, 0.0])
 tags['external'] = Reflective_boundary(domain) 
-tags['Downstream'] = Transmissive_boundary(domain)
+tags['Downstream'] = Dirichlet_boundary([downstream_depth, 2, 0.0])
 domain.set_boundary(tags)
 
@@ -104 +105 @@
 import time
 t0 = time.time()
-for t in domain.evolve(yieldstep = 50, finaltime = 5000):
+for t in domain.evolve(yieldstep = 5, finaltime = 500):
     domain.write_time()