Changeset 5575

Timestamp:

Jul 25, 2008, 4:17:26 PM (16 years ago)

Author:

kristy

Message:

Addition of Polyline information with sts files

Location:

anuga_work/production/busselton

Files:

• build_boundary.py (added)
• export_results.py (modified) (3 diffs)
• project.py (modified) (6 diffs)
• run_busselton.py (modified) (5 diffs)

anuga_work/production/busselton/export_results.py

@@ -17 +17 @@
 #time_dir = '20080707_164738_run_trial_0.6_exmouth_5120_kvanputt'
 #time_dir = '20080707_165101_run_trial_0.6_exmouth_5140_kvanputt'
-time_dir = '20080707_165619_run_trial_0.6_exmouth_5160_kvanputt'
+#time_dir = '20080707_165619_run_trial_0.6_exmouth_5160_kvanputt'
+#time_dir = '20080708_101451_run_trial_0.6_exmouth_4743_jsexton'
+#time_dir = '20080708_101731_run_trial_0.6_exmouth_4777_jsexton'
+#time_dir = '20080710_095410_run_trial_0.6_exmouth_4695_jsexton'
+time_dir = '20080717_092900_run_final_0_exmouth_newExtent_kvanputt'
 
 
 cellsize = 25
 #cellsize = 150
-#timestep = 0
+timestep = 0
 directory = project.output_dir
 name = directory+sep+'test'+sep+time_dir+sep+project.scenario_name
@@ -35 +39 @@
 print 'output dir:', name
 
-#var = [0,4] 
+var = [0,4] 
 #var = [2,3] # depth and Speed
