Changeset 3615 for anuga_work/production/hobart_2006/run_hobart.py

Timestamp:

Sep 18, 2006, 5:19:25 PM (18 years ago)

Author:

sexton

Message:

updates to Hobart script(data not yet received)

File:

• anuga_work/production/hobart_2006/run_hobart.py (modified) (9 diffs)

anuga_work/production/hobart_2006/run_hobart.py

@@ -6 +6 @@
 
 The scenario is defined by a triangular mesh created from project.polygon,
-the elevation data and a simulated submarine landslide.
+the elevation data and a tsunami wave generated by MOST.
 
 Ole Nielsen and Duncan Gray, GA - 2005 and Nick Bartzis, GA - 2006
 """
-
-
 #-------------------------------------------------------------------------------# Import necessary modules
 #-------------------------------------------------------------------------------
@@ -67 +65 @@
 # filenames
 onshore_dem_name = project.onshore_dem_name
-coast_points = project.coast_dem_name
-offshore_points = project.offshore_dem_name
 meshname = project.meshname+'.msh'
 source_dir = project.boundarydir
@@ -81 +77 @@
 
 #creates pts file for onshore DEM
-dem2pts(onshore_dem_name,
-        easting_min=project.eastingmin,
-        easting_max=project.eastingmax,
-        northing_min=project.northingmin,
-        northing_max= project.northingmax,
-        use_cache=True, 
-        verbose=True)
-
-convert_dem_from_ascii2netcdf(islands_dem_name, use_cache=True, verbose=True)
-
-#creates pts file for islands DEM
-dem2pts(islands_dem_name, use_cache=True, verbose=True)
+dem2pts(onshore_dem_name, use_cache=True, verbose=True)
 
 print'create G1'
-G1 = Geospatial_data(file_name = project.offshore_dem_name + '.xya')
+G1 = Geospatial_data(file_name = project.offshore_dem_name_local + '.xya')
 print'create G2'
-G2 = Geospatial_data(file_name = project.onshore_dem_name + '.pts')
+G2 = Geospatial_data(file_name = project.offshore_dem_name_aho + '.xya')
 print'create G3'
-G3 = Geospatial_data(file_name = project.coast_dem_name + '.xya')
-print'add G1+G2+G3'
-G = G1 + G2 + G3
+G3 = Geospatial_data(file_name = project.onshore_dem_name + '.pts')
+print'create G4'
+G4 = Geospatial_data(file_name = project.coast_dem_name + '.xya')
+print'add G1+G2+G3+G4'
+G = G1 + G2 + G3 + G4
 print'export G'
 G.export_points_file(project.combined_dem_name + '.pts')
@@ -113 +100 @@
 from anuga.pmesh.mesh_interface import create_mesh_from_regions
 
-#new
-region_res = 200000
-coast_res = 25000
-hobart_res = 5000
+# use 75 for onshore components (12.5m DEM)
+island_res = 10000
+hobart_res = 10000
+peninsular_res = 10000
 interior_regions = [[project.poly_hobart, hobart_res],
-                    [project.poly_coast, coast_res],
-                    [project.poly_region, region_res]]
+                    [project.poly_tasman_peninsula, peninsula_res],
+                    [project.poly_bruny, island_res]]
 
 print 'number of interior regions', len(interior_regions)
@@ -126 +113 @@
 _ = cache(create_mesh_from_regions,
           project.polyAll,
-          {'boundary_tags': {'top': [0], 'topleft': [1],
-                             'topleft1': [2], 'bottomleft': [3],
-                             'bottom': [4], 'bottomright': [5],
-                             'topright':[6]},
+          {'boundary_tags': {'bottom': [0], 'right': [1],
+                             'top': [2], 'left': [3]},
            'maximum_triangle_area': 100000,
            'filename': meshname,           
@@ -177 +162 @@
 print 'start ferret2sww'
 from anuga.pyvolution.data_manager import ferret2sww
-
+'''
 south = project.south 
 north = project.north
@@ -209 +194 @@
        dependencies = source_dir + project.boundary_basename + '.sww')
 
-
+'''
 print 'Available boundary tags', domain.get_boundary_tags()
 
-Bf = File_boundary(source_dir + project.boundary_basename + '.sww', 
+#Bf = File_boundary(source_dir + project.boundary_basename + '.sww', 
                     domain, verbose = True)
 Br = Reflective_boundary(domain)
@@ -222 +207 @@
                    f=lambda t: [(60<t<480)*6, 0, 0])
 
-domain.set_boundary( {'top': Bf, 'topleft': Bf,
-                      'topleft1': Bf, 'bottomleft': Bd,
-                      'bottom': Br, 'bottomright': Br, 'topright': Bd} )
+# for MOST BC
+#domain.set_boundary( {'top': Bd, 'left': Bd,
+#                      'bottom': Bf, 'right': Bf} )
+
+# for testing
+domain.set_boundary( {'top': Bd, 'left': Bd,
+                      'bottom': Bd, 'right': Bw} )
 
 #-------------------------------------------------------------------------------                                 
@@ -234 +223 @@
 for t in domain.evolve(yieldstep = 240, finaltime = 7200): 
     domain.write_time()
-    domain.write_boundary_statistics(tags = 'top')     
+    domain.write_boundary_statistics(tags = 'bottom')     
 
 for t in domain.evolve(yieldstep = 120, finaltime = 12600
                        ,skip_initial_step = True): 
     domain.write_time()
-    domain.write_boundary_statistics(tags = 'top')     
-
-for t in domain.evolve(yieldstep = 60, finaltime = 19800
-                       ,skip_initial_step = True): 
-    domain.write_time()
-    domain.write_boundary_statistics(tags = 'top')     
-    
-for t in domain.evolve(yieldstep = 120, finaltime = 25200
-                       ,skip_initial_step = True): 
-    domain.write_time()
-    domain.write_boundary_statistics(tags = 'top')     
-
-for t in domain.evolve(yieldstep = 240, finaltime = 36000
-                       ,skip_initial_step = True): 
-    domain.write_time()
-    domain.write_boundary_statistics(tags = 'top')
+    domain.write_boundary_statistics(tags = 'bottom')     
     
 print 'That took %.2f seconds' %(time.time()-t0)