function ecodata=calc_eco(file_loc,dname,drftdata,calfile,calind,flag)

[scfa,dava,lambda]=read_ecobb2fcal(file_loc,calfile,calind);

goodind=find(drftdata(7,:)~=-999);
drftdata=drftdata(:,goodind);

vsclo=drftdata(34,:);
vschi=drftdata(36,:);
dchl=drftdata(37,:);
sal=drftdata(40,:);

Betalo=(vsclo-dava(1)).*scfa(1);
Betahi=(vschi-dava(2)).*scfa(2);
Chl=(dchl-dava(3)).*scfa(3);

vsfwlo=calc_bs(Betalo,str2num(lambda(1,:)),sal);
vsfwhi=calc_bs(Betahi,str2num(lambda(2,:)),sal);

Betaplo =  Betalo-vsfwlo;
Betaphi = Betahi-vsfwhi;

Bbplo = (2*pi*1.1).*Betaplo;
Bbphi = (2*pi*1.1).*Betaphi;

Bbswlo = (0.0029308*(str2num(lambda(1,:))/500)^-4.24)/2;
Bbswhi = (0.0029308*(str2num(lambda(2,:))/500)^-4.24)/2;

bblo = Bbplo + Bbswlo;
bbhi = Bbphi + Bbswhi;

ecodata=[Betalo;Betaplo;Bbplo;bblo;Betahi;Betaphi;Bbphi;bbhi;Chl];

x=1;
function vsfw=calc_bs(beta,lambda,sal)
    %Raleigh scattering
    Rcor = 1.38*(lambda/500)^-4.32;
    
    %Salinity Correction
    Scor = (1+0.3.*sal./37).*10^-4;
    
    %solid angle
    Acor = (1+((1+cos(2*117))/2)*(1-0.09)/(1+0.09));
    
    vsfw = Rcor.*Scor.*Acor;