initialize c_phi and c_grad_phi as arrays