--- source/energy/energy.src.orig	2003-08-15 15:00:43.000000000 -0700
+++ source/energy/energy.src	2005-04-13 10:08:18.000000000 -0700
@@ -33,6 +33,9 @@
       LOGICAL QSECD
       REAL*8 EIMSLV                             !##ASPENER
 C
+##IF APBS
+##INCLUDE '~/charmm_fcm/apbs.fcm'
+##ENDIF
 ##INCLUDE '~/charmm_fcm/block.fcm'
 ##INCLUDE '~/charmm_fcm/lambda.fcm'
 ##INCLUDE '~/charmm_fcm/cnst.fcm'
@@ -1922,6 +1925,20 @@
       ENDIF                                                  !##MTS
 ##ENDIF
 ##ENDIF
+##IF APBS
+C    APBS calculated forces
+      IF(ENE3) THEN                                          !##MTS
+        IF(MYNOD.EQ.INODE(20)) THEN  !##PARALLEL
+          IF (QETERM(PBELEC) .AND. QETERM(PBNP) .AND. QFAPBS ) THEN
+            CALL apbsfrc(NATOM,X,Y,Z,CG,ETERM(PBELEC),ETERM(PBNP),
+     $                   DX,DY,DZ,ecalls, .false.)
+            IF (TIMER.GT.1)
+     $         CALL WRTTIM(' APBS Solvation energy times:')
+          ENDIF
+        ENDIF                        !##PARALLEL
+      ENDIF                                                  !##MTS
+
+##ENDIF
 C
 C-----------------------------------------------------------------------
 ##IF PRIMSH
