Changeset 5411 for anuga_work/development/convergence_okushiri_2008/run_okushiri_truescale.py

Timestamp:

Jun 23, 2008, 1:01:17 PM (16 years ago)

Author:

Leharne

Message:

Updates to okushiri convergence study

File:

• anuga_work/development/convergence_okushiri_2008/run_okushiri_truescale.py (modified) (5 diffs)

anuga_work/development/convergence_okushiri_2008/run_okushiri_truescale.py

@@ -1 @@
-"""Validation of the AnuGA implementation of the shallow water wave equation.
+"""Anuga convergence study using a true-scale version of the Okushiri
+Island tsunami wavetank experiment
 
-This script sets up a modified version of the Okushiri Island benchmark
+This script runs a modified version of the Okushiri Island benchmark
 
 as published at
@@ -12 +13 @@
 
 This version up-scales the original 1:400 scale wave-tank experiment, to
-"true-scale".
+"true-scale" with mesh defined using interior polygons.
 
-The original validation data was downloaded and made available in this directory
-for convenience but the original data is available at
+The original validation data is available at
 http://www.cee.cornell.edu/longwave/index.cfm?page=benchmark&problem=2
 where a detailed description of the problem is also available.
 
-
-Run create_okushiri_truescale.py to process the boundary condition and build a the
-mesh before running this script.
+Run create_okushiri_truescale.py to convert bathymetry and input boundary
+condition into NetCDF format.
 
 """
+#----------------------------------------
+# Import necessary modules
+#----------------------------------------
 
-# Module imports
+# Standard modules
+from os import sep,umask
+from os.path import dirname, basename
+from os import mkdir, access, F_OK
+from shutil import copy
+import time
+import sys
+
+# Related major packages
 from anuga.shallow_water import Domain
 from anuga.shallow_water import Reflective_boundary
 from anuga.shallow_water import Transmissive_Momentum_Set_Stage_boundary
 from anuga.abstract_2d_finite_volumes.util import file_function
-
-import project_truescale
+from anuga.shallow_water.data_manager import copy_code_files, start_screen_catcher
+from anuga.pmesh.mesh_interface import create_mesh_from_regions
 
 
-#-------------------------
+import project_truescale # Definition of filenames and interior polygons
+
+copy_code_files(project_truescale.output_run_time_dir,'run_okushiri_truescale.py',
+                'project_truescale.py' )
+myid = 0
+numprocs = 1
+start_screen_catcher(project_truescale.output_run_time_dir, myid, numprocs)
+
+#--------------------------------------------------------------------------
+# Create the triangular mesh based on overall bounding polygon with a
+# tagged boundary and interior regions defined in project_truescale.py
+# along with resolutions (maximal area of per triangle) for each polygon
+#--------------------------------------------------------------------------
+
+create_mesh_from_regions(project_truescale.poly_all,
+                         boundary_tags={'wall': [0, 1, 3],'wave': [2]},
+                         maximum_triangle_area=project_truescale.res_poly_all,
+                         #interior_regions=project_truescale.interior_regions,
+                         filename=project_truescale.mesh_name+'.msh',
+                         verbose=True)
+
+#------------------------------
 # Create Domain from mesh
-#-------------------------
-domain = Domain(project_truescale.mesh_filename, use_cache=True, verbose=True)
+#------------------------------
+domain = Domain(project_truescale.mesh_name+'.msh', use_cache=True, verbose=True)
 print domain.statistics()
 
@@ -43 +74 @@
 
 # Write vertex coordinates to file
-filename='okushiri_mesh_vertex_coordinates_truescale.txt'
+filename=project.vertex_filename
 fid=open(filename,'w')
 fid.write('x (m), y (m)\n')
@@ -68 +99 @@
 #-------------------------
 domain.set_name(project_truescale.output_filename)  # Name of output sww file
-domain.set_datadir(project_truescale.output_dir)
+domain.set_datadir(project_truescale.output_run_time_dir) # Name of output directory
 domain.set_default_order(2)               # Apply second order scheme 
 domain.set_all_limiters(0.9)              # Max second order scheme (old lim)
-domain.set_minimum_storable_height(0.4)   # Don't store h < 0.4m
+domain.set_minimum_storable_height(0.1)   # Don't store h < 0.1m
 domain.tight_slope_limiters = 1
 domain.beta_h = 0.0
 
-#Timings on AMD64-242 (beta_h=0)
-#  tight_slope_limiters = 0:
-#    3035s - 3110s
-#  tight_slope_limiters = 1:
-#    3000s - 3008s
-#
-# beta_h>0: In the order of 3200s
 
 #-------------------------
@@ -109 +133 @@
 t0 = time.time()
 
-for t in domain.evolve(yieldstep = 1, finaltime = 450):
+for t in domain.evolve(yieldstep=project_truescale.yieldstep, finaltime = 450):
     print domain.timestepping_statistics(track_speeds=False,
                                          triangle_id=triangle_id)