Changeset 6321

Timestamp:

Feb 12, 2009, 3:21:39 PM (16 years ago)

Author:

kristy

Message:

Adjusted for larger area, not running yet

Location:

anuga_work/production/patong

Files:

• build_patong2.py (deleted)
• project.py (modified) (4 diffs)
• run_patong.py (modified) (4 diffs)

anuga_work/production/patong/project.py

@@ -92 +92 @@
 
 # gauges - used in get_timeseries.py
-gauge_name = 'patong.csv'
+gauge_name = 'gauges.csv'
 
 # BOUNDING POLYGON - used in build_boundary.py and run_patong.py respectively
@@ -98 +98 @@
 # Check the run_patong.py for boundary_tags
 # thinned ordering file from Hazard Map: format is index,latitude,longitude (with title)
-order_filename = 'thinned_boundary_ordering.csv'
+order_filename = 'urs_order.csv'
 #landward bounding points
-landward = 'landward_bounding_polygon.csv'
+landward = 'landward_bounding_polygon2.csv'
 
 #------------------------------------------------------------------------------
@@ -107 +107 @@
 # Output filename for elevation
 # this is a combination of all the data (utilisied in build_boundary)
-combined_name ='patong_combined_elevation'
+combined_name ='patong_combined_elevation_large'
 combined_smaller_name = 'patong_combined_elevation_smaller'
 
@@ -174 +174 @@
 
 # Initial bounding polygon for data clipping 
-poly_all = read_polygon(polygons_dir+'poly_all.csv')
+poly_all = read_polygon(polygons_dir+'bounding_polygon.csv')
 res_poly_all = 150000*res_factor
 

anuga_work/production/patong/run_patong.py

@@ -33 +33 @@
 from anuga.interface import create_sts_boundary
 from anuga.interface import csv2building_polygons
+from file_length import file_length
 
 from anuga.caching import cache
@@ -80 +81 @@
 # Counting segments
 N = len(urs_bounding_polygon)-1
-
-# boundary tags refer to project.landward 4 points equals 5 segments start at N
-boundary_tags={'back': [N+1,N+2,N+3], 'side': [N,N+4], 'ocean': range(N)}
+# Number of landward_boundary points
+M = file_length(project.landward_dir)
+
+# Boundary tags refer to project.landward_boundary
+# 4 points equals 5 segments start at N
+boundary_tags={'back': range(N+1,N+M),
+               'side': [N,N+M],
+               'ocean': range(N)}
 
 #--------------------------------------------------------------------------
@@ -176 +182 @@
 
 Br = Reflective_boundary(domain)
-#Bd = Dirichlet_boundary([project.tide,0,0])
-Bs = Transmissive_stage_zero_momentum_boundary(domain)
+Bd = Dirichlet_boundary([project.tide,0,0])
+#Bs = Transmissive_stage_zero_momentum_boundary(domain)
 Bf = Field_boundary(urs_boundary_name+'.sts', 
                     domain,
                     mean_stage= project.tide,
                     time_thinning=project.time_thinning,
-                    default_boundary=Bs,
+                    default_boundary=Bd,
                     use_cache=True,
                     verbose=True,
@@ -188 +194 @@
 
 domain.set_boundary({'back': Br,
-                     'side': Bs,
+                     'side': Bd,
                      'ocean': Bf})