function [q,g,h,r]=geomb2(p,e,u,time)
%Boundary condition data

b=[
 1   1   1   1   1   1   1   1   1   1   1    1   1   1   1   1   1   1   1   1   1   1   1   1   1   1   1   1   1   1   %N
 1   1   1   1   1   1   1   1   1   1   1    0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   %M
 1   1   1   1   1   1   1   1   1   1   1    1   1   1   1   1   1   1   1   1   1   1   1   1   1   1   1   1   1   1   %Q
 1   1   1   1   1   1   1   1   1   1   1    1   1   1   1   1   1   1   1   1   1   1   1   1   1   1   1   1   1   1   %G
 1   1   1   1   1   1   1   1   1   1   1   48  48  48  48  48  48  48  48  48  48  48  48  48  48  48  48  48  48  48  %H
 1   1   1   1   1   1   1   1   1   1   1   48  48  48  48  48  48  48  48  48  48  48  48  48  48  48  48  48  48  48  %R
 48  48  48  48  48  48  48  48  48  48  48  0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   %Q
 48  48  48  48  48  48  48  48  48  48  48  0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   %G
 49  49  49  49  49  49  49  49  49  49  49  0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   %H
 48  48  48  49  49  49  48  48  48  48  48  0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   %R 
];

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