- Timestamp:
- Dec 5, 2011, 2:31:59 PM (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
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__":
Note: See TracChangeset
for help on using the changeset viewer.