function [a_phi] = phyto_Lee(wavelength) 

%phytoplankton absorption spectrum as parameterized by Lee et al. 1996
%note options to change F, a1/a2 sigma (s)
%scaled as specific absorption coefficient

%Lee basis
wave1=400:1:570;
%F=[1.6 4.2];,a2=0.86+0.16*log(1);,s=14.17+0.9*log(1);
F=[1.6 4.2];,a2=[0.85 0.21];,s=14.17+0.9*log(1);
aph1=(exp(-F(1)*(log((wave1-340)/100)).^2));

wave3=656:1:700;
aph3=a2(1)*exp(-((wave3-676).^2)./(2*s.^2));

wave2=571:1:655;
aph2=aph1(171)+(aph3(1)-aph1(171))/(656-570)*(wave2-570);
a=[aph1';aph2';aph3'];
w=[wave1';wave2';wave3'];
I=find(w==676);a=0.0145*a/(a(I))
a_phi=interp1(w,a,wavelength,'linear');