ksupg.f

Go to the documentation of this file.

!                   ****************
                    SUBROUTINE ksupg
!                   ****************
!
     &(kx,ky,xmul,u,v,mesh)
!
!***********************************************************************
! BIEF   V6P1                                   21/08/2010
!***********************************************************************
!
!brief    COMPUTES A VECTOR, USED BY THE METHOD:
!+                STREAMLINE UPWIND PETROV GALERKIN (SUPG)
!+                WITH AN OFF-CENTERING OF 1.
!code
!+                    DX   U
!+             KX = -----------
!+                   2 NORM(U)
!+
!+                    DY   V
!+             KY = -----------
!+                   2 NORM(U)
!
!history  J-M HERVOUET (LNH)
!+        08/12/94
!+        V5P6
!+
!
!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
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| KX             |-->| FIRST COMPONENT OF RESULTING VECTOR
!| KY             |-->| SECOND COMPONENT OF RESULTING VECTOR
!| MESH           |-->| MESH STRUCTURE
!| U              |-->| FIRST COMPONENT OF VELOCITY
!| V              |-->| SECOND COMPONENT OF VELOCITY
!| XMUL           |-->| MULTIPLICATION COEFICIENT
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief, ex_ksupg => ksupg
!
      USE declarations_special
      IMPLICIT NONE
!
      TYPE(bief_mesh), INTENT(IN)   :: MESH
!
      TYPE(bief_obj), INTENT(INOUT) :: KX,KY
      TYPE(bief_obj), INTENT(IN)    :: U,V
!
      DOUBLE PRECISION, INTENT(IN)  :: XMUL
!
!-----------------------------------------------------------------------
!
      IF(u%ELM.EQ.11.OR.u%ELM.EQ.12.OR.u%ELM.EQ.13) THEN
!
        CALL kspg11(kx%R,ky%R,mesh%XEL%R,mesh%YEL%R,u%R,v%R,
     &              mesh%IKLE%I,mesh%NELEM,mesh%NELMAX,xmul)
!
!  ELEMENT NOT IMPLEMENTED: ERROR
!
      ELSE
        WRITE(lu,101) u%ELM
101     FORMAT(1x,'KSUPG (BIEF): U%ELM = ',1i6,' ELEMENT NOT AVAILABLE')
        CALL plante(1)
        stop
      ENDIF
!
!-----------------------------------------------------------------------
!
      RETURN
      END