Changeset 3318

Timestamp:

Jul 12, 2006, 5:40:30 PM (18 years ago)

Author:

duncan

Message:

more damage stuff. still working on it

Location:

inundation/damage

Files:

• inundation_damage.py (modified) (7 diffs)
• test_inundation_damage.py (modified) (2 diffs)

inundation/damage/inundation_damage.py

@@ -6 +6 @@
 
 from math import sqrt
+from Scientific.Functions.Interpolation import InterpolatingFunction
+from Numeric import array
 
-from pyvolution.data_manager import Exposure_csv
-from pyvolution.util import file_function
-from geospatial_data.geospatial_data import ensure_absolute
-from utilities.numerical_tools import INF
+
 try:  
     import kinds  
 except ImportError:  
-    # Hand-built mockup of the things we need from the kinds package,
-    #since it  
+    # Hand-built mockup of the things we need from the kinds package, since it  
     # was recently removed from the standard Numeric distro.  Some users may  
     # not have it by default.  
@@ -27 +25 @@
      
     kinds = _kinds()
+    
+
+from pyvolution.data_manager import Exposure_csv
+from pyvolution.util import file_function
+from geospatial_data.geospatial_data import ensure_absolute
+from utilities.numerical_tools import INF
+from anuga_config import epsilon
 
 def inundation_damage(sww_file, exposure_file_in,
@@ -37 +42 @@
     csv = Exposure_csv(exposure_file_in)
     geospatial = csv.get_location()
-    max_depth, max_momentum = calc_max_depth_and_momentum(sww_file,
+    max_depths, max_momentums = calc_max_depth_and_momentum(sww_file,
                                                           geospatial,
                                                           verbose=verbose,
                                                           use_cache=use_cache)
-    csv.set_column("MAX INUNDATION DEPTH (m)",max_depth)
-    csv.set_column("MOMENTUM (m^2/s) ",max_momentum)
-    csv.save(exposure_file_out)
+    #Do the two damage curves
+    # Ole probably knows if ther is good code for this lying around
+    #Do the prob of struct collapse
+    #Do economics
+    # Save info back to csv file
 
     
@@ -55 +62 @@
     csv = Exposure_csv(exposure_file_in)
     geospatial = csv.get_location()
-    max_depth, max_momentum = calc_max_depth_and_momentum(sww_file,
+    max_depths, max_momentums = calc_max_depth_and_momentum(sww_file,
                                                           geospatial,
                                                           verbose=verbose,
                                                           use_cache=use_cache)
-    csv.set_column("MAX INUNDATION DEPTH (m)",max_depth, overwrite=overwrite)
-    csv.set_column("MOMENTUM (m^2/s) ",max_momentum, overwrite=overwrite)
+    csv.set_column("MAX INUNDATION DEPTH (m)",max_depths, overwrite=overwrite)
+    csv.set_column("MOMENTUM (m^2/s) ",max_momentums, overwrite=overwrite)
     csv.save(exposure_file_out)
     
@@ -68 +75 @@
                                  use_cache = True):
     """
-    Calculate the maximum depth above ground height
+    Calculate the maximum inundation height above ground floor (mihagf) for a list
+    of locations.
+
+    The mihagf value is in the range -ground_floor_height to overflow errors.
     """
     quantities =  ['stage', 'elevation', 'xmomentum', 'ymomentum']
@@ -79 +89 @@
                                  use_cache=use_cache)
     # initialise the max lists
-    max_depth = [kinds.default_float_kind.MIN]*point_count
-    max_momentum = [kinds.default_float_kind.MIN]*point_count
+    max_depths = [ground_floor_height]*point_count
+    max_momentums = [ground_floor_height]*point_count
     for point_i, point in enumerate(points):
         for time in callable_sww.get_time():
@@ -90 +100 @@
             uh = quantities[2]
             vh = quantities[3]
+
+            # This can easily give -ve values,
+            # but the array is initialised to 0.0,
+            # so 0.0 will be the minimum value.
             depth = w - z - ground_floor_height
               
-            if depth > max_depth[point_i]:
-                max_depth[point_i] = depth
+            if depth > max_depths[point_i]:
+                max_depths[point_i] = depth
             if w == INF or z == INF or uh == INF or vh == INF:
                 continue
             momentum = sqrt(uh*uh + vh*vh)
-            if momentum > max_momentum[point_i]:
-                max_momentum[point_i] = momentum
-    return max_depth, max_momentum       
+            if momentum > max_momentums[point_i]:
+                max_momentums[point_i] = momentum
+    return max_depths, max_momentums
+
+class EventDamageModel:
+    """
+    Initially I'm doing this as a class to group all of these functions/methods.
+
+
+    """
+    double_brick_damage_curve = array([[-kinds.default_float_kind.MAX, 0.0],
+                                    [0.0-epsilon, 0.0],
+                                    [0.0,1.6],
+                                    [2.5,70.0],
+                                    [kinds.default_float_kind.MAX,70]])
+     
+    def __init__(self,max_depths, shore_distances, walls):
+        self.max_depths = max_depths
+        self.shore_distances = shore_distances
+        self.walls = walls
+        assert len(self.max_depths) == len(self.shore_distances)
+        self.structure_count = len(self.max_depths)
+        structure_damage_curve = InterpolatingFunction((array([0,1]),),array([0,10]))
+
+    def calc_damage_percentages(self):
+
+        # the data being created
+        percent_struct_damage = [0.0]*self.structure_count
+        percent_contents_damage = [0.0]*self.structure_count
+        #for i_building,
+
+        for i,max_height,shore_distance,wall in map(None,
+                                                    range(self.structure_count),
+                                                    self.max_depths,
+                                                    self.shore_distances,
+                                                    self.walls):
+            print "i",i
+            print "max_height",max_height
+            print "shore_distance",shore_distance
+           
         
+#############################################################################
+if __name__ == "__main__":
+    from Scientific.Functions.Interpolation import InterpolatingFunction
+    from Numeric import array, ravel
+    a = array([[0,0],[1,10]])
+    c = array([0,1])
+
+    axis = ravel(a[:,0:1])
+    values = ravel(a[:,1:])
+    
+    #axis = array([0,1])
+    #values = array([0,10])
+    print "axis",axis
+    print "values",values 
+    i = InterpolatingFunction((axis,), values)
+    print "value",i(0.5)
+
+    
+    i = InterpolatingFunction((array([0,1]),),array([0,10]))
+    print "value",i(0.5)
+
+        

inundation/damage/test_inundation_damage.py

@@ -128 +128 @@
         fd = open(self.csv_file,'wb')
         writer = csv.writer(fd)
-        writer.writerow(['LONGITUDE','LATITUDE','STR_VALUE','C_VALUE','ROOF_TYPE','WALLS'])
-        writer.writerow(['151.5','-34','199770','130000','Metal','Timber'])
-        writer.writerow(['151','-34.5','150000','76000','Metal','Double Brick'])
+        writer.writerow(['LONGITUDE','LATITUDE','STR_VALUE','C_VALUE','ROOF_TYPE','WALLS', 'SHORE_DIST'])
+        writer.writerow(['151.5','-34','199770','130000','Metal','Timber',20])
+        writer.writerow(['151','-34.5','150000','76000','Metal','Double Brick',200])
 
         
@@ -153 +153 @@
                                    out_csv, verbose=False)
         
+    def test_calc_damage_percentages(self):
+        max_depths = [-0.3, 0.0, 1.0]
+        shore_distances = [100, 200, 300]
+        walls = ['Double Brick','Timber','Brick Veneer']
+
+        edm = EventDamageModel(max_depths, shore_distances, walls)
+        edm.calc_damage_percentages()
+        
 #-------------------------------------------------------------
 if __name__ == "__main__":