source: trunk/src/psit.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: 586 bytes
Line 
1;+
2;-
3;--------------------------------------------------------
4function psit, zet
5;--------------------------------------------------------
6
7x=(1.-15*zet)^.5
8psik=2*alog((1+x)/2)
9x=(1.-34.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+2./3*zet(ii))^1.5+.6667*(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.