Changeset 6322

Timestamp:

Feb 12, 2009, 3:40:04 PM (16 years ago)

Author:

ole

Message:

Added file length, changed resolution in time and space.

Location:

anuga_work/production/patong

Files:

• file_length.py (added)
• project.py (modified) (4 diffs)
• run_patong.py (modified) (7 diffs)

anuga_work/production/patong/project.py

@@ -48 +48 @@
 finaltime=15000         # final time for simulation 15000
 
-setup='trial'  # Final can be replaced with trial or basic.
+setup='final'  # Final can be replaced with trial or basic.
                # Either will result in a coarser mesh that will allow a
                # faster, but less accurate, simulation.
@@ -66 +66 @@
     res_factor=1
     time_thinning=4
-    yieldstep=10
+    yieldstep=5
 
 use_buildings = True
@@ -180 +180 @@
 # Inundation area
 poly_ia = read_polygon(polygons_dir+'inundation_area.csv')
-res_ia = 75*res_factor
+res_ia = 59*res_factor
 poly_saddle = read_polygon(polygons_dir+'saddle_10m.txt') 
 
 # Area of Interest 1 elevation from -10m to 20m
 poly_aoi1 = read_polygon(polygons_dir+'aoi.csv')
-res_aoi1 = 200*res_factor
+res_aoi1 = 150*res_factor
 
 # Area of Significance 1 elevation from -20m to a 200m buffer of the 20m contour
 poly_aos1 = read_polygon(polygons_dir+'aos.csv')
-res_aos1 = 900*res_factor
+res_aos1 = 700*res_factor
 
 
 # Area of Shallow water and coastal land that needs a finer res than 1000000
 poly_sw = read_polygon(polygons_dir+'sw.csv')
-res_sw = 6000*res_factor
+res_sw = 5000*res_factor
 
 
@@ -204 +204 @@
 building_saddle = read_polygon(polygons_dir+'building_saddle.csv')
 building_area_polygons = [building_main, building_main_small, building_main_south, building_saddle]
-bld_res = 25*res_factor
+bld_res = 20*res_factor
 
 

anuga_work/production/patong/run_patong.py

@@ -57 +57 @@
 print 'output_dir', output_dir
 
-copy_code_files(output_dir,__file__, 
-                dirname(project.__file__)+sep+ project.__name__+'.py' )
-
-start_screen_catcher(output_dir, myid, numprocs)
+#copy_code_files(output_dir,__file__, 
+#                dirname(project.__file__)+sep+ project.__name__+'.py' )
+#
+#start_screen_catcher(output_dir, myid, numprocs)
 
 
@@ -83 +83 @@
 # Number of landward_boundary points
 M = file_length(project.landward_dir)
+
 
 # Boundary tags refer to project.landward_boundary
@@ -105 +106 @@
 
 print domain.statistics()
-
 
 #-------------------------------------------------------------------------
@@ -122 +122 @@
 print 'friction'
 domain.set_quantity('friction', project.friction)
+
 print 'elevation'
 domain.set_quantity('elevation',
@@ -182 +183 @@
 
 Br = Reflective_boundary(domain)
-Bd = Dirichlet_boundary([project.tide,0,0])
+Bs = Dirichlet_boundary([project.tide,0,0])
 #Bs = Transmissive_stage_zero_momentum_boundary(domain)
 Bf = Field_boundary(urs_boundary_name+'.sts', 
@@ -188 +189 @@
                     mean_stage= project.tide,
                     time_thinning=project.time_thinning,
-                    default_boundary=Bd,
+                    default_boundary=Bs,
                     use_cache=True,
                     verbose=True,
@@ -194 +195 @@
 
 domain.set_boundary({'back': Br,
-                     'side': Bd,
+                     'side': Bs,
                      'ocean': Bf})