1 subroutine mdiffn(ib,ie,il,dlyr)
7 implicit real*8 (a-h,o-z)
10 real*8 fiit(4,2*nsz,nph), fiiz(4,2*nsz,jpart,nph)
16 fiit(is,kk,ll)=fiic(is,kk,ll)
20 fiit(1,kk,ll)=fiic(1,kk,nophi-ll+2)
21 fiit(2,kk,ll)=fiic(2,kk,nophi-ll+2)
22 fiit(3,kk,ll)=-fiic(3,kk,nophi-ll+2)
23 fiit(4,kk,ll)=-fiic(4,kk,nophi-ll+2)
31 mmp=nophi*(1+(mmp-1)/nophi)-mmp+2
34 fiiz(is,kk,ip,ll)=fiit(is,kk,mmp)
46 ppin(ij,ll,ii,kk)=ppin(ij,nophi-ll+2,ii,kk)
47 ppin(ij+16,ll,ii,kk)=ppin(ij+16,nophi-ll+2,ii
49 ppin(ij,ll,ii,kk)=-ppin(ij,nophi-ll+2,ii,kk)
50 ppin(ij+16,ll,ii,kk)=-ppin(ij+16,nophi-ll+2,ii
68 prod1=prod1+(c(kk)*tmsl*ppin(ij,ll,it,kk)+
69 1 trsl*ppin(ij+16,ll,it,kk))*fiiz(j,kk,ip,ll)
73 sumta=sumta+sumtb*dcmu(kk)
75 ftmp(is,it,ip)=sumta*ddphi