Shelf Response for two KdV Solitons

Shelf response for an identical pair of KdV solitons in the interface impinging upon the Shelf.

Information and Reference:
YouTube video: Changes in Shelf Response due to Increases in KdV Soliton Amplitude
http://www.youtube.com/watch?v=qywAEmS7_d4

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. http://oceanphysics.weebly.com/