Changeset 2441
- Timestamp:
- Feb 22, 2006, 6:08:30 PM (19 years ago)
- Location:
- production/onslow_2006
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
production/onslow_2006/project.py
r2436 r2441 67 67 polyAll = [d0, d1, d2, d3, d4, d5] 68 68 69 #Interior region - Onslow - take 1 70 i0 = [310000, 7610000] 71 i1 = [300000, 7610000] 72 i2 = [300000, 7600000] 73 i3 = [310000, 7600000] 69 #Interior region - Onslow town 74 70 75 polyInterior = [i0, i1, i2, i3] 71 i0 = [304000, 7608000] 72 i1 = [302000, 7605000] 73 i2 = [303000, 7602000] 74 i3 = [305000, 7601000] 75 i4 = [309000, 7603000] 76 i5 = [307000, 7606500] 77 78 poly_onslow = [i0, i1, i2, i3, i4, i5] 79 80 #Thevenard Island 81 j0 = [294000, 7629000] 82 j1 = [285000, 7625000] 83 j2 = [294000, 7621000] 84 j3 = [299000, 7625000] 85 86 poly_thevenard = [j0, j1, j2, j3] 87 88 # Direction Is 89 k0 = [309000, 7619000] 90 k1 = [304000, 7619000] 91 k2 = [304000, 7616500] 92 k3 = [309000, 7616500] 93 94 poly_direction = [k0, k1, k2, k3] 76 95 77 96 #!slump_origin = [385000.0, 6255000.0] #Absolute UTM -
production/onslow_2006/run_onslow.py
r2437 r2441 76 76 # original 77 77 interior_res = 5000 78 interior_regions = [[project.polyInterior, interior_res]] 78 interior_regions = [[project.poly_onslow, interior_res], 79 [project.poly_thevenard, interior_res], 80 [project.poly_direction, interior_res]] 79 81 80 82 #FIXME: Fix caching of this one once the mesh_interface is ready … … 150 152 f=lambda t: [(20<t<200)*6, 0, 0]) 151 153 152 domain.set_boundary( {'top': B w, 'topleft': Br,154 domain.set_boundary( {'top': Br, 'topleft': Br, 153 155 'left': Br, 'bottom': Br, 154 156 'bottomright': Br, 'topright': Br} ) … … 162 164 t0 = time.time() 163 165 164 for t in domain.evolve(yieldstep = 10, finaltime = 300):166 for t in domain.evolve(yieldstep = 10, finaltime = 20): 165 167 domain.write_time() 166 168 domain.write_boundary_statistics(tags = 'top')
Note: See TracChangeset
for help on using the changeset viewer.