Changeset 1392
- Timestamp:
- May 14, 2005, 5:57:37 PM (20 years ago)
- Location:
- inundation/ga/storm_surge
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
inundation/ga/storm_surge/pyvolution/mesh.py
r1387 r1392 130 130 c = sqrt((x2-x0)**2+(y2-y0)**2) 131 131 132 self.radii[i]= self.areas[i]/(a+b+c)132 self.radii[i]=2.0*self.areas[i]/(a+b+c) 133 133 134 134 -
inundation/ga/storm_surge/pyvolution/netherlands.py
r1387 r1392 65 65 66 66 #Wall 67 if x[i] > 0.995:68 z[i] = -x[i]/20+0.367 #if x[i] > 0.995: 68 # z[i] = -x[i]/20+0.3 69 69 70 70 return z/2 … … 80 80 N = 130 #size = 33800 81 81 N = 600 #Size = 720000 82 N = 5083 M = 40 82 N = 30 83 84 84 85 85 #N = 15 … … 87 87 print 'Creating domain' 88 88 #Create basic mesh 89 points, vertices, boundary = rectangular(N, M)89 points, vertices, boundary = rectangular(N, N) 90 90 91 91 #Create shallow water domain … … 156 156 t0 = time.time() 157 157 158 for t in domain.evolve(yieldstep = 0.01, finaltime = 5.0):158 for t in domain.evolve(yieldstep = 1.0, finaltime = 100.0): 159 159 domain.write_time() 160 160 #domain.visualiser.scale_z = 1.0 -
inundation/ga/storm_surge/pyvolution/test_mesh.py
r1366 r1392 370 370 [5, 2, 3]] 371 371 try: 372 mesh = Mesh(points, triangles) 372 mesh = Mesh(points, triangles) 373 373 except: 374 374 pass 375 375 else: 376 376 raise "triangle edge duplicates not caught" 377 377 378 378 def test_rectangular_mesh_basic(self): 379 379 M=1 … … 805 805 806 806 807
Note: See TracChangeset
for help on using the changeset viewer.