Changeset 2979 for development

Timestamp:

May 25, 2006, 5:30:42 PM (19 years ago)

Author:

ole

Message:

Arbitrary working dir

Location:

development/stochastic_study

Files:

• create_mesh.py (modified) (1 diff)
• create_realisations.py (modified) (2 diffs)
• extract_timeseries.py (modified) (1 diff)
• plot_spread.py (modified) (1 diff)
• project.py (modified) (3 diffs)
• run_model.py (modified) (5 diffs)

development/stochastic_study/create_mesh.py

@@ -193 +193 @@
 
     import project
-    m.export_mesh_file(project.mesh_filename)
+    m.export_mesh_file(project.working_dir + project.mesh_filename)
     

development/stochastic_study/create_realisations.py

@@ -15 +15 @@
 
 bathymetry_filename = project.bathymetry_filename
-mesh_filename = project.mesh_filename
+mesh_filename = project.working_dir + project.mesh_filename
 
 
@@ -80 +80 @@
               %block_number
 
-        filename = project.basename + '_block%d.pck'\
-                   %block_number
+        filename = project.working_dir + project.basename +\
+                   '_block%d.pck' %block_number
 
         block_number += 1    

development/stochastic_study/extract_timeseries.py

@@ -62 +62 @@
 
 #Read model output
-filename = project.basename + '_original.sww'
+#filename = project.basename + '_original.sww'
+filename = project.basename + '.sww'
 
 f = cache(file_function, filename,

development/stochastic_study/plot_spread.py

@@ -24 +24 @@
 
 j = 0 # Count realisations
-for filename in os.listdir('.'):
+for filename in os.listdir(project.working_dir):
     if filename.startswith(project.basename) and filename.endswith(gauge):
         if j < data.shape[1]:        

development/stochastic_study/project.py

@@ -2 +2 @@
 """
 
+# Working dir
+working_dir = '/mnt/store/stochastic/'
+import os
+try:
+    os.stat(working_dir)
+except:
+    working_dir = '.' + os.sep
+    
 # Inputs
 boundary_filename = 'input_wave.tms'
@@ -12 +20 @@
 
 # Gauges (3 timeseries (Ch 5-7-9))
-gauges = [[4.521, 1.196],  [4.521, 1.696],  [4.521, 2.196]] 
+gauges = [[4.521, 1.196],  [4.521, 1.696],  [4.521, 2.196]]
 gauge_names = ['ch5', 'ch7', 'ch9']
 
@@ -19 +27 @@
 
 # Stats (Suresh ?)
-number_of_realisations = 9
+number_of_realisations = 2
 #std_dev = 0.0026  #Range is 26.035 cm
-std_dev = 0.0013  #Range is 26.035 cm
+#std_dev = 0.0013  #Range is 26.035 cm
+std_dev = 0.0006  #Range is 26.035 cm
 mean = 0.0
 blocksize = 100 #How many realisations to fit at a time
 
 number_of_bins = 10
-
-
-
-    
-
-
-
-

development/stochastic_study/run_model.py

@@ -44 +44 @@
 print 'Creating domain from', project.mesh_filename
 
-domain = Domain(project.mesh_filename,
+domain = Domain(project.working_dir + project.mesh_filename,
                 use_cache=False,
                 verbose=True)                
@@ -53 +53 @@
 
 
-domain.set_datadir('.')
+domain.set_datadir(project.working_dir)
 domain.set_quantities_to_be_stored(['stage', 'xmomentum', 'ymomentum'])
 
@@ -82 +82 @@
 
 realisation = 0
-for filename in os.listdir('.'):
+for filename in os.listdir(project.working_dir):
     if filename.startswith(project.basename) and filename.endswith('.pck'):
         print 'Reading %s' %filename
@@ -121 +121 @@
             # Now extract the 3 timeseries (Ch 5-7-9) and store them
             # in three files for this realisation
-            f = file_function(domain.filename + '.sww',
+            f = file_function(project.working_dir + domain.filename + '.sww',
                               quantities='stage',
                               interpolation_points=project.gauges,
@@ -127 +127 @@
 
 
-            simulation_name = domain.filename + '_realisation_%d' %realisation
+            simulation_name = project.working_dir + \
+                              domain.filename + '_realisation_%d' %realisation
 
             for k, name in enumerate(project.gauge_names):