Changeset 1845

Timestamp:

Sep 21, 2005, 2:42:41 PM (19 years ago)

Author:

ole

Message:

Comments

File:

• production/karratha_2005/run_karratha.py (modified) (8 diffs)

production/karratha_2005/run_karratha.py

@@ -1 @@
-"""Convert from Arcview ASCII DEMs via native netcdf dem format
-to native pts netcdf format for use with least_squares fits
+"""Script for running a tsunami inundation scenario for Karratha, WA, Australia.
+
+Source data such as elevation and boundary data is assumed to be available in
+directories specified by project.py
+The output sww file is stored in project.outputdir 
+
+The scenario is defined by a triangular mesh created from project.polygon,
+the elevation data and boundary data obtained from a tsunami simulation done with MOST.
+
+Ole Nielsen, GA - 2005
 """
 
@@ -7 +15 @@
 
 
-from pyvolution.shallow_water import Domain, Reflective_boundary, File_boundary, Dirichlet_boundary, Time_boundary, Transmissive_boundary
+from pyvolution.shallow_water import Domain, Reflective_boundary, File_boundary,\
+     Dirichlet_boundary, Time_boundary, Transmissive_boundary
 from pyvolution.data_manager import convert_dem_from_ascii2netcdf,\
      dem2pts, ferret2sww
@@ -31 +40 @@
 source_dir = project.boundarydir
 
-
-#Data
 from pyvolution.data_manager import ferret2sww
-
 
 south = project.south 
@@ -42 +48 @@
 
 
-#Mesh
+
+#Create Triangular Mesh
 from create_mesh import create_mesh
-
 m = cache(create_mesh,
           project.polygon,
@@ -90 +96 @@
 
 
-#Setup IC
+
+#Setup Initial Conditions
 tide = 0
 domain.set_quantity('friction', 0)
@@ -99 +106 @@
                     verbose = True)
 
-#Setup BC
 
+
+#Setup Boundary Conditions
 print domain.get_boundary_tags()
 
 Bf = File_boundary(project.boundary_basename + '.sww', domain, verbose = True)
 domain.starttime = 3000  #Obtained from MOST
-
 
 Br = Reflective_boundary(domain)
@@ -113 +120 @@
                    f=lambda t: [(60<t<660)*4, 0, 0])
 
-#domain.set_boundary( {'wall': Bd, 'ocean': Bf} )
 domain.set_boundary( {'back': Br,'side': Bd, 'ocean': Bf} )
 
@@ -120 +126 @@
 
 #Run
-
 for t in domain.evolve(yieldstep = 60, finaltime = 40000):
     domain.write_time()