]>
Commit | Line | Data |
---|---|---|
21886bb6 | 1 | * |
2 | * $Id$ | |
3 | * | |
4 | * $Log$ | |
5 | * Revision 1.1.1.2 1996/10/30 08:27:32 cernlib | |
6 | * Version 7.04 | |
7 | * | |
8 | * Revision 1.1.1.1 1996/04/12 15:29:00 plothow | |
9 | * Version 7.01 | |
10 | * | |
11 | * | |
12 | #include "pdf/pilot.h" | |
13 | SUBROUTINE AFGIN2(I,NDRV,X,S,ANS) | |
14 | DIMENSION F1(32),F2(32),F3(32) | |
15 | DIMENSION AF(3),AS(3) | |
16 | DIMENSION CELCO(8,20,32) | |
17 | COMMON/W5051IB/CELCO | |
18 | DATA DELTA/0.8000E-01/ | |
19 | ANS=0. | |
20 | IF(X.GT.0.9985) RETURN | |
21 | N=3 | |
22 | IS=S/DELTA+1 | |
23 | C IF(IS.GE.17) IS=17 | |
24 | IS1=IS+1 | |
25 | IS2=IS1+1 | |
26 | DO 1 L=1,32 | |
27 | KL=L+32*NDRV | |
28 | F1(L)=CELCO(I,IS,KL) | |
29 | F2(L)=CELCO(I,IS1,KL) | |
30 | F3(L)=CELCO(I,IS2,KL) | |
31 | 1 CONTINUE | |
32 | AF(1)=AFGETFV(X,F1) | |
33 | AF(2)=AFGETFV(X,F2) | |
34 | AF(3)=AFGETFV(X,F3) | |
35 | AS(1)=(IS-1)*DELTA | |
36 | AS(2)=AS(1)+DELTA | |
37 | AS(3)=AS(2)+DELTA | |
38 | CALL AFPOLIN(AS,AF,N,S,AANS,DY) | |
39 | ANS=AANS | |
40 | RETURN | |
41 | END |