Esope source files

   1 : C RICOQ0    SOURCE    CHAT      05/01/13    02:56:14     5004           
   2 :       SUBROUTINE RICOQ0(RIG,X13,X23,Y13,Y23,FAC)                                
   3 :       IMPLICIT INTEGER(I-N)
   4 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
   5 :       DIMENSION RIG(18,18),VEC(18)                                              
   6 : *                                                                               
   7 :       S=(X13*Y23-X23*Y13)*2.D0                                                  
   8 :       FAC4=FAC*4.D0                                                             
   9 :       DO 100 I=1,18                                                             
  10 : 100   VEC(I)=0.                                                                 
  11 :       VEC(1)=-X23/S                                                             
  12 :       VEC(2)=-Y23/S                                                             
  13 :       VEC(7)=X13/S                                                              
  14 :       VEC(8)=Y13/S                                                              
  15 :       VEC(13)=(X23-X13)/S                                                       
  16 :       VEC(14)=(Y23-Y13)/S                                                       
  17 :       DO 102 I=1,3                                                              
  18 :       I6=6*I                                                                    
  19 :       VEC(6) =0.D0                                                              
  20 :       VEC(12)=0.D0                                                              
  21 :       VEC(18)=0.D0                                                              
  22 :       VEC(I6)=1.D0                                                              
  23 :       DO 103 J=1,18                                                             
  24 :       DO 103 K=1,18                                                             
  25 :       RIG(J,K)=RIG(J,K)+VEC(J)*VEC(K)*FAC4                                      
  26 :  103  CONTINUE                                                                  
  27 :  102  CONTINUE                                                                  
  28 : *                                                                               
  29 :       RETURN                                                                    
  30 :       END                                                                       
  31 :