! -*- F90 -*-


      subroutine DOGevolvep0(xin,qin,p2in,ip2in,pdf) 
      include 'parmsetup.inc' 
      real*8 xin,qin,q2in,p2in,pdf(-6:6),xval(45),qcdl4,qcdl5 
      real*8 upv,dnv,usea,dsea,str,chm,bot,top,glu 
      character*16 name(nmxset) 
      integer nmem(nmxset),ndef(nmxset),mmem 
      common/NAME/name,nmem,ndef,mmem 
      integer ns 
                                                                        
      save 
                                                                        
      call DOPHO1(xin,qin,upv,dnv,usea,dsea,str,chm,bot,glu) 
                                                                        
      pdf(-6)= 0.0d0 
      pdf(6)= 0.0d0 
      pdf(-5)= bot 
      pdf(5 )= bot 
      pdf(-4)= chm 
      pdf(4 )= chm 
      pdf(-3)= str 
      pdf(3 )= str 
      pdf(-2)= usea 
      pdf(2 )= upv 
      pdf(-1)= dsea 
      pdf(1 )= dnv 
      pdf(0 )= glu 
                                                                        
      return 
!cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
      entry DOGevolvep1(xin,qin,p2in,ip2in,pdf) 
                                                                        
      call DOPHO2(xin,qin,upv,dnv,usea,dsea,str,chm,bot,glu) 
                                                                        
      pdf(-6)= 0.0d0 
      pdf(6)= 0.0d0 
      pdf(-5)= bot 
      pdf(5 )= bot 
      pdf(-4)= chm 
      pdf(4 )= chm 
      pdf(-3)= str 
      pdf(3 )= str 
      pdf(-2)= usea 
      pdf(2 )= upv 
      pdf(-1)= dsea 
      pdf(1 )= dnv 
      pdf(0 )= glu 
                                                                        
      return 
!cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
      entry DOGread(nset) 
      read(1,*)nmem(nset),ndef(nset) 
      return 
!                                                                       
!cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
      entry DOGalfa(alfas,qalfa) 
        call getnset(iset) 
        call GetOrderAsM(iset,iord) 
        call Getlam4M(iset,imem,qcdl4) 
        call Getlam5M(iset,imem,qcdl5) 
        call aspdflib(alfas,Qalfa,iord,qcdl5) 
      return 
!                                                                       
!cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
      entry DOGinit(Eorder,Q2fit) 
      return 
!                                                                       
!cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
      entry DOGpdf(mem) 
      imem = mem 
      return 
!                                                                       
 1000 format(5e13.5) 
      END                                           
!cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
       SUBROUTINE DOPHO1(DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DGL) 
!********************************************************************   
!*                                                                  *   
!*    Parametrization of parton distribution functions              *   
!*    in the photon (LO analysis) - asymptotic solution of AP eq.!  *   
!*                                                                  *   
!* authors:  D.Duke and H.Owens (DO)                                *   
!*           /Phys.Rev. D26 (1982) 1600/                            *   
!*                                                                  *   
!* Prepared by:                                                     *   
!*             Krzysztof Charchula, DESY                            *   
!*             bitnet: F1PCHA@DHHDESY3                              *   
!*             decnet: 13313::CHARCHULA                             *   
!*                                                                  *   
!* Modified by:                                                     *   
!*             H. Plothow-Besch/CERN-PPE                            *   
!*                                                                  *   
!********************************************************************   
!                                                                       
      implicit real*8 (a-h,o-z) 
       double precision                                                 &
     &        CQ(5),                                                    &
     &        DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DGL                     
      PARAMETER (ALPEM=7.29927D-3, PI=3.141592D0) 
      PARAMETER (ALAM=0.2D0) 
      DATA CQ/0.33333D0,0.66666D0,0.33333D0,0.66666D0,0.33333D0/ 
!                                                                       
       Q2 = DQ*DQ 
       ALAM2=ALAM**2 
       FQ=ALPEM/(2.*PI)*LOG(Q2/ALAM2) 
!                                                                       
!...gluons                                                              
       POMG=0.194*(1.-DX)**1.03/(DX**0.97) 
       DGL=POMG*FQ 
!                                                                       
!...quarks                                                              
        POM1=(1.81-1.67*DX+2.16*DX**2) 
        POM2=DX**0.7/(1.-0.4*LOG(1.-DX)) 
        POM3=38.D-4*(1.-DX)**1.82/(DX**1.18) 
          DDB=(CQ(1)**2*POM1*POM2+POM3)*FQ 
          DDV=DDB 
          DUB=(CQ(2)**2*POM1*POM2+POM3)*FQ 
          DUV=DUB 
          DSB=(CQ(3)**2*POM1*POM2+POM3)*FQ 
          DCB=(CQ(4)**2*POM1*POM2+POM3)*FQ 
          DBB=(CQ(5)**2*POM1*POM2+POM3)*FQ 
       RETURN 
      END                                           
!ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
       SUBROUTINE DOPHO2(DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DGL) 
!********************************************************************   
!*                                                                  *   
!*    Parametrization of parton distribution functions              *   
!*    in the photon (LO analysis) - asymptotic solution of AP eq.!  *   
!*                                                                  *   
!* authors:  D.Duke and H.Owens (DO)                                *   
!*           /Phys.Rev. D26 (1982) 1600/                            *   
!*                                                                  *   
!* Prepared by:                                                     *   
!*             Krzysztof Charchula, DESY                            *   
!*             bitnet: F1PCHA@DHHDESY3                              *   
!*             decnet: 13313::CHARCHULA                             *   
!*                                                                  *   
!* Modified by:                                                     *   
!*             H. Plothow-Besch/CERN-PPE                            *   
!*                                                                  *   
!********************************************************************   
!                                                                       
      implicit real*8 (a-h,o-z) 
      double precision                                                  &
     &        CQ(5),                                                    &
     &        DX,DQ,DUV,DDV,DUB,DDB,DSB,DCB,DBB,DGL                     
      PARAMETER (ALPEM=7.29927D-3,PI=3.141592D0) 
      PARAMETER (ALAM=0.4D0) 
      DATA CQ/0.33333D0,0.66666D0,0.33333D0,0.66666D0,0.33333D0/ 
!                                                                       
       Q2 = DQ*DQ 
       ALAM2=ALAM**2 
       FQ=ALPEM/(2.*PI)*LOG(Q2/ALAM2) 
!                                                                       
!...gluons                                                              
       POMG=0.194*(1.-DX)**1.03/(DX**0.97) 
       DGL=POMG*FQ 
!                                                                       
!...quarks                                                              
        POM1=(1.81-1.67*DX+2.16*DX**2) 
        POM2=DX**0.7/(1.-0.4*LOG(1.-DX)) 
        POM3=38.D-4*(1.-DX)**1.82/(DX**1.18) 
          DDB=(CQ(1)**2*POM1*POM2+POM3)*FQ 
          DDV=DDB 
          DUB=(CQ(2)**2*POM1*POM2+POM3)*FQ 
          DUV=DUB 
          DSB=(CQ(3)**2*POM1*POM2+POM3)*FQ 
          DCB=(CQ(4)**2*POM1*POM2+POM3)*FQ 
          DBB=(CQ(5)**2*POM1*POM2+POM3)*FQ 
       RETURN 
      END