inideflo.F

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file inideflo.F

*72*********************************************************************
      SUBROUTINE inideflo
************************      
*HI    Speed optimized *
************************      
       IMPLICIT NONE

*
* to avoid variable conflictions, a second keep element is necessary
* with the same common block name (see LEPTO2)
*

      COMMON /leptou/ cut(14),lst(40),parl(30),
     &                xlp,ylp,w2lp,q2lp,ulp
      REAL cut,parl,xlp,ylp,w2lp,q2lp,ulp
      INTEGER lst
      SAVE /leptou/

      INTEGER imxpdf
      parameter(imxpdf=40)
      COMMON /pepadm/cpdfnam(2,imxpdf),ipdfnam(2,imxpdf),
     &       iplst(10),cunpol,cpol
      CHARACTER*256 cpdfnam,cunpol,cpol
      INTEGER iplst,ipdfnam
      SAVE /pepadm/

**************************************************************
*      
*   IPLST(1) =  0 (default) : number of PDF warnings
*   IPLST(2) = 11 (default) : unit -1- for pdf files
*   IPLST(3) = 12 (default) : unit -2- for pdf files      
*
**************************************************************      


      INTEGER k,j,lx1,lq
      DOUBLE PRECISION xx,qq
      common/ xarray / xarray
      common/ grid / gu,gd,gl,gs,gg,gp,gn
      DOUBLE PRECISION xarray(102), gu(76,26), gd(76,26), gl(76,26),
     &              gs(76,26), gg(76,26), gp(76,26), gn(76,26) 

      OPEN(iplst(3),file=cpol,status='OLD')
      DO  k = 1, 76 
         DO  j = 1, 26
            READ(iplst(3),40)gu(k,j),gd(k,j),gl(k,j),gs(k,j), 
     &                 gg(k,j),gp(k,j),gn(k,j)
         END DO
      END DO
      CLOSE(iplst(3))
 40   FORMAT (7(1pe15.7))

      DO lx1=0,25
         xx=10.d0**(-4.d0+lx1/25.d0*3.d0)
         xarray(lx1+1) = dlog(xx)
      ENDDO
      DO lx1=1,50
         xx=0.1d0+lx1/50.d0*0.9d0
         xarray(lx1+26) = dlog(xx)
      ENDDO
      DO lq=0,25
         qq=0.6d0*10**(5.d0*lq/25.d0)
         xarray(lq+77) = dlog(qq)
      ENDDO
 
      RETURN
      END