X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Fwham%2Fsrc%2Fparmread.F;h=1b66ce01b81dca697fbfe632a0f2e7c3c10d8a0a;hb=0f16e6a86d4e9cfd5c8583879b6b97a35a2f0eda;hp=d8d7817b6b333dc3078c1c669119e81397ac3b22;hpb=d07177b41ad32fa404c1c6cff69df2a83d9cbbc8;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"