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-10.15*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 c=min(50,.35*zet); 0016 psi(ii)=-((1+1.0*zet(ii)).^1.0+.667*(zet(ii)-14.28)./exp(c(ii))+8.525); 0017 end;