Changeset 5408 for anuga_work

Timestamp:

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

Author:

kristy

Message:

 

Location:

anuga_work/production

Files:

• busselton/project.py (modified) (1 diff)
• busselton/run_busselton.py (modified) (4 diffs)
• exmouth_2006/export_results.py (modified) (4 diffs)
• exmouth_2006/project.py (modified) (1 diff)
• exmouth_2006/run_exmouth.py (modified) (4 diffs)
• perth/project.py (modified) (2 diffs)
• perth/run_perth.py (modified) (1 diff)

anuga_work/production/busselton/project.py

@@ -37 +37 @@
 alpha = 0.1
 friction=0.01
-starttime=10000
+starttime=0
 midtime=21600
 #finaltime=25000
-finaltime=10000
-export_cellsize=50
+finaltime=86000
+export_cellsize=25
 setup='final'
 source='exmouth'

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__":
     

anuga_work/production/exmouth_2006/export_results.py

@@ -10 +10 @@
 #time_dir = '20070619_042140_run_final_0_onslow_nbartzis' # MSL 1 in Onslow 10000 yr
 #time_dir = '20070619_042542_run_final_1.4_exmouth_nbartzis' # HAT 1 in Exmouth 10000 yr
-time_dir = '20070619_042417_run_final_0_exmouth_nbartzis' # MSL 1 in Exmouth 10000 yr
+#time_dir = '20070619_042417_run_final_0_exmouth_nbartzis' # MSL 1 in Exmouth 10000 yr
+#time_dir = '20080610_041607_run_final_1.4_exmouth_kvanputt' # HAT 1 in Exmouth 10000 yr
+time_dir = '20080613_055032_run_final_1.4_exmouth_kvanputt' # HAT 1 in Exmouth 10000 yr slope 1
+
 cellsize = 25
 timestep = None
@@ -17 +20 @@
 name = directory+time_dir+sep+project.scenario_name
 
-is_parallel = True
+is_parallel = False
 if is_parallel == True: nodes = 4
 print 'output dir:', name
 
-var = [2,3] # depth and speed
+var = [2,3, 4] # depth and speed
 
 for which_var in var:
@@ -33 +36 @@
 
     if which_var == 2:  # Depth
-        outname = name + '_depth_north'
+        outname = name + '_depth'
         quantityname = 'stage-elevation'  
 
     if which_var == 3:  # Speed
-        outname = name + '_speed_north'
+        outname = name + '_speed'
         quantityname = '(xmomentum**2 + ymomentum**2)**0.5/(stage-elevation+1.e-6/(stage-elevation))'  #Speed
 
@@ -67 +70 @@
                     timestep = timestep,
                     cellsize = cellsize,      
-                    easting_min = project.e_min_area,
-                    easting_max = project.e_max_area,
-                    northing_min = project.n_min_area,
-                    northing_max = project.n_max_area,        
+                    easting_min = project.eastingmin,
+                    easting_max = project.eastingmax,
+                    northing_min = project.northingmin,
+                    northing_max = project.northingmax,        
                     reduction = max, 
                     verbose = True,

anuga_work/production/exmouth_2006/project.py

@@ -133 +133 @@
 output_dir = anuga_dir+'outputs'+sep
 output_build_time_dir = output_dir+build_time+sep
-#output_run_time_dir = output_dir +run_time+dir_comment+sep
-output_run_time_dir = output_dir +'20070619_042140_run_final_0_onslow_nbartzis'+sep
+output_run_time_dir = output_dir +run_time+dir_comment+sep
+#output_run_time_dir = output_dir +'20070619_042140_run_final_0_onslow_nbartzis'+sep
 output_run_time_dir_name = output_run_time_dir + scenario_name  #Used by post processing
 

anuga_work/production/exmouth_2006/run_exmouth.py

@@ -66 +66 @@
 #    kwargs['bathy_file']=project.combined_small_dir_name + '.txt'
     kwargs['boundary_file']=project.boundaries_in_dir_name + '.sww'
-    '''
+  
     print 'output_dir',kwargs['output_dir']
     if myid == 0:
@@ -108 +108 @@
                              use_cache=False,
                              verbose=True)
-        copy(project.meshes_dir_name+'.msh',project.output_dir+project.scenario_name+'.msh')
     barrier()
 
@@ -170 +169 @@
     domain.beta_h = 0 #sets the surface of the triangle to follow the bathy
     #domain.H0=0.01 #controls the flux limiter (limiter2007)
-    #domain.tight_slope_limiters = 1 #minimises creep
+    domain.tight_slope_limiters = 0 #minimises creep
 
     #-------------------------------------------------------------------------
@@ -227 +226 @@
     kwargs['completed']=str(time.time()-t0)
     
-    '''
+    
     if myid == 0:
         store_parameters(**kwargs)

anuga_work/production/perth/project.py

@@ -38 +38 @@
 friction=0.01
 starttime=10000
-midtime=21600
-finaltime=10000
+finaltime=80000
 export_cellsize=50
 setup='final'
@@ -213 +212 @@
 ###################################################################
 
-# exporting asc grid
-eastingmin = 379000
-eastingmax = 385000
-northingmin = 6450700
-northingmax = 6459800
-
-
-
+
+#Geordie Bay extract ascii grid
+xminGeordie = 359500
+xmaxGeordie = 360850
+yminGeordie = 6459400
+ymaxGeordie = 6460430
+
+#Sorrento extract ascii grid
+xminSorrento = 379560
+xmaxSorrento = 380920
+yminSorrento = 6477750
+ymaxSorrento = 6479030
+
+#Perth extract ascii grid
+xminPerth = 379106
+xmaxPerth = 383100
+yminPerth = 6450500
+ymaxPerth = 6461150
+
+#Rockingham extract ascii grid
+xminRockingham = 375500
+xmaxRockingham = 383100
+yminRockingham = 6427900
+ymaxRockingham = 6430100
+

anuga_work/production/perth/run_perth.py

@@ -225 +225 @@
     kwargs['starttime']=project.starttime
     kwargs['yieldstep']=project.yieldstep
-    kwargs['midtime']=project.midtime
     kwargs['finaltime']=project.finaltime