psiu_25

0001 
0002 function psi=psiuo(zet)
0003 
0004     x=(1-15*zet).^.25;
0005     psik=2*log((1+x)/2)+log((1+x.*x)/2)-2*atan(x)+2*atan(1);
0006     x=(1-12.87*zet).^.3333;
0007     psic=1.5*log((1+x+x.*x)/3)-sqrt(3)*atan((1+2*x)/sqrt(3))+4*atan(1)/sqrt(3);
0008     f=zet.*zet/(1+zet.*zet);
0009     psi=(1-f).*psik+f.*psic;                                               
0010    ii=find(zet>0);
0011    if ~isempty(ii);
0012 
0013     psi(ii)=-4.7*zet(ii);
0014     %c(ii)=min(50,.35*zet(ii));
0015     %psi(ii)=-((1+1.0*zet(ii))^1.0+.667*(zet(ii)-14.28)/exp(c(ii))+8.525);
0016     end;