Changeset 5992

Timestamp:

Nov 25, 2008, 9:51:52 AM (16 years ago)

Author:

Leharne

Message:

Phase 2 hazard map comparison plot

File:

• anuga_work/production/perth/ph2_compare.py (modified) (9 diffs)

anuga_work/production/perth/ph2_compare.py

@@ -7 +7 @@
 
 import project
-from pylab import plot, xlabel, ylabel, savefig, ion, close, axis, title, legend, grid
+from pylab import plot, xlabel, ylabel, savefig, ion, close, axis, title, legend, grid, figure
 from os import sep
 
@@ -15 +15 @@
 
 def get_max_boundary_data(filename):
+    from anuga.utilities.numerical_tools import mean
     fid = open(filename)
     lines = fid.readlines()
@@ -22 +23 @@
         fields = line.split(',')
         stage.append(float(fields[3]))
-    return max(stage), min(stage)
+    return mean(stage)
 
 def get_stage_data(filename,depth,no_models):
@@ -60 +61 @@
 
 event_number = 27283
-boundary_max, boundary_min = get_max_boundary_data(project.boundaries_dir+str(event_number)+sep+'max_sts_stage.csv')
+boundary_mean = get_max_boundary_data(project.boundaries_dir+str(event_number)+sep+'max_sts_stage.csv')
 
 ###################################################
@@ -70 +71 @@
 depth = [5.0,10.,20.,50.]
 no_models = 3
-filename = directory + sep + 'perth_ph2_compare_v2.csv'
+filename = directory + sep + 'perth_ph2_compare_v5.csv'
 mean_stages, max_stages = get_stage_data(filename, depth, no_models)
 
@@ -81 +82 @@
 d1 = 100.
 d2 = arange(d1,1,-0.1)
-h1 = boundary_max
+h1 = boundary_mean
 green = []
 for d in d2:
@@ -89 +90 @@
 ion()
 figure(1)
+plot(depth,mean_stages[:,2],'>g',d2,green,'-g')
+xlabel('depth (m)') 
+ylabel('stage (m)')
+title('ANUGA outputs (average stage) versus Green\'s approximation \n \
+for event 27283 at Perth')
+legend(['original','Green\'s law'])
+#axis([5,105,min(min(stages))*0.9,max(max(stages))*1.1])
+grid(True)
+figname = 'ph2compare_perth_mean_ORIG_' + str(event_number) + '_mean'
+savefig(figname)
+
+figure(2)
+plot(depth,mean_stages[:,1],'+r',depth,mean_stages[:,2],'>g',d2,green,'-g')
+xlabel('depth (m)') 
+ylabel('stage (m)')
+title('ANUGA outputs (average stage) versus Green\'s approximation \n \
+for event 27283 at Perth')
+legend(['250m poly','original','Green\'s law'])
+#axis([5,105,min(min(stages))*0.9,max(max(stages))*1.1])
+grid(True)
+figname = 'ph2compare_perth_mean_250AP_' + str(event_number) + '_mean'
+savefig(figname)
+
+figure(3)
 plot(depth,mean_stages[:,0],'ob',depth,mean_stages[:,1],'+r',depth,mean_stages[:,2],'>g',d2,green,'-g')
 xlabel('depth (m)') 
@@ -97 +122 @@
 #axis([5,105,min(min(stages))*0.9,max(max(stages))*1.1])
 grid(True)
-figname = 'ph2compare_perth_mean_' + str(event_number)
+figname = 'ph2compare_perth_mean_ALL_' + str(event_number) + '_mean'
 savefig(figname)
 
-figure(2)
+figure(4)
 plot(depth,max_stages[:,0],'ob',depth,max_stages[:,1],'+r',depth,max_stages[:,2],'>g',d2,green,'-g')
 xlabel('depth (m)') 
@@ -109 +134 @@
 #axis([5,105,min(min(stages))*0.9,max(max(stages))*1.1])
 grid(True)
-figname = 'ph2compare_perth_max_' + str(event_number)
+figname = 'ph2compare_perth_max_ALL_' + str(event_number) + '_mean'
+savefig(figname)
 close('all')