function m = ep_chl_morel(e)
% EP_CHL_MOREL - Calculate Morel's chl from different K(lambda)
%
% This is not a production function, its meant to be changed and used
% for debugging
%

% Brian Schlining
% 10 Aug 1999

ep_ini
ep = 1;
i = strmatch('Pressure',ColumnS);
Depth = e(ep).dat(:,i);
Depth = 20; % Depth in meters

%i = strmatch('K412',ColumnS);
%m = ((e(ep).dat(:,i) - 0.0196)./.1125).^(1/0.680);
i = strmatch('Ed[3m+]412',ColumnS);
j = strmatch('Ed[20m]412',ColumnS);
K412 = (log(0.98*e(ep).dat(:,i)./e(ep).dat(:,j)))./Depth;
m = ((K412 - 0.0196)./0.1125).^(1/0.680);

%i = strmatch('K443',ColumnS);
%m(:,2) = ((e(ep).dat(:,i) - 0.0166)./.0996).^(1/0.704);
i = strmatch('Ed[3m+]443',ColumnS);
j = strmatch('Ed[20m]443',ColumnS);
K443 = (log(0.98*e(ep).dat(:,i)./e(ep).dat(:,j)))./Depth;
m(:,2) = ((K443 - 0.0166)./.0996).^(1/0.704);

%i = strmatch('K490',ColumnS);
%m(:,3) = ((e(ep).dat(:,i) - 0.0217)./.0690).^(1/0.702);
i = strmatch('Ed[3m+]490',ColumnS);
j = strmatch('Ed[20m]490',ColumnS);
K490 = (log(0.98*e(ep).dat(:,i)./e(ep).dat(:,j)))./Depth;
m(:,3) = ((K490 - 0.0217)./.0690).^(1/0.702);

%i = strmatch('K510',ColumnS);
%m(:,4) = ((e(ep).dat(:,i) - 0.0384)./.0578).^(1/0.690);
i = strmatch('Ed[3m+]510',ColumnS);
j = strmatch('Ed[20m]510',ColumnS);
K510 = (log(0.98*e(ep).dat(:,i)./e(ep).dat(:,j)))./Depth;
m(:,4) = ((K510 - 0.0384)./.0578).^(1/0.690);

%i = strmatch('K555',ColumnS);
%m(:,5) = ((e(ep).dat(:,i) - 0.0683)./.0400).^(1/0.645);
i = strmatch('Ed[3m+]555',ColumnS);
j = strmatch('Ed[20m]555',ColumnS);
K555 = (log(0.98*e(ep).dat(:,i)./e(ep).dat(:,j)))./Depth;
m(:,5) = ((K555 - 0.0683)./.0400).^(1/0.645);

%i = strmatch('K670',ColumnS);
%m(:,6) = ((e(ep).dat(:,i) - 0.4300)./.0515).^(1/0.695);
%i = strmatch('Ed[3m+]670',ColumnS);
%j = strmatch('Ed[20m]670',ColumnS);
%K670= (log(0.98*e(ep).dat(:,i)./e(ep).dat(:,j)))./Depth;
%m(:,6) = ((K670 - 0.4300)./.0515).^(1/0.695);