C C $Id: getrforc.F,v 1.2 1998/07/16 16:39:57 jjv5 Exp arjan $ C C------------------------------------------------------------------------ subroutine getrforc implicit double precision(a-h,o-z) #include "divcon.dim" #include "divcon.h" do ir=1,nres ibeg = irpnt(ir) iend = irpnt(ir+1)-1 irx = 3*ir-2 iry = irx+1 irz = iry+1 rforcmc(irx) = 0.0 rforcmc(iry) = 0.0 rforcmc(irz) = 0.0 do i=ibeg,iend i3 = 3*i-2 rforcmc(irx) = rforcmc(irx) - grad(i3) rforcmc(iry) = rforcmc(iry) - grad(i3+1) rforcmc(irz) = rforcmc(irz) - grad(i3+2) enddo enddo end