Esope source files

   1 : C KDOM4B    SOURCE    CHAT      05/01/13    00:54:16     5004           
   2 :       SUBROUTINE KDOM4B(X1,Y1,X2,Y2,X3,Y3,VOLU,SURF,XCEN,YCEN,XC2D,YC2D)
   3 : C
   4 : C************************************************************************
   5 : C
   6 : C PROJET            :  CASTEM 2000
   7 : C
   8 : C NOM               :  KDOM4B
   9 : C                      
   9 : C DESCRIPTION       :  Subroutine called by KDOM4A
  10 : C                      Axial-symmetric case
  11 : C                      We compute 
  12 : C                      the surface of the 2D mesh
  13 : C                      the 3D volume of the 3D mesh
  14 : C                      its center of gravity
  15 : C     
  16 : C LANGAGE           :  FORTRAN 77 + ESOPE 2000 (avec estensions CISI)
  17 : C      
  18 : C AUTEUR            :  A. BECCANTINI, DRN/DMT/SEMT/LTMF                        
  19 : C
  20 : C************************************************************************
  21 : C
  22 : C INPUT             : (X1,Y1), (X2,Y2), (X3,Y3) coordinates of the TRI3
  23 : C 
  24 : C OUTPUT            : VOLU: volume
  25 : c                     SURF: surface
  26 : C                     (XCEN,YCEN) center of gravity
  27 : C                     (XC2D,YC2D) center of gravity of the surface
  28 : C
  29 : C************************************************************************
  30 : C
  31 : C Created the 24/02/04
  32 : C
  33 : C
  34 :       IMPLICIT INTEGER(I-N)
  35 :       REAL*8 X1,Y1,X2,Y2,X3,Y3,VOLU,SURF,XCEN,YCEN
  36 :      &     ,RJACO,XC2D,YC2D
  37 : C
  38 :       RJACO=(X2-X1)*(Y3-Y1)-(X3-X1)*(Y2-Y1)
  39 :       RJACO=ABS(RJACO)
  40 :       SURF=0.5D0*RJACO
  41 :       XC2D=(X1+X2+X3)/3.0D0
  42 :       YC2D=(Y1+Y2+Y3)/3.0D0
  43 :       VOLU=RJACO*(X1+X2+X3)/6.0D0
  44 : C
  45 :       XCEN=(X1+X2)**2
  46 :       XCEN=XCEN+((X1+X3)**2)
  47 :       XCEN=XCEN+((X2+X3)**2)
  48 :       XCEN=XCEN*RJACO/24.0D0
  49 :       XCEN=XCEN/VOLU
  50 :       YCEN=(X1+X2)*(Y1+Y2)
  51 :       YCEN=YCEN+((X1+X3)*(Y1+Y3))
  52 :       YCEN=YCEN+((X2+X3)*(Y2+Y3))
  53 :       YCEN=YCEN*RJACO/24.0D0
  54 :       YCEN=YCEN/VOLU
  55 : C      
  56 :       RETURN 
  57 :       END
  58 :       
  59 : 
  60 :  
  61 :