-var = [2,3,4] # elevation, depth and Speed
+#var = [2,3,4] # elevation, depth and Speed
 
 
@@ -82 +86 @@
         sww2dem(name, basename_out = outname,
                     quantity = quantityname,
-                    #timestep = timestep,
+                    timestep = timestep,
                     cellsize = cellsize,      
                     #easting_min = project_grad.e_min_area,

anuga_work/production/busselton/project.py

@@ -10 +10 @@
 #from anuga.coordinate_transforms.redfearn import degminsec2decimal_degrees, convert_points_from_latlon_to_utm
 from anuga.utilities.system_tools import get_user_name, get_host_name
+from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon
 
 # file and system info
@@ -16 +17 @@
 
 home = getenv('INUNDATIONHOME') + sep +'data'+sep #Sandpit's parent dir   
+muxhome = getenv('MUXHOME')
 user = get_user_name()
 host = get_host_name()
@@ -33 +35 @@
 scenario = 'busselton_tsunami_scenario'
 
-tide = 0
+tide = 0.6
 
 alpha = 0.1 
 friction=0.01
 starttime=0
-finaltime=50000 # changed for rev_number trial
+finaltime=80000 
 export_cellsize=25
 setup='final'
-source='exmouth'
+source='polyline'
 
 
@@ -135 +137 @@
 #boundaries_source = '1'
 
-if source=='exmouth':
-    boundaries_name = 'busselton_3103_30052008' #exmouth gun
-    boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'exmouth'+sep+'1_10000'+sep
-
-if source=='test':
-    boundaries_name = 'other' #exmouth gun
-    boundaries_in_dir = anuga_dir+'boundaries'+sep
-
+##if source=='exmouth':
+##    boundaries_name = 'busselton_3103_30052008' # exmouth gun
+##    boundaries_in_dir = anuga_dir+'boundaries'+sep+'urs'+sep+'exmouth'+sep+'1_10000'+sep
+##
+##if source=='test':
+##    boundaries_name = 'other' #exmouth gun
+##    boundaries_in_dir = anuga_dir+'boundaries'+sep
+##
 
 #boundaries locations
-boundaries_in_dir_name = boundaries_in_dir + boundaries_name
+#boundaries_in_dir_name = boundaries_in_dir + boundaries_name
 boundaries_dir = anuga_dir+'boundaries'+sep
 boundaries_dir_name = boundaries_dir + scenario_name
+boundaries_dir_mux = muxhome
 
 #output locations
 output_dir = anuga_dir+'outputs'+sep
-#output_build_time_dir = output_dir +'test'+sep + build_time + dir_comment + sep
 output_build_time_dir = output_dir +build_time + dir_comment + sep
-#output_run_time_dir = output_dir +run_time+dir_comment+sep
-output_run_time_dir = output_dir +'test'+sep + run_time + dir_comment +sep # take test out
+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
 
@@ -170 +171 @@
 
 ###############################
-# Domain definitions
+# Interior region definitions
 ###############################
 
-from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon
-
-# bounding polygon for study area
+# Initial bounding polygon for data clipping
 poly_all = read_polygon(polygons_dir+'poly_all_extend.csv')
 res_poly_all = 100000*res_factor
-
-
-###############################
-# Interior region definitions
-###############################
 
 #digitized polygons
@@ -209 +203 @@
 
 
-boundary_tags={'back': [3, 4, 5, 6, 7],
-               'side': [2, 8], 'ocean': [0, 1, 9]}
-
 trigs_min = number_mesh_triangles(interior_regions, poly_all, res_poly_all)
+print 'min number triangles', trigs_min
 
 poly_mainland=read_polygon(polygons_dir+'initial_condition.csv')
-
-print 'min number triangles', trigs_min
 
 

anuga_work/production/busselton/run_busselton.py

@@ -39 +39 @@
 from anuga.damage_modelling.inundation_damage import add_depth_and_momentum2csv, inundation_damage
 from anuga.fit_interpolate.benchmark_least_squares import mem_usage 
+from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon
 
 # Application specific imports
@@ -71 +72 @@
     print "Processor Name:",get_processor_name()
 
-    # filenames
-#    meshes_dir_name = project.meshes_dir_name+'.msh'
-
-    # creates copy of code in output dir
-    print 'min triangles', project.trigs_min,
-    print 'Note: This is generally about 20% less than the final amount'
-
+    #-----------------------------------------------------------------------
+    # Domain definitions
+    #-----------------------------------------------------------------------
+
+    # Read in boundary from ordered sts file
+    urs_bounding_polygon=create_sts_boundary(project.scenario_name)
+
+    # Reading the landward defined points, this incorporates the original clipping
+    # polygon minus the 100m contour
+    landward_bounding_polygon = read_polygon(project.polygons_dir+'landward_bounding_polygon.txt')
+
+    # Combine sts polyline with landward points
+    bounding_polygon = urs_bounding_polygon + landward_bounding_polygon
+    
+    # counting segments
+    N = len(urs_bounding_polygon)-1
+    boundary_tags={'back': [N+1,N+2,N+3,N+4, N+5], 'side': [N,N+6],'ocean': range(N)}
+
+    
     #--------------------------------------------------------------------------
     # Create the triangular mesh based on overall clipping polygon with a
@@ -106 +119 @@
     print 'Setup computational domain'
 
-    #domain = cache(Domain, (meshes_dir_name), {'use_cache':True, 'verbose':True}, verbose=True)
-    #above don't work
     domain = Domain(project.meshes_dir_name+'.msh', use_cache=False, verbose=True)
     print 'memory usage before del domain',mem_usage()
@@ -168 +179 @@
     print 'Available boundary tags', domain.get_boundary_tags()
     print 'domain id', id(domain)
-    #print 'Reading Boundary file',project.boundaries_dir_namea + '.sww'
-
+   
     print'set_boundary'
 
+    boundary_urs_out=project.boundaries_dir_name
+    
+    print 'Available boundary tags', domain.get_boundary_tags()
+    Bf = File_boundary(boundary_urs_out+'.sts', 
+                   domain, time_thinning=1,
+                   use_cache=True,
+                   verbose = True,
+                   boundary_polygon=bounding_polygon)
+    
     Br = Reflective_boundary(domain)
     Bd = Dirichlet_boundary([kwargs['tide'],0,0])
-    Bw = Dirichlet_boundary([kwargs['tide']+10.0,0,0]) # To be deleted for FESA runs
-
-    Bf = Field_boundary(kwargs['boundary_file'],
-                    domain, time_thinning=kwargs['time_thinning'], mean_stage=kwargs['tide'], 
-                    use_cache=True, verbose=True)
+    #Bw = Dirichlet_boundary([kwargs['tide']+10.0,0,0]) # To be deleted for FESA runs
+
     print 'finished reading boundary file'
+
     domain.set_boundary({'back': Bd,
-                             'side': Bd,
-                             'ocean': Bd}) #changed from Bf to Bd for large wave
+                         'side': Bd,
+                         'ocean': Bf}) #changed from Bf to Bd for large wave
 
     kwargs['input_start_time']=domain.starttime
@@ -198 +215 @@
         domain.write_boundary_statistics(tags = 'ocean')
 
-        # To be deleted for FESA runs
-        if allclose(t, 240):
-            domain.set_boundary({'back': Br, 'side': Bd, 'ocean': Bw})
-
-        if allclose(t, 1440):
-            domain.set_boundary({'back': Br, 'side': Bd, 'ocean': Bd})
+##        # To be deleted for FESA runs
+##        if allclose(t, 240):
+##            domain.set_boundary({'back': Br, 'side': Bd, 'ocean': Bw})
+##
+##        if allclose(t, 1440):
+##            domain.set_boundary({'back': Br, 'side': Bd, 'ocean': Bd})