if ((wcorr4.eq.0.0d0 .and. wcorr.gt.0.0d0) .and. ipot.lt.6) then
call multibody_hb(ecorr,ecorr5,ecorr6,n_corr,n_corr1)
endif
-C
+C
+C if (wliptran.gt.0) then
+C call Eliptransfer(eliptran)
+C else
+C eliptran=0.0d0
+C endif
C If performing constraint dynamics, call the constraint energy
C after the equilibration time
if(usampl.and.totT.gt.eq_time) then
energia(20)=Uconst+Uconst_back
call sum_energy(energia,.true.)
C call enerprint
- write (iout,*) "Exit ETOTAL_LONG"
- call flush(iout)
+C write (iout,*) "Exit ETOTAL_LONG"
+C call flush(iout)
return
end
c------------------------------------------------------------------------------
else
esccor=0.0d0
endif
+ if (wliptran.gt.0) then
+ call Eliptransfer(eliptran)
+ else
+ eliptran=0.0d0
+ endif
+C print *,eliptran,wliptran
C
C Put energy components into an array
C
energia(17)=estr
energia(19)=edihcnstr
energia(21)=esccor
- write (iout,*) "ETOTAL_SHORT before SUM_ENERGY"
+ energia(22)=eliptran
+C write (iout,*) "ETOTAL_SHORT before SUM_ENERGY"
call flush(iout)
call sum_energy(energia,.true.)
- write (iout,*) "Exit ETOTAL_SHORT"
- call flush(iout)
+C write (iout,*) "Exit ETOTAL_SHORT"
+C call flush(iout)
return
end