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