Changeset 3630


Ignore:
Timestamp:
Sep 20, 2006, 11:59:26 AM (18 years ago)
Author:
sexton
Message:

Incorporation of gauge timeseries into runup convergence example (comment out prior interrogation using 'gaugelines'. Require .xya file (added). Output images written to same directory as sww file.

Location:
anuga_core/source/anuga/examples
Files:
1 added
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_core/source/anuga/examples/runup_convergence.py

    r3618 r3630  
    3838normalized_amplitude = amplitude/offshore_depth
    3939simulation_name = 'runup_convergence'   
    40 
     40coastline_x = -highest_point/slope
    4141
    4242# Basin dimensions (m)
     
    155155#-----------------------------------------------------------------------------
    156156
    157 
     157'''
    158158# Define line of gauges through center of domain
    159159def gauge_line(west,east,north,south):
     
    213213print 'Coastline (meters form west):       ', coastline
    214214
    215 
    216 
     215'''
     216# Generate time series of "gauge" situated at right hand boundary
     217from anuga.abstract_2d_finite_volumes.util import sww2timeseries
     218production_dirs = {'.': 'test'}
     219swwfiles = {}
     220for label_id in production_dirs.keys():
     221   
     222    swwfile = simulation_name + '.sww'
     223    swwfiles[swwfile] = label_id
     224   
     225texname, elev_output = sww2timeseries(swwfiles,
     226                                      'boundary_gauge.xya',
     227                                      production_dirs,
     228                                      report = False,
     229                                      reportname = 'test',
     230                                      plot_quantity = ['stage', 'speed'],
     231                                      surface = False,
     232                                      time_min = None,
     233                                      time_max = None,
     234                                      title_on = True,
     235                                      verbose = True)
    217236
    218237# Stop here
Note: See TracChangeset for help on using the changeset viewer.