function [a_phi] = phyto_species_abs(wavelength) 

%set of four species specific phytoplankton absorption spectra
%measured by C. Roesler

a=load('phyto_species_abs.dat');
wave=400:1:700;
a_phi=interp1(wave,a,wavelength,'linear');%dinoflagellate diatom chlorophyte cryptomonad
