Changeset 5643


Ignore:
Timestamp:
Aug 13, 2008, 8:14:27 AM (16 years ago)
Author:
kristy
Message:

Updated script so that get_timeseries.py can still be acquired even though some of the gauges could be out of the domain.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • anuga_core/source/anuga/abstract_2d_finite_volumes/util.py

    r5482 r5643  
    22622262                            vel = 0.0
    22632263                        else:
    2264                             momentum = sqrt(point_quantities[2]**2 +\
    2265                                             point_quantities[3]**2)
    2266 #                            vel = momentum/depth             
    2267                             vel = momentum/(point_quantities[0] - point_quantities[1])
    2268 #                            vel = momentum/(depth + 1.e-6/depth)             
    2269                        
     2264                            if point_quantities[2] < 1.0e6:
     2265                                momentum = sqrt(point_quantities[2]**2 +\
     2266                                                point_quantities[3]**2)
     2267    #                            vel = momentum/depth             
     2268                                vel = momentum/(point_quantities[0] - point_quantities[1])
     2269    #                            vel = momentum/(depth + 1.e-6/depth)
     2270                            else:
     2271                                momentum = 0
     2272                                vel = 0
     2273                           
    22702274                        points_list.append(vel)
    22712275                       
Note: See TracChangeset for help on using the changeset viewer.