Ignore:
Timestamp:
Nov 4, 2010, 5:23:35 PM (14 years ago)
Author:
steve
Message:

Pulling together basic validate tests for anuga

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/anuga_validation/convergence_study/wave.py

    r5442 r8067  
    1111
    1212import sys
    13 from anuga.abstract_2d_finite_volumes.mesh_factory import rectangular_cross
    14 from anuga.shallow_water import Domain
    15 from anuga.shallow_water import Reflective_boundary
    16 from anuga.shallow_water import Dirichlet_boundary
    17 from anuga.shallow_water import Time_boundary
    18 from anuga.shallow_water import Transmissive_boundary
    19 from anuga.shallow_water import Transmissive_Momentum_Set_Stage_boundary
    20 from anuga.shallow_water.data_manager import start_screen_catcher, copy_code_files
     13from anuga import rectangular_cross
     14from anuga import Domain
     15from anuga import Reflective_boundary
     16from anuga import Dirichlet_boundary
     17from anuga import Time_boundary
     18from anuga import Transmissive_boundary
     19from anuga import Transmissive_momentum_set_stage_boundary
     20#from anuga.data_manager import start_screen_catcher, copy_code_files
    2121
    2222from math import cos
    23 from Numeric import zeros, Float
     23import numpy as num
    2424from time import localtime, strftime, gmtime
    2525from os import sep
     
    3636output_file = 'wave'
    3737
    38 copy_code_files(output_dir,__file__)
     38#copy_code_files(output_dir,__file__)
    3939#start_screen_catcher(output_dir+sep)
    4040
    41 interactive_visualisation = False
     41interactive_visualisation = True
    4242
    4343#------------------------------------------------------------------------------
     
    6060# Setup Algorithm
    6161#------------------------------------------------------------------------------
    62 domain.set_timestepping_method('rk2')
     62
     63domain.set_timestepping_method(2)
    6364domain.set_default_order(2)
    6465
    6566print domain.get_timestepping_method()
    6667
    67 domain.use_edge_limiter = True
    68 domain.tight_slope_limiters = False
    69 domain.use_centroid_velocities = False
     68#domain.use_edge_limiter = True
     69#domain.tight_slope_limiters = False
     70#domain.use_centroid_velocities = False
     71
     72domain.set_beta(1.9)
    7073
    7174domain.CFL = 1.0
    7275
    73 domain.beta_w      = 1.0
    74 domain.beta_w_dry  = 0.0
    75 domain.beta_uh     = 1.0
    76 domain.beta_uh_dry = 0.0
    77 domain.beta_vh     = 1.0
    78 domain.beta_vh_dry = 0.0
     76#domain.beta_w      = 1.0
     77#domain.beta_w_dry  = 0.0
     78#domain.beta_uh     = 1.0
     79#domain.beta_uh_dry = 0.0
     80#domain.beta_vh     = 1.0
     81#domain.beta_vh_dry = 0.0
    7982
    8083
Note: See TracChangeset for help on using the changeset viewer.