first introduction of valence constrains - not working yet
[unres.git] / source / wham / src-M / readrtns.F
index 33924b5..5dd092e 100644 (file)
@@ -91,6 +91,8 @@
       zscfile=index(controlcard,"ZSCFILE").gt.0
       with_dihed_constr = index(controlcard,"WITH_DIHED_CONSTR").gt.0
       write (iout,*) "with_dihed_constr ",with_dihed_constr
+      with_theta_constr = index(controlcard,"WITH_THETA_CONSTR").gt.0
+      write (iout,*) "with_theta_constr ",with_theta_constr
       call readi(controlcard,'CONSTR_DIST',constr_dist,0)
       dyn_ss=index(controlcard,"DYN_SS").gt.0
       return