Esope source files

   1 : C KDOM3B    SOURCE    CHAT      05/01/13    00:54:03     5004           
   2 :       SUBROUTINE KDOM3B(X1,Y1,X3,Y3,VOLU,X2,Y2)
   3 : C
   4 : C************************************************************************
   5 : C
   6 : C PROJET            :  CASTEM 2000
   7 : C
   8 : C NOM               :  KDOM3B
   9 : C                      
   9 : C DESCRIPTION       :  Subroutine called by KDOM3A and KDOM4C
  10 : C                      Axial-symmetric case
  11 : C                      We compute 
  12 : C                      the surface of the 2D line
  13 : C                      its center of gravity
  14 : C     
  15 : C LANGAGE           :  FORTRAN 77 + ESOPE 2000 (avec estensions CISI)
  16 : C      
  17 : C AUTEUR            :  A. BECCANTINI, DRN/DMT/SEMT/LTMF                        
  18 : C
  19 : C************************************************************************
  20 : C
  21 : C INPUT             : (X1,Y1), (X3,Y3) coordinates of the line extrema
  22 : C 
  23 : C OUTPUT            : VOLUME: surface
  24 : C                     (X2,Y2) center of gravity
  25 : C
  26 : C************************************************************************
  27 : C
  28 : C Created the 24/02/04
  29 : C
  30 : C
  31 :       IMPLICIT INTEGER(I-N)
  32 :       REAL*8 X1,Y1,X3,Y3,X2,Y2,RL13,VOLU,W1,W2
  33 : C
  34 :       W1=1.0D0-(1.0D0/(3.0D0**0.5D0))
  35 :       W2=1.0D0+(1.0D0/(3.0D0**0.5D0))
  36 : C
  37 :       RL13 = ((Y3-Y1)*(Y3-Y1)) + ((X3-X1)*(X3-X1))
  38 :       RL13 = RL13**0.5D0
  39 :       VOLU=0.5D0*(X1+X3)*RL13
  40 :       X2=0.125D0*RL13*(((X1*W1+X3*W2)**2)+((X1*W2+X3*W1)**2))
  41 :       X2=X2/VOLU
  42 :       Y2=0.125D0*RL13*(((X1*W1+X3*W2)*(Y1*W1+Y3*W2))+
  43 :      &     ((X3*W1+X1*W2)*(Y3*W1+Y1*W2)))
  44 :       Y2=Y2/VOLU
  45 :       RETURN 
  46 :       END
  47 :       
  48 : 
  49 :  
  50 :