source: trunk/src/psiu.pro @ 152

Last change on this file since 152 was 89, checked in by pinsard, 13 years ago

add missing tools for coare correction
  • Property svn:keywords set to URL
File size: 653 bytes
Line 
1;+
2;-
3;--------------------------------------------------------
4function psiu, zet
5;--------------------------------------------------------
6
7x=(1-15*zet)^.25
8psik=2*alog((1+x)/2)+alog((1+x*x)/2)-2*atan(x)+2*atan(1)
9x=(1-10.15*zet)^.3333
10psic=1.5*alog((1+x+x*x)/3)-sqrt(3)*atan((1+2*x)/sqrt(3))+4*atan(1)/sqrt(3)
11f=zet*zet/(1+zet*zet)
12psi=(1-f)*psik+f*psic                                               
13
14ii=where(zet ge 0.)
15
16if (ii(0) ne -1) then begin
17  c=min([50,.35*zet])
18  psi(ii)=-((1+1.*zet(ii))^1.+.667*(zet(ii)-14.28)/exp(c(ii))+8.525)
19endif
20
21return, psi
22
23end
24;--------------------------------------------------------
25
Note: See TracBrowser for help on using the repository browser.