Changeset 3008 for development/stochastic_study

Timestamp:

May 29, 2006, 2:09:54 PM (19 years ago)

Author:

ole

Message:

Allowed parallel version to run sequentially

Location:

development/stochastic_study

Files:

• project.py (modified) (2 diffs)
• run_model.py (modified) (4 diffs)

development/stochastic_study/project.py

@@ -3 +3 @@
 
 # Working dir
-working_dir = '/mnt/store/stochastic/'
+working_dir = '/mnt/store/stochastic3/'
 import os
 try:
@@ -27 +27 @@
 
 # Stats (Suresh ?)
-number_of_realisations = 8000
+number_of_realisations = 1000
+#number_of_realisations = 16000
+
 #std_dev = 0.0026  #Range is 26.035 cm
-#std_dev = 0.0013  #Range is 26.035 cm
-std_dev = 0.0006  #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 = 10 #How many realisations to fit at a time

development/stochastic_study/run_model.py

@@ -32 +32 @@
 from pyvolution.util import file_function
 from caching.caching import cache
-import pypar
 
 # Application specific imports
@@ -42 +41 @@
 #-----------------------------------------------------------------------------
 
-myid = pypar.rank()
-numprocs = pypar.size()
-processor_name = pypar.Get_processor_name()
+try:
+    import pypar
+except:
+    print 'Could not import pypar'
+    myid = 0
+    numprocs = 1
+    processor_name = 'nautilus'
+else:    
+    myid = pypar.rank()
+    numprocs = pypar.size()
+    processor_name = pypar.Get_processor_name()
 
 print 'I am process %d of %d running on %s' %(myid, numprocs, processor_name)
@@ -91 +98 @@
 
 
-realisation = 0
+realisation = 100
 for filename in os.listdir(project.working_dir):
     if filename.startswith(project.basename) and filename.endswith('.pck'):
-        #print 'P%d: Reading %s' %(myid, filename)
-        fid = open(filename)
+        print 'P%d: Reading %s' %(myid, filename)
+        fid = open(project.working_dir + filename)
         V = cPickle.load(fid)
         fid.close()
@@ -110 +117 @@
                 name = project.basename + '_P%d' %myid    
                 domain.set_name(name)                     #Output name
+                print 'V', V.shape
                 domain.set_quantity('elevation', V[:,i])  #Assign bathymetry