Changeset 6225 for anuga_work/production/patong/run_patong.py

Timestamp:

Jan 21, 2009, 4:57:58 PM (15 years ago)

Author:

ole

Message:

Work on Patong w buildings

File:

• anuga_work/production/patong/run_patong.py (modified) (2 diffs)

anuga_work/production/patong/run_patong.py

@@ -116 +116 @@
 domain.set_quantity('friction', project.friction)
 print 'elevation'
-domain.set_quantity('elevation', 
+domain.set_quantity('elevation',
                     filename=project.combined_dir_name+'.pts',
+                    alpha=project.alpha,                    
                     use_cache=True,
-                    verbose=True,
-                    alpha=project.alpha)
+                    verbose=True)
+
 
 # Add buildings from file
 print 'Reading building polygons'    
 building_polygons, building_heights = csv2building_polygons(project.building_polygon_file)
+                                                            #clipping_polygons=project.building_area_polygons)
 
 print 'Creating %d building polygons' % len(building_polygons)
-def create_polygon_function(building_polygons):
+def create_polygon_function(building_polygons, geo_reference=None):
     L = []
     for i, key in enumerate(building_polygons):
@@ -133 +135 @@
         poly = building_polygons[key]
         elev = building_heights[key]
-        L.append((poly, elev))#
+        L.append((poly, elev))
 
         buildings = Polygon_function(L, default=0.0,
-                                     geo_reference=domain.geo_reference)
+                                     geo_reference=geo_reference)
     return buildings
 
 buildings = cache(create_polygon_function,
                   building_polygons,
+                  {'geo_reference': domain.geo_reference},
                   verbose=True)
 
-print 'Add buildings'
-domain.add_quantity('elevation', buildings)
+print 'Adding buildings'
+domain.add_quantity('elevation',
+                    buildings,
+                    use_cache=False, # FIXME(OLE): This seems to pickup the stage IC??
+                    verbose=True)
+