Changeset 4063 for anuga_work/production/newcastle_2006/run_newcastle_slide.py

Timestamp:

Dec 6, 2006, 5:35:00 PM (17 years ago)

Author:

sexton

Message:

updates for nsw slide modelling and cairns demo

File:

• anuga_work/production/newcastle_2006/run_newcastle_slide.py (modified) (9 diffs)

anuga_work/production/newcastle_2006/run_newcastle_slide.py

@@ -2 +2 @@
 
 Source data such as elevation and boundary data is assumed to be available in
-directories specified by project.py
-The output sww file is stored in project.outputtimedir 
+directories specified by project_slide.py
+The output sww file is stored in project_slide.outputtimedir 
 
-The scenario is defined by a triangular mesh created from project.polygon,
+The scenario is defined by a triangular mesh created from project_slide.polygon,
 the elevation data and a tsunami wave generated by s submarine mass failure.
 
@@ -38 +38 @@
 
 # creates copy of code in output dir 
-copy_code_files(project.outputtimedir,__file__,dirname(project.__file__)+sep+ project.__name__+'.py' )
+copy_code_files(project_slide.outputtimedir,__file__,dirname(project_slide.__file__)+sep+ project_slide.__name__+'.py' )
 myid = 0
 numprocs = 1
-start_screen_catcher(project.outputtimedir, myid, numprocs)
+start_screen_catcher(project_slide.outputtimedir, myid, numprocs)
 
-print 'USER:    ', project.user
+print 'USER:    ', project_slide.user
 
 #-------------------------------------------------------------------------------
@@ -79 +79 @@
 print 'add'
 G = G11.clip(Geospatial_data(project_slide.poly_surveyclip)) +\
-    G12.clip(Geospatial_data(project_slide.polyAll))
-    (G4.clip(Geospatial_data(project_slide.polyAll))).clip_outside(Geospatial_data(project_slide.poly_surveyclip))
+    G12.clip(Geospatial_data(project_slide.polyAll)) +\
+    (G4.clip(Geospatial_data(project_slide.polyAll)).clip_outside(Geospatial_data(project_slide.poly_surveyclip)))
 print 'export points'
 G.export_points_file(project_slide.combined_dem_name + '.pts')
@@ -88 +88 @@
 #----------------------------------------------------------------------------
 # Create the triangular mesh based on overall clipping polygon with a tagged
-# boundary and interior regions defined in project.py along with
+# boundary and interior regions defined in project_slide.py along with
 # resolutions (maximal area of per triangle) for each polygon
 #-------------------------------------------------------------------------------
@@ -94 +94 @@
 from anuga.pmesh.mesh_interface import create_mesh_from_regions
 remainder_res = 500000
-local_res = 25000
-newcastle_res = 5000
-coast_res = 500
-interior_regions = [[project.poly_newcastle1, local_res],
-                    [project.poly_newcastle2, newcastle_res],
-                    [project.poly_newcastle3, coast_res]]
+local_res = 50000
+newcastle_res = 1000
+interior_regions = [[project_slide.poly_local, local_res],
+                    [project_slide.poly_newcastle, newcastle_res]]
 
 from caching import cache
 _ = cache(create_mesh_from_regions,
-          project.polyAll,
+          project_slide.polyAll,
            {'boundary_tags': {'e0': [0], 'e1': [1], 'e2': [2],
-                              'e3': [3], 'e4':[4], 'e5': [5],
-                              'e6': [6]},
+                              'e3': [3], 'e4':[4]},
            'maximum_triangle_area': remainder_res,
            'filename': meshname,
@@ -122 +119 @@
 print domain.statistics()
 
-domain.set_name(project.basename)
-domain.set_datadir(project.outputtimedir)
+domain.set_name(project_slide.basename)
+domain.set_datadir(project_slide.outputtimedir)
 domain.set_quantities_to_be_stored(['stage', 'xmomentum', 'ymomentum'])
 domain.set_minimum_storable_height(0.01)
@@ -135 +132 @@
 domain.set_quantity('friction', 0.0) 
 domain.set_quantity('elevation', 
-                    filename = project.combined_dem_name + '.pts',
+                    filename = project_slide.combined_dem_name + '.pts',
                     use_cache = True,
                     verbose = True,
@@ -146 +143 @@
 from smf import slide_tsunami
 
-tsunami_source = slide_tsunami(length=30000.0,
-                               depth=400.0,
-                               slope=6.0,
-                               thickness=176.0, 
-                               radius=3330,
-                               dphi=0.23,
-                               x0=project.slump_origin[0], 
-                               y0=project.slump_origin[1], 
-                               alpha=0.0, 
+tsunami_source = slide_tsunami(length=project_slide.bulli_length,
+                               width=project_slide.bulli_width,
+                               depth=project_slide.bulli_depth,
+                               slope=project_slide.bulli_slope,
+                               thickness=project_slide.bulli_thickness, 
+                               x0=project_slide.slide_origin_c[0], 
+                               y0=project_slide.slide_origin_c[1], 
+                               alpha=project_slide.bulli_alpha, 
                                domain=domain)
 
@@ -165 +161 @@
 Bd = Dirichlet_boundary([tide,0,0])
 
-domain.set_boundary( {'e0': Bd,  'e1': Bd, 'e2': Bd, 'e3': Bd, 'e4': Bd,
-                      'e5': Bd,  'e6': Bd} )
+domain.set_boundary( {'e0': Bd,  'e1': Bd, 'e2': Bd, 'e3': Bd, 'e4': Bd} )