c lprn=.true.
c write (iout,*) "EBACK_SC_COR",iphi_start,iphi_end,nterm_sccor
esccor=0.0D0
- do i=iphi_start,iphi_end
+ do i=iphi_start-1,iphi_end+1
esccor_ii=0.0D0
isccori=isccortyp(itype(i-2))
isccori1=isccortyp(itype(i-1))
& (itype(i-1).eq.10))
& .or. ((intertyp.eq.1).and.(itype(i-2).ne.10))
& .or. ((intertyp.eq.2).and.(itype(i-1).ne.10))) cycle
+ if ((intertyp.eq.2).and.(i.eq.iphi_start-1)) cycle
+ if ((intertyp.eq.1).and.(i.eq.iphi_end+1)) cycle
do j=1,nterm_sccor(isccori,isccori1)
v1ij=v1sccor(j,intertyp,isccori,isccori1)
v2ij=v2sccor(j,intertyp,isccori,isccori1)
esccor=esccor+v1ij*cosphi+v2ij*sinphi
gloci=gloci+j*(v2ij*cosphi-v1ij*sinphi)
enddo
- gloc_sc(intertyp,i-3,icg)=gloc_sc(i-3,icg)+wtor*gloci
+ gloc_sc(intertyp,i-3,icg)=gloc_sc(intertyp,i-3,icg)+wsccor*gloci
if (lprn)
& write (iout,'(2(a3,2x,i3,2x),2i3,6f8.3/26x,6f8.3/)')
& restyp(itype(i-2)),i-2,restyp(itype(i-1)),i-1,itori,itori1,