Changeset 5411 for anuga_work/development/convergence_okushiri_2008/project_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/project_truescale.py (modified) (4 diffs)

anuga_work/development/convergence_okushiri_2008/project_truescale.py

@@ -1 @@
-"""Common filenames for truescale Okushiri Island validation
+"""Common filenames for truescale Okushiri Island convergence study
 Formats are given as ANUGA native netCDF where applicable.
 
@@ -11 +11 @@
 from anuga.utilities.system_tools import get_user_name, get_host_name
 
-codename = 'project_grad.py' # FIXME can be obtained automatically as __file__
-
-home = join(getenv('INUNDATIONHOME'), 'data', 'anuga_validation',
-            'convergence_okushiri_2008')# Location of Data   
+home = join(getenv('INUNDATIONHOME'),'data', 'anuga_validation',
+            'convergence_okushiri_2008') # Location of Data
 user = get_user_name()
 host = get_host_name()
@@ -20 +18 @@
 umask(002)
 
-#time stuff
-time = strftime('%Y%m%d_%H%M%S',gmtime()) #gets time for new dir
-run_time = time+'_run'
-
-#Set anuga directories
-anuga_dir = join(home, 'anuga')
-
-dir_comment='_'+setup+'_'+str(tide)+'_'+\
-             str(user)+'_'+boundary_event+'_'+which_data
-
-mesh_dir = join(anuga_dir, 'meshes')
-mesh_name = join(mesh_dir, 'okushiri_truescale')
-
-polygons_dir = join(anuga_dir, 'polygons') # Created with ArcGIS (csv files)
-tide_dir = join(anuga_dir, 'tide_data')
-
-#output locations
-output_dir = join(anuga_dir, 'outputs')+sep
-output_run_time_dir = output_dir +run_time+dir_comment+sep
-output_run_time_dir_name = output_run_time_dir + scenario_name  #Used by post processing
-
-#gauges
-gauges_dir = join(anuga_dir,'gauges')
-gauge_name = 'gauge_location_onslow.csv'
-gauges_dir_name = gauges_dir + gauge_name
+#-------------------
+# Input file names 
+#-------------------
 
 # Given boundary wave
@@ -55 +31 @@
 bathymetry_filename = 'okushiri_truescale_bathymetry.pts'
 
-# Triangular mesh
-mesh_filename = 'okushiri_truescale.msh'
+
+#------------------------------------
+# Output file names and directories 
+#------------------------------------
 
 # Model output
 output_filename = 'okushiri_truescale.sww'
 
+# Time stuff
+time = strftime('%Y%m%d_%H%M%S',gmtime()) #gets time for new dir
+run_time = time+'_run'
+
+# Run parameters
+finaltime=450
+setup='original'
+
+if setup =='original':
+    print 'original resolution'
+    base_resolution=1
+    yieldstep=1
+if setup =='double':
+    print 'double original resolution'
+    base_resolution=0.5
+    yieldstep=1
+if setup =='half':
+    print 'half original resolution'
+    base_resolution=2
+    yieldstep=1
+if setup =='no polygons':
+    print 'half original resolution'
+    base_resolution=2
+    yieldstep=1
+  
+    
+# Set anuga directories
+anuga_dir = join(home,'anuga')+sep
+
+dir_comment='_'+setup+'_'+str(user)
+
+mesh_dir = join(anuga_dir, 'meshes')+sep
+mesh_name = join(mesh_dir, 'okushiri_truescale')
+
+polygons_dir = join(anuga_dir, 'polygons')+sep # Created with ArcGIS (csv files)
+
+# Output locations
+output_dir = join(anuga_dir, 'outputs')+sep
+output_run_time_dir = output_dir+run_time+dir_comment+sep
+output_run_time_dir_name = output_run_time_dir + output_filename #Used by post processing
+
+# Gauges
+gauges_dir = join(anuga_dir,'gauges')+sep
+gauge_name = 'gauge_location_okushiri.csv'
+gauges_dir_name = gauges_dir+gauge_name
+
+# Vertex coordinates
+vertex_filename = output_run_time_dir+setup+'vertex_coordinates.txt'
+
+#------------------------------
+# Polygon definitions
+#------------------------------
+
+poly_all = read_polygon(polygons_dir+'bounding_polygon.csv')
+res_poly_all = 16000*base_resolution
+
+#print 'Area of bounding polygon', polygon_area(poly_all)/1000000.0
+
+poly_gulleys = read_polygon(polygons_dir+'gulleys_polygon.csv')
+res_gulleys = 3.2*base_resolution
+
+poly_island = read_polygon(polygons_dir+'island_polygon.csv')
+res_island = 32*base_resolution
+
+poly_rhs = read_polygon(polygons_dir+'rhs_polygon.csv')
+res_rhs = 80*base_resolution
+
+interior_regions = [[poly_gulleys,res_gulleys],[poly_island,res_island],
+                    [poly_rhs,res_rhs]]
+                    
+trigs_min = number_mesh_triangles(interior_regions, poly_all, res_poly_all)
+
+print 'min number triangles', trigs_min
 
 
+
+
+
+