Changeset 7807


Ignore:
Timestamp:
Jun 8, 2010, 9:48:54 AM (10 years ago)
Author:
hudson
Message:

Cairns demo works with new API.

Location:
trunk/anuga_core/documentation/user_manual/demos/cairns
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/anuga_core/documentation/user_manual/demos/cairns/project.py

    r7801 r7807  
    1 """Common filenames and locations for topographic data, meshes and outputs."""
     1""" Common filenames and locations for topographic data, meshes and outputs.
     2    This file defines the parameters of the scenario you wish to run.
     3"""
    24
    35from anuga.geometry.polygon import read_polygon, plot_polygons, \
     
    4547
    4648# Define resolutions (max area per triangle) for each polygon
    47 default_res = 1000000000    # Background resolution
    48 islands_res = 10000000
    49 cairns_res = 10000000
    50 shallow_res = 50000000
     49# Make these numbers larger to reduce the number of triangles in the model,
     50# and hence speed up the simulation
     51default_res = 10000000    # Background resolution
     52islands_res = 100000
     53cairns_res = 100000
     54shallow_res = 500000
    5155
    5256# Define list of interior regions with associated resolutions
  • trunk/anuga_core/documentation/user_manual/demos/cairns/runcairns.py

    r7801 r7807  
    8181if project.scenario == 'slide':
    8282    # Function for submarine slide
    83     tsunami_source = slide_tsunami(length=35000.0,
     83    tsunami_source = anuga.slide_tsunami(length=35000.0,
    8484                                   depth=project.slide_depth,
    8585                                   slope=6.0,
     
    9696print 'Available boundary tags', domain.get_boundary_tags()
    9797
    98 Bd = Dirichlet_boundary([tide, 0, 0]) # Mean water level
    99 Bs = Transmissive_stage_zero_momentum_boundary(domain) # Neutral boundary
     98Bd = anuga.Dirichlet_boundary([tide, 0, 0]) # Mean water level
     99Bs = anuga.Transmissive_stage_zero_momentum_boundary(domain) # Neutral boundary
    100100
    101101if project.scenario == 'fixed_wave':
    102102    # Huge 50m wave starting after 60 seconds and lasting 1 hour.
    103     Bw = Time_boundary(domain=domain,
     103    Bw = anuga.Time_boundary(domain=domain,
    104104                       function=lambda t: [(60<t<3660)*50, 0, 0])
    105105    domain.set_boundary({'ocean_east': Bw,
     
    121121t0 = time.time()
    122122
    123 from Numeric import allclose
     123from numpy import allclose
    124124
    125125if project.scenario == 'slide':
Note: See TracChangeset for help on using the changeset viewer.