X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc_MD-M%2Finitialize_p.F;h=15f1e197831f698be3ec6c2d470d823029bd0af3;hb=03a99f9b2b61db5e2be5f20b3d5d912d5394af6c;hp=6c5f4eb343d5dd789cc46ab3001851490edd043a;hpb=42b30ef61fe3a8077aee8292bd43e4de4e6b8522;p=unres.git

diff --git a/source/unres/src_MD-M/initialize_p.F b/source/unres/src_MD-M/initialize_p.F
index 6c5f4eb..15f1e19 100644
--- a/source/unres/src_MD-M/initialize_p.F
+++ b/source/unres/src_MD-M/initialize_p.F
@@ -630,12 +630,13 @@ c     &  " ivec_start",ivec_start," ivec_end",ivec_end
         call int_bounds(ndih_constr,idihconstr_start,idihconstr_end)
       endif
       if (ntheta_constr.eq.0) then
-        idihconstr_start=1
-        idihconstr_end=0
+        ithetaconstr_start=1
+        ithetaconstr_end=0
       else
         call int_bounds
      &  (ntheta_constr,ithetaconstr_start,ithetaconstr_end)
       endif
+C        print *,ntheta_constr,ithetaconstr_start,ithetaconstr_end
 c      nsumgrad=(nres-nnt)*(nres-nnt+1)/2
 c      nlen=nres-nnt+1
       nsumgrad=(nres-nnt)*(nres-nnt+1)/2
@@ -670,7 +671,7 @@ c      nlen=nres-nnt+1
      & ' ivec_start',ivec_start,' ivec_end',ivec_end,
      & ' iset_start',iset_start,' iset_end',iset_end,
      & ' idihconstr_start',idihconstr_start,' idihconstr_end',
-     &   idihconstr_end
+     &   idihconstr_end,
      & ' ithetaconstr_start',ithetaconstr_start,' ithetaconstr_end',
      &   ithetaconstr_end