]>
Commit | Line | Data |
---|---|---|
e74335a4 | 1 | * $Id$ |
2 | C | |
3 | C | |
4 | C | |
5 | C | |
6 | FUNCTION WDSAX1(X) | |
7 | C ********THREE PARAMETER WOOD SAXON | |
8 | C FOR PROJECTILE | |
9 | C HINT1(72)=R, HINT1(73)=D, HINT1(74)=W, HINT1(75)=FNORM | |
10 | C | |
11 | #include "hiparnt.inc" | |
12 | WDSAX1=HINT1(75)*(1.+HINT1(74)*(X/HINT1(72))**2)/ | |
13 | & (1+EXP((X-HINT1(72))/HINT1(73))) | |
14 | IF (HINT1(74).LT.0.) THEN | |
15 | IF (X.GE.HINT1(72)/SQRT(ABS(HINT1(74)))) WDSAX1=0. | |
16 | ENDIF | |
17 | RETURN | |
18 | END |