Changeset 8331 for trunk/anuga_work/development/mem_time_tests/parameters/timelen/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/parameters/timelen/ex1.py (modified) (5 diffs)

trunk/anuga_work/development/mem_time_tests/parameters/timelen/ex1.py

@@ -11 +11 @@
 #set up variables to store the results and the log files of each experiment
 f = sys.argv[1]
+g = float(f)/100.0
 home = os.getenv('INUNDATIONHOME')
 scenariodirV = add_directories(home, ["data","mem_time_test", "parameters",
@@ -17 +18 @@
 log._setup = False
 
-log.resource_usage_timing(prefix = 'BeforeSimulation') #get memory usage here
+log.timingInfo(msg=('variable1,'+str(f))) #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 here
 
 #------------------------------------------------------------------------------
@@ -27 +32 @@
 domain.set_datadir(scenariodirV)  
 
-log.resource_usage_timing(prefix = 'AfterMesh') #get memory usage here
+log.resource_usage_timing(prefix = 'aftermesh') #get memory usage here
 #------------------------------------------------------------------------------
 # Setup initial conditions
 #------------------------------------------------------------------------------
 def topography(x, y):
-    return -x/10 # linear bed slope
+    return 0 # linear bed slope
 domain.set_quantity('elevation', topography) # Use function for elevation
 domain.set_quantity('friction', 0.01) # Constant friction
-domain.set_quantity('stage',expression='elevation')
+domain.set_quantity('stage', 10)
 
 log.resource_usage_timing(prefix='afterinitialconditions')#get memory usage here
@@ -42 +47 @@
 # Setup boundary conditions
 #------------------------------------------------------------------------------
-Bi = anuga.Dirichlet_boundary([0.4, 0, 0]) # Inflow
 Br = anuga.Reflective_boundary(domain) # Solid reflective wall
-domain.set_boundary({'left': Bi, 'right': Br, 'top': Br, 'bottom': Br})
+domain.set_boundary({'left': Br, 'right': Br, 'top': Br, 'bottom': Br})
 
 log.resource_usage_timing(prefix='afterboundary')#get memory usage here
@@ -51 +55 @@
 # Evolve system through time
 #------------------------------------------------------------------------------
-for t in domain.evolve(yieldstep=0.2, finaltime=f):
+for t in domain.evolve(yieldstep=g, finaltime=f):
     print domain.timestepping_statistics()
 
 
 log.resource_usage_timing(prefix='aftersimulation') #get memory usage here
-
+log.timingInfo(msg=('aftertime,'+str(log.TimeStamp()))) #get the time at the end of the simulation