Changeset 5695 for anuga_validation/Hinwood_2008/calc_rmsd.py

Timestamp:

Aug 27, 2008, 10:04:24 AM (16 years ago)

Author:

duncan

Message:

Get the output directory going

File:

• anuga_validation/Hinwood_2008/calc_rmsd.py (modified) (6 diffs)

anuga_validation/Hinwood_2008/calc_rmsd.py

@@ -1 +1 @@
 """
-
-All err results are going into the same dir, and it can't really be changed.
+Functions used to calculate the root mean square deviation.
+
+Duncan Gray, GA - 2007 
+
 """
 
@@ -13 +15 @@
 from csv import writer
 from time import localtime, strftime
+from os.path import join
 
 # Related major packages
 from Numeric import zeros, Float, where, greater, less, compress, sqrt, sum
-
 from anuga.shallow_water.data_manager import csv2dict
 from anuga.utilities.numerical_tools import ensure_numeric, err, norm
 from anuga.utilities.interp import interp
 
+# Scenario specific imports
+import project
+
 def get_max_min_condition_array(min, max, vector):
+    """
+    Given a vector of values, and minimum and maximum values, return a
+    vector of 0/1's that can be used to cut arrays so only the times
+    in the min max range are used.
+    
+    precondition: The vector values are ascending.
+    
+    """
     
     SMALL_MIN = -1e10  # Not that small, but small enough
@@ -48 +61 @@
         id = run_data['scenario_id']
         outputdir_name = id + outputdir_tag
-        file_sim = outputdir_name + '_' +  quantity + ".csv"
+        file_sim = join(project.output_dir,outputdir_name + '_' +  \
+                        quantity + ".csv")
         file_exp = id + '_exp_' + quantity + '.csv'
-        file_err = outputdir_name + "_" + quantity + "_err.csv"
+        file_err = join(project.output_dir,outputdir_name + "_" + \
+                        quantity + "_err.csv")
         
 
@@ -60 +75 @@
         time_sim = ensure_numeric(time_sim)
         time_exp = ensure_numeric(time_exp)
-        
         condition = get_max_min_condition_array(run_data['wave_times'][0],
                                                 run_data['wave_times'][1],
@@ -145 +159 @@
         id = scenario['scenario_id']
         outputdir_name = id + outputdir_tag
-        file_err =  outputdir_name + "_" + quantity + "_err.csv"
+        file_err =  join(project.output_dir,outputdir_name + "_" + \
+                         quantity + "_err.csv")
         file_errs.append(file_err)
     return file_errs
@@ -196 +211 @@
     #scenarios = [scenarios[0]] # !!!!!!!!!!!!!!!!!!!!!!
 
-    outputdir_tag = "_nolmts_wdth_0.1_z_0.0_ys_0.01_mta_0.01"
+    outputdir_tag = "_nolmts_wdth_0.1_z_0.0_ys_0.5_mta_0.01_A"
     calc_norms = True
     #calc_norms = False