Ignore:
Timestamp:
Jun 13, 2008, 10:29:14 AM (16 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_polyline.py

    r5393 r5401  
    4040# Application specific imports
    4141import project                 # Definition of file names and polygons
    42 import create_boundary
    4342
    4443
     
    5453
    5554timestamp = strftime('%Y%m%d_%H%M%S',localtime())
    56 basename = scenario + '_polyline'
     55basename = scenario[:4] + '_polyline'
    5756
    5857
     
    7170# resolutions (maximal area of per triangle) for each polygon
    7271#------------------------------------------------------------------------------
    73 extent_res     = 10000000.0
     72extent_res     = 1000000.0
    7473contour20m_res = 50000.0
    7574island_res     = 5000.0
    7675bay_res        = 2000.0
    77 patong_res     = 100.0
     76patong_res     = 400.0
    7877
    7978
     
    8281interior_regions = [[project.patong, patong_res],
    8382                    [project.bay, bay_res],
    84                     [project.contour20m, contour20m_res],
    85                     [project.island_north, island_res],
    86                     [project.island_south, island_res],
    87                     [project.island_south2, island_res]]
     83                    [project.contour20m, contour20m_res]]#,
     84                    #[project.island_north, island_res],
     85                    #[project.island_south, island_res],
     86                    #[project.island_south2, island_res]]
    8887
    8988
    9089#for coarse run to test gauges
    91 extent_res = 10000000.0
    92 interior_regions=None
     90#extent_res = 10000000.0
     91#interior_regions=None
    9392
    9493from Numeric import arange,allclose
    95 boundary_tags={'ocean': arange(1,85).tolist(), 'otherocean': [0,89], 'land': [86, 87], 'both': [85,88]}
    96 
    97 #trigs_min = number_mesh_triangles(interior_regions, create_boundary.bounding_polygon,extent_res)
     94boundary_tags={'ocean': arange(0,41).tolist(), 'otherocean': [41,44], 'land': [42,43]}
     95
     96#trigs_min = number_mesh_triangles(interior_regions, project.bounding_polygon,extent_res)
    9897#print 'Minimum number of traingles ', trigs_min
    9998
     
    105104
    106105_ = cache(create_mesh_from_regions,
    107           create_boundary.bounding_polygon,
     106          project.bounding_polygon,
    108107          {'boundary_tags': boundary_tags,
    109108           'maximum_triangle_area': extent_res,
     
    111110           'interior_regions': interior_regions},
    112111          verbose = True,
    113           dependencies = ['create_boundary.py']
     112          dependencies = ['project.py']
    114113          #, evaluate=True
    115114          )
     
    130129print domain.statistics()
    131130
    132 domain.set_name(basename+'friction'+timestamp)
     131domain.set_name(basename+timestamp)
    133132#domain.set_name(basename)
    134133domain.set_datadir(scenario)
     
    168167#------------------------------------------------------------------------------
    169168boundary_urs_in='data/boxing'
    170 boundary_urs_out=create_boundary.base_name
     169boundary_urs_out=project.base_name
    171170
    172171print 'Available boundary tags', domain.get_boundary_tags()
    173172Bf = File_boundary(boundary_urs_out+'.sts',
    174173                   domain, time_thinning=1,
    175                    use_cache=False,verbose = True,boundary_polygon=create_boundary.bounding_polygon)
     174                   use_cache=True,verbose = True,boundary_polygon=project.bounding_polygon)
    176175
    177176Br = Reflective_boundary(domain)
     
    196195# start when boundary stage first becomes non-zero.
    197196
    198 for t in domain.evolve(yieldstep = 5.0, finaltime = 18000.0,
     197for t in domain.evolve(yieldstep = 20.0, finaltime = 18000.0,
    199198                       skip_initial_step = False):
    200199    if allclose(t,10800.0):
Note: See TracChangeset for help on using the changeset viewer.