function [q,g,h,r]=scatterb(p,e,u,time)
%SCATTERB       Boundary condition data
%
%
%
% Copyright (c) 1994-98 by The MathWorks, Inc.
% $Revision: 1.7 $
bl=[
1 1 1 1 1 1 1 1
1 1 1 1 0 0 0 0
1 1 1 1 5 5 5 5
1 1 1 1 1 1 1 1
1 1 1 1 45 45 45 45
13 13 13 13 54 54 54 54
48 48 48 48 48 48 48 48
48 48 48 48 42 42 42 42
49 49 49 49 105 105 105 105
45 45 45 45 48 48 48 48
101 101 101 101 0 0 0 0
120 120 120 120 0 0 0 0
112 112 112 112 0 0 0 0
40 40 40 40 0 0 0 0
45 45 45 45 0 0 0 0
105 105 105 105 0 0 0 0
42 42 42 42 0 0 0 0
54 54 54 54 0 0 0 0
48 48 48 48 0 0 0 0
42 42 42 42 0 0 0 0
120 120 120 120 0 0 0 0
41 41 41 41 0 0 0 0
];

if any(size(u))
  [q,g,h,r]=pdeexpd(p,e,u,time,bl);
else
  [q,g,h,r]=pdeexpd(p,e,time,bl);
end

