Changeset 6095

Timestamp:

Dec 19, 2008, 3:26:16 PM (16 years ago)

Author:

jgriffin

Message:

 

File:

• anuga_work/production/carnarvon/ph2_compare_table.py (modified) (4 diffs)

anuga_work/production/carnarvon/ph2_compare_table.py

@@ -65 +65 @@
     rel_error = zeros(len(mean_stages[1,:]),Float)
     for i in range(len(mean_stages[1,:])):
-        abs_error[i] = (mean_stages[1,i]-mean_stages[1,2])
+        abs_error[i] = (mean_stages[1,i]-mean_stages[1,no_models-1])
         rel_error[i]= (abs_error[i]/mean_stages[1,no_models-1])*100
     return abs_error,rel_error
@@ -73 +73 @@
 ###################################################
 
-event_number = 27283
+event_number = 27255
 boundary_mean = get_max_boundary_data(project.boundaries_dir+str(event_number)+sep+'max_sts_stage.csv')
 
@@ -83 +83 @@
 directory = project.home+project.state+sep+project.scenario+sep+'map_work'
 depth = [5.0,10.,20.,50.]
-no_models = 3
+no_models = 2
 #plot?
 plot_flag = False
 #Table??
 table_flag = True
-filename = directory + sep + 'Carnarvon_Ph2_compare.csv'
+filename = directory + sep + 'Carnarvon_Ph2_compare_'+str(event_number)+'.csv'
+
 mean_stages, max_stages = get_stage_data(filename, depth, no_models)
-
+print mean_stages
 ###################################################
 # Compare with Green's function and plot
@@ -170 +171 @@
 if table_flag==True:
     import csv
-    table=[mean_stages[:,0],mean_stages[:,1],mean_stages[:,2]]
     csvfile = directory+sep+'ph2compare_Carnarvon_max_ALL_' + str(event_number) + '_mean.csv'
     fid1 = open(csvfile,'w')
     writer = csv.writer(fid1,delimiter=',',)
-    writer.writerow(('depth (m)','250m no poly','250m poly','original'))
-    for i in range(len(depth)):
-        writer.writerow((depth[i],mean_stages[i,0],mean_stages[i,1],mean_stages[i,2]))
-    writer.writerow(('100',boundary_mean,boundary_mean,boundary_mean))
-    writer.writerow(('M-value',M_value[0],M_value[1],M_value[2]))
-    writer.writerow(('Absolute Error', abs_error[0],abs_error[1],abs_error[2]))
-    writer.writerow(('Relative Error %', rel_error[0],rel_error[1],rel_error[2]))
-
+
+    if no_models==1:
+        writer.writerow(('depth (m)','250m no poly'))
+        for i in range(len(depth)):
+            writer.writerow((depth[i],mean_stages[i,0]))
+        writer.writerow(('100',boundary_mean,))
+        writer.writerow(('M-value',M_value[0]))
+        #writer.writerow(('Absolute Error', abs_error[0],abs_error[1]))
+        #writer.writerow(('Relative Error %', rel_error[0],rel_error[1]))
+
+    if no_models==2:
+        writer.writerow(('depth (m)','250m no poly','original'))
+        for i in range(len(depth)):
+            writer.writerow((depth[i],mean_stages[i,0],mean_stages[i,1]))
+        writer.writerow(('100',boundary_mean,boundary_mean))
+        writer.writerow(('M-value',M_value[0],M_value[1]))
+        writer.writerow(('Absolute Error', abs_error[0],abs_error[1]))
+        writer.writerow(('Relative Error %', rel_error[0],rel_error[1]))
+        
+    if no_models==3:
+        writer.writerow(('depth (m)','250m no poly','250m poly','original'))
+        for i in range(len(depth)):
+            writer.writerow((depth[i],mean_stages[i,0],mean_stages[i,1],mean_stages[i,2]))
+        writer.writerow(('100',boundary_mean,boundary_mean,boundary_mean))
+        writer.writerow(('M-value',M_value[0],M_value[1],M_value[2]))
+        writer.writerow(('Absolute Error', abs_error[0],abs_error[1],abs_error[2]))
+        writer.writerow(('Relative Error %', rel_error[0],rel_error[1],rel_error[2]))
 
     fid1.close()