function  [new_oxy,Soc,Boc,tcor,wt,pcor,tau,gama]=...
    oxy_fit_press_beckman(curr,oxy_temp,OXSAT,temp,press,bot_oxy,docdt,...
    Soc,Boc,tcor,pcor,wt,tau,gama);
%OXY_FIT_PRESS_BECKMAN - pressure dependent fit for oxygen coefficients
% 

%	$Revision: 1.3 $
%	$Date: 2008/04/09 14:31:51 $
%	$Id: oxy_fit_press_beckman.m,v 1.3 2008/04/09 14:31:51 duncombe Exp duncombe $
% 
%

x=[curr OXSAT temp press docdt oxy_temp];
beta(1,1)=Soc;
beta(1,2)=Boc;
beta(1,3)=tcor;
beta(1,4)=wt;
beta(1,5)=pcor;
beta(1,6)=tau;
beta(1,7)=gama;
beta0=non_lin_fit(x,bot_oxy,'oxygen_fit_press_beckman',beta);
Soc=beta0(1);
Boc=beta0(2);
tcor=beta0(3);
wt=beta0(4);
pcor=beta0(5);
tau=beta0(6);
gama=beta0(7);

new_oxy=((Soc.*(curr+(tau.*docdt)))+Boc+(gama.*press)).* ...
	exp(tcor.*(temp+(wt*(oxy_temp-temp)))).*OXSAT.*exp(pcor.*press);

return;

