function dt = adr_dt(t, ut, params, AD)
% ODE for advection diffusion reaction via finite differences
A=params.A; E=params.E; p=params.p; 
R=A*(1-ut).*exp(-E./(ut-p)); % reaction term
dt=AD*ut +R;              

