working cutoff

[unres4.git] / source / wham / conform_compar.F90

diff --git a/source/wham/conform_compar.F90 b/source/wham/conform_compar.F90
index 1720af2..1d3a632 100644 (file)
--- a/source/wham/conform_compar.F90
+++ b/source/wham/conform_compar.F90
@@ -25,7 +25,7 @@
 !-----------------------------------------------------------------------------
 ! conf_compar.F
 !-----------------------------------------------------------------------------
-      subroutine conf_compar(jcon,lprn,print_class)
+      subroutine conf_compar(jcon)
 !      implicit real*8 (a-h,o-z)
       use energy_data, only:icont,ncont,nnt,nct,maxcont!,&
 !                      nsccont_frag_ref,isccont_frag_ref
@@ -65,9 +65,14 @@
                  nc_match,ncon_match,iclass_rms,ishifft_rms,ishiff,ishif
       integer :: k,kk,iclass_con,iscor,ik,ishifft_con,idig,iex,im
 !      print *,"Enter conf_compar",jcon
-      call angnorm12(rmsang)
+       lprn=.false.
+       print_class=.false.
+        write (iout,*) "before anything"
+        call flush(iout)
+
+!      call angnorm12(rmsang)
 ! Level 1: check secondary and supersecondary structure
-      call elecont(lprn,ncont,icont,nnt,nct)
+!      call elecont(lprn,ncont,icont,nnt,nct)
       if (lprn) then
         write (iout,*) "elecont finished"
         call flush(iout)
@@ -1435,7 +1440,7 @@
         ymedi=yi+0.5*dyi
         zmedi=zi+0.5*dzi
         do 4 j=i+2,ien-1
-          ind=ind+1
+!          ind=ind+1
           iteli=itel(i)
           itelj=itel(j)
           if (j.eq.i+2 .and. itelj.eq.2) iteli=2