Changeset 8005
 Timestamp:
 Sep 4, 2010, 9:47:43 PM (12 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

trunk/anuga_core/source/anuga/structures/testing_wide_bridge.py
r7995 r8005 166 166 domain.forcing_terms.append(culvert_energy) 167 167 """ 168 from anuga.structures.boyd_box_operator import Boyd_box_operator168 #from anuga.structures.boyd_box_operator import Boyd_box_operator 169 169 #culvert0 = Culvert_operator(domain, 170 170 # end_point0=[40.0, 75.0], … … 180 180 # 181 181 culverts = [] 182 number_of_culverts = 1 182 number_of_culverts = 1 183 183 for i in range(number_of_culverts): 184 184 culvert_width = 50.0/number_of_culverts … … 186 186 ep0 = [40.0, y] 187 187 ep1 = [50.0, y] 188 culverts.append( Boyd_box_operator(domain,188 culverts.append(anuga.Boyd_box_operator(domain, 189 189 end_point0=ep0, 190 190 end_point1=ep1, 191 width=culvert_width, 192 height=10.0, 191 losses=1.5, 192 width=3.658, 193 height=3.658, 193 194 apron=5.0, 194 195 use_momentum_jet=True, 195 use_velocity_head= False,196 use_velocity_head=True, 196 197 manning=0.013, 198 description='bridge culvert', 197 199 verbose=False)) 198 200 … … 244 246 #Btds = Time_boundary(domain, lambda t: [0.0+ 0.0075*(1+num.sin(2*pi*(t4)/20)), 0.0, 0.0]) 245 247 246 Btus = anuga.Dirichlet_boundary([ 18.0, 0, 0]) # Outflow water at 5.0247 Btds = anuga.Dirichlet_boundary([ 0.0, 0, 0]) # Outflow water at 1.0248 Btus = anuga.Dirichlet_boundary([20.0, 0, 0]) # Outflow water at 5.0 249 Btds = anuga.Dirichlet_boundary([19.0, 0, 0]) # Outflow water at 1.0 248 250 domain.set_boundary({'left': Btus, 'right': Btds, 'top': Br, 'bottom': Br}) 249 251 … … 253 255 # 254 256 255 for t in domain.evolve(yieldstep = 1, finaltime = 10 0):257 for t in domain.evolve(yieldstep = 1, finaltime = 10): 256 258 print domain.timestepping_statistics() 257 259 print domain.volumetric_balance_statistics() 258 for culvert in culverts: 260 for i, culvert in enumerate(culverts): 261 print 'culvert: ', i 259 262 print culvert.structure_statistics() 260 263
Note: See TracChangeset
for help on using the changeset viewer.