N = 12 log.critical('Creating domain') #Create basic mesh points, vertices, boundary = rectangular(N, N / 2, len1=1.2, len2=0.6, origin=(-0.07, 0)) log.critical('Number of elements=%d' % len(vertices)) #Create shallow water domain domain = Domain(points, vertices, boundary) domain.smooth = False domain.default_order = 2 domain.set_name('show_balanced_limiters') domain.store = True domain.format = 'sww' #Native netcdf visualisation format #Set bed-slope and friction inflow_stage = 0.1 manning = 0.1 Z = Weir(inflow_stage) log.critical('Field values') domain.set_quantity('elevation', Z) domain.set_quantity('friction', manning) ###################### # Boundary conditions #
print 'Setting up domain' length = 40. width = 5. dx = dy = 1 # Resolution: Length of subdivisions on both axes #dx = dy = .5 # Resolution: Length of subdivisions on both axes #dx = dy = .5 # Resolution: Length of subdivisions on both axes #dx = dy = .1 # Resolution: Length of subdivisions on both axes points, vertices, boundary = rectangular_cross(int(length / dx), int(width / dy), len1=length, len2=width) domain = Domain(points, vertices, boundary) domain.set_name('Test_Culv_Flat_WL') # Output name domain.set_default_order(2) domain.H0 = 0.01 domain.tight_slope_limiters = 1 print 'Size', len(domain) #------------------------------------------------------------------------------ # Setup initial conditions #------------------------------------------------------------------------------ def topography(x, y): """Set up a weir A culvert will connect either side
# Domain # N = 12 log.critical('Creating domain') #Create basic mesh points, vertices, boundary = rectangular(N, N/2, len1=1.2,len2=0.6, origin=(-0.07, 0)) log.critical('Number of elements=%d' % len(vertices)) #Create shallow water domain domain = Domain(points, vertices, boundary) domain.smooth = False domain.default_order = 2 domain.set_name('show_balanced_limiters') domain.store = True domain.format = 'sww' #Native netcdf visualisation format #Set bed-slope and friction inflow_stage = 0.1 manning = 0.1 Z = Weir(inflow_stage) log.critical('Field values') domain.set_quantity('elevation', Z) domain.set_quantity('friction', manning) ###################### # Boundary conditions
# Setup computational domain #------------------------------------------------------------------------------ print 'Setting up domain' length = 40. width = 5. dx = dy = 1 # Resolution: Length of subdivisions on both axes #dx = dy = .5 # Resolution: Length of subdivisions on both axes #dx = dy = .5 # Resolution: Length of subdivisions on both axes #dx = dy = .1 # Resolution: Length of subdivisions on both axes points, vertices, boundary = rectangular_cross(int(length/dx), int(width/dy), len1=length, len2=width) domain = Domain(points, vertices, boundary) domain.set_name('Test_Culv_Flat_WL') # Output name domain.set_default_order(2) domain.H0 = 0.01 domain.tight_slope_limiters = 1 print 'Size', len(domain) #------------------------------------------------------------------------------ # Setup initial conditions #------------------------------------------------------------------------------ def topography(x, y): """Set up a weir A culvert will connect either side """