X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Fwham%2Fsrc%2Fparmread.F;fp=source%2Fwham%2Fsrc%2Fparmread.F;h=1c1c818a7bc44f745d42a3565beb233491d7c7a7;hb=5b0e27e99ac1e61b8b155f3cc6351fac58c391ca;hp=29f1021d01a44fcdc924334117883b9d7774051d;hpb=3dad03d6cf13c486395188f78daff2121cf8391c;p=unres.git

diff --git a/source/wham/src/parmread.F b/source/wham/src/parmread.F
index 29f1021..1c1c818 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)
 
@@ -596,8 +596,8 @@ C
         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"