include 'COMMON.FREE'
character*1 t1,t2,t3
character*1 onelett(4) /"G","A","P","D"/
- character*1 toronelet(-2:2)/"p","a","G","A","P"/
logical lprint
dimension blower(3,3,maxlob)
character*800 controlcard
wtor=ww(13)
wtor_d=ww(14)
wvdwpp=ww(16)
+ wstrain=ww(15)
wbond=ww(18)
wsccor=ww(19)
do j=-ntortyp+1,ntortyp-1
do k=-ntortyp+1,ntortyp-1
read (itordp,'(3a1)') 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
write (iout,*) "Error in double torsional parameter file",
& i,j,k,t1,t2,t3
stop "Error in double torsional parameter file"