Changeset 7647 for anuga_work/production/bunbury_storm_surge_2009/run_model.py

Timestamp:

Mar 2, 2010, 2:31:30 PM (13 years ago)

Author:

fountain

Message:

updates to mandurah and bunbury storm surge models

File:

• anuga_work/production/bunbury_storm_surge_2009/run_model.py (modified) (2 diffs)

anuga_work/production/bunbury_storm_surge_2009/run_model.py

@@ -159 +159 @@
                      'ocean': Bf}) 
 
-#-----------------------------------------------------------------------------
-#  Setup dynamic topography
-#-----------------------------------------------------------------------------
-from project import floodgate_height, floodgate_thickness,floodgate_LHS_x,floodgate_LHS_y,floodgate_RHS_x,floodgate_RHS_y,floodgate_start,floodgate_duration
-
-floodgate_end = floodgate_start + floodgate_duration #the time in seconds at which the floodgates are fully closed
-
-def floodgate(x,y):
-    z = 0*x
-    N = len(x)
-    for i in range(N):
-        ymin = ((floodgate_RHS_y - floodgate_LHS_y)/(floodgate_RHS_x - floodgate_LHS_x))*(x[i]-floodgate_LHS_x)+ floodgate_LHS_y - 0.5*floodgate_thickness
-        ymax = ((floodgate_RHS_y - floodgate_LHS_y)/(floodgate_RHS_x - floodgate_LHS_x))*(x[i]-floodgate_LHS_x)+ floodgate_LHS_y + 0.5*floodgate_thickness
-        xmin = floodgate_LHS_x
-        xmax = floodgate_RHS_x
-        if ymin < y[i] < ymax and xmin < x[i] < xmax:
-            z[i] += (floodgate_height/((floodgate_duration)/project.yieldstep))
-
-    return z
-    
 #-------------------------------------------------------------------------------
 # Evolve system through time
 #-------------------------------------------------------------------------------
-
-raising_floodgates = False
-lowering_floodgates = False
 
 t0 = time.time()
@@ -193 +170 @@
     print domain.boundary_statistics(tags='ocean')
 
-#raise floodgate at a given time interval (to make floodgate lower again repeat this code below
-#for a later time interval and change rising to false and falling to true
-    if floodgate_start < t < floodgate_end:
-        rising = True
-        falling = False
-
-    if rising:
-        domain.add_quantity('elevation',floodgate)
-
-    if falling:
-        domain.add_quantity('elevation', lambda x,y: -1*floodgate(x,y))
 
 print 'Simulation took %.2f seconds' % (time.time()-t0)