#ifdef MPL
c endif
#endif
-#define DEBUG
#ifdef DEBUG
call enerprint(energia,fact)
#endif
-#undef DEBUG
if (calc_grad) then
C
C Sum up the components of the Cartesian gradient.
include 'COMMON.NAMES'
include 'COMMON.FFIELD'
include 'COMMON.CONTROL'
+ include 'COMMON.TORCNSTR'
double precision coskt(mmaxtheterm),sinkt(mmaxtheterm),
& cosph1(maxsingle),sinph1(maxsingle),cosph2(maxsingle),
& sinph2(maxsingle),cosph1ph2(maxdouble,maxdouble),