Changeset 5617

Timestamp:

Aug 6, 2008, 9:49:29 AM (16 years ago)

Author:

kristy

Message:

Inserted new mesh triangle over dredged area

File:

• anuga_work/production/perth/run_perth.py (modified) (4 diffs)

anuga_work/production/perth/run_perth.py

@@ -115 +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       
-
+##        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       
 
     #-------------------------------------------------------------------------
@@ -196 +195 @@
     domain.set_quantities_to_be_stored(['stage', 'xmomentum', 'ymomentum'])
     domain.tight_slope_limiters = 1
-    #domain.set_maximum_allowed_speed(0.1) # Allow a little runoff (0.1 is OK)
     print 'domain id', id(domain)
-
 
     #-------------------------------------------------------------------------
@@ -240 +237 @@
     t0 = time.time()
 
-    using_sts_boundary = True 
-    
     for t in domain.evolve(yieldstep = project.yieldstep, finaltime = kwargs['finaltime']
                        ,skip_initial_step = False): 
         domain.write_time()
         domain.write_boundary_statistics(tags = 'ocean')
-        if using_sts_boundary is True:
-            if Bf.values >= 99:
-                domain.set_boundary({'ocean': Bd})
-                using_sts_boundary = False
-
+            
     x, y = domain.get_maximum_inundation_location()
     q = domain.get_maximum_inundation_elevation()
@@ -282 +273 @@
     kwargs['output_dir']=project.output_run_time_dir
     kwargs['elevation_file']=project.combined_dir_name+'.pts'
-#    kwargs['elevation_file']=project.combined_small_dir_name + '.pts'
     kwargs['boundary_file']=project.boundaries_in_dir_name + '.sww'
     kwargs['file_name']=project.home+'detail.csv'