Changeset 1353 for inundation/ga/storm_surge/analytical solutions

Timestamp:

May 9, 2005, 12:33:09 PM (20 years ago)

Author:

chris

Message:

 

Location:

inundation/ga/storm_surge/analytical solutions

Files:

• Analytical solution_wave_runup.py (modified) (1 diff)
• Analytical_solution_contracting_channel_import_mesh.py (modified) (5 diffs)
• MacDonald_77541.tsh (modified) (1 diff)
• MacDonald_first_order_size77541.sww (added)
• Merimbula_matlab_77668.fig (added)
• Sydney.py (modified) (2 diffs)
• Traditional_asymetrical_dam_break_second_order_dry_size5093.sww (added)
• sampson_strang_mesh.pt (added)

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

@@ -55 +55 @@
 domain.store = True
 domain.format = 'sww'
-domain.filename = 'run-up_second_order'
+domain.filename = 'run-up_second_order_again'
 
 print "Number of triangles = ", len(domain)

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

@@ -6 +6 @@
 for use with the generic finite volume framework
 
-   Copyright 2004
-   Christopher Zoppou, Stephen Roberts, Ole Nielsen, Duncan Gray
-   Geoscience Australia
+   Copyright 2005
+   Christopher Zoppou, Stephen Roberts
+   ANU
    
 Specific methods pertaining to the 2D shallow water equation
@@ -18 +18 @@
 """
 
-#-----------------------------
+#---------------
 # Module imports 
 import sys
@@ -29 +29 @@
 from pmesh2domain import pmesh_to_domain_instance
 
-
-#------------------------------
+#-------
 # Domain 
 filename = 'converging_channel_30846.tsh'
@@ -37 +36 @@
 print 'Number of triangles = ', len(domain)
 
-domain.default_order = 1
+#----------------
+# Order of scheme
+domain.default_order = 2
 domain.smooth = True
 
+#-------------------------------------
 # Provide file name for storing output
 domain.store = True     #Store for visualisation purposes
 domain.format = 'sww'   #Native netcdf visualisation format
-domain.filename = 'contracting_channel_first-order'
+domain.filename = 'contracting_channel_second-order'
 
-#Decide which quantities are to be stored at each timestep
+#----------------------------------------------------------
+# Decide which quantities are to be stored at each timestep
 domain.quantities_to_be_stored = ['stage', 'xmomentum', 'ymomentum']
 
-#Reduction operation for get_vertex_values              
+#------------------------------------------
+# Reduction operation for get_vertex_values              
 from util import mean
 domain.reduction = mean 
 
-#------------------------------
-# Boundary Conditions
+#------------------------
+# Set boundary Conditions
 tags = {}
 tags['upstream'] = Dirichlet_boundary([0.2, 1.2, 0.0])
@@ -60 +64 @@
 domain.set_boundary(tags)
 
-#-----------------
-#Initial condition
+#----------------------
+# Set initial condition
 domain.set_quantity('elevation', 0.0)
 domain.set_quantity('stage', 0.2)
     
-
-#-----------------
-#Evolution
+#----------
+# Evolution
 import time
 t0 = time.time()

inundation/ga/storm_surge/analytical solutions/MacDonald_77541.tsh

@@ -122152 +122152 @@
 0 # <Hole #> <x> <y> ...Mesh Holes...
 0 # <Region #> <x> <y> <tag>...Mesh Regions...
+0 # <Region #> [Max Area] ...Mesh Regions...

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

@@ -17 +17 @@
      Dirichlet_boundary, Transmissive_boundary
 from pmesh2domain import pmesh_to_domain_instance
-from util import Polygon_function
+from util import Polygon_function,read_polygon
+
 
 ######################
@@ -33 +34 @@
 domain.store = True
 domain.format = 'sww'
-domain.filename = 'sydney_ubd'
+domain.filename = 'sydney_ubd_again'
 
 #Reduction operation for get_vertex_values