Ignore:
Timestamp:
Mar 19, 2008, 8:02:05 PM (16 years ago)
Author:
steve
Message:

Ole and I seem to have got edge limiteri working.
Set domain.use_edge_limiter = True
to test the method

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_core/source/anuga/shallow_water/test_data_manager.py

    r5162 r5175  
    244244        # Get the variables
    245245        range = fid.variables['stage_range'][:]
     246        #print range
    246247        assert allclose(range,[-0.93519, 0.15]) or\
    247                allclose(range,[-0.9352743, 0.15]) # Old slope limiters
     248               allclose(range,[-0.9352743, 0.15]) or\
     249               allclose(range,[-0.93522203, 0.15000001]) # Old slope limiters
    248250       
    249251        range = fid.variables['xmomentum_range'][:]
    250 
    251 
    252 
     252        #print range
    253253        assert allclose(range,[0,0.4695096]) or \
    254254               allclose(range,[0,0.47790655]) or\
    255                allclose(range,[0,0.46941957])
     255               allclose(range,[0,0.46941957]) or\
     256               allclose(range,[0,0.47769409])
     257
    256258       
    257259        range = fid.variables['ymomentum_range'][:]
    258         #assert allclose(range,[0,0.02174380])
    259        
    260         assert allclose(range,[0,0.02174439]) or\
     260        #print range
     261        assert allclose(range,[0,0.02174380]) or\
     262              allclose(range,[0,0.02174439]) or\
    261263               allclose(range,[0,0.02283983]) or\
    262                allclose(range,[0,0.0217342]) # Old slope limiters
     264               allclose(range,[0,0.0217342]) or\
     265               allclose(range,[0,0.0227564]) # Old slope limiters
    263266       
    264267        fid.close()
     
    330333
    331334        extrema = fid.variables['xmomentum.extrema'][:]
    332         assert allclose(extrema,[-0.06062178, 0.47873023])
     335        assert allclose(extrema,[-0.06062178, 0.47873023]) or allclose(extrema, [-0.06062178, 0.47847986])
    333336       
    334337        extrema = fid.variables['ymomentum.extrema'][:]
    335         assert allclose(extrema,[0.00, 0.0625786])        
     338        assert allclose(extrema,[0.00, 0.0625786]) or allclose(extrema,[0.00, 0.06062178])
    336339
    337340        time_interval = fid.variables['extrema.time_interval'][:]
Note: See TracChangeset for help on using the changeset viewer.