v_inr 2*v_inlet*(1-r^2/r_compartment^2) "r-dependent inlet velocity factor"
v_in v_inr*profile(t) "Inlet velocity"
m_d intop1(phils*(z>l_dr_limit)*2*pi*r) "Droplet volume"