Ignore:
Timestamp:
Apr 9, 2008, 9:09:18 AM (16 years ago)
Author:
herve
Message:

run scripts for tonga

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/tonga/tongatapu.py

    r5139 r5194  
    2424
    2525# Related major packages
     26from anuga.shallow_water import Time_boundary
    2627from anuga.shallow_water import Domain
    2728from anuga.shallow_water import Dirichlet_boundary
     
    3132from Numeric import allclose
    3233from anuga.shallow_water.data_manager import export_grid
    33 
     34from anuga.shallow_water.data_manager import convert_dem_from_ascii2netcdf, dem2pts
    3435from anuga.pmesh.mesh_interface import create_mesh_from_regions
    3536from anuga.shallow_water.data_manager import start_screen_catcher, copy_code_files,store_parameters
     
    8182    # resolutions (maximal area of per triangle) for each polygon
    8283    #--------------------------------------------------------------------------
    83 
    84     #IMPORTANT don't cache create_mesh_from_region and Domain(mesh....) as it
    85     # causes problems with the ability to cache set quantity which takes alot of times
    8684   
    8785    if myid == 0:
     
    9896    barrier()
    9997
     98    scenario='fixed_wave'
     99   
    100100    #-------------------------------------------------------------------------
    101101    # Setup computational domain
     
    133133        domain.set_quantity('elevation',
    134134                            filename = kwargs['bathy_file'],
    135                             use_cache = True,
     135                            use_cache = False,
    136136                            verbose = True,
    137137                            alpha = kwargs['alpha'])
     
    170170        tsunami_source = slide_tsunami(length=35000.0,
    171171                                       depth=project.slide_depth,
    172                                        slope=6.0,
    173                                        thickness=500.0,
     172                                       slope=10.0,
     173                                       thickness=800.0,
    174174                                       x0=project.slide_origin[0],
    175175                                       y0=project.slide_origin[1],
     
    192192    # FIXME (Ole): Change this to Transmissive Momentum as
    193193    # suggested by Rajaraman (ticket:208)
     194        from math import sin, pi
    194195        Bw = Time_boundary(domain = domain,
    195                            f=lambda t: [(60<t<3660)*100, 0, 0])
    196         domain.set_boundary({'ocean_east': Bw,
    197                              'ocean_southeast': Bd,
    198                              'land': Bd,
    199                              'ocean_north': Bd,
    200                              'ocean_west':})
     196                           f=lambda t: [5.0*sin(t*pi*1/15), 0, 0])
     197        domain.set_boundary({'ocean_east':Bw,
     198                             'ocean_north':Bd,
     199                             'ocean_west':Bd,
     200                             'land':Bd,
     201                             'ocean_southeast':Bd})
    201202
    202203    # boundary conditions for slide scenario
    203204    if scenario == 'slide':
    204         domain.set_boundary({'ocean_east': Bd,
    205                              'ocean_southeast': Bd,
    206                              'land': Bd,
    207                              'ocean_north': Bd,
    208                              'ocean_west':})
     205        domain.set_boundary({'ocean_east':Bd,
     206                             'ocean_north':Bd,
     207                             'ocean_west':Bd,
     208                             'land':Bd,
     209                             'ocean_southeast':Bd})
    209210
    210211    #----------------------------------------------------------------------------
     
    213214    t0 = time.time()
    214215
    215     for t in domain.evolve(yieldstep = 240, finaltime = kwargs['finaltime']):
     216    for t in domain.evolve(yieldstep = 1, finaltime = kwargs['finaltime']):
    216217        domain.write_time()
    217218        domain.write_boundary_statistics(tags = 'ocean_east')     
     
    240241    print 'memory usage after del domain',mem_usage()
    241242   
    242     swwfile = kwargs['output_dir']+kwargs['scenario_name']
     243    swwfile = kwargs['output_dir']+kwargs['aa_scenario_name']
    243244    print'swwfile',swwfile
    244245   
    245     export_grid(swwfile, extra_name_out = 'town',
     246    export_grid(swwfile, extra_name_out = 'tonga',
    246247            quantities = ['speed','depth','elevation','stage'], # '(xmomentum**2 + ymomentum**2)**0.5' defaults to elevation
    247248            #quantities = ['speed','depth'], # '(xmomentum**2 + ymomentum**2)**0.5' defaults to elevation
    248249            timestep = None,
    249250            reduction = max,
    250             cellsize = kwargs['export_cellsize'],
     251            cellsize = 25,
    251252            NODATA_value = -1E-030,
    252253            easting_min = project.eastingmin,
     
    271272    kwargs['midtime']=project.midtime
    272273    kwargs['finaltime']=project.finaltime
    273    
    274274    kwargs['output_dir']=project.output_run_time_dir
    275275    kwargs['bathy_file']=project.combined_dir_name+'.txt'
Note: See TracChangeset for help on using the changeset viewer.