Changeset 5439


Ignore:
Timestamp:
Jun 26, 2008, 10:58:08 AM (16 years ago)
Author:
ole
Message:

Updated 'two_levels' example to show that tight_slope_limiters take care of upwards momentum or 'creep'. This used to be treated with the h-limiter with beta_h > 0 but this is no longer needed.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/development/demos/twolevels.py

    r4539 r5439  
    99#
    1010from os import sep, path
    11 from mesh_factory import rectangular
    12 from shallow_water import Domain, Reflective_boundary, Dirichlet_boundary,\
    13      Constant_height
     11from anuga.abstract_2d_finite_volumes.mesh_factory import rectangular
     12from shallow_water import Domain, Reflective_boundary, Dirichlet_boundary
    1413from Numeric import array
    15 from anuga.pyvolution.util import Polygon_function, read_polygon
     14from anuga.utilities.polygon import Polygon_function, read_polygon
    1615
    1716
     
    3332#PLAY WITH THIS [0;1]:
    3433#
    35 # beta_h == 0.0 reveals the problem
    36 # beta_h > 0.2 alleviates it
    37 domain.beta_h = 0.2   
     34# To reveal the problem set
     35# domain.tight_slope_limiters = False
     36domain.tight_slope_limiters = False
     37# To fix it set
     38#domain.tight_slope_limiters = True
     39domain.tight_slope_limiters = True
     40
    3841
    3942#IC
Note: See TracChangeset for help on using the changeset viewer.