function [T_esti,Vesti] = front_speed(css,tht, rhob, Kf, epsilon, q,fac)
%Estimate speed of the front
%feed q in m/s please
%   Detailed explanation goes here
epsilon =0;
R = (tht+rhob*(Kf/fac)*(1-epsilon)*mean(css).^(-epsilon));
Vesti = norm(q/R);
T_esti = 1/Vesti/60/60/24/365.25;
end