Changeset 8331 for trunk/anuga_work/development/mem_time_tests/scenarios/channelflow/ex1.py

Timestamp:

Feb 10, 2012, 10:34:08 AM (12 years ago)

Author:

pittj

Message:

update for the new experiments

File:

• trunk/anuga_work/development/mem_time_tests/scenarios/channelflow/ex1.py (modified) (3 diffs)

trunk/anuga_work/development/mem_time_tests/scenarios/channelflow/ex1.py

@@ -19 +19 @@
 log._setup = False
 
-log.resource_usage_timing(prefix = 'BeforeSimulation')#get memory usage
+log.timingInfo(msg=('variable1,'+str(length))) #write the variable to be measured to file
+log.timingInfo(msg=('variable2,'+str(lower))) #write the variable to be measured to file
+log.timingInfo(msg=('beforetime,'+str(log.TimeStamp()))) #get the time at the beginning of the simulation
+
+log.resource_usage_timing(prefix = 'beforesimulation')#get memory usage
 #------------------------------------------------------------------------------
 # Setup computational domain
 #------------------------------------------------------------------------------
-points, vertices, boundary = anuga.rectangular_cross(50,150,len1=length, len2=length) # Mesh
+points, vertices, boundary = anuga.rectangular_cross(100,300,len1=length, len2=length) # Mesh
 domain = anuga.Domain(points, vertices, boundary) # Create domain
 domain.set_datadir(scenariodirV) 
 domain.set_name('channel1.sww') # Output name
 
-log.resource_usage_timing(prefix = 'AfterMesh') #get memory usage
+log.resource_usage_timing(prefix = 'aftermesh') #get memory usage
 #--------------------------------------------------------------------------
 # Setup Initial Conditions
 #--------------------------------------------------------------------------
 def topography(x, y):
-    z = -x/10
+    z = 0
     N = len(x)
     
@@ -52 +56 @@
 # Setup boundary conditions
 #------------------------------------------------------------------------------
-Bi = anuga.Dirichlet_boundary([0.5, 1, 0]) # Inflow
-Bo = anuga.Dirichlet_boundary([-100,0,0]) 
+Bi = anuga.Dirichlet_boundary([0.5, 10, 0]) # Inflow
+Bo = anuga.Dirichlet_boundary([-0.5,0,0]) #outflow
 Br = anuga.Reflective_boundary(domain) # Solid reflective wall
 domain.set_boundary({'left': Bi, 'right': Bo, 'top': Br, 'bottom': Br})
@@ -61 +65 @@
 # Evolve system through time
 #------------------------------------------------------------------------------
-for t in domain.evolve(yieldstep=0.2, finaltime=40.0):
+for t in domain.evolve(yieldstep=120, finaltime=3600):
     domain.write_time()
 
 log.resource_usage_timing(prefix='aftersimulation') #get memory usage
+log.timingInfo(msg=('aftertime,'+str(log.TimeStamp()))) #get the time at the end of the simulation