Changeset 7572

Timestamp:

Nov 26, 2009, 2:02:35 PM (15 years ago)

Author:

jgriffin

Message:

Now calculates differences in area and writes to file

File:

• anuga_work/production/new_south_wales/batemans_bay/compare_inundation_areas.py (modified) (3 diffs)

anuga_work/production/new_south_wales/batemans_bay/compare_inundation_areas.py

@@ -13 +13 @@
 from os.path import join
 from collections import defaultdict
+import csv
+
+
+
 # Create the Geoprocessor object
 gp = arcgisscripting.create(9.3)
@@ -30 +34 @@
 boundary_path = "N:\\georisk_models\\inundation\\data\\new_south_wales\\batemans_bay_tsunami_scenario_2009\\anuga\\boundaries\\"
 boundary_file_name = 'wave_energy_summary.csv'
+out_file_name = 'area_comparisons.csv'
+#figure = 'area_comparisons.png'
+
 boundary_file = join(boundary_path,boundary_file_name)
+out_file = join(boundary_path, out_file_name)
+#figure_path = join(boundary_path, 'figures',figure)
+plot = True
 
 event_dict = {'Event1_MSL':58346,
-              'Event1_HAT':58346,
-##              'Event2_MSL':51204,
+##              'Event1_HAT':58346,
+              'Event2_MSL':51204,
 ##              'Event2_HAT':51204,
-              'Event3_MSL':58284}#,
+              'Event3_MSL':58284,#,
 ##              'Event3_HAT':58284,
-##              'Puysegur_200yr':58129,
-##              'Puysegur_500yr':58115,
-##              'Puysegur_1000yr':58226,
-##              'Puysegur_5000yr':58286,
-##              'Puysegur_15000yr':58272,
-##              'New_Hebrides_200yr':51077,
-##              'New_Hebrides_500yr':51378,
-##              'New_Hebrides_1000yr':51347,
-##              'New_Hebrides_2000yr':51292,
-##              'New_Hebrides_5000yr':51424,}
-
-# Dictionary for storing area comparison
+              'Puysegur_200yr':58129,
+              'Puysegur_500yr':58115,
+              'Puysegur_1000yr':58226,
+              'Puysegur_5000yr':58286,
+              'Puysegur_15000yr':58272,
+              'New_Hebrides_200yr':51077,
+              'New_Hebrides_500yr':51378,
+              'New_Hebrides_1000yr':51347,
+              'New_Hebrides_2000yr':51292,
+              'New_Hebrides_5000yr':51424,}
+
+# Dictionaries for storing area comparison
 area_comp = defaultdict(list)
+area_percent = defaultdict(list)
 
 boundaries = file(boundary_file,'r')
@@ -236 +247 @@
 
     for key in anuga_area_dict:
+        if not key in analytical_area_dict:
+            analytical_area_dict[key] = 0.
         diff = anuga_area_dict[key]-analytical_area_dict[key]
+        percent_diff = 100*diff/anuga_area_dict[key]
         area_comp[key].append(diff)
-    print 'area_comp', area_comp
-    
+        area_percent[key].append(percent_diff)
+print 'area_comp', area_comp
+print 'area_percent', area_percent
+
+##csv_dict = csv.DictWriter(open(out_file,'w'),[1,2,3,4,5,6,7,8,9])
+##csv_dict.writerow(dict(zip([1,2,3,4,5,6,7,8,9],[1,2,3,4,5,6,7,8,9])))
+##csv_dict.writerow(area_percent)
+
+csv_out = csv.writer(open(out_file,'w'))
+csv_out.writerow([1,2,3,4,5,6,7,8,9])
+for i in range(len(area_percent[1])):
+    value_list =[]
+    for key in area_percent:
+        value_list.append(area_percent[key][i])
+    csv_out.writerow(value_list)
+     
+out_file.close()    
+##for key, value in area_percent:
+##    key_list = []
+##    for val in value:
+##        key_list.append(key)
+##    pylab.scatter(key_list,value)
+        
+
+
+
+
 boundaries.close()
+
+