1 subroutine SMRSPevolve(xin,qin,pdf)
2 include 'parmsetup.inc'
5 real*8 xin,qin,pdf(-6:6),xval(45),qcdl4,qcdl5
6 real*8 upv,dnv,usea,dsea,str,chm,bot,top,glu
10 common /SMRSP/ F(7,NX,NQ,3)
11 character*16 name(nmxset)
12 integer nmem(nmxset),ndef(nmxset),mmem
13 common/NAME/name,nmem,ndef,mmem
18 call getnmem(iset,imem)
21 if(iimem.eq.0) iimem = 2
23 call SMRSPxx(iimem,xin,qin,upv,dnv,usea,str,chm,bot,glu)
42 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
44 read(1,*)nmem(nset),ndef(nset)
48 read(1,*)(F(m,k,l,j),m=1,5),F(7,k,l,j),F(6,k,l,j)
54 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
55 entry SMRSPalfa(alfas,qalfa)
57 call getnmem(iset,imem)
58 call GetOrderAsM(iset,iord)
59 call Getlam4M(iset,imem,qcdl4)
60 call Getlam5M(iset,imem,qcdl5)
61 call aspdflib(alfas,Qalfa,iord,qcdl5)
64 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
65 entry SMRSPinit(Eorder,Q2fit)
68 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
72 call setnmem(iset,mem)
77 ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
79 SUBROUTINE SMRSPxx(iset,X,SCALE,UPV,DNV,SEA,STR,CHM,BOT,GLU)
81 C :::::::::::: PION STRUCTURE FUNCTION :: 10% SEA :::::::::::::::::
83 implicit real*8 (a-h,o-z)
87 common /SMRSP/ F(7,NX,NQ,3)
88 DIMENSION G(7),XX(NX),N0(7)
89 c DIMENSION F(7,NX,NQ),G(7),XX(NX),N0(7)
90 DATA XX/1.D-5,2.D-5,4.D-5,6.D-5,8.D-5,
91 . 1.D-4,2.D-4,4.D-4,6.D-4,8.D-4,
92 . 1.D-3,2.D-3,4.D-3,6.D-3,8.D-3,
93 . 1.D-2,2.D-2,4.D-2,6.D-2,8.D-2,
94 . .1D0,.125D0,.15D0,.175D0,.2D0,.225D0,.25D0,.275D0,
95 . .3D0,.325D0,.35D0,.375D0,.4D0,.425D0,.45D0,.475D0,
96 . .5D0,.525D0,.55D0,.575D0,.6D0,.65D0,.7D0,.75D0,
97 . .8D0,.85D0,.9D0,.95D0,.975D0,1.D0/
98 DATA XMIN,XMAX,QSQMIN,QSQMAX/1.D-5,1.D0,5.D0,1310720.D0/
99 DATA N0/0,0,3,5,0,5,0/
100 DATA ZEROD/0.D0/,PONED/0.1D0/,ONED/1.D0/,ONEDO/1.1D0/,TWOD/2.D0/
103 IF(INIT.NE.0) GOTO 10
108 25 F(I,N,M,iset)=F(I,N,M,iset)/(ONED-XX(N))**N0(I)
111 XX(J)= LOG10(XX(J))+ONEDO
114 30 F(I,J,K,iset)= LOG(F(I,J,K,iset))
115 + *F(I,NTENTH,K,iset)/ LOG(F(I,NTENTH,K,iset))
119 40 F(I,NX,M,iset)=ZEROD
124 IF(QSQ.LT.QSQMIN) QSQ=QSQMIN
125 IF(QSQ.GT.QSQMAX) QSQ=QSQMAX
127 IF(X.LT.PONED) XXX= LOG10(X)+ONEDO
130 IF(XXX.GT.XX(N+1)) GOTO 70
131 A=(XXX-XX(N))/(XX(N+1)-XX(N))
132 RM= LOG(QSQ/QSQMIN)/ LOG(TWOD)
136 G(I)= (ONED-A)*(ONED-B)*F(I,N,M,iset)+(ONED-A)*B*F(I,N,M+1,iset)
137 . + A*(ONED-B)*F(I,N+1,M,iset) + A*B*F(I,N+1,M+1,iset)
138 IF(N.GE.NTENTH) GOTO 65
139 IF(I.EQ.7.OR.I.EQ.1) GOTO 65
140 FAC=(ONED-B)*F(I,NTENTH,M,iset)+B*F(I,NTENTH,M+1,iset)
143 G(I)=G(I)*(ONED-X)**N0(I)
145 C UPBAR DISTRIBUTION = D DISTRIBUTION
148 C THIS SEA IS (UBAR+DBAR)/2