X-Git-Url: http://mmka.chem.univ.gda.pl/gitweb/?a=blobdiff_plain;f=source%2Funres%2Fsrc_MD%2Fenergy_p_new_barrier.F;h=a2eb0f9288190efdc4d89a714dd53ea4ad1afcad;hb=2d33f82339bda8759795452c7e56f5c4ecb31edc;hp=8a9a3eb2bd20a40fa2312d5ad22af33c15d83c0e;hpb=e123f68beb7954771018d0311468c68250b5f2d1;p=unres.git diff --git a/source/unres/src_MD/energy_p_new_barrier.F b/source/unres/src_MD/energy_p_new_barrier.F index 8a9a3eb..a2eb0f9 100644 --- a/source/unres/src_MD/energy_p_new_barrier.F +++ b/source/unres/src_MD/energy_p_new_barrier.F @@ -134,7 +134,7 @@ C cmc cmc Sep-06: egb takes care of dynamic ss bonds too cmc - if (dyn_ss) call dyn_set_nss +c if (dyn_ss) call dyn_set_nss c print *,"Processor",myrank," computed USCSC" #ifdef TIMING @@ -331,6 +331,7 @@ C energia(23)=evdw_m c print *," Processor",myrank," calls SUM_ENERGY" call sum_energy(energia,.true.) + if (dyn_ss) call dyn_set_nss c print *," Processor",myrank," left SUM_ENERGY" #ifdef TIMING #ifdef MPI @@ -1589,6 +1590,8 @@ C IF (dyn_ss_mask(i).and.dyn_ss_mask(j)) THEN call dyn_ssbond_ene(i,j,evdwij) evdw=evdw+evdwij + if (energy_dec) write (iout,'(a6,2i5,0pf7.3,a3)') + & 'evdw',i,j,evdwij,' ss' ELSE ind=ind+1 itypj=itype(j)