Changeset 6334

Timestamp:

Feb 13, 2009, 1:48:31 PM (16 years ago)

Author:

ole

Message:

Played with resolution and timestepping

Location:

anuga_work/production/patong

Files:

• project.py (modified) (1 diff)
• run_patong.py (modified) (1 diff)

anuga_work/production/patong/project.py

@@ -175 +175 @@
 # Initial bounding polygon for data clipping 
 poly_all = read_polygon(polygons_dir+'bounding_polygon.csv')
-res_poly_all = 150000*res_factor
+res_poly_all = 160000*res_factor
 
 
 # Inundation area
 poly_ia = read_polygon(polygons_dir+'inundation_area.csv')
-res_ia = 59*res_factor
+res_ia = 75*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 = 150*res_factor
+res_aoi1 = 200*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 = 700*res_factor
+res_aos1 = 900*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 = 5000*res_factor
+res_sw = 7000*res_factor
 
 

anuga_work/production/patong/run_patong.py

@@ -204 +204 @@
 t0 = time.time()
 
-# Skip over the first 5000 seconds
-for t in domain.evolve(yieldstep=500,
-                       finaltime=5000):
+# Skip over the first 6000 seconds
+for t in domain.evolve(yieldstep=2000,
+                       finaltime=6000):
     print domain.timestepping_statistics()
     print domain.boundary_statistics(tags='ocean')