Changeset 7799


Ignore:
Timestamp:
Jun 7, 2010, 11:41:49 AM (10 years ago)
Author:
habili
Message:

changed the order of assigning "no values" to stage, current_x and current_y.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/anuga_1_1/anuga_work/production/bunbury_storm_surge_2009/process_gems_grids.py

    r7794 r7799  
    122122        #    print 'third block'
    123123
    124         if i > 13 + 3*nrows and i < (14 + 4*nrows):   
     124        if i > (13 + 3*nrows) and i < (14 + 4*nrows):   
    125125            d += L.strip().split()
    126126           
     
    158158no_value_index = numpy.where(((depth < -9000) + (current_x < -9000) + (current_y < -9000)) == True)[0]
    159159
    160 numpy.put(stage, no_value_index, -9999)
    161 numpy.put(current_x, no_value_index, -9999)
    162 numpy.put(current_y, no_value_index, -9999)
    163 numpy.put(depth, no_value_index, 0)
     160#numpy.put(stage, no_value_index, -9999)
     161#numpy.put(current_x, no_value_index, -9999)
     162#numpy.put(current_y, no_value_index, -9999)
     163#numpy.put(depth, no_value_index, 0)
    164164   
    165165# Taking absolute value is to account for -ve depths obtained when depth-elevation
     
    174174ymomentum = current_y*depth #momentum*numpy.cos(numpy.radians(current_y))
    175175
     176numpy.put(stage, no_value_index, 0)
    176177numpy.put(xmomentum, no_value_index, 0)
    177178numpy.put(ymomentum, no_value_index, 0)
Note: See TracChangeset for help on using the changeset viewer.