Changeset 6368

Timestamp:

Feb 19, 2009, 12:06:56 PM (15 years ago)

Author:

kristy

Message:

took out build_urs_boundary.py from the top of the run script

File:

• anuga_work/production/busselton/standardised_version/run_model.py (modified) (3 diffs)

anuga_work/production/busselton/standardised_version/run_model.py

@@ -37 +37 @@
 from anuga.interface import create_sts_boundary
 from anuga.interface import csv2building_polygons
-<<<<<<< .mine
-=======
 from file_length import file_length
 
->>>>>>> .r6341
 from anuga.shallow_water.data_manager import start_screen_catcher
 from anuga.shallow_water.data_manager import copy_code_files
@@ -50 +47 @@
 from setup_model import project
 
-
-#-------------------------------------------------------------------------------
-# Get gauges (timeseries of index points)
-#-------------------------------------------------------------------------------
-def get_sts_gauge_data(filename, verbose=False):
-    fid = NetCDFFile(filename+'.sts', 'r')      #Open existing file for read
-    permutation = fid.variables['permutation'][:]
-    x = fid.variables['x'][:] + fid.xllcorner   #x-coordinates of vertices
-    y = fid.variables['y'][:] + fid.yllcorner   #y-coordinates of vertices
-    points = num.transpose(num.asarray([x.tolist(), y.tolist()]))
-    time = fid.variables['time'][:] + fid.starttime
-    elevation = fid.variables['elevation'][:]
-       
-    basename = 'sts_gauge'
-    quantity_names = ['stage', 'xmomentum', 'ymomentum']
-    quantities = {}
-    for i, name in enumerate(quantity_names):
-        quantities[name] = fid.variables[name][:]
-
-    #---------------------------------------------------------------------------
-    # Get maximum wave height throughout timeseries at each index point
-    #---------------------------------------------------------------------------
-
-    maxname = 'max_sts_stage.csv'
-    fid_max = open(os.path.join(project.event_sts, maxname), 'w')
-    fid_max.write('index, x, y, max_stage \n')    
-    for j in range(len(x)):
-        index = permutation[j]
-        stage = quantities['stage'][:,j]
-        xmomentum = quantities['xmomentum'][:,j]
-        ymomentum = quantities['ymomentum'][:,j]
-
-        fid_max.write('%d, %.6f, %.6f, %.6f\n' % (index, x[j], y[j], max(stage)))
-      
-    #---------------------------------------------------------------------------
-    # Get minimum wave height throughout timeseries at each index point
-    #---------------------------------------------------------------------------
-
-    minname = 'min_sts_stage.csv'
-    fid_min = open(os.path.join(project.event_sts, minname), 'w')
-    fid_min.write('index, x, y, max_stage \n')    
-    for j in range(len(x)):
-        index = permutation[j]
-        stage = quantities['stage'][:,j]
-        xmomentum = quantities['xmomentum'][:,j]
-        ymomentum = quantities['ymomentum'][:,j]
-
-        fid_min.write('%d, %.6f, %.6f, %.6f\n' %(index, x[j], y[j], min(stage)))
-
-        out_file = os.path.join(project.event_sts,
-                                basename+'_'+str(index)+'.csv')
-        fid_sts = open(out_file, 'w')
-        fid_sts.write('time, stage, xmomentum, ymomentum \n')
-
-        #-----------------------------------------------------------------------
-        # End of the get gauges
-        #-----------------------------------------------------------------------
-        for k in range(len(time)-1):
-            fid_sts.write('%.6f, %.6f, %.6f, %.6f\n'
-                          % (time[k], stage[k], xmomentum[k], ymomentum[k]))
-
-        fid_sts.close()      
-
-    fid.close()
-
-    return quantities,elevation,time
 
 
@@ -226 +157 @@
 
 # Create the STS file
-<<<<<<< .mine
-build_urs_boundary(project.mux_data_folder,
-                   project.mux_input,
-                   os.path.join(project.event_folder,
-                                project.scenario_name))
-=======
 print 'project.mux_data_folder=%s' % project.mux_data_folder
 bub.build_urs_boundary(project.mux_input_filename,
                        os.path.join(project.event_folder,
                                     project.scenario_name))
->>>>>>> .r6341
 
 # Read in boundary from ordered sts file