Changeset 5458

Timestamp:

Jul 2, 2008, 11:51:19 AM (17 years ago)

Author:

jakeman

Message:

updates of scripts used for boxingday validation to include boundary ordering

Location:

anuga_work/development/boxingday08

Files:

• project.py (modified) (3 diffs)
• run_boxingday_most.py (modified) (3 diffs)
• run_boxingday_polyline.py (modified) (4 diffs)

anuga_work/development/boxingday08/project.py ¶

@@ -77 +77 @@
 #read in polyline boundary
 from anuga.shallow_water.data_manager import urs2sts,create_sts_boundary
-base_name='polyline'
+urs_name='polyline'
+#base_name=urs_name
+tide=0.35
+base_name='tide_polyline'
+#base_name='most_polyline'
 import os
 if os.path.exists(base_name+'.sts'):
-        print 'sts boundary file already created.'
+    print 'sts boundary file already created.'
 else:
-        urs2sts(base_name,mean_stage=0.0,verbose=False)
+    print 'creatin sts file'
+    urs2sts(urs_name,base_name,mean_stage=tide,verbose=False)
 
 order_filename='sts_order.txt'
@@ -118 +123 @@
 mesh_elevname = 'phuket_mesh_elev'
 boundary_most_in = dir+sep+'out'
-boundary_most_out = dir+sep+'most_boundary_condition'
+#boundary_most_out = dir+sep+'most_boundary_condition'
+boundary_most_out = dir+sep+'tide_most_boundary_condition'
 
 dir='mesh_polygons'
@@ -140 +146 @@
     from pylab import plot,show,hold
     bounding_polygon=ensure_numeric(bounding_polygon)#for plotting only
+    #boundary=ensure_numeric(boundary)
+    #plot(boundary[:,0],boundary[:,1],'o')
     plot(bounding_polygon[:,0],bounding_polygon[:,1],'o')
     plot(bounding_polygon[:,0],bounding_polygon[:,1])

anuga_work/development/boxingday08/run_boxingday_most.py ¶

@@ -143 +143 @@
 #------------------------------------------------------------------------------
      
-tide = 0.0
+tide = 0.35
 domain.set_quantity('stage', tide)
 domain.set_quantity('friction', 0.01) 
@@ -187 +187 @@
            'origin': domain.geo_reference.get_origin(),
            'mean_stage': tide,
-           'zscale': 1,                 #Enhance tsunami
+           'zscale': 1.0,                 #Enhance tsunami
            'fail_on_NaN': False,
            'inverted_bathymetry': True},
@@ -195 +195 @@
 
 print 'Available boundary tags', domain.get_boundary_tags()
-#Bf = File_boundary(project.boundary_most_out+'.sww', 
-#                    domain, time_thinning=1, use_cache=True,verbose = True)
+Bf = File_boundary(project.boundary_most_out+'.sww', 
+                    domain, time_thinning=1, use_cache=True,verbose = True)
 Br = Reflective_boundary(domain)
 Bd = Dirichlet_boundary([tide,0.0,0.0])
 
-domain.set_boundary({'ocean': Bd,
+domain.set_boundary({'ocean': Bf,
                      'otherocean': Bd,
                      'land': Br,

anuga_work/development/boxingday08/run_boxingday_polyline.py ¶

@@ -46 +46 @@
 #------------------------------------------------------------------------------
 
-#scenario = 'poor_simulation'
-scenario = 'good_simulation'
+scenario = 'poor_simulation'
+#scenario = 'good_simulation'
 
 if os.access(scenario, os.F_OK) == 0:
@@ -92 +92 @@
 
 from Numeric import arange,allclose
+#This needs to be done by the user. Not easy at the moment
 boundary_tags={'ocean': arange(0,41).tolist(), 'otherocean': [41,44], 'land': [42,43]}
 
@@ -145 +146 @@
 #------------------------------------------------------------------------------
      
-tide = 0.0
+tide = 0.35
 domain.set_quantity('stage', tide)
 domain.set_quantity('friction', 0.01) 
@@ -165 +166 @@
 # Setup boundary conditions
 #------------------------------------------------------------------------------
-boundary_urs_in='data/boxing'
 boundary_urs_out=project.base_name