Changeset 4049 for anuga_work/production/dampier_2006/run_dampier.py
- Timestamp:
- Dec 4, 2006, 11:16:58 AM (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
anuga_work/production/dampier_2006/run_dampier.py
r3998 r4049 29 29 from anuga.shallow_water import File_boundary 30 30 from anuga.shallow_water import Reflective_boundary 31 from Numeric import allclose 31 32 32 33 from anuga.pmesh.mesh_interface import create_mesh_from_regions … … 84 85 print 'start create mesh from regions' 85 86 interior_regions = [#[project.karratha_polygon, 25000], 86 #[project.poly_coast, 10000],87 #[project.poly_pipeline, 2000],88 #[project.poly_facility, 500]]89 [project.poly_interior, 1000]]87 [project.poly_coast, 10000], 88 [project.poly_pipeline, 2000], 89 [project.poly_facility, 500]] 90 # [project.poly_interior, 1000]] 90 91 91 from anuga.utilities.polygon import plot_polygons92 figname = project.output_run_time_dir + 'poly_pic'93 plot_polygons([interior_regions,94 project.bounding_polygon],95 figname,96 verbose = True)97 # if access(project.meshes_time_dir,F_OK) == 0:98 # mkdir(project.meshes_time_dir)99 if access(project.meshes_dir,F_OK) == 0:100 mkdir(project.meshes_dir)101 print 'start create mesh from regions'102 interior_regions = [#[project.karratha_polygon, 25000],103 # [project.cipma_polygon, 1000],104 # [project.poly_pipeline, 5000],105 # [project.poly_facility, 500]]106 [project.poly_interior, 1000]]107 # meshes_dir_name = project.meshes_dir_name + '.msh'108 92 109 93 create_mesh_from_regions(project.bounding_polygon, … … 147 131 from caching import cache 148 132 if myid == 0: 149 cache( urs2sww,133 cache(ferret2sww, 150 134 (boundaries_in_dir_name, 151 135 # boundaries_time_dir_name), … … 160 144 # 'minlon': project.west, 161 145 # 'maxlon': project.east, 162 'mint': 0, 'maxt': 35 000,146 'mint': 0, 'maxt': 35100, 163 147 'origin': domain.geo_reference.get_origin(), 164 148 'mean_stage': project.tide, … … 173 157 # Setup initial conditions 174 158 #------------------------------------------------------------------------- 175 print 'Setup initial conditions' 176 177 domain.set_quantity('stage', tide) 178 domain.set_quantity('friction', 0.0) 179 #combined_time_dir_name = project.topographies_dir+build_time+project.combined_name 180 print 'Start Set quantity' 181 182 domain.set_quantity('elevation', 159 if myid == 0: 160 161 print 'Setup initial conditions' 162 163 domain.set_quantity('stage', tide) 164 domain.set_quantity('friction', 0.0) 165 #combined_time_dir_name = project.topographies_dir+build_time+project.combined_name 166 print 'Start Set quantity' 167 168 169 domain.set_quantity('elevation', 183 170 filename = project.topographies_dir + build_time + sep + project.combined_final_name + '.pts', 184 171 use_cache = True, 185 172 verbose = True, 186 173 alpha = 0.1) 187 #domain.set_quantity('elevation', -50)188 print 'Finished Set quantity'189 174 #domain.set_quantity('elevation', -50) 175 print 'Finished Set quantity' 176 barrier() 190 177 191 178 … … 247 234 248 235 for t in domain.evolve(yieldstep = 120, finaltime = 9000): 249 #for t in domain.evolve(yieldstep = 120, finaltime = 130):250 236 domain.write_time() 251 237 domain.write_boundary_statistics(tags = 'ocean') 252 print 'time: ',time.time()238 # print 'time: ',time.time() 253 239 254 240 if allclose(t, 6000): 255 241 domain.set_quantity('xmomentum', 0) 256 242 domain.set_quantity('ymomentum', 0) 257 import sys; sys.exit()243 #import sys; sys.exit() 258 244 259 245 for t in domain.evolve(yieldstep = 60, finaltime = 28800
Note: See TracChangeset
for help on using the changeset viewer.