Ignore:
Timestamp:
Aug 1, 2008, 1:33:26 PM (14 years ago)
Author:
kristy
Message:
 
File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_work/production/perth/get_gauges.py

    r5581 r5592  
    3636    for i, name in enumerate(quantity_names):
    3737        quantities[name] = fid.variables[name][:]
    38         if inside_points is not None:
    39             quantities[name] = take(quantities[name],indices,1)
    40         else:
    41             msg = 'No gauges found inside polygon'
    42             raise msg
    43            
     38##        if inside_points is not None:
     39##            quantities[name] = take(quantities[name],indices,1)
     40##        else:
     41##            msg = 'No gauges found inside polygon'
     42##            raise msg
     43##           
    4444    for j in range(len(x)-1):
    4545        locx=int(x[j])
     
    6363polygon=project.poly_all
    6464points,quantities,elevation,time=get_sts_gauge_data(os.path.join(project.boundaries_dir,project.scenario_name),polygon,verbose=False)
    65 
     65print quantities['stage'][0,:]
     66print len(points), len(elevation), len(quantities['stage'][0,:])
    6667assert len(points)==len(elevation)==len(quantities['stage'][0,:])
Note: See TracChangeset for help on using the changeset viewer.