source: trunk/src/psit.pro @ 169

Last change on this file since 169 was 164, checked in by pinsard, 12 years ago

fix for developper documentation
  • Property svn:keywords set to URL
File size: 614 bytes
Line 
1;+
2;
3; ========
4; psit.pro
5; ========
6;
7; .. function:: psit(zet)
8;
9; DESCRITPION
10; ===========
11;
12; EVOLUTIONS
13; ==========
14;
15; $Id$
16;
17; $URL$
18;
19; - fplod 20120319
20;
21;   * minimal header
22;
23;-
24function psit, zet
25
26x=(1.-15*zet)^.5
27psik=2*alog((1+x)/2)
28x=(1.-34.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+2./3*zet(ii))^1.5+.6667*(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.