- Timestamp:
- Jun 20, 2005, 4:53:05 PM (19 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
inundation/ga/storm_surge/parallel/run_parallel_advection.py
r1471 r1520 22 22 M = 50 23 23 24 points, vertices, boundary, full_send_dict, ghost_recv_dict = parallel_rectangular(N, M, len1=1.0) 24 points, vertices, boundary, full_send_dict, ghost_recv_dict = \ 25 parallel_rectangular(N, M, len1_g=1.0) 25 26 26 27 #Create advection domain with direction (1,-1) … … 31 32 32 33 #turn on the visualisation 33 domain.initialise_visualiser() 34 rect = [0.0, 0.0, 1.0, 1.0] 35 domain.initialise_visualiser(rect=rect) 34 36 35 37 #Boundaries … … 55 57 return self.h*((x>self.x0)&(x<self.x1)) 56 58 59 domain.set_quantity('stage', Set_Stage(0.2,0.4,1.0)) 60 57 61 if myid == 0: 58 domain.set_quantity('stage', Set_Stage(0.2,0.4,1.0))59 62 import time 60 63 t0 = time.time() 61 64 62 65 #Check that the boundary value gets propagated to all elements 63 for t in domain.evolve(yieldstep = 0. 5, finaltime = 3.0):66 for t in domain.evolve(yieldstep = 0.1, finaltime = 3.0): 64 67 if myid == 0: 65 68 domain.write_time()
Note: See TracChangeset
for help on using the changeset viewer.