Changeset 8263
- Timestamp:
- Dec 5, 2011, 2:31:59 PM (13 years ago)
- Location:
- trunk
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/anuga_core/source/anuga_parallel/test_all.py
r8256 r8263 20 20 #List files that should be excluded from the testing process. 21 21 #E.g. if they are known to fail and under development 22 exclude_files = ['test_failure.py']23 22 24 23 # Directories that should not be searched for test files. -
trunk/anuga_work/development/2010-projects/anuga_1d/base/generic_domain.py
r8251 r8263 1395 1395 1396 1396 1397 Stage = self.quantities['stage']1398 Xmom = self.quantities['xmomentum']1399 Bed = self.quantities['elevation']1400 Height = self.quantities['height']1401 Velocity = self.quantities['velocity']1402 1403 #Arrays1404 w_C = Stage.centroid_values1405 uh_C = Xmom.centroid_values1406 z_C = Bed.centroid_values1407 h_C = Height.centroid_values1408 u_C = Velocity.centroid_values1409 1410 import numpy as np1411 print '== before forcing ======='1412 print 'w_C', np.any(np.isnan(w_C))1413 print 'uh_C', np.any(np.isnan(uh_C))1414 print 'z_C', np.any(np.isnan(z_C))1415 print 'h_C', np.any(np.isnan(h_C))1416 print 'u_C', np.any(np.isnan(u_C))1417 print 'w_ex_update', np.any(np.isnan(Stage.explicit_update))1397 # Stage = self.quantities['stage'] 1398 # Xmom = self.quantities['xmomentum'] 1399 # Bed = self.quantities['elevation'] 1400 # Height = self.quantities['height'] 1401 # Velocity = self.quantities['velocity'] 1402 # 1403 # #Arrays 1404 # w_C = Stage.centroid_values 1405 # uh_C = Xmom.centroid_values 1406 # z_C = Bed.centroid_values 1407 # h_C = Height.centroid_values 1408 # u_C = Velocity.centroid_values 1409 # 1410 # import numpy as np 1411 # print '== before forcing =======' 1412 # print 'w_C', np.any(np.isnan(w_C)) 1413 # print 'uh_C', np.any(np.isnan(uh_C)) 1414 # print 'z_C', np.any(np.isnan(z_C)) 1415 # print 'h_C', np.any(np.isnan(h_C)) 1416 # print 'u_C', np.any(np.isnan(u_C)) 1417 # print 'w_ex_update', np.any(np.isnan(Stage.explicit_update)) 1418 1418 1419 1419 timestep = self.timestep … … 1424 1424 self.compute_forcing_terms() 1425 1425 1426 print '==after forcing ======='1427 print 'w_C', np.any(np.isnan(w_C))1428 print 'uh_C', np.any(np.isnan(uh_C))1429 print 'z_C', np.any(np.isnan(z_C))1430 print 'h_C', np.any(np.isnan(h_C))1431 print 'u_C', np.any(np.isnan(u_C))1432 print 'w_ex_update', np.any(np.isnan(Stage.explicit_update))1426 # print '==after forcing =======' 1427 # print 'w_C', np.any(np.isnan(w_C)) 1428 # print 'uh_C', np.any(np.isnan(uh_C)) 1429 # print 'z_C', np.any(np.isnan(z_C)) 1430 # print 'h_C', np.any(np.isnan(h_C)) 1431 # print 'u_C', np.any(np.isnan(u_C)) 1432 # print 'w_ex_update', np.any(np.isnan(Stage.explicit_update)) 1433 1433 1434 1434 #Update conserved_quantities … … 1437 1437 Q.update(timestep) 1438 1438 1439 print '==after quantity update ======='1440 print 'w_C', np.any(np.isnan(w_C))1441 print 'uh_C', np.any(np.isnan(uh_C))1442 print 'z_C', np.any(np.isnan(z_C))1443 print 'h_C', np.any(np.isnan(h_C))1444 print 'u_C', np.any(np.isnan(u_C))1445 print 'w_ex_update', np.any(np.isnan(Stage.explicit_update))1439 # print '==after quantity update =======' 1440 # print 'w_C', np.any(np.isnan(w_C)) 1441 # print 'uh_C', np.any(np.isnan(uh_C)) 1442 # print 'z_C', np.any(np.isnan(z_C)) 1443 # print 'h_C', np.any(np.isnan(h_C)) 1444 # print 'u_C', np.any(np.isnan(u_C)) 1445 # print 'w_ex_update', np.any(np.isnan(Stage.explicit_update)) 1446 1446 1447 1447 if __name__ == "__main__": -
trunk/anuga_work/development/2010-projects/anuga_1d/sqpipe/parabolic_canal.py
r8254 r8263 69 69 N = 100 70 70 print "Evaluating domain with %d cells" %N 71 domain = dom.Domain(*uniform_mesh(N, x_0 = -2.0*L_x, x_1 = 2.0*L_x), bulk_modulus = 75.0)71 domain = dom.Domain(*uniform_mesh(N, x_0 = -2.0*L_x, x_1 = 2.0*L_x), bulk_modulus = 100.0) 72 72 73 73 domain.set_spatial_order(2) -
trunk/anuga_work/development/2010-projects/anuga_1d/sqpipe/sqpipe_domain.py
r8254 r8263 170 170 171 171 self.state = new_state 172 #self.state = numpy.where(h>=t, 1, 0)172 self.state = numpy.where(h>=t, 1, 0) 173 173 174 174 def distribute_to_vertices_and_edges(self): -
trunk/anuga_work/development/2010-projects/anuga_1d/sqpipe/test_sqpipe_domain.py
r8254 r8263 6 6 domain = anuga_1d.sqpipe.parabolic_canal.get_domain(dom) 7 7 finaltime = 10000.0 8 yieldstep = 50.08 yieldstep = 10.0 9 9 anuga_1d.sqpipe.parabolic_canal.animate_domain(domain, yieldstep, finaltime) 10 10 anuga_1d.sqpipe.parabolic_canal.plot_domain(domain)
Note: See TracChangeset
for help on using the changeset viewer.