Changeset 5408 for anuga_work/production/busselton/run_busselton.py

Timestamp:

Jun 17, 2008, 9:03:30 AM (16 years ago)

Author:

kristy

Message:

 

File:

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

anuga_work/production/busselton/run_busselton.py

@@ -123 +123 @@
         from polygon import Polygon_function
         #following sets the stage/water to be offcoast only
-        IC = Polygon_function( [(project.poly_mainland, -1.0)], default = kwargs['tide'],
+        IC = Polygon_function( [(project.poly_mainland, 0)], default = kwargs['tide'],
                                  geo_reference = domain.geo_reference)
         domain.set_quantity('stage', IC)
@@ -173 +173 @@
     Br = Reflective_boundary(domain)
     Bd = Dirichlet_boundary([kwargs['tide'],0,0])
-    Bo = Dirichlet_boundary([kwargs['tide']+10.0,0,0])
-
-    if project.source != 'test':
-        print 'start reading boundary file'
-        
-        Bf = Field_boundary(kwargs['boundary_file'],
+   
+
+    Bf = Field_boundary(kwargs['boundary_file'],
                     domain, time_thinning=kwargs['time_thinning'], mean_stage=kwargs['tide'], 
                     use_cache=True, verbose=True)
-        print 'finished reading boundary file'
-        domain.set_boundary({'back': Br,
+    print 'finished reading boundary file'
+    domain.set_boundary({'back': Bd,
                              'side': Bd,
                              'ocean': Bf})
-    else: 
-        print 'set ocean'               
-        domain.set_boundary({'back': Br,
-                             'side': Bd,
-                             'ocean': Bd}) 
 
     kwargs['input_start_time']=domain.starttime
@@ -200 +192 @@
     t0 = time.time()
 
-    for t in domain.evolve(yieldstep = 240, finaltime = kwargs['starttime']): 
+    for t in domain.evolve(yieldstep = 240, finaltime = kwargs['finaltime']): 
         domain.write_time()
         domain.write_boundary_statistics(tags = 'ocean')     
-
-        if allclose(t, 240):
-            domain.set_boundary({'back': Br, 'side': Bd, 'ocean': Bo})
-
-        if allclose(t, 1440):
-            domain.set_boundary({'back': Br, 'side': Bd, 'ocean': Bd})
-
-#    for t in domain.evolve(yieldstep = kwargs['yieldstep'], finaltime = kwargs['midtime']
-#                       ,skip_initial_step = True): 
-#        domain.write_time()
-#        domain.write_boundary_statistics(tags = 'ocean')     
-#    
-#    for t in domain.evolve(yieldstep = 240, finaltime = kwargs['finaltime']
-#                       ,skip_initial_step = True): 
-#        domain.write_time()
-#        domain.write_boundary_statistics(tags = 'ocean')   
 
     x, y = domain.get_maximum_inundation_location()
@@ -235 +211 @@
     
     print 'memory usage before del domain1',mem_usage()
-    
-def export_model(**kwargs):    
-    #store_parameters(**kwargs)
-    
-#    print 'memory usage before del domain',mem_usage()
-    #del domain
-    print 'memory usage after del domain',mem_usage()
-    
-    swwfile = kwargs['output_dir']+kwargs['aa_scenario_name']
-    print'swwfile',swwfile
-    
-    export_grid(swwfile, extra_name_out = 'town',
-            quantities = ['speed','depth','elevation','stage'], # '(xmomentum**2 + ymomentum**2)**0.5' defaults to elevation
-            #quantities = ['speed','depth'], # '(xmomentum**2 + ymomentum**2)**0.5' defaults to elevation
-            timestep = None,
-            reduction = max,
-            cellsize = kwargs['export_cellsize'],
-            NODATA_value = -9999,
-            easting_min = project.eastingmin,
-            easting_max = project.eastingmax,
-            northing_min = project.northingmin,
-            northing_max = project.northingmax,
-            verbose = False,
-            origin = None,
-            datum = 'GDA94',
-            format = 'asc')
-            
-    inundation_damage(swwfile+'.sww', project.buildings_filename, 
-                      project.buildings_filename_out)
-    
-#-------------------------------------------------------------
+
+ #-------------------------------------------------------------
 if __name__ == "__main__":