do j=1,ntortyp
do k=1,ntortyp
read (itordp,'(3a1)',end=114,err=114) t1,t2,t3
-<<<<<<< HEAD
c write (iout,*) "OK onelett",
c & i,j,k,t1,t2,t3
if (t1.ne.toronelet(i) .or. t2.ne.toronelet(j)
& .or. t3.ne.toronelet(k)) then
-=======
- if (t1.ne.onelett(i) .or. t2.ne.onelett(j)
- & .or. t3.ne.onelett(k)) then
->>>>>>> devel
write (iout,*) "Error in double torsional parameter file",
& i,j,k,t1,t2,t3
#ifdef MPI
if (lprint) then
write (iout,*)
write (iout,*) 'Constants for double torsionals'
-<<<<<<< HEAD
do iblock=1,2
do i=0,ntortyp-1
do j=-ntortyp+1,ntortyp-1
do k=-ntortyp+1,ntortyp-1
-=======
- do i=1,ntortyp
- do j=1,ntortyp
- do k=1,ntortyp
->>>>>>> devel
write (iout,*) 'ityp',i,' jtyp',j,' ktyp',k,
& ' nsingle',ntermd_1(i,j,k),' ndouble',ntermd_2(i,j,k)
write (iout,*)
C Modified 11 May 2012 by Adasko
CCC
C
-<<<<<<< HEAD
- read (isccor,*,end=119,err=119) nsccortyp
+ read (isccor,*,end=1113,err=1113) nsccortyp
#ifdef SCCORPDB
- read (isccor,*,end=119,err=119) (isccortyp(i),i=1,ntyp)
+ read (isccor,*,end=1113,err=1113) (isccortyp(i),i=1,ntyp)
do i=-ntyp,-1
isccortyp(i)=-isccortyp(-i)
enddo
iscprol=isccortyp(20)
-=======
- read (isccor,*,end=1113,err=1113) nsccortyp
- read (isccor,*,end=1113,err=1113) (isccortyp(i),i=1,ntyp)
->>>>>>> devel
c write (iout,*) 'ntortyp',ntortyp
maxinter=3
cc maxinter is maximum interaction sites
do l=1,maxinter
do i=1,nsccortyp
do j=1,nsccortyp
-<<<<<<< HEAD
- read (isccor,*,end=119,err=119) nterm_sccor(i,j),nlor_sccor(i,j)
-=======
read (isccor,*,end=1113,err=1113) nterm_sccor(i,j),
& nlor_sccor(i,j)
->>>>>>> devel
v0ijsccor=0.0d0
v0ijsccor1=0.0d0
v0ijsccor2=0.0d0
nterm_sccor(-i,-j)=nterm_sccor(i,j)
nterm_sccor(i,-j)=nterm_sccor(i,j)
do k=1,nterm_sccor(i,j)
-<<<<<<< HEAD
read (isccor,*,end=119,err=119) kk,v1sccor(k,l,i,j)
& ,v2sccor(k,l,i,j)
if (j.eq.iscprol) then
endif
endif
endif
-=======
read (isccor,*,end=1113,err=1113) kk,v1sccor(k,l,i,j)
& ,v2sccor(k,l,i,j)
->>>>>>> devel
v0ijsccor=v0ijsccor+si*v1sccor(k,l,i,j)
v0ijsccor1=v0ijsccor+si*v1sccor(k,l,-i,j)
v0ijsccor2=v0ijsccor+si*v1sccor(k,l,i,-j)
si=-si
enddo
do k=1,nlor_sccor(i,j)
-<<<<<<< HEAD
- read (isccor,*,end=119,err=119) kk,vlor1sccor(k,i,j),
-=======
read (isccor,*,end=1113,err=1113) kk,vlor1sccor(k,i,j),
->>>>>>> devel
& vlor2sccor(k,i,j),vlor3sccor(k,i,j)
v0ijsccor=v0ijsccor+vlor1sccor(k,i,j)/
&(1+vlor3sccor(k,i,j)**2)
c b1(1,i)=0.0d0
c b1(2,i)=0.0d0
B1tilde(1,i) = b(3)
-<<<<<<< HEAD
B1tilde(2,i) =-b(5)
B1tilde(1,-i) =-b(3)
B1tilde(2,-i) =b(5)
-=======
- B1tilde(2,i) =-b(5)
->>>>>>> devel
c b1tilde(1,i)=0.0d0
c b1tilde(2,i)=0.0d0
B2(1,i) = b(2)