function [q,g,h,r]=geomb_cap(p,e,u,time)
%Boundary condition data
b=[
 1   1   1   1   1   1   %N
 1   1   1   1   1   1   %M
 1   1   1   1   1   1   %Q
 1   1   1   1   1   1   %G
 1   1   1   1   1   1   %H
 1   1   1   1   1   1   %R
 48  48  48  48  48  48  %Q
 48  48  48  48  48  48  %G
 49  49  49  49  49  49  %H
 48  48  48  48  49  48 %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
