Changeset 6592

Timestamp:

Mar 24, 2009, 12:13:42 PM (15 years ago)

Author:

ole

Message:

Made changes to Patong so that STS files are stored at mean sea level.
The tidal correction should happen in Field_boundary only.
Also added notes about the danger of not recomputing the STS file.

Location:

anuga_work/production/patong/new_version

Files:

• build_urs_boundary.py (modified) (2 diffs)
• run_model.py (modified) (3 diffs)

anuga_work/production/patong/new_version/build_urs_boundary.py

@@ -140 +140 @@
                 weights=mux_weights,
                 zone=project.zone,
-                mean_stage=project.tide,
+                mean_stage=0.0, # NOTE: This has to be 0.0 if we use mean_stage in Field_boundary
                 verbose=True)
     else:                           # a single mux stem file, assume 1.0 weight
@@ -158 +158 @@
                 ordering_filename=order_filename,
                 weights=mux_weights,
-                mean_stage=project.tide,
+                mean_stage=0.0,  # NOTE: This has to be 0.0 if we use mean_stage in Field_boundary
                 verbose=True)
 

anuga_work/production/patong/new_version/run_model.py

@@ -68 +68 @@
 
 # Create the STS file
+# FIXME (Ole): This is deadly dangerous if buildcode changes (as was the case 24th March 2009)
+# We need to use caching instead!
+
 print 'project.mux_data_folder=%s' % project.mux_data_folder
 if not os.path.exists(project.event_sts + '.sts'):
@@ -140 +143 @@
     #clipping_polygons=project.building_area_polygons)
 
-    print 'Creating %d building polygons' % len(building_polygons)
     def create_polygon_function(building_polygons, geo_reference=None):
         L = []
@@ -153 +155 @@
         return buildings
 
+    print 'Creating %d building polygons' % len(building_polygons)
     buildings = cache(create_polygon_function,
                       building_polygons,