source: trunk/src/psiu.pro

Last change on this file was 204, checked in by pinsard, 10 years ago

fix thanks to coding rules; typo

  • Property svn:keywords set to Id URL
File size: 643 bytes
Line 
1;+
2;
3; ========
4; psiu.pro
5; ========
6;
7; .. function:: psiu(zet)
8;
9; DESCRIPTION
10; ===========
11;
12; EVOLUTIONS
13; ==========
14;
15; $Id$
16;
17; $URL$
18;
19; - fplod 20120319
20;
21;   * minimal header
22;
23;-
24function psiu, zet
25
26x=(1-15*zet)^.25
27psik=2*alog((1+x)/2)+alog((1+x*x)/2)-2*atan(x)+2*atan(1)
28x=(1-10.15*zet)^.3333
29psic=1.5*alog((1+x+x*x)/3)-sqrt(3)*atan((1+2*x)/sqrt(3))+4*atan(1)/sqrt(3)
30f=zet*zet/(1+zet*zet)
31psi=(1-f)*psik+f*psic                                               
32
33ii=where(zet ge 0.)
34
35if (ii(0) ne -1) then begin
36    c=min([50,.35*zet])
37    psi(ii)=-((1+1.*zet(ii))^1.+.667*(zet(ii)-14.28)/exp(c(ii))+8.525)
38endif
39
40return, psi
41
42end
Note: See TracBrowser for help on using the repository browser.