cbavarnl.F File Reference

#include "implicit_f.inc"
#include "vectorize.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	cbavarnl (jft, jlt, ng, ixc, nloc_dmg, varnl, nddl, nc1, nc2, nc3, nc4, nel)

Function/Subroutine Documentation

cbavarnl()

subroutine cbavarnl	(	integer	jft,
		integer	jlt,
		integer	ng,
		integer, dimension(nixc,*)	ixc,
		type(nlocal_str_), target	nloc_dmg,
		dimension(nel,nddl), intent(inout)	varnl,
		integer	nddl,
		integer, dimension(*), intent(inout)	nc1,
		integer, dimension(*), intent(inout)	nc2,
		integer, dimension(*), intent(inout)	nc3,
		integer, dimension(*), intent(inout)	nc4,
		integer	nel )

Definition at line 30 of file cbavarnl.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------      
      USE nlocal_reg_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER JFT,JLT,NG,IXC(NIXC,*),NDDL,NEL
      INTEGER, INTENT(INOUT) :: 
     .   NC1(*), NC2(*), NC3(*), NC4(*)
      my_real, INTENT(INOUT) ::
     .   varnl(nel,nddl)
      TYPE(NLOCAL_STR_), TARGET :: NLOC_DMG
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,  K, 
     .        INOD(4), IPOS(4),L_NLOC
      my_real 
     .   vpg(2,4),pg1,pg,ksi,eta,sf1,sf2,sf3,sf4
       parameter(pg=.577350269189626)
       parameter(pg1=-.577350269189626)
      my_real, 
     .   DIMENSION(:), POINTER :: dnl 
C--------------------------
C     INITIALISATION
C-------------------------
      DATA vpg/pg1,pg1,pg,pg1,pg,pg,pg1,pg/
C---------------------------------------------------------------------------------------------------------------
C --- Non-local variable VARNL in Gauss point NG
C    
C --- Position of the Gauss point
       ksi = vpg(1,ng)
       eta = vpg(2,ng)
C --- Shape function value
       sf1 = (1-ksi)*(1-eta)
       sf2 = (1+ksi)*(1-eta)
       sf3 = (1+ksi)*(1+eta)
       sf4 = (1-ksi)*(1+eta)
C --- Length of the non-local vector       
       l_nloc = nloc_dmg%L_NLOC
       dnl  => nloc_dmg%DNL(1:l_nloc) ! DNL = non local variable increment
C --- Loop over 
       DO i=jft,jlt       
         nc1(i)  = ixc(2,i)
         nc2(i)  = ixc(3,i)
         nc3(i)  = ixc(4,i)
         nc4(i)  = ixc(5,i)
       ENDDO
       DO k = 1,nddl
#include "vectorize.inc"
         DO i=jft,jlt
           inod(1) = nloc_dmg%IDXI(nc1(i))
           inod(2) = nloc_dmg%IDXI(nc2(i))
           inod(3) = nloc_dmg%IDXI(nc3(i))
           inod(4) = nloc_dmg%IDXI(nc4(i)) 
           ipos(1) = nloc_dmg%POSI(inod(1)) 
           ipos(2) = nloc_dmg%POSI(inod(2)) 
           ipos(3) = nloc_dmg%POSI(inod(3)) 
           ipos(4) = nloc_dmg%POSI(inod(4)) 
           varnl(i,k) = zero
           varnl(i,k) = fourth*(sf1*dnl(ipos(1)+k-1) + sf2*dnl(ipos(2)+k-1) 
     .                       + sf3*dnl(ipos(3)+k-1) + sf4*dnl(ipos(4)+k-1))
         ENDDO
       ENDDO
C