File Exchange

image thumbnail

Shelf Response for two KdV Solitons

version (4 KB) by Edwin Alfonso
Shelf response for an identical pair of KdV solitons in the interface impinging upon the Shelf.


Updated 09 Dec 2013

View Version History

View License

Information and Reference:
YouTube video: Changes in Shelf Response due to Increases in KdV Soliton Amplitude

Calculates the Shelf Response for an identical pair of incident KdV solitons. Based on the article of Chapman, D. C. and G. S. Giese, 1990. J. Phys. Oceanogr. 20 1459-1467 pp.:
"The incident internal pulse at the interface generates a surface pulse at the shelf break with amplitude Q0/2 which travels across the shelf toward the coast. It's amplitude doubles
to Q0 at the coast where it reflects from the coastal wall and travels back across the shelf toward the shelf break. Upon reaching the shelf break, part of the pulse is reflected
shoreward while part continues into the deep ocean leaking energy to deep-ocean surface and internal waves." The incident pulses are negative amplitude KdV solitons.The pulse
amplitude can be set. The time shift between pulses can be set. The shelf depth and length, the deep-ocean upper layer depth, lower layer depth, total depth and respective densities
can be changed to typical measured values. Restriction: the shelf depth is less than the deep-ocean upper layer depth: S < H1. To get the dimensional values of eta_s and zeta_pulse
multiply each by eta_i; for u_s multiply by sqrt(g/H)*eta_i.

The m-file was written for Matlab R12.1 by Edwin Alfonso-Sosa, Ocean Physics Education Inc., 2013.

Cite As

Edwin Alfonso (2020). Shelf Response for two KdV Solitons (, MATLAB Central File Exchange. Retrieved .

Comments and Ratings (2)


Zong Hua Chen


MATLAB Release Compatibility
Created with R12.1
Compatible with any release
Platform Compatibility
Windows macOS Linux

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!