VERB_code_2.3
RadialDiffCoeff.m
1 function [DLL] = RadialDiffCoeff(Kp,L)
2 %RadialDiffCoeff - This function returns the radial diffusion coefficients
3 %for L between 2 and 8 and for a provided value of Kp.
4 % This function enables one to reproduce the first row of figures in this paper
5 
6 DLL = 10^(0.506*Kp - 9.325)*(L.^10);
7 
8 %To see plot, uncomment following lines.
9 
10 %L = linspace(2,8);
11 %plot(L,1./DLL(2,L),'.-k');
12 %set(gca,'yscale','log','ylim',[1e-2 1e6],'ylabel','1/Dll','xlabel,'L');
13 end
14