C
read (itorp,*,end=113,err=113) ntortyp
read (itorp,*,end=113,err=113) (itortyp(i),i=1,ntyp)
+ itortyp(ntyp1)=ntortyp
do iblock=1,2
do i=-ntyp,-1
itortyp(i)=-itortyp(-i)
c ee(1,2,i)=0.0d0
c ee(2,1,i)=ee(1,2,i)
enddo
-c lprint=.true.
+ lprint=.true.
if (lprint) then
do i=1,nloctyp
write (iout,*) 'Type',i
enddo
enddo
endif
-c lprint=.false.
+ lprint=.false.
C
C Read electrostatic-interaction parameters
endif
enddo
enddo
+ write(iout,*) "tube param"
+ read(itube,*) epspeptube,sigmapeptube,tubetranenepep
+ sigmapeptube=sigmapeptube**6
+ sigeps=dsign(1.0D0,epspeptube)
+ epspeptube=dabs(epspeptube)
+ pep_aa_tube=4.0d0*epspeptube*sigmapeptube**2
+ pep_bb_tube=-sigeps*4.0d0*epspeptube*sigmapeptube
+ write(iout,*) pep_aa_tube,pep_bb_tube,tubetranenepep
+ do i=1,ntyp
+ read(itube,*) epssctube,sigmasctube,acavtub(i),bcavtub(i),
+ & ccavtub(i),dcavtub(i),tubetranene(i)
+ sigmasctube=sigmasctube**6
+ sigeps=dsign(1.0D0,epssctube)
+ epssctube=dabs(epssctube)
+ sc_aa_tube_par(i)=4.0d0*epssctube*sigmasctube**2
+ sc_bb_tube_par(i)=-sigeps*4.0d0*epssctube*sigmasctube
+ write(iout,*) sc_aa_tube_par(i), sc_bb_tube_par(i),tubetranene(i)
+ enddo
+
#ifdef OLDSCP
C
C Define the SC-p interaction constants (hard-coded; old style)