Changeset 5751 for anuga_work/production/geraldton/run_geraldton.py
- Timestamp:
- Sep 9, 2008, 2:14:12 PM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/geraldton/run_geraldton.py
r5669 r5751 78 78 # Domain definitions 79 79 #----------------------------------------------------------------------- 80 ## 81 ## # Read in boundary from ordered sts file 82 ## urs_bounding_polygon=create_sts_boundary(os.path.join(project.boundaries_dir,project.scenario_name)) 83 ## 84 ## # Reading the landward defined points, this incorporates the original clipping 85 ## # polygon minus the 100m contour 86 ## landward_bounding_polygon = read_polygon(project.boundaries_dir+'landward_boundary.txt') 87 ## 88 ## # Combine sts polyline with landward points 89 ## bounding_polygon = urs_bounding_polygon + landward_bounding_polygon 90 ## 91 ## # counting segments 92 ## N = len(urs_bounding_polygon)-1 93 ## boundary_tags={'back': [N+2,N+3], 'side': [N,N+1,N+4],'ocean': range(N)} 94 95 bounding_polygon = project.poly_all 96 97 98 80 81 # Read in boundary from ordered sts file 82 urs_bounding_polygon=create_sts_boundary(os.path.join(project.boundaries_dir,project.scenario_name)) 83 84 # Reading the landward defined points, this incorporates the original clipping 85 # polygon minus the 100m contour 86 landward_bounding_polygon = read_polygon(project.boundaries_dir+'landward_boundary_polygon.txt') 87 88 # Combine sts polyline with landward points 89 bounding_polygon = urs_bounding_polygon + landward_bounding_polygon 90 91 # counting segments 92 N = len(urs_bounding_polygon)-1 93 boundary_tags={'back': [N+2,N+3], 'side': [N,N+1,N+4],'ocean': range(N)} 94 95 print 'boundary tags',boundary_tags 96 99 97 #-------------------------------------------------------------------------- 100 98 # Create the triangular mesh based on overall clipping polygon with a … … 111 109 112 110 create_mesh_from_regions(bounding_polygon, 113 boundary_tags= project.boundary_tags,111 boundary_tags=boundary_tags, 114 112 maximum_triangle_area=project.res_poly_all, 115 113 interior_regions=project.interior_regions, … … 208 206 209 207 boundary_urs_out=project.boundaries_dir_name 208 209 Br = Reflective_boundary(domain) 210 Bd = Dirichlet_boundary([kwargs['tide'],0,0]) 210 211 211 212 print 'Available boundary tags', domain.get_boundary_tags() … … 213 214 domain, mean_stage= project.tide, 214 215 time_thinning=1, 216 default_boundary=Bd, 215 217 use_cache=True, 216 218 verbose = True, 217 219 boundary_polygon=bounding_polygon) 218 220 219 220 Br = Reflective_boundary(domain)221 Bd = Dirichlet_boundary([kwargs['tide'],0,0])222 223 ## Bf = Field_boundary(kwargs['boundary_file'],224 ## domain, time_thinning=kwargs['time_thinning'], mean_stage=kwargs['tide'],225 ## use_cache=False, verbose=True)226 227 221 domain.set_boundary({'back': Br, 228 222 'side': Bd, 229 'ocean': B d}) # change baxk to Bf when running properly223 'ocean': Bf}) 230 224 231 225 kwargs['input_start_time']=domain.starttime … … 290 284 run_model(**kwargs) 291 285 292 if myid==0: 293 export_model(**kwargs) 294 #barrier 286 #barrier
Note: See TracChangeset
for help on using the changeset viewer.