c
#ifdef CRYST_BOND
read (ibond,*) vbldp0,vbldpdum,akp,mp,ip,pstok
c
#ifdef CRYST_BOND
read (ibond,*) vbldp0,vbldpdum,akp,mp,ip,pstok
& vbldsc0(j,i),aksc(j,i),abond0(j,i)
enddo
enddo
& vbldsc0(j,i),aksc(j,i),abond0(j,i)
enddo
enddo
+ if (lprint) then
+ write (iout,*) "Lipid transfer parameters"
+ write (iout,'(a5,f10.5)') "pept",pepliptran
+ do i=1,ntyp
+ write (iout,'(a5,f10.5)') restyp(i),liptranene(i)
+ enddo
+#ifdef AIX
+ call flush_(iout)
+#else
+ call flush(iout)
+#endif
+ endif
write (iout,'(a3,i4,2x,5(1pe14.5))') restyp(i),i,theta0(i),
& sig0(i),(gthet(j,i),j=1,3)
enddo
write (iout,'(a3,i4,2x,5(1pe14.5))') restyp(i),i,theta0(i),
& sig0(i),(gthet(j,i),j=1,3)
enddo
write (iout,'(a3,1h&,2x,5(f8.3,1h&))') restyp(i),theta0(i),
& 100*sig0(i),gthet(1,i)*0.1D0,gthet(2,i),gthet(3,i)*10.0D0
enddo
write (iout,'(a3,1h&,2x,5(f8.3,1h&))') restyp(i),theta0(i),
& 100*sig0(i),gthet(1,i)*0.1D0,gthet(2,i),gthet(3,i)*10.0D0
enddo
do i=1,ntyp
c write (2,*) 'i=',i
read (ithep_pdb,*,err=111,end=111)
do i=1,ntyp
c write (2,*) 'i=',i
read (ithep_pdb,*,err=111,end=111)
do i=1,ntyp
read (irotam_pdb,'(3x,i3,f8.3)',end=112,err=112) nlob(i),dsc(i)
if (i.eq.10) then
do i=1,ntyp
read (irotam_pdb,'(3x,i3,f8.3)',end=112,err=112) nlob(i),dsc(i)
if (i.eq.10) then
write (iout,'(6x,6f10.5)') (v2(k,i,j),k=1,nterm_old)
enddo
enddo
write (iout,'(6x,6f10.5)') (v2(k,i,j),k=1,nterm_old)
enddo
enddo
do i=1,ntortyp
do j=1,ntortyp
write (iout,*) 'ityp',i,' jtyp',j
do i=1,ntortyp
do j=1,ntortyp
write (iout,*) 'ityp',i,' jtyp',j
CCC
C
read (isccor,*,end=119,err=119) nsccortyp
CCC
C
read (isccor,*,end=119,err=119) nsccortyp
#ifdef SCCORPDB
read (isccor,*,end=119,err=119) (isccortyp(i),i=1,ntyp)
do i=-ntyp,-1
#ifdef SCCORPDB
read (isccor,*,end=119,err=119) (isccortyp(i),i=1,ntyp)
do i=-ntyp,-1
v0ijsccor=v0ijsccor+vlor1sccor(k,i,j)/
&(1+vlor3sccor(k,i,j)**2)
enddo
v0ijsccor=v0ijsccor+vlor1sccor(k,i,j)/
&(1+vlor3sccor(k,i,j)**2)
enddo
- write (iout,'(/a/)') 'Torsional constants:'
+ write (iout,'(/a/)') 'SC-torsional constants:'
+ do l=1,maxinter
+ write (iout,*) "Torsional type",l
do i=1,nsccortyp
do j=1,nsccortyp
write (iout,*) 'ityp',i,' jtyp',j
do i=1,nsccortyp
do j=1,nsccortyp
write (iout,*) 'ityp',i,' jtyp',j
write (iout,'(2f10.5)') EEold(j,1,i),EEold(j,2,i)
enddo
enddo
write (iout,'(2f10.5)') EEold(j,1,i),EEold(j,2,i)
enddo
enddo
write (iout,'(4f8.3,4e12.4)') eps_scp(i,1),rscp(i,1),
& eps_scp(i,2),rscp(i,2),aad(i,1),bad(i,1),aad(i,2),bad(i,2)
enddo
write (iout,'(4f8.3,4e12.4)') eps_scp(i,1),rscp(i,1),
& eps_scp(i,2),rscp(i,2),aad(i,1),bad(i,1),aad(i,2),bad(i,2)
enddo