Procedures | Cast3M

   1 : * @COUTOR1  PROCEDUR  AM        92/10/23    21:15:12     731            
   2 : *-------------------------------------------------                              
   3 : ******          PROCEDURE @COUTOR1          ******
   4 : *-------------------------------------------------                              
   5 : *
   6 : *
   7 : *    CETTE PROCEDURE A ETE MISE GRACIEUSEMENT
   8 : *   A DISPOSOTION DE LA COMMUNAUTE  CASTEM2000
   9 : *    PAR M. LIBEYRE ( CEA/DSM/DRFC )
  10 : *
  11 : *     TEL : ( 33 1 ) 42 25 46 03
  12 : *
  13 : *-------------------------------------------------                              
  14 : DEBPROC @COUTOR1 IEL*MAILLAGE
  15 :                 CHT*CHPOINT CHN*CHPOINT CHB*CHPOINT ;
  16 : PI = 3.14159 ;
  17 : P2 = IEL POIN INITIAL ; 
  18 : P3 = IEL POIN FINAL ;
  19 : DIMGEO1 = VALEUR DIME ;
  20 : SI (DIMGEO1 > 2);
  21 :  NXI2 = EXTR CHN NX P2 ;
  22 :  NYI2 = EXTR CHN NY P2 ;  
  23 :  NZI2 = EXTR CHN NZ P2 ;
  24 :  VN2 = NXI2 NYI2 NZI2 ;  
  25 :  TXI2 = EXTR CHT TX P2 ;
  26 :  TYI2 = EXTR CHT TY P2 ;  
  27 :  TZI2 = EXTR CHT TZ P2 ;
  28 :  VT2 = TXI2 TYI2 TZI2 ;  
  29 :  BXI2 = EXTR CHB BX P2 ;
  30 :  BYI2 = EXTR CHB BY P2 ;  
  31 :  BZI2 = EXTR CHB BZ P2 ;
  32 :  VB2 = BXI2 BYI2 BZI2 ;  
  33 :  NXI3 = EXTR CHN NX P3 ;
  34 :  NYI3 = EXTR CHN NY P3 ;  
  35 :  NZI3 = EXTR CHN NZ P3 ;
  36 :  VN3 = NXI3 NYI3 NZI3 ;  
  37 :  TXI3 = EXTR CHT TX P3 ;
  38 :  TYI3 = EXTR CHT TY P3 ;  
  39 :  TZI3 = EXTR CHT TZ P3 ;
  40 :  VT3 = TXI3 TYI3 TZI3 ;
  41 :  BXI3 = EXTR CHB BX P3 ;
  42 :  BYI3 = EXTR CHB BY P3 ;  
  43 :  BZI3 = EXTR CHB BZ P3 ;
  44 :  VB3 = BXI3 BYI3 BZI3 ;
  45 :  XR = VT3 PSCA VT2 ;  
  46 :  YR = VT3 PSCA VN2 ;
  47 : * MESS 'XR =' XR   'YR =' YR ;
  48 :  ALPHA = ATG YR XR ;
  49 : * MESS 'ALPHA =' ALPHA ;
  50 :  DS2 = NORM (MOIN P2 P3) ;
  51 :  RR = (ALPHA*PI/180.)/DS2 ;
  52 :  SI (RR NEG 0.) ;
  53 :   R = 1./RR ;
  54 :  SINON ;
  55 :   R = 1.E99 ;
  56 :  FINSI ;
  57 : * MESS 'R =' R ;
  58 :  XT = VB2 PSCA VB3 ;
  59 :  YT = VB2 PSCA VN3 ;
  60 : * MESS 'XT =' XT   'YT =' YT ;
  61 :  BETA = ATG YT XT ;  
  62 : * MESS 'BETA =' BETA ;
  63 :  TT = -1*(BETA*PI/180.)/DS2 ;
  64 :  SI (TT NEG 0.) ;
  65 :   T = 1./TT ;
  66 :  SINON ;
  67 :   T = 1.E99 ;
  68 :  FINSI ;
  69 : * MESS 'T =' T ;
  70 : 
  71 : SINON;  
  72 :  NXI2 = EXTR CHN NX P2 ;
  73 :  NYI2 = EXTR CHN NY P2 ;  
  74 :  VN2 = NXI2 NYI2 ;  
  75 :  TXI2 = EXTR CHT TX P2 ;
  76 :  TYI2 = EXTR CHT TY P2 ;  
  77 :  VT2 = TXI2 TYI2 ;  
  78 :  BXI2 = EXTR CHB BX P2 ;
  79 :  BYI2 = EXTR CHB BY P2 ;  
  80 :  VB2 = BXI2 BYI2 ;  
  81 :  NXI3 = EXTR CHN NX P3 ;
  82 :  NYI3 = EXTR CHN NY P3 ;  
  83 :  VN3 = NXI3 NYI3 ;  
  84 :  TXI3 = EXTR CHT TX P3 ;
  85 :  TYI3 = EXTR CHT TY P3 ;  
  86 :  VT3 = TXI3 TYI3 ;
  87 :  BXI3 = EXTR CHB BX P3 ;
  88 :  BYI3 = EXTR CHB BY P3 ;  
  89 :  VB3 = BXI3 BYI3 ;
  90 :  XR = VT3 PSCA VT2 ;  
  91 :  YR = VT3 PSCA VN2 ;
  92 : * MESS 'XR =' XR   'YR =' YR ;
  93 :  ALPHA = ATG YR XR ;
  94 : * MESS 'ALPHA =' ALPHA ;
  95 :  DS2 = NORM (MOIN P2 P3) ;
  96 :  RR = (ALPHA*PI/180.)/DS2 ;
  97 :  SI (RR NEG 0.) ;
  98 :   R = 1./RR ;
  99 :  SINON ;
 100 :   R = 1.E99 ;
 101 :  FINSI ;
 102 : * MESS 'R =' R ;
 103 :  BETA = 0.;
 104 :  T = 0.;
 105 : FINSI;
 106 : FINPROC DS2 R T ALPHA BETA ;