X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc_MD-M%2Fenergy_p_new_barrier.F;h=4a47935f900c0056959556a27240aafcf00eb80d;hb=9ab6a6c94ffa93357b777d7a01615bc2cb964622;hp=aafba26456434bd64e88cdd7600098857971afde;hpb=1102b0fb298602dd1f056bfe4de1d42e350f16f3;p=unres.git diff --git a/source/unres/src_MD-M/energy_p_new_barrier.F b/source/unres/src_MD-M/energy_p_new_barrier.F index aafba26..4a47935 100644 --- a/source/unres/src_MD-M/energy_p_new_barrier.F +++ b/source/unres/src_MD-M/energy_p_new_barrier.F @@ -3593,7 +3593,7 @@ C I TU KURWA do j=ielstart(i),ielend(i) C do j=16,17 C write (iout,*) i,j - if (j.le.1) cycle +C if (j.le.1) cycle if (itype(j).eq.ntyp1.or. itype(j+1).eq.ntyp1 C changes suggested by Ana to avoid out of bounds c & .or.((j+2).gt.nres) @@ -3801,6 +3801,7 @@ cd & xmedi,ymedi,zmedi,xj,yj,zj write (iout,'(a6,2i5,0pf7.3,2i5,2e11.3)') &'evdw1',i,j,evdwij &,iteli,itelj,aaa,evdw1 + write (iout,*) sss write (iout,'(a6,2i5,0pf7.3,2f8.3)') 'ees',i,j,eesij, &fac_shield(i),fac_shield(j) endif