Ignore:
Timestamp:
Jun 16, 2010, 8:23:42 AM (13 years ago)
Author:
hudson
Message:

Fixed failing validations

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_validation/convergence_study/dam_break.py

    r5442 r7846  
    1010
    1111import sys
    12 from anuga.abstract_2d_finite_volumes.mesh_factory import rectangular_cross
    13 from anuga.shallow_water import Domain
    14 from anuga.shallow_water import Reflective_boundary
    15 from anuga.shallow_water import Dirichlet_boundary
    16 from anuga.shallow_water import Time_boundary
    17 from anuga.shallow_water import Transmissive_boundary
    18 from anuga.shallow_water import Transmissive_Momentum_Set_Stage_boundary
    19 from anuga.shallow_water.data_manager import start_screen_catcher, copy_code_files
     12import anuga
    2013
    2114from math import cos
    22 from Numeric import zeros, Float
     15from numpy import zeros, float
    2316from time import localtime, strftime, gmtime
    2417
     
    3427output_file = 'dam_break'
    3528
    36 copy_code_files(output_dir,__file__)
     29anuga.copy_code_files(output_dir,__file__)
    3730#start_screen_catcher(output_dir+'_')
    3831
     
    4740
    4841# structured mesh
    49 points, vertices, boundary = rectangular_cross(int(L/dx), int(W/dy), L, W, (0.0, -W/2))
     42points, vertices, boundary = anuga.rectangular_cross(int(L/dx), int(W/dy), L, W, (0.0, -W/2))
    5043
    51 domain = Domain(points, vertices, boundary)
     44domain = anuga.Domain(points, vertices, boundary)
    5245
    5346domain.set_name(output_file)               
     
    9790#------------------------------------------------------------------------------
    9891from math import sin, pi, exp
    99 Br = Reflective_boundary(domain)      # Solid reflective wall
    100 Bt = Transmissive_boundary(domain)    # Continue all values on boundary
    101 Bd = Dirichlet_boundary([1,0.,0.]) # Constant boundary values
     92Br = anuga.Reflective_boundary(domain)      # Solid reflective wall
     93Bt = anuga.Transmissive_boundary(domain)    # Continue all values on boundary
     94Bd = anuga.Dirichlet_boundary([1,0.,0.]) # Constant boundary values
    10295
    10396
Note: See TracChangeset for help on using the changeset viewer.