Spill packageEfficient Algorithms for Finding Sills in Digital Topographic Maps |
Thomas D. Otto, Andreas Thurnherr, Efficient Algorithms for Finding Sills in Digital Topographic Maps, submitted 2006
In stratified geophysical flows, the energetically optimal exchange of dense fluid across a topographic barrier generally takes place at the deepest unblocked connection, which is typically a saddle point (sill). The flow at or near a sill is often hydraulically controlled, in which case is is called a controlling sill. Oceanographic examples include overflows of newly formed dense water at high latitudes as well as sills in channels connecting major ocean basins, such as the Strait of Gibraltar. Controlling sills are often associated with strong flows, making them ideal sites for monitoring transport variability and detecting change. The locations and depths of controlling sills also provide strong constraints for the downstream hydrographic properties below sill depth. Here, two algorithms for finding sills in digital topographic maps are presented. The first approximates the sill height to arbitrary accuracy in O(k) steps, where k is the number of data points in the map. The second algorithm, which requires O(k \log k) steps, additionally returns the sill location. Several Tests carried out with realistic problems from physical oceanography reveal that the second algorithm runs faster in practice, even though its worst-case behavior is worse.
binary Linux 32 bits (kernel 2.4)
binary Linux 32 bits (kernel 2.6)
Notice:
An example map (RJd.grd) is included in the tar-ball.
The GMT software package is needed
using the source code.
The stack size must be adaptet, if using big maps. (ulimit -s 65000 for mac
OS X).
Please send comments or critics to Thomas D. Otto - January 2006.
Fundação Oswaldo Cruz |