Esope source files

   1 : C CALJXY    SOURCE    CHAT      05/01/12    21:47:03     5004           
   2 :       SUBROUTINE CALJXY(XYZ,P,IES,XY,ND,NP)                                     
   3 : C************************************************************************       
   4 : C                                                                               
   5 : C     CE SP CALCULE XY = P * XYZ (EN GROS)                                      
   6 : C                                                                               
   7 : C                                                                               
   8 : C                                                                               
   9 : C                                                                               
  10 : C************************************************************************       
  11 :       IMPLICIT INTEGER(I-N)
  12 :       IMPLICIT REAL*8 (A-H,O-Z)                                                 
  13 :       REAL*8 XYZ(IES,NP),P(IES,IES),XY(ND,NP)                                   
  14 :                                                                                 
  15 : C***                                                                            
  16 :                                                                                 
  17 :                                                                                 
  18 :       DO 43 I=1,NP                                                              
  19 :       DO 43 N=1,ND                                                              
  20 :       XY(N,I)=0.                                                                
  21 :       DO 43 M=1,IES                                                             
  22 :       XY(N,I)=XY(N,I)+P(N,M)*(XYZ(M,I)-XYZ(M,1))                                
  23 :  43   CONTINUE                                                                  
  24 :                                                                                 
  25 :                                                                                 
  26 : C     WRITE(6,*)'SUB CALJXY XYZ'                                                
  27 : C     WRITE(6,1983)XYZ                                                          
  28 : C1983 FORMAT(3(2X,1PE11.4))                                                     
  29 : C     WRITE(6,*)'SUB CALJCC XY'                                                 
  30 : C     WRITE(6,1982)XY                                                           
  31 : C1982 FORMAT(2(2X,1PE11.4))                                                     
  32 :                                                                                 
  33 :                                                                                 
  34 :       RETURN                                                                    
  35 :       END                                                                       
  36 :