function [ ts ] = elliptic_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.
% 
%-----

% $$$ disp( ["t(1)=",num2str(t(1))] );
% $$$ disp( ["x(1)=",num2str(x(1))] )
 
inds = find( x==0.0 );     
   
ts = zeros(size(x)); 
if( ~isempty(inds) )
  ts(inds)=1.0;
end

% $$$ disp( ["ts(1)=",num2str(ts(1))] )
% $$$ pause