Ignore:
Timestamp:
Jun 13, 2008, 10:29:14 AM (17 years ago)
Author:
ole
Message:

John reducing domain extent and coarsening of some regions.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/development/boxingday08/run_boxingday_most.py

    r5393 r5401  
    4040# Application specific imports
    4141import project                 # Definition of file names and polygons
    42 import create_boundary
    4342
    4443#------------------------------------------------------------------------------
     
    5352
    5453timestamp = strftime('%Y%m%d_%H%M%S',localtime())
    55 basename = scenario + 'source'
     54basename = scenario[:4] + '_most'
    5655
    5756
     
    7069# resolutions (maximal area of per triangle) for each polygon
    7170#------------------------------------------------------------------------------
    72 extent_res     = 10000000.0
     71extent_res     = 1000000.0
    7372contour20m_res = 50000.0
    7473island_res     = 5000.0
    7574bay_res        = 2000.0
    76 patong_res     = 100.0
     75patong_res     = 400.0
     76
     77
    7778
    7879# make polygon that contains land that does not affect result.
    7980interior_regions = [[project.patong, patong_res],
    8081                    [project.bay, bay_res],
    81                     [project.contour20m, contour20m_res],
    82                     [project.island_north, island_res],
    83                     [project.island_south, island_res],
    84                     [project.island_south2, island_res]]
     82                    [project.contour20m, contour20m_res]]#,
     83                    #[project.island_north, island_res],
     84                    #[project.island_south, island_res],
     85                    #[project.island_south2, island_res]]
    8586
    8687#for coarse run to test gauges
    87 extent_res = 10000000.0
    88 interior_regions=None
     88#extent_res = 1000000.0 #causes inaccurate stage
     89#interior_regions=None
    8990
    9091from Numeric import arange,allclose
    91 boundary_tags={'ocean': arange(1,85).tolist(), 'otherocean': [0,89], 'land': [86, 87], 'both': [85,88]}
    92 
    93 
    94 #trigs_min = number_mesh_triangles(interior_regions, create_boundary.bounding_polygon,extent_res)
     92boundary_tags={'ocean': arange(0,41).tolist(), 'otherocean': [41,44], 'land': [42,43]}
     93
     94
     95#trigs_min = number_mesh_triangles(interior_regions, project.bounding_polygon,extent_res)
    9596#print 'Minimum number of traingles ', trigs_min
    9697
     
    102103
    103104_ = cache(create_mesh_from_regions,
    104           create_boundary.bounding_polygon,
     105          project.bounding_polygon,
    105106          {'boundary_tags': boundary_tags,
    106107           'maximum_triangle_area': extent_res,
    107108           'filename': meshname,           
    108109           'interior_regions': interior_regions},
    109           verbose = True,
    110           dependencies = ['run_boxingday_most.py']
     110          verbose = True#,
     111          #dependencies = ['run_boxingday_most.py']
    111112          #, evaluate=True
    112113          )
     
    128129print domain.statistics()
    129130
    130 domain.set_name(basename+'friction'+timestamp)
     131domain.set_name(basename+timestamp)
    131132domain.set_datadir(scenario)
    132133domain.set_quantities_to_be_stored(['stage', 'xmomentum', 'ymomentum'])
     
    197198print 'Available boundary tags', domain.get_boundary_tags()
    198199Bf = File_boundary(project.boundary_most_out+'.sww',
    199                     domain, time_thinning=1, verbose = True)
     200                    domain, time_thinning=1, use_cache=True,verbose = True)
    200201Br = Reflective_boundary(domain)
    201202Bd = Dirichlet_boundary([tide,0.0,0.0])
     
    219220# start when boundary stage first becomes non-zero.
    220221
    221 for t in domain.evolve(yieldstep = 5.0, finaltime = 18000.0,
     222for t in domain.evolve(yieldstep = 20.0, finaltime = 18000.0,
    222223                       skip_initial_step = False):
    223224    if allclose(t,10800.0):
Note: See TracChangeset for help on using the changeset viewer.