- Timestamp:
- May 25, 2009, 3:02:58 PM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_core/source/anuga/culvert_flows/test_culvert_routines_box_1pct.py
r7073 r7074 69 69 sum_loss) 70 70 71 print ('%s,%.2f,%.2f,%.2f' %('ANUGAcalcsTEST01 Q-v-d',Q,v,d))72 print('%s,%.2f,%.2f,%.2f' %('Spreadsheet_Boydcalcs', 0.5526, 1.146, 0.1339))71 #print ('%s,%.2f,%.2f,%.2f' %('ANUGAcalcsTEST01 Q-v-d',Q,v,d)) 72 #print('%s,%.2f,%.2f,%.2f' %('Spreadsheet_Boydcalcs', 0.5526, 1.146, 0.1339)) 73 73 assert num.allclose(Q, 0.5526, rtol=1.0e-1) #inflow 74 74 assert num.allclose(v, 1.146, rtol=1.0e-1) #outflow velocity … … 119 119 sum_loss) 120 120 121 print ('%s,%.2f,%.2f,%.2f' %('ANUGAcalcsTEST02 Q-v-d',Q,v,d))122 print ('%s,%.2f,%.2f,%.2f' %('Spreadsheet_Boydcalcs', 0.224, 0.152, 0.409))121 #print ('%s,%.2f,%.2f,%.2f' %('ANUGAcalcsTEST02 Q-v-d',Q,v,d)) 122 #print ('%s,%.2f,%.2f,%.2f' %('Spreadsheet_Boydcalcs', 0.224, 0.152, 0.409)) 123 123 assert num.allclose(Q, 0.224, rtol=1.0e-1) #inflow 124 124 assert num.allclose(v, 0.152, rtol=1.0e-1) #outflow velocity … … 168 168 manning, 169 169 sum_loss) 170 print ('%s,%.2f'%('SPEC_E = ',inlet_specific_energy))171 print ('%s,%.2f'%('Delta E = ',delta_total_energy))172 print ('%s,%.2f,%.2f,%.2f' %('ANUGAcalcsTEST03 Q-v-d',Q,v,d))173 print ('%s,%.2f,%.2f,%.2f' %('Spreadsheet_Boydcalcs', 13.554, 3.329, 1.131))170 #print ('%s,%.2f'%('SPEC_E = ',inlet_specific_energy)) 171 #print ('%s,%.2f'%('Delta E = ',delta_total_energy)) 172 #print ('%s,%.2f,%.2f,%.2f' %('ANUGAcalcsTEST03 Q-v-d',Q,v,d)) 173 #print ('%s,%.2f,%.2f,%.2f' %('Spreadsheet_Boydcalcs', 13.554, 3.329, 1.131)) 174 174 assert num.allclose(Q, 13.554, rtol=1.0e-2) #inflow 175 175 assert num.allclose(v, 3.329, rtol=1.0e-2) #outflow velocity … … 221 221 manning, 222 222 sum_loss) 223 print ('%s,%.2f'%('SPEC_E = ',inlet_specific_energy))224 print ('%s,%.2f'%('Delta E = ',delta_total_energy))225 print ('%s,%.2f,%.2f,%.2f' %('ANUGAcalcsTEST04 Q-v-d',Q,v,d))226 print ('%s,%.2f,%.2f,%.2f' %('Spreadsheet_Boydcalcs', 5.164, 2.047, 0.70))223 #print ('%s,%.2f'%('SPEC_E = ',inlet_specific_energy)) 224 #print ('%s,%.2f'%('Delta E = ',delta_total_energy)) 225 #print ('%s,%.2f,%.2f,%.2f' %('ANUGAcalcsTEST04 Q-v-d',Q,v,d)) 226 #print ('%s,%.2f,%.2f,%.2f' %('Spreadsheet_Boydcalcs', 5.164, 2.047, 0.70)) 227 227 assert num.allclose(Q, 5.164, rtol=1.0e-2) #inflow 228 228 assert num.allclose(v, 2.047, rtol=1.0e-2) #outflow velocity … … 273 273 manning, 274 274 sum_loss) 275 print ('%s,%.3f'%('SPEC_E = ',inlet_specific_energy))276 print ('%s,%.3f'%('Delta E = ',delta_total_energy))277 278 print ('%s,%.3f,%.3f,%.3f' %('ANUGAcalcsTEST05Q-v-d',Q,v,d))279 print ('%s,%.3f,%.3f,%.3f' %('Spreadsheet_Boydcalcs',8.808, 2.039, 1.20))275 #print ('%s,%.3f'%('SPEC_E = ',inlet_specific_energy)) 276 #print ('%s,%.3f'%('Delta E = ',delta_total_energy)) 277 278 #print ('%s,%.3f,%.3f,%.3f' %('ANUGAcalcsTEST05Q-v-d',Q,v,d)) 279 #print ('%s,%.3f,%.3f,%.3f' %('Spreadsheet_Boydcalcs',8.808, 2.039, 1.20)) 280 280 assert num.allclose(Q, 8.808, rtol=1.0e-2) #inflow 281 281 assert num.allclose(v, 2.039, rtol=1.0e-2) #outflow velocity … … 327 327 manning, 328 328 sum_loss) 329 print ('%s,%.3f'%('SPEC_E = ',inlet_specific_energy))330 print ('%s,%.3f'%('Delta E = ',delta_total_energy))331 332 print ('%s,%.3f,%.3f,%.3f' %('ANUGAcalcsTEST06 Q-v-d',Q,v,d))333 print ('%s,%.3f,%.3f,%.3f' %('Spreadsheet_Boydcalcs',13.546, 3.136, 1.20))329 #print ('%s,%.3f'%('SPEC_E = ',inlet_specific_energy)) 330 #print ('%s,%.3f'%('Delta E = ',delta_total_energy)) 331 332 #print ('%s,%.3f,%.3f,%.3f' %('ANUGAcalcsTEST06 Q-v-d',Q,v,d)) 333 #print ('%s,%.3f,%.3f,%.3f' %('Spreadsheet_Boydcalcs',13.546, 3.136, 1.20)) 334 334 assert num.allclose(Q, 13.546, rtol=1.0e-2) #inflow 335 335 assert num.allclose(v, 3.136, rtol=1.0e-2) #outflow velocity
Note: See TracChangeset
for help on using the changeset viewer.