PDF/spdf/aurpi1.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.2  1996/10/30 08:27:36  cernlib
   6 * Version 7.04
   7 *
   8 * Revision 1.1.1.1  1996/04/12 15:29:03  plothow
   9 * Version 7.01
  10 *
  11 *
  12 #include "pdf/pilot.h"
  13 C
  14       SUBROUTINE AURPI1(I,NDRV,X,S,ANS)
  15 #include "pdf/expdp.inc"
  16      +       CALCPIO(8,20,25)
  17       COMMON/W5051I4/CALCPIO
  18       REAL   F1(25),F2(25)
  19       DATA DELTA/.10/
  20       ANS=0.
  21       IF(X.GT.0.9985) RETURN
  22       IF(I.EQ.3.AND.X.GT.0.95) RETURN
  23       IF(I.EQ.8.AND.X.GT.0.95) RETURN
  24       IS=S/DELTA+1
  25       IS1=IS+1
  26       DO 1 L=1,25
  27       KL=L+NDRV*25
  28       F1(L)=CALCPIO(I,IS,KL)
  29       F2(L)=CALCPIO(I,IS1,KL)
  30     1 CONTINUE
  31       A1=AUGETFV(X,F1)
  32       A2=AUGETFV(X,F2)
  33       S1=(IS-1)*DELTA
  34       S2=S1+DELTA
  35       ANS=A1*(S-S2)/(S1-S2)+A2*(S-S1)/(S2-S1)
  36       RETURN
  37       END