|
Last change
on this file since 86 was
3,
checked in by slebonnois, 15 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
|
| Rev | Line | |
|---|
| [3] | 1 | FUNCTION GAMMB(W,T,TAU1,TAU2) |
|---|
| 2 | IMPLICIT REAL (A-H,O-Z) |
|---|
| 3 | DATA HK/7.638967E-12/,PI/3.141592654/ |
|---|
| 4 | cfix ,TAU10/0./,TAU20/0./ 367. |
|---|
| 5 | C NOTE. HK = 1.05459e-27/1.38054e-16 |
|---|
| 6 | cfix IF (TAU1.NE.TAU10) GO TO 10 |
|---|
| 7 | cfix IF (TAU2.EQ.TAU20) GO TO 20 |
|---|
| 8 | cfix 10 CONTINUE |
|---|
| 9 | TAU12=TAU1*TAU1 |
|---|
| 10 | TAU22=TAU2*TAU2 |
|---|
| 11 | HBH = 0.5*HK/T |
|---|
| 12 | Z2=SQRT ( TAU22 + HBH**2) / TAU1 |
|---|
| 13 | cfix TAU10 = TAU1 |
|---|
| 14 | cfix TAU20=TAU2 |
|---|
| 15 | cfix 20 CONTINUE |
|---|
| 16 | WSQR=W*W |
|---|
| 17 | Z = SQRT (1.+ WSQR*TAU12) * Z2 |
|---|
| 18 | C COMPUTE THE MODIFIED BESSEL FUNCTION OF THE SECOND KIND (K1) |
|---|
| 19 | C 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.
