Changeset 5194 for anuga_work/production/tonga/tongatapu.py

Timestamp:

Apr 9, 2008, 9:09:18 AM (17 years ago)

Author:

herve

Message:

run scripts for tonga

File:

• anuga_work/production/tonga/tongatapu.py (modified) (10 diffs)

anuga_work/production/tonga/tongatapu.py

@@ -24 +24 @@
 
 # Related major packages
+from anuga.shallow_water import Time_boundary
 from anuga.shallow_water import Domain
 from anuga.shallow_water import Dirichlet_boundary
@@ -31 +32 @@
 from Numeric import allclose
 from anuga.shallow_water.data_manager import export_grid
-
+from anuga.shallow_water.data_manager import convert_dem_from_ascii2netcdf, dem2pts
 from anuga.pmesh.mesh_interface import create_mesh_from_regions
 from anuga.shallow_water.data_manager import start_screen_catcher, copy_code_files,store_parameters
@@ -81 +82 @@
     # resolutions (maximal area of per triangle) for each polygon
     #--------------------------------------------------------------------------
-
-    #IMPORTANT don't cache create_mesh_from_region and Domain(mesh....) as it
-    # causes problems with the ability to cache set quantity which takes alot of times
     
     if myid == 0:
@@ -98 +96 @@
     barrier()
 
+    scenario='fixed_wave'
+    
     #-------------------------------------------------------------------------
     # Setup computational domain
@@ -133 +133 @@
         domain.set_quantity('elevation', 
                             filename = kwargs['bathy_file'],
-                            use_cache = True,
+                            use_cache = False,
                             verbose = True,
                             alpha = kwargs['alpha'])
@@ -170 +170 @@
         tsunami_source = slide_tsunami(length=35000.0,
                                        depth=project.slide_depth,
-                                       slope=6.0,
-                                       thickness=500.0, 
+                                       slope=10.0,
+                                       thickness=800.0, 
                                        x0=project.slide_origin[0], 
                                        y0=project.slide_origin[1], 
@@ -192 +192 @@
     # FIXME (Ole): Change this to Transmissive Momentum as
     # suggested by Rajaraman (ticket:208)
+        from math import sin, pi
         Bw = Time_boundary(domain = domain,
-                           f=lambda t: [(60<t<3660)*100, 0, 0])
-        domain.set_boundary({'ocean_east': Bw,
-                             'ocean_southeast': Bd,
-                             'land': Bd,
-                             'ocean_north': Bd,
-                             'ocean_west':})
+                           f=lambda t: [5.0*sin(t*pi*1/15), 0, 0])
+        domain.set_boundary({'ocean_east':Bw,
+                             'ocean_north':Bd,
+                             'ocean_west':Bd,
+                             'land':Bd,
+                             'ocean_southeast':Bd})
 
     # boundary conditions for slide scenario
     if scenario == 'slide':
-        domain.set_boundary({'ocean_east': Bd,
-                             'ocean_southeast': Bd,
-                             'land': Bd,
-                             'ocean_north': Bd,
-                             'ocean_west':})
+        domain.set_boundary({'ocean_east':Bd,
+                             'ocean_north':Bd,
+                             'ocean_west':Bd,
+                             'land':Bd,
+                             'ocean_southeast':Bd})
 
     #----------------------------------------------------------------------------
@@ -213 +214 @@
     t0 = time.time()
 
-    for t in domain.evolve(yieldstep = 240, finaltime = kwargs['finaltime']): 
+    for t in domain.evolve(yieldstep = 1, finaltime = kwargs['finaltime']): 
         domain.write_time()
         domain.write_boundary_statistics(tags = 'ocean_east')     
@@ -240 +241 @@
     print 'memory usage after del domain',mem_usage()
     
-    swwfile = kwargs['output_dir']+kwargs['scenario_name']
+    swwfile = kwargs['output_dir']+kwargs['aa_scenario_name']
     print'swwfile',swwfile
     
-    export_grid(swwfile, extra_name_out = 'town',
+    export_grid(swwfile, extra_name_out = 'tonga',
             quantities = ['speed','depth','elevation','stage'], # '(xmomentum**2 + ymomentum**2)**0.5' defaults to elevation
             #quantities = ['speed','depth'], # '(xmomentum**2 + ymomentum**2)**0.5' defaults to elevation
             timestep = None,
             reduction = max,
-            cellsize = kwargs['export_cellsize'],
+            cellsize = 25,
             NODATA_value = -1E-030,
             easting_min = project.eastingmin,
@@ -271 +272 @@
     kwargs['midtime']=project.midtime
     kwargs['finaltime']=project.finaltime
-    
     kwargs['output_dir']=project.output_run_time_dir
     kwargs['bathy_file']=project.combined_dir_name+'.txt'