Changeset 6171
- Timestamp:
- Jan 14, 2009, 4:43:43 PM (16 years ago)
- Location:
- anuga_core/source/anuga/abstract_2d_finite_volumes
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_core/source/anuga/abstract_2d_finite_volumes/test_util.py
r6145 r6171 186 186 d_uh = num.reshape(num.take(xmomentum[last_time_index, :], [0,5,10,15]), (4,1)) 187 187 d_vh = num.reshape(num.take(ymomentum[last_time_index, :], [0,5,10,15]), (4,1)) 188 D = num.concatenate( 188 D = num.concatenate((d_stage, d_uh, d_vh), axis=1) 189 189 190 190 #Reference interpolated values at midpoints on diagonal at … … 222 222 d_uh = num.reshape(num.take(xmomentum[timestep, :], [0,5,10,15]), (4,1)) 223 223 d_vh = num.reshape(num.take(ymomentum[timestep, :], [0,5,10,15]), (4,1)) 224 D = num.concatenate( 224 D = num.concatenate((d_stage, d_uh, d_vh), axis=1) 225 225 226 226 #Reference interpolated values at midpoints on diagonal at … … 244 244 d_uh = num.reshape(num.take(xmomentum[timestep, :], [0,5,10,15]), (4,1)) 245 245 d_vh = num.reshape(num.take(ymomentum[timestep, :], [0,5,10,15]), (4,1)) 246 D = num.concatenate( 246 D = num.concatenate((d_stage, d_uh, d_vh), axis=1) 247 247 248 248 #Reference interpolated values at midpoints on diagonal at … … 265 265 d_uh = num.reshape(num.take(xmomentum[timestep, :], [0,5,10,15]), (4,1)) 266 266 d_vh = num.reshape(num.take(ymomentum[timestep, :], [0,5,10,15]), (4,1)) 267 D = num.concatenate( 267 D = num.concatenate((d_stage, d_uh, d_vh), axis=1) 268 268 269 269 #Reference interpolated values at midpoints on diagonal at … … 278 278 d_uh = num.reshape(num.take(xmomentum[timestep, :], [0,5,10,15]), (4,1)) 279 279 d_vh = num.reshape(num.take(ymomentum[timestep, :], [0,5,10,15]), (4,1)) 280 D = num.concatenate( 280 D = num.concatenate((d_stage, d_uh, d_vh), axis=1) 281 281 282 282 #Reference interpolated values at midpoints on diagonal at … … 389 389 d_uh = num.reshape(num.take(xmomentum[last_time_index, :], [0,5,10,15]), (4,1)) 390 390 d_vh = num.reshape(num.take(ymomentum[last_time_index, :], [0,5,10,15]), (4,1)) 391 D = num.concatenate( 391 D = num.concatenate((d_stage, d_uh, d_vh), axis=1) 392 392 393 393 #Reference interpolated values at midpoints on diagonal at -
anuga_core/source/anuga/abstract_2d_finite_volumes/util.py
r6169 r6171 2215 2215 for k, quantity in enumerate(quantities): 2216 2216 if quantity != 'time' and quantity != 'elevation': 2217 num = int(k*100+j) 2218 pylab.figure(num) 2217 pylab.figure(int(k*100+j)) 2219 2218 pylab.ylabel(quantities_label[quantity]) 2220 2219 pylab.plot(t,
Note: See TracChangeset
for help on using the changeset viewer.