Ignore:
Timestamp:
Aug 10, 2011, 9:12:16 AM (14 years ago)
Author:
steve
Message:

setting up inlet testing script

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/anuga_work/development/2010-projects/anuga_1d/sww/sww_domain.py

    r8193 r8194  
    282282    #print 'any' ,numpy.any( h_V[:,1] < 0.0)
    283283
    284     #h00 = 1e-12
     284    h00 = 1e-12
    285285    #print h00
    286286
    287     h_V[:,:] = numpy.where (h_V <= h0, 0.0, h_V)
    288     u_V[:,:] = numpy.where (h_V <= h0, 0.0, u_V)
     287    h_V[:,:] = numpy.where (h_V <= h00, 0.0, h_V)
     288    u_V[:,:] = numpy.where (h_V <= h00, 0.0, u_V)
    289289
    290290#    # protect from edge values going negative
Note: See TracChangeset for help on using the changeset viewer.