Debug PDE MATLAB code

I am using a MATLAB code to solve for a 1D diffusion problem in packed bed column. The code should solve a system of 14 nonlinear PDE equations. The code gives the error below. I tried to check that my vectors are consistent and could not find a problem. I will appreciate it if someone can debug the code and get it to work. Error using pdepe (line 231) Invalid output of ICFUN. ICFUN must return a column vector. Error in pdebfb_newTestcorrect (line 160) sol = pdepe(m,@pdebfbpde,@pdebfbic,@pdebfbbc,x,t);