Changeset 1989
- Timestamp:
- Nov 1, 2005, 5:18:08 PM (18 years ago)
- Location:
- production/karratha_2005
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
production/karratha_2005/create_mesh.py
r1868 r1989 14 14 points = [] 15 15 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 17 18 assert zone == refzone 18 19 points.append([easting, northing]) … … 64 65 65 66 66 def create_mesh(bounding_polygon, boundary_tags, resolution, 67 filename = None, interior_regions = None): 67 def create_mesh_from_regions(bounding_polygon, 68 boundary_tags, 69 resolution, 70 filename = None, 71 interior_regions = None): 68 72 """Create mesh from bounding polygon, tags for all segments and resolution. 69 73 -
production/karratha_2005/project.py
r1987 r1989 106 106 107 107 refinery_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 113 point_south = degminsec2decimal_degrees(-20,39,46) 114 point_north = degminsec2decimal_degrees(-20,39,42) 115 point_west = degminsec2decimal_degrees(116,42,0) 116 point_east = degminsec2decimal_degrees(116,42,10) 117 118 d0 = [point_south, point_west] 119 d1 = [point_south, point_east] 120 d2 = [point_north, point_east] 121 d3 = [point_north, point_west] 122 123 point_polygon = [d0, d1, d2, d3] 124 125 126 #Neils areas around interesting points 127 neil1_point1 = [degminsec2decimal_degrees(-20,35,34), 128 degminsec2decimal_degrees(116,45,18)] 129 neil1_point2 = [degminsec2decimal_degrees(-20,36,15), 130 degminsec2decimal_degrees(116,46,18)] 131 neil1_point3 = [degminsec2decimal_degrees(-20,35,9), 132 degminsec2decimal_degrees(116,47,17)] 133 neil1_point4 = [degminsec2decimal_degrees(-20,34,26), 134 degminsec2decimal_degrees(116,46,17)] 135 136 neil1_polygon = [neil1_point1, neil1_point2, neil1_point3, neil1_point4] 137 138 139 140 neil2_point1 = [degminsec2decimal_degrees(-20,39,36), 141 degminsec2decimal_degrees(116,41,33)] 142 neil2_point2 = [degminsec2decimal_degrees(-20,40,10), 143 degminsec2decimal_degrees(116,42,13)] 144 neil2_point3 = [degminsec2decimal_degrees(-20,38,39), 145 degminsec2decimal_degrees(116,43,49)] 146 neil2_point4 = [degminsec2decimal_degrees(-20,38,5), 147 degminsec2decimal_degrees(116,43,9)] 148 149 neil2_polygon = [neil2_point1, neil2_point2, neil2_point3, neil2_point4] -
production/karratha_2005/run_karratha.py
r1920 r1989 83 83 84 84 #Create Triangular Mesh 85 from create_mesh import create_mesh 85 from create_mesh import create_mesh_from_regions 86 86 87 interior_regions = [[project.karratha_polygon, 25000], 88 [project.dampier_polygon, 8000], 89 [project.refinery_polygon, 8000]] 87 interior_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 90 94 91 m = cache(create_mesh ,95 m = cache(create_mesh_from_regions, 92 96 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]}, 94 99 'resolution': 100000, 95 100 'filename': project.meshname + '.msh', … … 107 112 108 113 114 #domain.set_name(project.basename + '_notsunami') 109 115 domain.set_name(project.basename) 110 116 domain.set_datadir(project.outputdir) 111 117 domain.store = True 118 #domain.smooth = False 119 112 120 113 121 domain.quantities_to_be_stored = ['stage', 'xmomentum', 'ymomentum'] … … 145 153 domain.set_boundary( {'back': Br,'side': Bd, 'ocean': Bf} ) #MOST tsunami 146 154 147 #domain.set_boundary( {'back': B r,'side': Bd, 'ocean': Bd} ) #Nothing155 #domain.set_boundary( {'back': Bd,'side': Bd, 'ocean': Bd} ) #Nothing 148 156 149 157
Note: See TracChangeset
for help on using the changeset viewer.