double precision etot,evdw,evdw_t,evdw2,ees,evdw1,ebe,etors,
& escloc,ehpb,ecorr,ecorr5,ecorr6,eello_turn4,eello_turn3,
& eturn6,eel_loc,edihcnstr,etors_d,estr,evdw2_14,esccor,
- & eliptran
+ & eliptran,etube
integer ind_point(maxpoint),upindE,indE
character*16 plik
esccor=enetb(19,i,iparm)
edihcnstr=enetb(20,i,iparm)
eliptran=enetb(22,i,iparm)
+ etube=enetb(25,i,iparm)
+
#ifdef DEBUG
write (iout,'(3i5,6f5.2,14f12.3)') i,ib,iparm,(ft(l),l=1,6),
& +ft(2)*wturn3*eello_turn3
& +ft(5)*wturn6*eturn6+ft(2)*wel_loc*eel_loc
& +edihcnstr+ft(2)*wtor_d*etors_d+ft(1)*wsccor*esccor
- & +wbond*estr+wliptran*eliptran
+ & +wbond*estr+wliptran*eliptran+wtube*Etube
else
etot=wsc*(evdw+ft(6)*evdw_t)+wscp*evdw2+ft(1)*welec*ees
& +wvdwpp*evdw1
& +ft(2)*wturn3*eello_turn3
& +ft(5)*wturn6*eturn6+ft(2)*wel_loc*eel_loc
& +edihcnstr+ft(2)*wtor_d*etors_d+ft(1)*wsccor*esccor
- & +wbond*estr+wliptran*eliptran
+ & +wbond*estr+wliptran*eliptran+wtube*Etube
endif
#else
if (shield_mode.gt.0) then
& +ft(2)*wturn3*eello_turn3
& +ft(5)*wturn6*eturn6+ft(2)*wel_loc*eel_loc+edihcnstr
& +ft(2)*wtor_d*etors_d+ft(1)*wsccor*esccor
- & +wbond*estr+wliptran*eliptran
+ & +wbond*estr+wliptran*eliptran+wtube*Etube
else
etot=wsc*(evdw+ft(6)*evdw_t)+wscp*evdw2
& +ft(1)*welec*(ees+evdw1)
& +ft(2)*wturn3*eello_turn3
& +ft(5)*wturn6*eturn6+ft(2)*wel_loc*eel_loc+edihcnstr
& +ft(2)*wtor_d*etors_d+ft(1)*wsccor*esccor
- & +wbond*estr+wliptran*eliptran
+ & +wbond*estr+wliptran*eliptran+wtube*Etube
endif
#endif
edihcnstr=enetb(20,i,iparm)
C edihcnstr=0.0d0
eliptran=enetb(22,i,iparm)
+ etube=enetb(25,i,iparm)
+
#ifdef SPLITELE
if (shield_mode.gt.0) then
etot=ft(1)*wsc*(evdw+ft(6)*evdw_t)+ft(1)*wscp*evdw2
& +ft(2)*wturn3*eello_turn3
& +ft(5)*wturn6*eturn6+ft(2)*wel_loc*eel_loc
& +edihcnstr+ft(2)*wtor_d*etors_d+ft(1)*wsccor*esccor
- & +wbond*estr+wliptran*eliptran
+ & +wbond*estr+wliptran*eliptran+wtube*Etube
else
etot=wsc*(evdw+ft(6)*evdw_t)+wscp*evdw2+ft(1)*welec*ees
& +wvdwpp*evdw1
& +ft(2)*wturn3*eello_turn3
& +ft(5)*wturn6*eturn6+ft(2)*wel_loc*eel_loc
& +edihcnstr+ft(2)*wtor_d*etors_d+ft(1)*wsccor*esccor
- & +wbond*estr+wliptran*eliptran
+ & +wbond*estr+wliptran*eliptran+wtube*Etube
endif
#else
if (shield_mode.gt.0) then
& +ft(2)*wturn3*eello_turn3
& +ft(5)*wturn6*eturn6+ft(2)*wel_loc*eel_loc+edihcnstr
& +ft(2)*wtor_d*etors_d+ft(1)*wsccor*esccor
- & +wbond*estr+wliptran*eliptran
+ & +wbond*estr+wliptran*eliptran+wtube*Etube
else
etot=wsc*(evdw+ft(6)*evdw_t)+wscp*evdw2
& +ft(1)*welec*(ees+evdw1)
& +ft(2)*wturn3*eello_turn3
& +ft(5)*wturn6*eturn6+ft(2)*wel_loc*eel_loc+edihcnstr
& +ft(2)*wtor_d*etors_d+ft(1)*wsccor*esccor
- & +wbond*estr+wliptran*eliptran
+ & +wbond*estr+wliptran*eliptran+wtube*Etube
endif
#endif
edihcnstr=enetb(20,t,iparm)
C edihcnstr=0.0d0
eliptran=enetb(22,i,iparm)
+ etube=enetb(25,i,iparm)
do k=0,nGridT
betaT=startGridT+k*delta_T
& +ft(2)*wturn3*eello_turn3
& +ft(5)*wturn6*eturn6+ft(2)*wel_loc*eel_loc
& +edihcnstr+ft(2)*wtor_d*etors_d+ft(1)*wsccor*esccor
- & +wbond*estr+wliptran*eliptran
+ & +wbond*estr+wliptran*eliptran+wtube*Etube
eprim=ftprim(1)*(ft(6)*evdw_t+evdw)
C & +ftprim(6)*evdw_t
& +ftprim(1)*wscp*evdw2
& +ft(2)*wturn3*eello_turn3
& +ft(5)*wturn6*eturn6+ft(2)*wel_loc*eel_loc
& +edihcnstr+ft(2)*wtor_d*etors_d+ft(1)*wsccor*esccor
- & +wbond*estr+wliptran*eliptran
+ & +wbond*estr+wliptran*eliptran+wtube*Etube
eprim=ftprim(6)*evdw_t+ftprim(1)*welec*ees
& +ftprim(1)*wtor*etors+
& ftprim(3)*wcorr*ecorr+ftprim(4)*wcorr5*ecorr5+
& +ft(2)*wturn3*eello_turn3
& +ft(5)*wturn6*eturn6+ft(2)*wel_loc*eel_loc+edihcnstr
& +ft(2)*wtor_d*etors_d+ft(1)*wsccor*esccor
- & +wbond*estr+wliptran*eliptran
+ & +wbond*estr+wliptran*eliptran+wtube*Etube
eprim=ftprim(1)*(evdw+ft(6)*evdw_t)
& +ftprim(1)*welec*(ees+evdw1)
& +ftprim(1)*wtor*etors+
& +ft(2)*wturn3*eello_turn3
& +ft(5)*wturn6*eturn6+ft(2)*wel_loc*eel_loc+edihcnstr
& +ft(2)*wtor_d*etors_d+ft(1)*wsccor*esccor
- & +wbond*estr+wliptran*eliptran
+ & +wbond*estr+wliptran*eliptran+wtube*Etube
eprim=ftprim(6)*evdw_t+ftprim(1)*welec*(ees+evdw1)
& +ftprim(1)*wtor*etors+
& ftprim(3)*wcorr*ecorr+ftprim(4)*wcorr5*ecorr5+