epsilon_p_i "0.95 0.8 0.8" []
cF_i "0.14 0.095 0.095" []
Da_i "1.9e-4 1.04e-4 1.04e-4" []