Changeset 2989

Timestamp:

May 27, 2006, 11:28:10 AM (19 years ago)

Author:

ole

Message:

Ready to run parallel stochastic study (I hope)

Location:

development/stochastic_study

Files:

• project.py (modified) (1 diff)
• run_model.py (modified) (4 diffs)

development/stochastic_study/project.py

@@ -27 +27 @@
 
 # Stats (Suresh ?)
-number_of_realisations = 16000
+number_of_realisations = 8000
 #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
 mean = 0.0
-blocksize = 100 #How many realisations to fit at a time
+blocksize = 10 #How many realisations to fit at a time
 
 number_of_bins = 10

development/stochastic_study/run_model.py

@@ -94 +94 @@
 for filename in os.listdir(project.working_dir):
     if filename.startswith(project.basename) and filename.endswith('.pck'):
-        #print 'Reading %s' %filename
+        #print 'P%d: Reading %s' %(myid, filename)
         fid = open(filename)
         V = cPickle.load(fid)
         fid.close()
 
+        #if myid == 0:
+        #    print 'V', V[6:7,:]
+            
         # For each column (each realisation)
         for i in range(V.shape[1]):
@@ -108 +111 @@
                 domain.set_name(name)                     #Output name
                 domain.set_quantity('elevation', V[:,i])  #Assign bathymetry
+
+                print 'P%d: Setting quantity %d: %s' %(myid, i, str(V[:4,i]))
+                
                 domain.set_time(0.0)                      #Reset time
 
@@ -132 +138 @@
                 # Now extract the 3 timeseries (Ch 5-7-9) and store them
                 # in three files for this realisation
+
+                print 'P%d: Extracting time series for realisation %d from file %s'\
+                      %(myid, realisation, project.working_dir + domain.filename + '.sww')
                 f = file_function(project.working_dir + domain.filename + '.sww',
                                   quantities='stage',
@@ -140 +149 @@
                 simulation_name = project.working_dir + \
                                   project.basename + '_realisation_%d' %realisation
+
+                print 'P%d: Writing to file %s'\
+                      %(myid, simulation_name + '_' + name + '.txt')
 
                 for k, name in enumerate(project.gauge_names):