Changeset 5281
 Timestamp:
 May 6, 2008, 4:18:43 PM (15 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

anuga_core/source/anuga/shallow_water/test_tsunami_okada.py
r5280 r5281 249 249 tsunami = Quantity(domain) 250 250 tsunami.set_values(Ts) 251 251 <<<<<<< .mine 252 interpolation_points=[] 253 ======= 254 255 #k=0.0 256 #for i in range(0,6): 257 # for j in range(0,6): 258 # p=j*4000 259 # Yt=p 260 # Xt=k 261 # Z=tsunami.get_values(interpolation_points=[[Xt,Yt]] 262 # ,location='edges') 263 # stage.append(Z[0]) 264 # k=k+4000 265 # 266 #assert allclose(stage,tmp,atol=1.e3) 267 268 # Here's a faster way  try that in the first test 269 interpolation_points=[] 270 >>>>>>> .r5278 252 271 k=0.0 253 272 for i in range(0,6): … … 256 275 Yt=p 257 276 Xt=k 258 Z=tsunami.get_values(interpolation_points=[[Xt,Yt]] 259 ,location='edges') 260 stage.append(Z[0]) 277 interpolation_points.append([Xt, Yt]) 278 261 279 k=k+4000 280 281 <<<<<<< .mine 282 Z=tsunami.get_values(interpolation_points=interpolation_points, 283 location='edges') 284 285 stage = Z # FIXME(Ole): Why the sign flip? 286 # Displacement in fortran code is looking downward 287 #print 'c est fini' 288 #print tmp 289 #print 'hello',stage 290 assert allclose(stage,tmp,atol=1.e3) 291 ======= 292 Z=tsunami.get_values(interpolation_points=interpolation_points, 293 location='edges') 294 >>>>>>> .r5278 295 296 stage = Z # FIXME(Ole): Why the sign flip? 262 297 263 298 assert allclose(stage, tmp, atol=1.e3) 264 print 'c est fini'265 299 266 300 #
Note: See TracChangeset
for help on using the changeset viewer.