source: trunk/LMDZ.TITAN/libf/phytitan/gammb.F @ 306

Last change on this file since 306 was 3, checked in by slebonnois, 14 years ago

Creation de repertoires:

  • chantiers : pour communiquer sur nos projets de modifs
  • documentation : pour stocker les docs

Ajout de:

  • libf/phytitan : physique de Titan
  • libf/chimtitan: chimie de Titan
  • libf/phyvenus : physique de Venus
File size: 806 bytes
Line 
1      FUNCTION GAMMB(W,T,TAU1,TAU2)
2      IMPLICIT REAL (A-H,O-Z)
3      DATA HK/7.638967E-12/,PI/3.141592654/
4cfix ,TAU10/0./,TAU20/0./    367.
5C        NOTE. HK = 1.05459e-27/1.38054e-16
6cfix  IF (TAU1.NE.TAU10) GO TO 10
7cfix  IF (TAU2.EQ.TAU20) GO TO 20
8cfix 10   CONTINUE
9      TAU12=TAU1*TAU1
10      TAU22=TAU2*TAU2
11      HBH = 0.5*HK/T
12      Z2=SQRT ( TAU22 + HBH**2) / TAU1
13cfix  TAU10 = TAU1
14cfix  TAU20=TAU2
15cfix 20   CONTINUE
16      WSQR=W*W
17      Z = SQRT (1.+ WSQR*TAU12) * Z2
18C        COMPUTE THE MODIFIED BESSEL FUNCTION OF THE SECOND KIND (K1)
19C        USING AN UNPUBLISHED APPROXIMATION GIVEN BY  COHEN.
20      F=1.5707963*(Z+0.5616)/(Z+0.4619)
21      BK1=SQRT(1.+Z*F)
22      BK1=BK1*EXP(-Z)
23      GAMMA=TAU1/PI*EXP(TAU2/TAU1+HBH*W)*BK1 / (1.+WSQR*TAU12)
24      GAMMB=GAMMA
25      RETURN
26      END
Note: See TracBrowser for help on using the repository browser.