function [ ts ] = shallow_init(t,x)
% 
% Written by:
% -- 
% John L. Weatherwax                2005-06-15
% 
% email: wax@alum.mit.edu
% 
% Please send comments and especially bug reports to the
% above email address.
% 
%-----

ts(1,:) = repmat( 40.0, 1, length(x) ); 
ts(2,:) = 20.0 - bottom_H(x);