kp=0.4;
sc=0.6;%.91*.91;
lam=13.3;
A=1.85;
ul=0:.1:2.9;
u=exp(ul);
[n m]=size(u);
t=22;
P=1010;
Rgas=287.1;
tdk=273.15;
rx=.75;
qs=24;
q=qs*rx;
rhoa=P*100/(Rgas*(t+tdk)*(1+0.61*q));
dt=0;
rhow=1022;

xo=[ 0 t+dt t qs q 141 419 0 600 1010 10 10 10];
x=[5 0 15 14 10 8.5 141 419 0 600 1010 15 15 15];
p=zeros(m,14);
p(:,1)=u';
sca=.83;
for i=1:m
   p(i,2:14)=xo;
   y=coul_lkb(p(i,:));
   dw(i)=y(13);
   us(i)=y(8);
   zo(i)=y(4);
   cd(i)=y(16);
   clkb(i)=y(21);
   jam=cor26_gs(p(i,:));
   s=solco2(t);
   sol(i)=s(1);%mole/m^3/atm
   al(i)=s(2);%dimensionless
   scw(i)=s(3);
   kc(i)=jam(22)*(660/scw(i));
end;
rr=zo.*us/1.5e-5;
f=3.8e-6*u.^3.4;
cdnh=kp./log(10./zo);
ha=lam/A;
wa=25;
gw=.03;
vs1=us./kc+sqrt(rhow/rhoa)*(log(.5./dw)/.4);
vs2=ha*sqrt(sca)+1./sqrt(cd)-5+.5*log(sc)/.4;
vs1m=(1-f)./vs1+f*gw*wa/sqrt(rhow/rhoa);
v=us./(vs1+al.*vs2);
v2=us./(1./vs1m+al.*vs2);
wann=0.31*u.^2./(scw/660).^.5;
kb=850*f;
fd=al.*(1+1./(14*al./scw.^.5).^(1/1.2)).^1.2;%Woolf equation
kbb=2450*f./fd;
u1=[5.5 12.5 ];
v1=[11.5 35];
u2=[3.6 5.7 6.42 6.7 6.8 7.9 10.7 11.45 12.4 13 15.7];
v2=[21 21.5 33.9 50 61 79 77.7 139 95 153 166];
u3=[4.67 5.35 6.5 7.5 8.6 9.5 10.9];
v3=[-13 -26 17 12 26 25 44];
v3x=[30 31 40 47 60 60 85];
u7=[0 1.1 2.5 5 8.6 8.7 9.84 12.7];
v7=[6 6.7 13 18 33.7 37 32 42.7];
u4=[0 1 2 3 3.6 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18];
v4=[0 .17 .34 .51 .61 1.75 4.6 7.45 10.3 13.15 16 18.85 21.7 24.55 27.4 33.3 39.2 45.1 51 56.9];
v5=6.5e-5*u4*100*3600;
figure;plot(u1,v1,'pentagram',u2,v2,'*',u3,v3,'o',u3,v3x,'^',u,wann,'-',u4,v4,'--',u4,v5,':',u7,v7,'diamond');
axis([0 18 0 160]);xlabel('\fontsize{12} U_{10} (m/s)');ylabel('\fontsize{12} V_t (cm/hr)');set(gca,'fontsize',12);
%Title('Air Side Bulk Transfer Coeff')
