X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Fwham%2Fsrc%2Fparmread.F;h=1b66ce01b81dca697fbfe632a0f2e7c3c10d8a0a;hb=c5370e6422f68c03f5801e5e8b127a5ba478695b;hp=d8d7817b6b333dc3078c1c669119e81397ac3b22;hpb=d209eee3dec438628ee9d5a9c990ad9ca9dec1c5;p=unres.git

diff --git a/source/wham/src/parmread.F b/source/wham/src/parmread.F
index d8d7817..1b66ce0 100644
--- a/source/wham/src/parmread.F
+++ b/source/wham/src/parmread.F
@@ -23,7 +23,6 @@ C
       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
@@ -76,6 +75,7 @@ c
       wtor=ww(13)
       wtor_d=ww(14)
       wvdwpp=ww(16)
+      wstrain=ww(15)
       wbond=ww(18)
       wsccor=ww(19)
 
@@ -529,12 +529,12 @@ C
 C
 C 6/23/01 Read parameters for double torsionals
 C
-      do i=0,ntortyp-1
-        do j=-ntortyp,ntortyp
+      do i=1,ntortyp
+        do j=1,ntortyp
           do k=1,ntortyp
             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"