Changeset 6290

Timestamp:

Feb 6, 2009, 9:41:31 AM (16 years ago)

Author:

jgriffin

Message:

New format and naming conventions

File:

• anuga_work/production/australia_ph2/sydney/build_boundary_7875.py (modified) (7 diffs)

anuga_work/production/australia_ph2/sydney/build_boundary_7875.py

@@ -3 +3 @@
 WA, Australia.  The boundary is based on the National Hazard Map.
 
-Input: order_filename from project_250m.py
-       event_number needs to be reflected in the project_250m file 
-Output: creates a sts file and csv files stored in project_250m.boundaries_dir_event.
+Input: order_filename from project.py
+       event_number needs to be reflected in the project file 
+Output: creates a sts file and csv files stored in project.boundaries_dir_event.
 The run_sydney.py is reliant on the output of this script.
 """
 import os
 from os import sep
-import project_250m
+from os.path import join
+import project
 
 from anuga.utilities.numerical_tools import ensure_numeric
@@ -29 +30 @@
 #--------------------------------------------------------------------------
 
-dir=os.path.join(project_250m.muxhome,'mux')
+dir=os.path.join(project.muxhome,'mux')
 
 # Refer to event_07875.list in home+state+sep+event08 for event details
@@ -63 +64 @@
     weights=weight_factor*ones(len(urs_filenames),Float)
 
-    scenario_name=project_250m.scenario_name
-    order_filename=os.path.join(project_250m.urs_boundary_name)
+    scenario_name=project.scenario_name
+    urs_order = join(project.urs_order)
 
-    print 'reading', order_filename
-    # Create ordered sts file
-    print 'creating sts file'
+    print 'reading', urs_order
 
     urs2sts(urs_filenames,
-            basename_out=os.path.join(project_250m.boundaries_dir_event,scenario_name),
-            ordering_filename=order_filename,
-            weights=weights,
-            mean_stage=project_250m.tide,
-            verbose=True)
+            basename_out = join(project.event_sts),
+            ordering_filename = urs_order,
+            weights = weights,
+            mean_stage = project.tide,
+            verbose = True)
 else:
     print 'number of sources do not match event.list'
@@ -99 +98 @@
         quantities[name] = fid.variables[name][:]
 
-#------------------------------------------------------------------------------
-# Get Maxium wave height throughout timeseries at each index point
-#------------------------------------------------------------------------------
+    #------------------------------------------------------------------------------
+    # Get Maxium wave height throughout timeseries at each index point
+    #------------------------------------------------------------------------------
 
     maxname = 'max_sts_stage.csv'
-    fid_max = open(project_250m.boundaries_dir_event+sep+maxname,'w')
+    fid_max = open(project.boundaries_dir_event+sep+maxname,'w')
     s = 'index, x, y, max_stage \n'
     fid_max.write(s)    
@@ -116 +115 @@
         fid_max.write(s)
       
-#------------------------------------------------------------------------------
-# Get Minium wave height throughout timeseries at each index point
-#------------------------------------------------------------------------------
+    #------------------------------------------------------------------------------
+    # Get Minium wave height throughout timeseries at each index point
+    #------------------------------------------------------------------------------
 
     minname = 'min_sts_stage.csv'
-    fid_min = open(project_250m.boundaries_dir_event+sep+minname,'w')
+    fid_min = open(project.boundaries_dir_event+sep+minname,'w')
     s = 'index, x, y, max_stage \n'
     fid_min.write(s)    
@@ -135 +134 @@
 
 
-        fid_sts = open(project_250m.boundaries_dir_event+sep+basename+'_'+ str(index)+'.csv', 'w')
+        fid_sts = open(project.boundaries_dir_event+sep+basename+'_'+ str(index)+'.csv', 'w')
         s = 'time, stage, xmomentum, ymomentum \n'
         fid_sts.write(s)
 
-#------------------------------------------------------------------------------
-# End of the get gauges
-#------------------------------------------------------------------------------
+    #------------------------------------------------------------------------------
+    # End of the get gauges
+    #------------------------------------------------------------------------------
         for k in range(len(time)-1):
             s = '%.6f, %.6f, %.6f, %.6f\n' %(time[k], stage[k], xmomentum[k], ymomentum[k])
@@ -153 +152 @@
     return quantities,elevation,time
 
-quantities,elevation,time=get_sts_gauge_data(os.path.join(project_250m.boundaries_dir_event,project_250m.scenario_name),verbose=False)
+quantities,elevation,time=get_sts_gauge_data(os.path.join(project.boundaries_dir_event,project.scenario_name),verbose=False)
 
 print len(elevation), len(quantities['stage'][0,:])