Changeset 5280
 Timestamp:
 May 6, 2008, 4:09:51 PM (16 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

anuga_core/source/anuga/shallow_water/test_tsunami_okada.py
r5279 r5280 249 249 tsunami = Quantity(domain) 250 250 tsunami.set_values(Ts) 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 251 271 252 k=0.0 272 253 for i in range(0,6): … … 275 256 Yt=p 276 257 Xt=k 277 interpolation_points.append([Xt, Yt]) 278 258 Z=tsunami.get_values(interpolation_points=[[Xt,Yt]] 259 ,location='edges') 260 stage.append(Z[0]) 279 261 k=k+4000 280 281 <<<<<<< .mine282 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 downward287 #print 'c est fini'288 #print tmp289 #print 'hello',stage290 assert allclose(stage,tmp,atol=1.e3)291 =======292 Z=tsunami.get_values(interpolation_points=interpolation_points,293 location='edges')294 >>>>>>> .r5278295 296 stage = Z # FIXME(Ole): Why the sign flip?297 262 298 263 assert allclose(stage, tmp, atol=1.e3) 264 print 'c est fini' 299 265 300 266 #
Note: See TracChangeset
for help on using the changeset viewer.