Esope source files

   1 : C VPAST2    SOURCE    CHAT      05/01/13    04:09:42     5004           
   2 :       SUBROUTINE VPAST2(VCORE,B)                                                
   3 : C=======================================================================        
   4 : C                                                                               
   5 : C     CONSTRUCTION DE LA MATRICE DE PASSAGE                                     
   6 : C    VCORE(3,2)=COORD DE LA LIGNE  (1 COL = 1 PT)                               
   7 : C SORTIES                                                                       
   8 : C    B(3,3)    =MATRICE DE PASSAGE                                              
   9 : C=======================================================================        
  10 :       IMPLICIT INTEGER(I-N)
  11 :       IMPLICIT REAL*8(A-H,O-Z)                                                  
  12 :       PARAMETER(XZER=0.D0,UN=1.D0)                                              
  13 :       DIMENSION VCORE(3,*),B(3,*)                                               
  14 : C                                                                               
  15 : C                                                                               
  16 :       X21=VCORE(1,2)-VCORE(1,1)                                                 
  17 :       Y21=VCORE(2,2)-VCORE(2,1)                                                 
  18 :       VL12=SQRT(X21*X21+Y21*Y21)                                                
  19 : C                                                                               
  20 :       VL12=UN/VL12                                                              
  21 :       B(1,1)=X21*VL12                                                           
  22 :       B(1,2)=Y21*VL12                                                           
  23 :       B(1,3)=XZER                                                               
  24 : C                                                                               
  25 :       B(3,1)= B(1,2)                                                            
  26 :       B(3,2)=-B(1,1)                                                            
  27 :       B(3,3)=XZER                                                               
  28 : C                                                                               
  29 :       B(2,1)=XZER                                                               
  30 :       B(2,2)=XZER                                                               
  31 :       B(2,3)=UN                                                                 
  32 : C                                                                               
  33 :       RETURN                                                                    
  34 :       END                                                                       
  35 :