Changeset 1989


Ignore:
Timestamp:
Nov 1, 2005, 5:18:08 PM (18 years ago)
Author:
ole
Message:

Renaming of mesh generation and new refinements

Location:
production/karratha_2005
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • production/karratha_2005/create_mesh.py

    r1868 r1989  
    1414    points = []
    1515    for point in polygon:
    16         zone, easting, northing  = redfearn(point[0], point[1]) #FIXME: Use point.latitude etc
     16        zone, easting, northing  = redfearn(point[0], point[1])
     17        #FIXME: Use point.latitude etc once we have a proper point set
    1718        assert zone == refzone
    1819        points.append([easting, northing])
     
    6465
    6566
    66 def create_mesh(bounding_polygon, boundary_tags, resolution,
    67                 filename = None, interior_regions = None):
     67def create_mesh_from_regions(bounding_polygon,
     68                             boundary_tags,
     69                             resolution,
     70                             filename = None,
     71                             interior_regions = None):
    6872    """Create mesh from bounding polygon, tags for all segments and resolution.
    6973
  • production/karratha_2005/project.py

    r1987 r1989  
    106106
    107107refinery_polygon = [d0, d1, d2, d3]
     108
     109
     110#Interior region around 468899, 7715177:
     111#lat (-20, 39, 44.93753), lon (116, 42, 5.09106)
     112
     113point_south = degminsec2decimal_degrees(-20,39,46)
     114point_north = degminsec2decimal_degrees(-20,39,42)
     115point_west = degminsec2decimal_degrees(116,42,0)
     116point_east = degminsec2decimal_degrees(116,42,10)
     117
     118d0 = [point_south, point_west]
     119d1 = [point_south, point_east]
     120d2 = [point_north, point_east]
     121d3 = [point_north, point_west]   
     122
     123point_polygon = [d0, d1, d2, d3]
     124
     125
     126#Neils areas around interesting points
     127neil1_point1 = [degminsec2decimal_degrees(-20,35,34),
     128                degminsec2decimal_degrees(116,45,18)]
     129neil1_point2 = [degminsec2decimal_degrees(-20,36,15),
     130                degminsec2decimal_degrees(116,46,18)]
     131neil1_point3 = [degminsec2decimal_degrees(-20,35,9),
     132                degminsec2decimal_degrees(116,47,17)]
     133neil1_point4 = [degminsec2decimal_degrees(-20,34,26),
     134                degminsec2decimal_degrees(116,46,17)]
     135
     136neil1_polygon = [neil1_point1, neil1_point2, neil1_point3, neil1_point4]
     137
     138
     139
     140neil2_point1 = [degminsec2decimal_degrees(-20,39,36),
     141                degminsec2decimal_degrees(116,41,33)]
     142neil2_point2 = [degminsec2decimal_degrees(-20,40,10),
     143                degminsec2decimal_degrees(116,42,13)]
     144neil2_point3 = [degminsec2decimal_degrees(-20,38,39),
     145                degminsec2decimal_degrees(116,43,49)]
     146neil2_point4 = [degminsec2decimal_degrees(-20,38,5),
     147                degminsec2decimal_degrees(116,43,9)]
     148
     149neil2_polygon = [neil2_point1, neil2_point2, neil2_point3, neil2_point4]
  • production/karratha_2005/run_karratha.py

    r1920 r1989  
    8383
    8484#Create Triangular Mesh
    85 from create_mesh import create_mesh
     85from create_mesh import create_mesh_from_regions
    8686
    87 interior_regions = [[project.karratha_polygon, 25000],
    88                     [project.dampier_polygon, 8000],
    89                     [project.refinery_polygon, 8000]]
     87interior_regions = [#[project.karratha_polygon, 25000],
     88                    #[project.dampier_polygon, 2000],
     89                    #[project.refinery_polygon, 2000],
     90                    #[project.point_polygon, 2000]]
     91                    [project.neil1_polygon, 1000],
     92                    [project.neil2_polygon, 1000]]
     93                   
    9094
    91 m = cache(create_mesh,
     95m = cache(create_mesh_from_regions,
    9296          project.polygon,
    93           {'boundary_tags': {'back': [7, 8], 'side': [0, 6], 'ocean': [1, 2, 3, 4, 5]},
     97          {'boundary_tags': {'back': [7, 8], 'side': [0, 6],
     98                             'ocean': [1, 2, 3, 4, 5]},
    9499          'resolution': 100000,
    95100          'filename': project.meshname + '.msh',
     
    107112
    108113
     114#domain.set_name(project.basename + '_notsunami')
    109115domain.set_name(project.basename)
    110116domain.set_datadir(project.outputdir)
    111117domain.store = True
     118#domain.smooth = False
     119
    112120
    113121domain.quantities_to_be_stored = ['stage', 'xmomentum', 'ymomentum']
     
    145153domain.set_boundary( {'back': Br,'side': Bd, 'ocean': Bf} ) #MOST tsunami
    146154
    147 #domain.set_boundary( {'back': Br,'side': Bd, 'ocean': Bd} )  #Nothing
     155#domain.set_boundary( {'back': Bd,'side': Bd, 'ocean': Bd} )  #Nothing
    148156
    149157
Note: See TracChangeset for help on using the changeset viewer.