]>
Commit | Line | Data |
---|---|---|
e74335a4 | 1 | * $Id$ |
2 | C | |
3 | C | |
4 | C THIS FUNCTION IS TO CALCULATE THE NUCLEAR PROFILE FUNCTION | |
5 | C OF THE COLLIDERING SYSTEM (IN UNITS OF 1/mb) | |
6 | C | |
7 | FUNCTION PROFILE(XB) | |
bc676b8e | 8 | #define BLANKET_SAVE |
e74335a4 | 9 | #include "pact_hijing.inc" |
10 | #include "hiparnt.inc" | |
11 | SAVE | |
12 | EXTERNAL FLAP, FLAP1, FLAP2 | |
13 | C | |
14 | BB=XB | |
15 | PROFILE=1.0 | |
16 | IF(IHNT2(1).GT.1 .AND. IHNT2(3).GT.1) THEN | |
17 | PROFILE=float(IHNT2(1))*float(IHNT2(3))*0.1* | |
18 | & GAUSS1(FLAP,0.0,HIPR1(34),0.01) | |
19 | ELSE IF(IHNT2(1).EQ.1 .AND. IHNT2(3).GT.1) THEN | |
20 | PROFILE=0.2*float(IHNT2(3))* | |
21 | & GAUSS1(FLAP2,0.0,HIPR1(35),0.001) | |
22 | ELSE IF(IHNT2(1).GT.1 .AND. IHNT2(3).EQ.1) THEN | |
23 | PROFILE=0.2*float(IHNT2(1))* | |
24 | & GAUSS1(FLAP1,0.0,HIPR1(34),0.001) | |
25 | ENDIF | |
26 | RETURN | |
27 | END |