Ignore:
Timestamp:
Mar 11, 2008, 8:43:22 PM (17 years ago)
Author:
steve
Message:

Updated some methods for quantity. Looks like we can use old
limiting system with larger values of beta.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_validation/convergence_study/animatesww2d_alt.py

    r4838 r5162  
    1 def animatesww2d_alt(sww_filename, movie_filename, range):
    2     """Plot cross section of model output
    3     """
     1## def animatesww2d_alt(sww_filename, movie_filename, range):
     2##     """Plot cross section of model output
     3##     """
    44
    5     # Read SWW file     
    6     from Scientific.IO.NetCDF import NetCDFFile
    7     fid = NetCDFFile(sww_filename, 'r')
     5##     # Read SWW file 
     6##     from Scientific.IO.NetCDF import NetCDFFile
     7##     fid = NetCDFFile(sww_filename, 'r')
    88   
    9     x = fid.variables['x']
    10     y = fid.variables['y']
    11     volumes = fid.variables['volumes']
     9##     x = fid.variables['x']
     10##     y = fid.variables['y']
     11##     volumes = fid.variables['volumes']
    1212   
    13     elevation = fid.variables['elevation']
    14     time = fid.variables['time']
    15     stage = fid.variables['stage']
    16     xmomentum = fid.variables['xmomentum']
    17     ymomentum = fid.variables['ymomentum']
     13##     elevation = fid.variables['elevation']
     14##     time = fid.variables['time']
     15##     stage = fid.variables['stage']
     16##     xmomentum = fid.variables['xmomentum']
     17##     ymomentum = fid.variables['ymomentum']
    1818   
    1919def animatesww2d(swwfile):
     
    3434    # interpolation points are for y = 0
    3535    # number of increments in x
    36     m = 100
     36    m = 1000
    3737    max_x = 100000.
    3838    x = 0
Note: See TracChangeset for help on using the changeset viewer.