Changeset 3704


Ignore:
Timestamp:
Oct 5, 2006, 8:50:08 PM (17 years ago)
Author:
steve
Message:

Changed beta defaults to allow more limiting at wet/dry interface

Location:
anuga_core/source/anuga
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • anuga_core/source/anuga/abstract_2d_finite_volumes/test_quantity.py

    r3689 r3704  
    10121012        """
    10131013        quantity = Conserved_quantity(self.mesh4)
    1014 
     1014        quantity.domain.beta_w = 0.9
     1015       
    10151016        #Test centroids
    10161017        quantity.set_values([2.,4.,8.,2.], location = 'centroids')
     
    10271028
    10281029        # limited value for beta_w = 0.9
     1030       
    10291031        assert allclose(quantity.vertex_values[1,:], [2.2, 4.9, 4.9])
    10301032        # limited values for beta_w = 0.5
  • anuga_core/source/anuga/config.py

    r3689 r3704  
    7070#
    7171#Good values are:
    72 beta_w      = 0.9
    73 beta_w_dry  = 0.9
    74 beta_uh     = 0.9
    75 beta_uh_dry = 0.9
    76 beta_vh     = 0.9
    77 beta_vh_dry = 0.9
    78 beta_h      = 0.2
     72
    7973
    8074# I think these are better SR but they conflict with the unit tests!
    81 # beta_w      = 1.0
    82 # beta_w_dry  = 0.2
    83 # beta_uh     = 1.0
    84 # beta_uh_dry = 0.2
    85 # beta_vh     = 1.0
    86 # beta_vh_dry = 0.2
    87 # beta_h      = 0.2
     75beta_w      = 1.0
     76beta_w_dry  = 0.2
     77beta_uh     = 1.0
     78beta_uh_dry = 0.2
     79beta_vh     = 1.0
     80beta_vh_dry = 0.2
     81beta_h      = 0.2
    8882
    8983
  • anuga_core/source/anuga/shallow_water/test_shallow_water_domain.py

    r3703 r3704  
    15171517        #Second order
    15181518        domain._order_ = 2
     1519        domain.beta_w      = 0.9
     1520        domain.beta_w_dry  = 0.9
     1521        domain.beta_uh     = 0.9
     1522        domain.beta_uh_dry = 0.9
     1523        domain.beta_vh     = 0.9
     1524        domain.beta_vh_dry = 0.9
    15191525        domain.distribute_to_vertices_and_edges()
    15201526        assert allclose(L[1], [2.2, 4.9, 4.9])
     
    15541560
    15551561        domain._order_ = 2
     1562        domain.beta_w      = 0.9
     1563        domain.beta_w_dry  = 0.9
     1564        domain.beta_uh     = 0.9
     1565        domain.beta_uh_dry = 0.9
     1566        domain.beta_vh     = 0.9
     1567        domain.beta_vh_dry = 0.9
    15561568        domain.distribute_to_vertices_and_edges()
    15571569        assert allclose(L[1], [0.57777777, 2.37777778, 2.37777778])
     
    15921604
    15931605        domain._order_ = 2
     1606        domain.beta_w      = 0.9
     1607        domain.beta_w_dry  = 0.9
     1608        domain.beta_uh     = 0.9
     1609        domain.beta_uh_dry = 0.9
     1610        domain.beta_vh     = 0.9
     1611        domain.beta_vh_dry = 0.9
    15941612        domain.distribute_to_vertices_and_edges()
    15951613        assert allclose(L[1], [1.07160494, 6.46058131, 7.28262855])
     
    17311749        #print E
    17321750        domain._order_ = 2
     1751        domain.beta_w      = 0.9
     1752        domain.beta_w_dry  = 0.9
     1753        domain.beta_uh     = 0.9
     1754        domain.beta_uh_dry = 0.9
     1755        domain.beta_vh     = 0.9
     1756        domain.beta_vh_dry = 0.9
    17331757        domain.beta_h = 0.0 #Use first order in h-limiter
    17341758        domain.distribute_to_vertices_and_edges()
     
    22262250        domain.smooth = False
    22272251        domain.default_order=2
    2228 
     2252        domain.beta_w      = 0.9
     2253        domain.beta_w_dry  = 0.9
     2254        domain.beta_uh     = 0.9
     2255        domain.beta_uh_dry = 0.9
     2256        domain.beta_vh     = 0.9
     2257        domain.beta_vh_dry = 0.9
     2258       
    22292259        # Boundary conditions
    22302260        Br = Reflective_boundary(domain)
     
    23712401        domain.visualise = False
    23722402        domain.default_order=2
     2403        domain.beta_w      = 0.9
     2404        domain.beta_w_dry  = 0.9
     2405        domain.beta_uh     = 0.9
     2406        domain.beta_uh_dry = 0.9
     2407        domain.beta_vh     = 0.9
     2408        domain.beta_vh_dry = 0.9       
    23732409        #domain.minimum_allowed_height = 0.0 #Makes it like the 'oldstyle' balance
    23742410
     
    24292465        domain.visualise = False
    24302466        domain.default_order=2
     2467        domain.beta_w      = 0.9
     2468        domain.beta_w_dry  = 0.9
     2469        domain.beta_uh     = 0.9
     2470        domain.beta_uh_dry = 0.9
     2471        domain.beta_vh     = 0.9
     2472        domain.beta_vh_dry = 0.9       
    24312473        domain.maximum_allowed_speed = 0.0 #Makes it like the 'oldstyle'
    24322474
     
    24752517        domain.visualise = False
    24762518        domain.default_order=domain._order_=2
     2519        domain.beta_w      = 0.9
     2520        domain.beta_w_dry  = 0.9
     2521        domain.beta_uh     = 0.9
     2522        domain.beta_uh_dry = 0.9
     2523        domain.beta_vh     = 0.9
     2524        domain.beta_vh_dry = 0.9       
    24772525
    24782526        # Boundary conditions
     
    27202768        domain.smooth = False
    27212769        domain.default_order=2
     2770        domain.beta_w      = 0.9
     2771        domain.beta_w_dry  = 0.9
     2772        domain.beta_uh     = 0.9
     2773        domain.beta_uh_dry = 0.9
     2774        domain.beta_vh     = 0.9
     2775        domain.beta_vh_dry = 0.9
    27222776
    27232777        #Bed-slope and friction at vertices (and interpolated elsewhere)
     
    28062860        domain.smooth = False
    28072861        domain.default_order=2
     2862        domain.beta_w      = 0.9
     2863        domain.beta_w_dry  = 0.9
     2864        domain.beta_uh     = 0.9
     2865        domain.beta_uh_dry = 0.9
     2866        domain.beta_vh     = 0.9
     2867        domain.beta_vh_dry = 0.9
    28082868        domain.beta_h = 0.0 #Use first order in h-limiter
    28092869
     
    28982958        domain.smooth = False
    28992959        domain.default_order=2
     2960        domain.beta_w      = 0.9
     2961        domain.beta_w_dry  = 0.9
     2962        domain.beta_uh     = 0.9
     2963        domain.beta_uh_dry = 0.9
     2964        domain.beta_vh     = 0.9
     2965        domain.beta_vh_dry = 0.9
    29002966        domain.beta_h = 0.0 #Use first order in h-limiter
    29012967
     
    29863052        domain.smooth = False
    29873053        domain.default_order=2
     3054        domain.beta_w      = 0.9
     3055        domain.beta_w_dry  = 0.9
     3056        domain.beta_uh     = 0.9
     3057        domain.beta_uh_dry = 0.9
     3058        domain.beta_vh     = 0.9
     3059        domain.beta_vh_dry = 0.9
    29883060        domain.beta_h = 0.0 #Use first order in h-limiter
    29893061
     
    31023174        domain.smooth = False
    31033175        domain.default_order=2
     3176        domain.beta_w      = 0.9
     3177        domain.beta_w_dry  = 0.9
     3178        domain.beta_uh     = 0.9
     3179        domain.beta_uh_dry = 0.9
     3180        domain.beta_vh     = 0.9
     3181        domain.beta_vh_dry = 0.9
    31043182        domain.beta_h = 0.0 #Use first order in h-limiter
    31053183
Note: See TracChangeset for help on using the changeset viewer.