Changeset 4537 for anuga_work

Timestamp:

Jun 7, 2007, 8:56:23 AM (17 years ago)

Author:

nick

Message:

 

Location:

anuga_work/development/alpha_validation

Files:

• find_alpha.py (modified) (6 diffs)
• pt_hedland_export.txt (added)

anuga_work/development/alpha_validation/find_alpha.py

@@ -7 +7 @@
 from anuga.shallow_water import Reflective_boundary
 from anuga.utilities.numerical_tools import cov
-from Numeric import array, resize,shape,Float,zeros
+from Numeric import array, resize,shape,Float,zeros,take,argsort
+from pylab import plot, ion, hold,savefig,semilogx,plotting
 
 
@@ -17 +18 @@
                          boundary_tags={'back': [2], 'side': [1,3],
                                         'ocean': [0]},
-                         maximum_triangle_area=100000,
+                         maximum_triangle_area=100,
                          filename=mesh_dir_name,
                          use_cache=True,
@@ -29 +30 @@
 G = Geospatial_data(file_name = topo_dir_name)
 print 'start split'
-G_small, G_other = G.split(0.1,True)
+
+G_sample,G_loss= G.split(0.1, True)
+
+G_small, G_other = G_sample.split(0.1,True)
 
 
 
-#alphas = [0.0001, 0.001, 0.01, 0.1, 1.0, 10.0, 100.0,1000.0,10000.0]
-alphas = [0.001,1,100]
+alphas = [0.000001, 0.00001, 0.0001, 0.001, 0.01, 0.1, 1.0, 10.0, 100.0,1000.0,10000.0]
+#alphas = [0.001,1,100]
 domains = {}
 #domains = []
@@ -66 +70 @@
 normal_cov=array(zeros([len(alphas),2]),typecode=Float)
 i=0
-#print 'domains',domains.keys()
-#keys_list=domains.keys()
-#print 'domains',keys_list.sort()
+
 for domain in domains:
 #    print'domain',domain
@@ -82 +84 @@
     sample_cov= cov(elevation_sample)
 
-    ele_cov= cov(elevation_sample,elevation_predicted)
+    ele_cov= cov(elevation_sample-elevation_predicted)
     normal_cov[i,:]= [domain,ele_cov/sample_cov]
 
@@ -91 +93 @@
 print'data',data[100:200]
 print'normal cov',normal_cov
+normal_cov0=normal_cov[:,0]
+normal_cov_new=take(normal_cov,argsort(normal_cov0))
+semilogx(normal_cov_new[:,0],normal_cov_new[:,1])
+savefig("alphas",dpi=300)
     
-
-# from file load x,y,z data
-# use file xy data to find predicted z data from domain
-# determine error and repeat