diff --git a/src/solvation/cosmo_thomson.F b/src/solvation/cosmo_thomson.F
index 6ac2a80066..a10f5452c2 100644
--- a/src/solvation/cosmo_thomson.F
+++ b/src/solvation/cosmo_thomson.F
@@ -255,7 +255,7 @@ c     of a unit sphere. Spherical coordinates are used.
       enddo
 
       call ga_dgop(0,energy,1,'+')
-      call ga_dgop(0,grad,2*n,'+')
+      call ga_dgop(0,grad(1),2*n,'+')
       call dscal(2*n,0.5d0,grad,1)
 
       end subroutine