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=[2.3 3.8 4.6 5.6 6.4 7.5 8.5 9.5 10.5 11.3 12.6 13.3 14.5 ];
v1=[5.67 7.05 5.45 7.09 9.15 13.57 23.72 21.63 43.1 55.3 60.7 69.3 104.9];
sg=[.27 .97 6.6 10.5 11.3 16.3 25 24 35 50 33 29 68];
nn=[2 5 59 94 127 91 99 96 69 26 19 7 3];
u2p=[7 12 7 8 5.8 8.1 4.8 7.2 7 6.8 7 7.5 7.7 9 9 9.5 11 11.5 8 8.8];
v2m=[24.9 28.6 7.5 14.7 10.1 27 4.8 3.6 97 16.5 10.2 12.4 10.7 15 15.5 13.5 19.2 19.2 25.5 34.4 ];
u3=[4.7 7.4 7.4];
v3=[8.8 22 21.2];
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,'o',u,wann,'-',u4,v4,'--',u4,v5,':');set(gca,'fontsize',12);
axis([0 18 0 120]);xlabel('\fontsize{12} U_{10} (m/s)');ylabel('\fontsize{12} V_t (cm/hr)');%Title('Air Side Bulk Transfer Coeff')
sga=sg./sqrt(nn-1);
uu=zeros(2,13);
uu(1,:)=u1;
uu(2,:)=u1;
vv=zeros(2,13);
vv(1,:)=v1-sga;
vv(2,:)=v1+sga;
hold on;
for i=1:13
   plot(uu(:,i),vv(:,i),'.-');
end;
jam=sg./v1;
sgb=sg./sqrt(nn-2);

figure;plot(u1(3:12),jam(3:12),'o');
set(gca,'fontsize',12);
axis([0 18 0 1.75]);
xlabel('\fontsize{12} U_{10} (m/s)');ylabel('\fontsize{12} \sigma_{V_t}/V_t ');
vv2(1,:)=jam-sgb./v1;
vv2(2,:)=jam+sgb./v1;
hold on;
for i=3:12
   plot(uu(:,i),vv2(:,i),'.-');
end;


