#else
read (isccor,*,end=119,err=119) (isccortyp(i),i=1,ntyp)
c write (iout,*) 'ntortyp',ntortyp
+ do i=-ntyp,-1
+ isccortyp(i)=-isccortyp(-i)
+ enddo
maxinter=3
cc maxinter is maximum interaction sites
+C NOW READING FOR LL aminoacid chirality
do l=1,maxinter
do i=1,nsccortyp
do j=1,nsccortyp
read (isccor,*,end=119,err=119)
& nterm_sccor(i,j),nlor_sccor(i,j)
+ nterm_sccor(-i,-j)=nterm_sccor(i,j)
+ v0ijsccor=0.0d0
+ si=-1.0d0
+
+ do k=1,nterm_sccor(i,j)
+ read (isccor,*,end=119,err=119) kk,v1sccor(k,l,i,j)
+ & ,v2sccor(k,l,i,j)
+ v0ijsccor=v0ijsccor+si*v1sccor(k,l,i,j)
+ si=-si
+ v1sccor(k,l,-i,-j)=v1sccor(k,l,i,j)
+ v2sccor(k,l,-i,-j)=-v2sccor(k,l,i,j)
+
+ enddo
+ do k=1,nlor_sccor(i,j)
+ read (isccor,*,end=119,err=119) kk,vlor1sccor(k,i,j),
+ & vlor2sccor(k,i,j),vlor3sccor(k,i,j)
+ v0ijsccor=v0ijsccor+vlor1sccor(k,i,j)/
+ &(1+vlor3sccor(k,i,j)**2)
+ enddo
+ v0sccor(l,i,j)=v0ijsccor
+ v0sccor(l,-i,-j)=v0ijsccor
+ enddo
+ enddo
+ enddo
+C NOW READING FOR LD aminoacid chirality
+ do l=1,maxinter
+ do i=1,nsccortyp
+ do j=-nsccortyp,-1
+ read (isccor,*,end=119,err=119)
+ & nterm_sccor(i,j),nlor_sccor(i,j)
+ nterm_sccor(-i,-j)=nterm_sccor(i,j)
v0ijsccor=0.0d0
si=-1.0d0
& ,v2sccor(k,l,i,j)
v0ijsccor=v0ijsccor+si*v1sccor(k,l,i,j)
si=-si
+ v1sccor(k,l,-i,-j)=v1sccor(k,l,i,j)
+ v2sccor(k,l,-i,-j)=-v2sccor(k,l,i,j)
+
enddo
do k=1,nlor_sccor(i,j)
read (isccor,*,end=119,err=119) kk,vlor1sccor(k,i,j),
&(1+vlor3sccor(k,i,j)**2)
enddo
v0sccor(l,i,j)=v0ijsccor
+ v0sccor(l,-i,-j)=v0ijsccor
enddo
enddo
enddo
+C symetry for GLY
+ i=10
+ do l=1,maxinter
+ do j=-nsccortyp,-1
+ do k=1,nterm_sccor(i,j)
+ v1sccor(k,l,10,j)=v1sccor(k,l,10,-j)
+ v2sccor(k,l,10,j)=-v2sccor(k,l,10,-j)
+ v1sccor(k,l,j,10)=v1sccor(k,l,-j,10)
+ v2sccor(k,l,j,10)=-v2sccor(k,l,-j,10)
+ enddo
+ enddo
+ enddo
+
close (isccor)
#endif
write (iout,*) 'Type',i
write (iout,'(a,i2,a,f10.5)') ('b(',ii,')=',b(ii,i),ii=1,13)
endif
-c B1(1,i) = b(3)
-c B1(2,i) = b(5)
-c B1(1,-i) = b(3)
-c B1(2,-i) = -b(5)
+ b(3,-i)= b(3,i)
+ b(5,-i)=-b(5,i)
+ b(4,-i)=-b(4,i)
+ b(2,-i)= b(2,i)
+ B1(1,i) = b(3,i)
+ B1(2,i) = b(5,i)
+ B1(1,-i) = b(3,i)
+ B1(2,-i) = -b(5,i)
c b1(1,i)=0.0d0
c b1(2,i)=0.0d0
-c B1tilde(1,i) = b(3)
-c B1tilde(2,i) =-b(5)
-c B1tilde(1,-i) =-b(3)
-c B1tilde(2,-i) =b(5)
+ B1tilde(1,i) = b(3,i)
+ B1tilde(2,i) =-b(5,i)
+ B1tilde(1,-i) =-b(3,i)
+ B1tilde(2,-i) =b(5,i)
c b1tilde(1,i)=0.0d0
c b1tilde(2,i)=0.0d0
-c B2(1,i) = b(2)
-c B2(2,i) = b(4)
-c B2(1,-i) =b(2)
-c B2(2,-i) =-b(4)
+ B2(1,i) = b(2,i)
+ B2(2,i) = b(4,i)
+ B2(1,-i) =b(2,i)
+ B2(2,-i) =-b(4,i)
cc B1tilde(1,i) = b(3,i)
cc B1tilde(2,i) =-b(5,i)
C B1tilde(1,-i) =-b(3,i)