agrad0.F File Reference

#include "implicit_f.inc"
#include "com01_c.inc"
#include "com04_c.inc"
#include "vect01_c.inc"
#include "param_c.inc"
#include "task_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	agrad0 (iparg, ixs, ixq, x, ale_connect, grad, veul, elbuf_str, itask)

Function/Subroutine Documentation

agrad0()

subroutine agrad0	(	integer, dimension(nparg,ngroup)	iparg,
		integer, dimension(nixs,numels)	ixs,
		integer, dimension(nixq,numelq)	ixq,
			x,
		type(t_ale_connectivity), intent(in)	ale_connect,
			grad,
			veul,
		type (elbuf_struct_), dimension (ngroup), target	elbuf_str,
		integer	itask )

Definition at line 40 of file agrad0.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE elbufdef_mod            
      USE ale_connectivity_mod
      use element_mod , only : nixs,nixq
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com01_c.inc"
#include      "com04_c.inc"
#include      "vect01_c.inc"
#include      "param_c.inc"
#include      "task_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IXS(NIXS,NUMELS), IXQ(NIXQ,NUMELQ)
      my_real x(3,numnod),grad(*), veul(*)
      INTEGER IPARG(NPARG,NGROUP),ITASK
      TYPE (ELBUF_STRUCT_), DIMENSION (NGROUP), TARGET :: ELBUF_STR
      TYPE(t_ale_connectivity), INTENT(IN) :: ALE_CONNECT
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER NG,NEL
      TYPE(G_BUFEL_)  ,POINTER :: GBUF     
C-----------------------------------------------
C   S o u r c e   L i n e s
C-----------------------------------------------
      DO ng=itask+1,ngroup,nthread
        gbuf => elbuf_str(ng)%GBUF
        jtur=iparg(12,ng)
        jthe=iparg(13,ng)
        ity =iparg(5,ng)
        IF((jtur == 1 .OR. jthe == 1) .AND. ity /= 51 .AND. iparg(8,ng) /= 1) THEN
          mtn=iparg(1,ng)
          jale=iparg(7,ng)
          jeul=iparg(11,ng)
          jlag=iparg(14,ng)
          llt=iparg(2,ng)
          nft=iparg(3,ng)
          lft=1
          nel = llt
          IF(n2d == 0)THEN
            IF(jale+jlag /= 0)THEN
              CALL agrad3(ixs,x,ale_connect,grad(6*nft+1))
            ELSEIF(jeul == 1)THEN
              CALL egrad3(ixs,x,ale_connect,grad(6*nft+1),veul)
            ELSE
              CALL tgrad3(gbuf%SIG,grad(6*nft+1),nel)
            ENDIF
          ELSE
            IF(jlag+jale /= 0)THEN
              CALL agrad2(ixq,x,ale_connect,grad(4*nft+1))
            ELSEIF(jeul == 1)THEN
              CALL egrad2(ixq,x,ale_connect,grad(4*nft+1))
            ELSE
              CALL tgrad2(gbuf%SIG,grad(4*nft+1),nel)
            ENDIF
          ENDIF
        ENDIF
      enddo!next NG
C
      CALL my_barrier
C
      RETURN