1      subroutine argos_cafe_mcti_w(isat,isga,vs,nums,wgt,dera)
2c
3      implicit none
4c
5#include "argos_cafe_common.fh"
6c
7      integer nums
8      integer isat(msa),isga(msa)
9      real*8 vs(msa,3),wgt(mat,mset)
10      real*8 dera(6,nsatot)
11c
12      integer i
13      real*8 vsq
14c
15      do 1 i=1,nums
16      vsq=half*(vs(i,1)*vs(i,1)+vs(i,2)*vs(i,2)+vs(i,3)*vs(i,3))
17      deriv(13,1)=deriv(13,1)-wgt(isat(i),4)*vsq
18      if(npgdec.gt.1) dera(6,isga(i))=dera(6,isga(i))-wgt(isat(i),4)*vsq
19      if(ipert2.and.ip2(13)) ep2(3)=ep2(3)+
20     + vsq*wgt(isat(i),1)*(wgt(isat(i),1)/wgt(isat(i),2)-one)
21      if(ipert3.and.ip3(13)) ep3(3)=ep3(3)+
22     + vsq*wgt(isat(i),1)*(wgt(isat(i),1)/wgt(isat(i),3)-one)
23    1 continue
24c
25      return
26      end
27c $Id$
28