Changeset 7074


Ignore:
Timestamp:
May 25, 2009, 3:02:58 PM (15 years ago)
Author:
ole
Message:

Removed verbosity from culvert unit tests

Location:
anuga_core/source/anuga/culvert_flows
Files:
5 edited

Legend:

Unmodified
Added
Removed
  • anuga_core/source/anuga/culvert_flows/culvert_routines.py

    r7073 r7074  
    8282    from anuga.config import velocity_protection
    8383    from anuga.utilities.numerical_tools import safe_acos as acos
    84     print "STARTING..BOYD................."
     84
    8585    local_debug ='false'
    8686    if inlet_depth > 0.1: #this value was 0.01:
  • anuga_core/source/anuga/culvert_flows/test_culvert_routines_box_10pct.py

    r7073 r7074  
    6969                                                 sum_loss)
    7070       
    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))
    7373        assert num.allclose(Q, 0.5526, rtol=1.0e-1) #inflow
    7474        assert num.allclose(v, 1.146, rtol=1.0e-1) #outflow velocity
     
    119119                                                 sum_loss)
    120120       
    121         print ('%s,%.2f,%.2f,%.2f' %('ANUGAcalcsTEST02 Q-v-d',Q,v,d))
    122         print ('%s,%.2f,%.2f,%.2f' %('Spreadsheet_Boydcalcs', 2.508, 1.897, 0.367))
     121        #print ('%s,%.2f,%.2f,%.2f' %('ANUGAcalcsTEST02 Q-v-d',Q,v,d))
     122        #print ('%s,%.2f,%.2f,%.2f' %('Spreadsheet_Boydcalcs', 2.508, 1.897, 0.367))
    123123        assert num.allclose(Q, 2.508, rtol=1.0e-1) #inflow
    124124        assert num.allclose(v, 1.897, rtol=1.0e-1) #outflow velocity
     
    168168                                                 manning,
    169169                                                 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))
    174174        assert num.allclose(Q, 13.554, rtol=1.0e-2) #inflow
    175175        assert num.allclose(v, 3.329, rtol=1.0e-2) #outflow velocity
     
    221221                                                 manning,
    222222                                                 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', 6.609, 2.621, 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', 6.609, 2.621, 0.70))
    227227        assert num.allclose(Q, 6.609, rtol=1.0e-2) #inflow
    228228        assert num.allclose(v, 2.621, rtol=1.0e-2) #outflow velocity
     
    273273                                                 manning,
    274274                                                 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' %('ANUGAcalcsTEST05 Q-v-d',Q,v,d))
    279         print ('%s,%.3f,%.3f,%.3f' %('Spreadsheet_Boydcalcs',2.961, 0.685, 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' %('ANUGAcalcsTEST05 Q-v-d',Q,v,d))
     279        #print ('%s,%.3f,%.3f,%.3f' %('Spreadsheet_Boydcalcs',2.961, 0.685, 1.20))
    280280        assert num.allclose(Q, 2.961, rtol=1.0e-2) #inflow
    281281        assert num.allclose(v, 0.685, rtol=1.0e-2) #outflow velocity
     
    326326                                                 manning,
    327327                                                 sum_loss)       
    328         print ('%s,%.3f'%('SPEC_E = ',inlet_specific_energy))
    329         print ('%s,%.3f'%('Delta E = ',delta_total_energy))
    330        
    331         print ('%s,%.3f,%.3f,%.3f' %('ANUGAcalcsTEST06 Q-v-d',Q,v,d))
    332         print ('%s,%.3f,%.3f,%.3f' %('Spreadsheet_Boydcalcs',15.537, 3.597, 1.20))
     328        #print ('%s,%.3f'%('SPEC_E = ',inlet_specific_energy))
     329        #print ('%s,%.3f'%('Delta E = ',delta_total_energy))
     330       
     331        #print ('%s,%.3f,%.3f,%.3f' %('ANUGAcalcsTEST06 Q-v-d',Q,v,d))
     332        #print ('%s,%.3f,%.3f,%.3f' %('Spreadsheet_Boydcalcs',15.537, 3.597, 1.20))
    333333        assert num.allclose(Q, 15.537, rtol=1.0e-2) #inflow
    334334        assert num.allclose(v, 3.597, rtol=1.0e-2) #outflow velocity
  • anuga_core/source/anuga/culvert_flows/test_culvert_routines_box_1pct.py

    r7073 r7074  
    6969                                                 sum_loss)
    7070       
    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))
    7373        assert num.allclose(Q, 0.5526, rtol=1.0e-1) #inflow
    7474        assert num.allclose(v, 1.146, rtol=1.0e-1) #outflow velocity
     
    119119                                                 sum_loss)
    120120       
    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))
    123123        assert num.allclose(Q, 0.224, rtol=1.0e-1) #inflow
    124124        assert num.allclose(v, 0.152, rtol=1.0e-1) #outflow velocity
     
    168168                                                 manning,
    169169                                                 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))
    174174        assert num.allclose(Q, 13.554, rtol=1.0e-2) #inflow
    175175        assert num.allclose(v, 3.329, rtol=1.0e-2) #outflow velocity
     
    221221                                                 manning,
    222222                                                 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))
    227227        assert num.allclose(Q, 5.164, rtol=1.0e-2) #inflow
    228228        assert num.allclose(v, 2.047, rtol=1.0e-2) #outflow velocity
     
    273273                                                 manning,
    274274                                                 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))
    280280        assert num.allclose(Q, 8.808, rtol=1.0e-2) #inflow
    281281        assert num.allclose(v, 2.039, rtol=1.0e-2) #outflow velocity
     
    327327                                                 manning,
    328328                                                 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))
    334334        assert num.allclose(Q, 13.546, rtol=1.0e-2) #inflow
    335335        assert num.allclose(v, 3.136, rtol=1.0e-2) #outflow velocity
  • anuga_core/source/anuga/culvert_flows/test_culvert_routines_pipe_10pct.py

    r7073 r7074  
    6969                                                 sum_loss)
    7070       
    71         print ('%s,%.2f,%.2f,%.2f' %('ANUGAcalcsTEST01 Q-v-d',Q,v,d))
    72         print('%s,%.2f,%.2f,%.2f' %('Spreadsheet_Boydcalcs', 0.113, 0.297, 0.443))
     71        #print ('%s,%.2f,%.2f,%.2f' %('ANUGAcalcsTEST01 Q-v-d',Q,v,d))
     72        #print('%s,%.2f,%.2f,%.2f' %('Spreadsheet_Boydcalcs', 0.113, 0.297, 0.443))
    7373        assert num.allclose(Q, 0.113, rtol=1.0e-1) #inflow
    7474        assert num.allclose(v, 0.297, rtol=1.0e-1) #outflow velocity
     
    119119                                                 sum_loss)
    120120       
    121         print ('%s,%.2f,%.2f,%.2f' %('ANUGAcalcsTEST02 Q-v-d',Q,v,d))
    122         print ('%s,%.2f,%.2f,%.2f' %('Spreadsheet_Boydcalcs', 0.585, 0.873, 0.687))
     121        #print ('%s,%.2f,%.2f,%.2f' %('ANUGAcalcsTEST02 Q-v-d',Q,v,d))
     122        #print ('%s,%.2f,%.2f,%.2f' %('Spreadsheet_Boydcalcs', 0.585, 0.873, 0.687))
    123123        assert num.allclose(Q, 0.585, rtol=1.0e-1) #inflow
    124124        assert num.allclose(v, 0.873, rtol=1.0e-1) #outflow velocity
     
    168168                                                 manning,
    169169                                                 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', 3.441, 3.042, 1.2))
     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', 3.441, 3.042, 1.2))
    174174        assert num.allclose(Q, 3.441, rtol=1.0e-2) #inflow
    175175        assert num.allclose(v, 3.042, rtol=1.0e-2) #outflow velocity
     
    221221                                                 manning,
    222222                                                 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', 1.676, 1.821, 0.91))
     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', 1.676, 1.821, 0.91))
    227227        assert num.allclose(Q, 1.676, rtol=1.0e-2) #inflow
    228228        assert num.allclose(v, 1.821, rtol=1.0e-2) #outflow velocity
     
    273273                                                 manning,
    274274                                                 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' %('ANUGAcalcsTEST05 Q-v-d',Q,v,d))
    279         print ('%s,%.3f,%.3f,%.3f' %('Spreadsheet_Boydcalcs',0.759, 0.671, 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' %('ANUGAcalcsTEST05 Q-v-d',Q,v,d))
     279        #print ('%s,%.3f,%.3f,%.3f' %('Spreadsheet_Boydcalcs',0.759, 0.671, 1.20))
    280280        assert num.allclose(Q, 0.759, rtol=1.0e-2) #inflow
    281281        assert num.allclose(v, 0.671, rtol=1.0e-2) #outflow velocity
     
    325325                                                 manning,
    326326                                                 sum_loss)       
    327         print ('%s,%.3f,%.3f,%.3f' %('ANUGAcalcsTEST06 Q-v-d',Q,v,d))
    328         print ('%s,%.3f,%.3f,%.3f' %('Spreadsheet_Boydcalcs',3.962, 3.503, 1.20))
     327        #print ('%s,%.3f,%.3f,%.3f' %('ANUGAcalcsTEST06 Q-v-d',Q,v,d))
     328        #print ('%s,%.3f,%.3f,%.3f' %('Spreadsheet_Boydcalcs',3.962, 3.503, 1.20))
    329329        assert num.allclose(Q, 3.962, rtol=1.0e-2) #inflow
    330330        assert num.allclose(v, 3.503, rtol=1.0e-2) #outflow velocity
  • anuga_core/source/anuga/culvert_flows/test_culvert_routines_pipe_1pct.py

    r7073 r7074  
    6969                                                 sum_loss)
    7070       
    71         print ('%s,%.2f,%.2f,%.2f' %('ANUGAcalcsTEST01 Q-v-d',Q,v,d))
    72         print('%s,%.2f,%.2f,%.2f' %('Spreadsheet_Boydcalcs', 0.113, 0.297, 0.443))
     71        #print ('%s,%.2f,%.2f,%.2f' %('ANUGAcalcsTEST01 Q-v-d',Q,v,d))
     72        #print('%s,%.2f,%.2f,%.2f' %('Spreadsheet_Boydcalcs', 0.113, 0.297, 0.443))
    7373        assert num.allclose(Q, 0.113, rtol=1.0e-1) #inflow
    7474        assert num.allclose(v, 0.297, rtol=1.0e-1) #outflow velocity
     
    119119                                                 sum_loss)
    120120       
    121         print ('%s,%.3f,%.3f,%.3f' %('ANUGAcalcsTEST02 Q-v-d',Q,v,d))
    122         print ('%s,%.3f,%.3f,%.3f' %('Spreadsheet_Boydcalcs', 0.108, 0.152, 0.721))
     121        #print ('%s,%.3f,%.3f,%.3f' %('ANUGAcalcsTEST02 Q-v-d',Q,v,d))
     122        #print ('%s,%.3f,%.3f,%.3f' %('Spreadsheet_Boydcalcs', 0.108, 0.152, 0.721))
    123123        assert num.allclose(Q, 0.108, rtol=1.0e-1) #inflow
    124124        assert num.allclose(v, 0.152, rtol=1.0e-1) #outflow velocity
     
    168168                                                 manning,
    169169                                                 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', 3.441, 3.042, 1.2))
     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', 3.441, 3.042, 1.2))
    174174        assert num.allclose(Q, 3.441, rtol=1.0e-2) #inflow
    175175        assert num.allclose(v, 3.042, rtol=1.0e-2) #outflow velocity
     
    221221                                                 manning,
    222222                                                 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', 1.559, 1.694, 0.91))
     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', 1.559, 1.694, 0.91))
    227227        assert num.allclose(Q, 1.559, rtol=1.0e-2) #inflow
    228228        assert num.allclose(v, 1.694, rtol=1.0e-2) #outflow velocity
     
    273273                                                 manning,
    274274                                                 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' %('ANUGAcalcsTEST05 Q-v-d',Q,v,d))
    279         print ('%s,%.3f,%.3f,%.3f' %('Spreadsheet_Boydcalcs',2.258, 1.996, 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' %('ANUGAcalcsTEST05 Q-v-d',Q,v,d))
     279        #print ('%s,%.3f,%.3f,%.3f' %('Spreadsheet_Boydcalcs',2.258, 1.996, 1.20))
    280280        assert num.allclose(Q, 2.258, rtol=1.0e-2) #inflow
    281281        assert num.allclose(v, 1.996, rtol=1.0e-2) #outflow velocity
     
    325325                                                 manning,
    326326                                                 sum_loss)       
    327         print ('%s,%.3f,%.3f,%.3f' %('ANUGAcalcsTEST06 Q-v-d',Q,v,d))
    328         print ('%s,%.3f,%.3f,%.3f' %('Spreadsheet_Boydcalcs',3.472, 3.070, 1.20))
     327        #print ('%s,%.3f,%.3f,%.3f' %('ANUGAcalcsTEST06 Q-v-d',Q,v,d))
     328        #print ('%s,%.3f,%.3f,%.3f' %('Spreadsheet_Boydcalcs',3.472, 3.070, 1.20))
    329329        assert num.allclose(Q, 3.472, rtol=1.0e-2) #inflow
    330330        assert num.allclose(v, 3.070, rtol=1.0e-2) #outflow velocity
Note: See TracChangeset for help on using the changeset viewer.