Changeset 4003 for anuga_core/documentation/user_manual/demos/cairns/runcairns.py

Timestamp:

Nov 17, 2006, 10:24:43 AM (17 years ago)

Author:

sexton

Message:

including data and supporting files for Cairns demo + updating scripts

File:

• anuga_core/documentation/user_manual/demos/cairns/runcairns.py (modified) (5 diffs)

anuga_core/documentation/user_manual/demos/cairns/runcairns.py

@@ -28 +28 @@
 from anuga.abstract_2d_finite_volumes.combine_pts import combine_rectangular_points_files 
 from anuga.geospatial_data.geospatial_data import *
-from anuga.abstract_2d_finite_volumes.util import Screen_Catcher
 
 # Application specific imports
@@ -36 +35 @@
 # Define scenario as either slump or fixed_wave.
 #-------------------------------------------------------------------------------
-scenario = 'slump' # 'fixedwave'
-mkdir (scenario)
-basename = scenario + project.basename
-
-#-------------------------------------------------------------------------------
-# Copy scripts to time stamped output directory and capture screen
-# output to file
-#-------------------------------------------------------------------------------
-
-# creates copy of code in output dir if dir doesn't exist
-if access(project.outputtimedir,F_OK) == 0 :
-    mkdir (project.outputtimedir)
-copy (project.codedirname, project.outputtimedir + project.codename)
-copy (project.codedir + 'runcairns.py', project.outputtimedir + 'runcairns2.py')
-print'output dir', project.outputtimedir
-
-#normal screen output is stored in 
-screen_output_name = project.outputtimedir + "screen_output.txt"
-screen_error_name = project.outputtimedir + "screen_error.txt"
-
-#used to catch screen output to file
-sys.stdout = Screen_Catcher(screen_output_name)
-sys.stderr = Screen_Catcher(screen_error_name)
-
-print 'USER:    ', project.user
+scenario = 'slump' # OR 'fixed_wave'
+if access(scenario,F_OK) == 0:
+    mkdir (scenario)
+basename = scenario + 'source'
 
 #-------------------------------------------------------------------------------
@@ -69 +47 @@
 
 # filenames
-dem_name = project.dem_name
-meshname = project.meshname+'.msh'
+dem_name = 'cairns' 
+meshname = 'cairns.msh'
 
 # creates DEM from asc data
@@ -116 +94 @@
 print domain.statistics()
 
-domain.set_name(basename) #domain.set_name(project.basename)
-domain.set_datadir(scenario) #domain.set_datadir(project.outputtimedir)
+domain.set_name(basename) 
+domain.set_datadir(scenario)
 domain.set_quantities_to_be_stored(['stage', 'xmomentum', 'ymomentum'])
 domain.set_minimum_storable_height(0.01)
@@ -129 +107 @@
 domain.set_quantity('friction', 0.0) 
 domain.set_quantity('elevation', 
-                    filename = project.dem_name + '.pts',
+                    filename = dem_name + '.pts',
                     use_cache = True,
                     verbose = True,