Ignore:
Timestamp:
Jun 13, 2008, 10:29:14 AM (15 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_urs.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 + '_urs'
     54basename = scenario[:4] + '_urs'
    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
    7776
    7877# make polygon that contains land that does not affect result.
     
    8079interior_regions = [[project.patong, patong_res],
    8180                    [project.bay, bay_res],
    82                     [project.contour20m, contour20m_res],
    83                     [project.island_north, island_res],
    84                     [project.island_south, island_res],
    85                     [project.island_south2, island_res]]
     81                    [project.contour20m, contour20m_res]]#,
     82                    #[project.island_north, island_res],
     83                    #[project.island_south, island_res],
     84                    #[project.island_south2, island_res]]
    8685
    8786#for coarse run to test gauges
    88 extent_res = 10000000.0
    89 interior_regions=None
     87#extent_res = 10000000.0
     88#interior_regions=None
    9089
    9190from Numeric import arange,allclose
    92 print 'Beware not all boundary segments that are tagged as ocean in polyline are tagged here. missing are 79,80,81,82,83,84 which are set to other ocean'
    93 boundary_tags={'ocean': arange(1,79).tolist(), 'otherocean': [0,89,84,83,82,81,80,79], 'land': [86, 87], 'both': [85,88]}
    94 
    95 
    96 #trigs_min = number_mesh_triangles(interior_regions, create_boundary.bounding_polygon,extent_res)
     91boundary_tags={'ocean': arange(0,41).tolist(), 'otherocean': [41,44], 'land': [42,43]}
     92
     93
     94#trigs_min = number_mesh_triangles(interior_regions, project.bounding_polygon,extent_res)
    9795#print 'Minimum number of traingles ', trigs_min
    9896
     
    104102
    105103_ = cache(create_mesh_from_regions,
    106           create_boundary.bounding_polygon,
     104          project.bounding_polygon,
    107105          {'boundary_tags': boundary_tags,
    108106           'maximum_triangle_area': extent_res,
     
    110108           'interior_regions': interior_regions},
    111109          verbose = True,
    112           #dependencies = ['create_boundary.py']
     110          #dependencies = ['project.py']
    113111          dependencies = ['run_boxingday_urs.py']
    114112          #, evaluate=True
     
    131129print domain.statistics()
    132130
    133 domain.set_name(basename+'friction'+timestamp)
     131domain.set_name(basename+timestamp)
    134132#domain.set_name(basename)
    135133domain.set_datadir(scenario)
     
    208206# start when boundary stage first becomes non-zero.
    209207
    210 for t in domain.evolve(yieldstep = 5.0, finaltime = 18000.0,
     208for t in domain.evolve(yieldstep = 20.0, finaltime = 18000.0,
    211209                       skip_initial_step = False):
    212210    if allclose(t,10800.0):
Note: See TracChangeset for help on using the changeset viewer.