p=log(2)/37;
q=1/24;
alpha=3500; %virus rilasciati da una cellula
beta=7*1e-10*alpha/24; %beta Friedmanetal06 * alpha (vedi paper per il motivo)
Du=((6.2-2.6)/2/40/24)^2/p; %ogni unita' spaziale vale 1 mm
Di=Du;

K=1e4;
Ru=2.6; %raggio iniziale non infette in mm
Ri=1; %raggio iniziale infette in mm

beta=beta*1e6/K;

Dc=0.8/24; %MCK04
qc=2/24; %Emma
% qc=log(2)/0.5; %too much
cstar=5^(1/3); %mug/mm^2
alphac=2.5/K; %valore massimo di Adrianne diviso...; 1/h/cell
gammac=alphac/50;

qT=0.18/24; %OK
alphaz=qT/50; %arbitrario

zU=10/K/10;
zI=zU;

DT=0.0023; %from microscopic considerations
chiT=1.65;
