Changeset 6256
- Timestamp:
- Feb 2, 2009, 12:13:18 AM (16 years ago)
- Location:
- anuga_work/production/busselton/standardised_version
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/busselton/standardised_version/project.py
r6254 r6256 256 256 257 257 # Combined all regions, must check that all are included! 258 interior_regions = [[poly_aoi1,res_aoi1],[poly_aoi2,res_aoi2] 259 ,[poly_aos1,res_aos1],[poly_aos2,res_aos2]260 ,[poly_aos3,res_aos3],[poly_aos4,res_aos4]261 ,[poly_sw1,res_sw1]]258 interior_regions = [[poly_aoi1,res_aoi1],[poly_aoi2,res_aoi2], 259 [poly_aos1,res_aos1],[poly_aos2,res_aos2], 260 [poly_aos3,res_aos3],[poly_aos4,res_aos4], 261 [poly_sw1,res_sw1]] 262 262 263 263 -
anuga_work/production/busselton/standardised_version/run_busselton.py
r6255 r6256 15 15 Ole Nielsen, Jane Sexton and Kristy Van Putten - 2008 16 16 """ 17 18 # Note repeated use of domain id 19 # Note naming of fundamental datasets and directories 20 17 21 18 22 #------------------------------------------------------------------------------ … … 41 45 from anuga_parallel.parallel_abstraction import get_processor_name 42 46 from anuga.caching import myhash 43 from anuga.damage_modelling.inundation_damage import add_depth_and_momentum2csv, inundation_damage47 #from anuga.damage_modelling.inundation_damage import add_depth_and_momentum2csv, inundation_damage 44 48 from anuga.fit_interpolate.benchmark_least_squares import mem_usage 45 49 from anuga.utilities.polygon import read_polygon, plot_polygons, polygon_area, is_inside_polygon … … 104 108 105 109 # FIXME(Ole): Introduce create_domain_from_regions: Simpler and caches well 106 print ' start create mesh from regions'110 print 'Start create mesh from regions' 107 111 108 112 create_mesh_from_regions(bounding_polygon, … … 120 124 121 125 domain = Domain(project.meshes_dir_name, use_cache=False, verbose=True) 122 print 'memory usage before del domain',mem_usage()123 124 126 print domain.statistics() 125 print 'triangles',len(domain)126 127 127 128 128 … … 141 141 domain.set_quantity('friction', project.friction) 142 142 143 print 'Start Set quantity',kwargs['elevation_file']143 print 'Start set quantity', kwargs['elevation_file'] 144 144 145 145 domain.set_quantity('elevation', … … 160 160 # Set domain parameters 161 161 #------------------------------------------------------ 162 print 'domain id', id(domain)163 162 domain.set_name(project.scenario_name) 164 domain.set_datadir(kwargs['output_dir']) 163 domain.set_datadir(kwargs['output_dir']) # FIXME: Use project.output_run_time_dir 165 164 domain.set_default_order(2) # Apply second order scheme 166 165 domain.set_minimum_storable_height(0.01) # Don't store anything less than 1cm 167 domain.set_store_vertices_uniquely(False) 168 domain.set_quantities_to_be_stored(['stage', 'xmomentum', 'ymomentum']) 169 domain.tight_slope_limiters = 1 170 print 'domain id', id(domain) 166 #domain.set_store_vertices_uniquely(False) 167 #domain.set_quantities_to_be_stored(['stage', 'xmomentum', 'ymomentum']) 168 #domain.tight_slope_limiters = 1 171 169 172 170 #------------------------------------------------------------------------- … … 174 172 #------------------------------------------------------------------------- 175 173 print 'Available boundary tags', domain.get_boundary_tags() 176 print 'domain id', id(domain)177 174 178 175 boundary_urs_out=project.boundaries_dir_event + sep + project.scenario_name … … 182 179 183 180 print 'Available boundary tags', domain.get_boundary_tags() 184 Bf = Field_boundary(boundary_urs_out+'.sts', # Change from file_boundary185 domain, mean_stage= project.tide,186 time_thinning=1,187 default_boundary=Bd,188 use_cache=True,189 verbose = True,190 boundary_polygon=bounding_polygon)181 Bf = Field_boundary(boundary_urs_out+'.sts', 182 domain, mean_stage= project.tide, 183 time_thinning=1, 184 default_boundary=Bd, 185 use_cache=True, 186 verbose = True, 187 boundary_polygon=bounding_polygon) 191 188 192 189 domain.set_boundary({'back': Br, … … 195 192 196 193 197 print 'finish set boundary'194 print 'Finish set boundary' 198 195 199 196 #---------------------------------------------------------------------------- … … 216 213 217 214 218 kwargs['completed']=str(time.time()-t0)219 215 220 store_parameters(**kwargs) 216 221 217 222 218
Note: See TracChangeset
for help on using the changeset viewer.