Changeset 5459 for anuga_work/development/Hinwood_2008/prepare_time_boundary.py

Timestamp:

Jul 2, 2008, 4:34:48 PM (15 years ago)

Author:

duncan

Message:

Current Hinwood scenario

File:

• anuga_work/development/Hinwood_2008/prepare_time_boundary.py (modified) (10 diffs)

anuga_work/development/Hinwood_2008/prepare_time_boundary.py

@@ -15 +15 @@
 from Scientific.IO.NetCDF import NetCDFFile
 from Numeric import array, zeros, Float
+from os.path import join
   
 from anuga.utilities.numerical_tools import ensure_numeric
@@ -20 +21 @@
 from interp import interp
 
-from os.path import join
+
+import project
 
 # from os import getenv
@@ -115 +117 @@
     
     # Read velocity file
+    #print "*********************"
+    #print "velocity_file", velocity_file
     vfid = open(velocity_file)
     lines = vfid.readlines()
@@ -133 +137 @@
 
     # Read the depth file
+    #print "depth_file", depth_file
     dfid = open(depth_file)
     lines = dfid.readlines()
@@ -159 +164 @@
     #print "len(dtimes)", len(vtimes)
     #print "len(depths_at_vtimes)", len(depths_at_vtimes)
-#     for i in range(len(depths_at_vtimes)):
+    #print "metadata_dic['scenario_id']", metadata_dic['scenario_id']
+ #    for i in range(len(depths_at_vtimes)):
 #         print "i", i
 #         print "vtimes[i]", vtimes[i]
 #         print "depths_at_vtimes[i]", depths_at_vtimes[i]
-#     print "depths_at_vtimes",  depths_at_vtimes
+
+    #print "depths_at_vtimes",  depths_at_vtimes
     depths_at_vtimes = depths_at_vtimes/1000.00 # convert from mm to m
     missing=missing/1000.00 # Do to missing what is done to depths_at_vtimes
@@ -192 +199 @@
 
     fid.close()
+    print "The start time for is", start_time
+    if not start_time == metadata_dic['ANUGA_start_time']:
+        raise ValueError, "The calc'ed and recorded start times are different"
+    # Modify the sensor array to reflect the new start time
+    sensors[:,0] -= start_time
     
     # Since there is a new time reference save the depth info using this
     # new reference.
     fid = open(depth_file[:-4] + '_exp_depth.csv','w')
-    sensors[:,0] -= start_time
     #print "depth_file", depth_file
     #print "start_time", start_time
@@ -226 +237 @@
     # new reference.
     fid = open(depth_file[:-4] + '_exp_stage.csv','w')
-    sensors[:,0] -= start_time
     
     # Write a header
@@ -248 +258 @@
 def prepare_time_boundary(metadata_dic, raw_data_dir, output_dir):
     """
+    Use this if a project instance has already been created.
     """
     scenario_id = metadata_dic['scenario_id']
@@ -256 +267 @@
     
     final_time = combine_velocity_depth(velocity_file, depth_file, out_file,
-                                        metadata_dic)  
+                                        metadata_dic)
+    #print "metadata_dic['xleft'][1]", metadata_dic['xleft'][1]
     if metadata_dic['xleft'][1] >= 0.0:
         # This should be a -ve value, since the still water level is the
@@ -264 +276 @@
     csv2tms(tsm_file, metadata_dic['xleft'][1])
     return final_time
+
+# Don't do this, since run-dam changes the  metadata_dic['xleft'][1],
+# which is used by this function
+def prepare_time_boundary_for_scenarios_warning(scenarios,
+                                        outputdir_tag):
+
+    for run_data in scenarios:       
+        id = run_data['scenario_id']
+        outputdir_name = id + outputdir_tag
+        pro_instance = project.Project(['data','flumes','Hinwood_2008'],
+                                       outputdir_name=outputdir_name)
+        prepare_time_boundary(run_data,
+                              pro_instance.raw_data_dir,
+                              pro_instance.boundarydir)
+        
     
 #-------------------------------------------------------------------
 if __name__ == "__main__":
-    pass
+    
+    import sys; sys.exit()
+    
+    from scenarios import scenarios
+    outputdir_tag = "XXXX"
+    prepare_time_boundary_for_scenarios(scenarios,
+                                        outputdir_tag)