Changeset 5609

Timestamp:

Aug 5, 2008, 1:47:42 PM (15 years ago)

Author:

kristy

Message:

updated run file with optimal alpha

File:

• anuga_work/production/busselton/run_busselton.py (modified) (5 diffs)

anuga_work/production/busselton/run_busselton.py

@@ -41 +41 @@
 from anuga.fit_interpolate.benchmark_least_squares import mem_usage 
 from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon
+from anuga.geospatial_data.geospatial_data import find_optimal_smoothing_parameter
 
 # Application specific imports
@@ -114 +115 @@
     #barrier()
     
+         covariance_value,alpha = find_optimal_smoothing_parameter (data_file= kwargs['elevation_file'],
+                                alpha_list=[0.001, 0.01, 0.1, 0.15, 0.2, 0.25, 0.3, 0.4, 0.5],
+                                mesh_file = project.meshes_dir_name+'.msh')
+        print 'optimal alpha', covariance_value,alpha       
 
     #-------------------------------------------------------------------------
@@ -181 +186 @@
     print 'domain id', id(domain)
    
-    print'set_boundary'
-
     boundary_urs_out=project.boundaries_dir_name
     
     print 'Available boundary tags', domain.get_boundary_tags()
-    Bf = File_boundary(boundary_urs_out+'.sts', 
-                   domain, time_thinning=1,
+    Bf = File_boundary(boundary_urs_out+'.sts', # Change from file_boundary
+                   domain, mean_stage= project.tide,
+                   time_thinning=1,
                    use_cache=True,
-                   verbose = True,
-                   boundary_polygon=bounding_polygon)
+                   verbose = True)
+                   #boundary_polygon=bounding_polygon)
     
     Br = Reflective_boundary(domain)
     Bd = Dirichlet_boundary([kwargs['tide'],0,0])
-    #Bw = Dirichlet_boundary([kwargs['tide']+10.0,0,0]) # To be deleted for FESA runs
-
+
+    print dir(Bf)
     print 'finished reading boundary file'
 
@@ -215 +219 @@
         domain.write_time()
         domain.write_boundary_statistics(tags = 'ocean')
-
-##        # To be deleted for FESA runs
-##        if allclose(t, 240):
-##            domain.set_boundary({'back': Br, 'side': Bd, 'ocean': Bw})
-##
-##        if allclose(t, 1440):
-##            domain.set_boundary({'back': Br, 'side': Bd, 'ocean': Bd})
-
 
 
@@ -255 +251 @@
     kwargs['output_dir']=project.output_run_time_dir
     kwargs['bathy_file']=project.combined_dir_name+'.txt'
-#    kwargs['bathy_file']=project.combined_small_dir_name + '.pts'
- 
-    kwargs['file_name']=project.home+'detail.csv'
+  kwargs['file_name']=project.home+'detail.csv'
     kwargs['aa_scenario_name']=project.scenario_name
     kwargs['ab_time']=project.time