C
read (itorp,*,end=113,err=113) ntortyp
read (itorp,*,end=113,err=113) (itortyp(i),i=1,ntyp)
+ itortyp(ntyp1)=ntortyp+1
c write (iout,*) 'ntortyp',ntortyp
do i=1,ntortyp
do j=1,ntortyp
v0ijsccor=v0ijsccor+vlor1sccor(k,i,j)/
&(1+vlor3sccor(k,i,j)**2)
enddo
- v0sccor(i,j)=v0ijsccor
+ v0sccor(l,i,j)=v0ijsccor
enddo
enddo
enddo
c ee(1,2,i)=0.0d0
c ee(2,1,i)=ee(1,2,i)
enddo
+C write(iout,*) "parm", B1(1,nloctyp+1),B1(2,nloctyp+1)
if (lprint) then
- do i=1,nloctyp
+ do i=1,nloctyp+1
write (iout,*) 'Type',i
write (iout,*) 'B1'
write(iout,*) B1(1,i),B1(2,i)