Changeset 4856 for anuga_work/production/wollongong_2006/run_kembla.py

Timestamp:

Nov 27, 2007, 4:18:39 PM (17 years ago)

Author:

sexton

Message:

rename plot_polygons_points to original name of plot_polygons

File:

• anuga_work/production/wollongong_2006/run_kembla.py (modified) (3 diffs)

anuga_work/production/wollongong_2006/run_kembla.py

@@ -41 +41 @@
 meshname = project_kembla.meshname+'.msh'
 
-# creates DEM from asc data
-convert_dem_from_ascii2netcdf(dem_name, use_cache=True, verbose=True)
-
-#creates pts file for onshore DEM
-dem2pts(dem_name, use_cache=True, verbose=True)
-
-print 'create onshore'
-G = Geospatial_data(file_name = project_kembla.dem_name + '.pts')
-
-print 'export points'
-G.export_points_file(project_kembla.combined_dem_name + '.pts')
+### creates DEM from asc data
+##convert_dem_from_ascii2netcdf(dem_name, use_cache=True, verbose=True)
+##
+###creates pts file for onshore DEM
+##dem2pts(dem_name, use_cache=True, verbose=True)
+##
+##print 'create onshore'
+##G = Geospatial_data(file_name = project_kembla.dem_name + '.pts')
+##
+##print 'export points'
+##G.export_points_file(project_kembla.combined_dem_name + '.pts')
 #G.export_points_file(project_kembla.combined_dem_name + '.xya')
 
@@ -114 +114 @@
                    f=lambda t: [(60<t<660)*6, 0, 0])
 
-domain.set_boundary( {'top': Bd,  'right': Bw, 'bottom': Bd, 'left': Bd} )
+domain.set_boundary( {'top': Bd,'right': Bw, 'bottom': Bd, 'left': Bd} )
 
 
@@ -125 +125 @@
 from anuga.abstract_2d_finite_volumes.quantity import Quantity
 
-for t in domain.evolve(yieldstep = 30, finaltime = 5000): 
+for t in domain.evolve(yieldstep = 10, finaltime = 1000): 
     domain.write_time()
-    domain.write_boundary_statistics(tags = '')
+    domain.write_boundary_statistics(tags = 'right')
     stagestep = domain.get_quantity('stage')