Ignore:
Timestamp:
Jul 25, 2008, 2:29:44 PM (16 years ago)
Author:
bridgette
Message:

update script for bridgettes use

Location:
anuga_work/development/landslide_4
Files:
3 added
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/development/landslide_4/benchmark_problem4.py

    r5568 r5572  
    1515from os import mkdir, access, F_OK, sep
    1616import sys
     17import project
    1718
    1819# Related major packages
    1920from anuga.shallow_water import Domain, Reflective_boundary, Dirichlet_boundary
    2021from anuga.abstract_2d_finite_volumes.util import start_screen_catcher, copy_code_files
     22from anuga.shallow_water.data_manager import start_screen_catcher, copy_code_files,store_parameters
     23
     24myid=0
     25numprocs=1
     26print 'output_dir',project.output_run_time_dir
     27copy_code_files(project.output_run_time_dir,'benchmark_problem4.py',
     28             'project.py' )
     29start_screen_catcher(project.output_run_time_dir, myid, numprocs)
    2130
    2231#----------------------------------------------------------------------------
     
    3544bounding_polygon = [[0,0],[length,0],[length,width],[0,width]]
    3645meshname = 'slide.msh'
    37 remainder_res = 1.
     46
    3847
    3948from caching import cache
     
    4251           {'boundary_tags': {'bottom': [0], 'right': [1],
    4352                              'top': [2], 'left': [3]},
    44            'maximum_triangle_area': remainder_res,
     53           'maximum_triangle_area': project.remainder_res,
    4554           'filename': meshname},
    4655           #'interior_regions': interior_regions},
     
    5766print domain.statistics()
    5867
    59 domain.set_name('slide')
     68domain.set_datadir(project.output_run_time_dir)
     69domain.set_name(project.sww_name)
    6070domain.set_quantities_to_be_stored(['stage', 'xmomentum', 'ymomentum'])
    6171domain.set_minimum_storable_height(0.01)
     
    93103tsunami_source = slide_tsunami(length=0.91,
    94104                               width=0.61,
    95                                depth=.1, # delta from Synolakis paper
     105                               depth=.01, # delta from Synolakis paper
    96106                               slope=degrees(atan(0.5)),
    97107                               thickness=a,
     
    122132
    123133t0 = time.time()
    124 
    125 for t in domain.evolve(yieldstep = 1, finaltime = 500):
     134finaltime = 500
     135for t in domain.evolve(yieldstep = 1, finaltime = project.finaltime):
    126136    domain.write_time()
    127137    stagestep = domain.get_quantity('stage')
Note: See TracChangeset for help on using the changeset viewer.