vc00ft.f

Go to the documentation of this file.

!                   *****************
                    SUBROUTINE vc00ft
!                   *****************
!
     &(xmul,x,y,z,ikle1,ikle2,ikle3,nbor,nelem,nelmax,w1,w2,w3)
!
!***********************************************************************
! BIEF   V6P1                                   21/08/2010
!***********************************************************************
!
!brief    COMPUTES THE FOLLOWING VECTOR IN FINITE ELEMENTS:
!code
!+                    /
!+    VEC(I) = XMUL  /    PSI(I)  D(OMEGA)
!+                  /OMEGA
!+
!+    PSI(I) IS A BASE OF TYPE P1 SEGMENT BUT ON A MESH WITH
!+    VERTICAL TRIANGLES IN THE X,Y,Z SPACE
!+
!+    F IS A VECTOR OF TYPE IELMF
!
!warning  THE JACOBIAN MUST BE POSITIVE
!warning  THE RESULT IS IN W IN NOT ASSEMBLED FORM
!
!history  J-M HERVOUET (LNH)    ; F LEPEINTRE (LNH)
!+        05/02/91
!+        V5P4
!+
!
!history  N.DURAND (HRW), S.E.BOURBAN (HRW)
!+        13/07/2010
!+        V6P0
!+   Translation of French comments within the FORTRAN sources into
!+   English comments
!
!history  N.DURAND (HRW), S.E.BOURBAN (HRW)
!+        21/08/2010
!+        V6P0
!+   Creation of DOXYGEN tags for automated documentation and
!+   cross-referencing of the FORTRAN sources
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| IKLE1          |-->| FIRST POINT OF TRIANGLES
!| IKLE2          |-->| SECOND POINT OF TRIANGLES
!| IKLE3          |-->| THIRD POINT OF TRIANGLES
!| NBOR           |-->| GLOBAL NUMBER OF BOUNDARY POINTS
!| NELEM          |-->| NUMBER OF ELEMENTS
!| NELMAX         |-->| MAXIMUM NUMBER OF ELEMENTS
!| W1             |<--| RESULT IN NON ASSEMBLED FORM
!| W2             |<--| RESULT IN NON ASSEMBLED FORM
!| W3             |<--| RESULT IN NON ASSEMBLED FORM
!| X              |-->| ABSCISSAE OF POINTS IN THE MESH
!| XMUL           |-->| MULTIPLICATION COEFFICIENT
!| Y              |-->| ORDINATES OF POINTS IN THE MESH
!| Z              |-->| ELEVATIONS OF POINTS IN THE MESH
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE declarations_special
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER, INTENT(IN) :: NELEM,NELMAX
      INTEGER, INTENT(IN) :: NBOR(*)
      INTEGER, INTENT(IN) :: IKLE1(nelmax),IKLE2(nelmax),IKLE3(nelmax)
!
      DOUBLE PRECISION, INTENT(IN)    :: X(*),Y(*),Z(*)
      DOUBLE PRECISION, INTENT(INOUT) :: W1(nelmax)
      DOUBLE PRECISION, INTENT(INOUT) :: W2(nelmax)
      DOUBLE PRECISION, INTENT(INOUT) :: W3(nelmax)
      DOUBLE PRECISION, INTENT(IN)    :: XMUL
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER IELEM,I1,I2,I3
      DOUBLE PRECISION XSUR3,COEF,X1,X2,X3,Y1,Y2
      DOUBLE PRECISION Y3,Z1,Z2,Z3,S
!
      INTRINSIC sqrt
!
!***********************************************************************
!
      xsur3 = xmul/3.d0
!
!   LOOP ON THE BOUNDARY SIDES
!
      DO ielem = 1,nelem
!
!       GLOBAL NUMBERING OF THE SIDE NODES
!
        i1 = nbor(ikle1(ielem))
        i2 = nbor(ikle2(ielem))
        i3 = nbor(ikle3(ielem))
!
        x1 = x(i1)
        y1 = y(i1)
        z1 = z(i1)
!
        x2 = x(i2)-x1
        x3 = x(i3)-x1
        y2 = y(i2)-y1
        y3 = y(i3)-y1
        z2 = z(i2)-z1
        z3 = z(i3)-z1
!
!       COMPUTES THE AREA OF THE TRIANGLE (BY VECTOR PRODUCT)
!
        s=0.5d0*sqrt(  (y2*z3-y3*z2)**2
     &                +(x3*z2-x2*z3)**2  )
!    &                +(X2*Y3-X3*Y2)**2  )  THIS TERM IS 0
!
        coef=xsur3*s
!
        w1(ielem) = coef
        w2(ielem) = coef
        w3(ielem) = coef
!
      ENDDO ! IELEM
!
!-----------------------------------------------------------------------
!
!
!-----------------------------------------------------------------------
!
!     NOTE: FOR A PLANE TRIANGLE (VC00AA):
!
!     XSUR3 = XMUL / 3.D0
!
!     DO IELEM = 1 , NELEM
!
!
!       COEF = XSUR3 * SURFAC(IELEM)
!
!       W1(IELEM) = COEF
!       W2(IELEM) = COEF
!       W3(IELEM) = COEF
!
!     ENDDO ! IELEM
!
!-----------------------------------------------------------------------
!
      RETURN
      END