+ if (shield_mode.gt.0) then
+ etot=ft(1)*wsc*(evdw+ft(6)*evdw_t)+ft(1)*wscp*evdw2
+ & +welec*ft(1)*(ees+evdw1)
+ & +wang*ebe+wtor*ft(1)*etors+wscloc*escloc
+ & +wstrain*ehpb+wcorr*ft(3)*ecorr+wcorr5*ft(4)*ecorr5
+ & +wcorr6*ft(5)*ecorr6+wturn4*ft(3)*eello_turn4
+ & +wturn3*ft(2)*eello_turn3+wturn6*ft(5)*eturn6
+ & +wel_loc*ft(2)*eel_loc+edihcnstr+wtor_d*ft(2)*etors_d
+ & +wbond*estr+wsccor*ft(1)*esccor+ehomology_constr
+ & +wliptran*eliptran+wsaxs*esaxs
+ else
+ etot=wsc*(evdw+ft(6)*evdw_t)+wscp*evdw2
+ & +welec*ft(1)*(ees+evdw1)
+ & +wang*ebe+wtor*ft(1)*etors+wscloc*escloc
+ & +wstrain*ehpb+wcorr*ft(3)*ecorr+wcorr5*ft(4)*ecorr5
+ & +wcorr6*ft(5)*ecorr6+wturn4*ft(3)*eello_turn4
+ & +wturn3*ft(2)*eello_turn3+wturn6*ft(5)*eturn6
+ & +wel_loc*ft(2)*eel_loc+edihcnstr+wtor_d*ft(2)*etors_d
+ & +wbond*estr+wsccor*ft(1)*esccor+!ethetacnstr
+ & +wliptran*eliptran+wsaxs*esaxs
+ endif