ig3dderishap.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	ig3dderishap (i, nctrl, r, drdxi, detjac, n, b, det)

Function/Subroutine Documentation

ig3dderishap()

subroutine ig3dderishap	(	integer	i,
		integer	nctrl,
		dimension(nctrl)	r,
		dimension(nctrl,3)	drdxi,
			detjac,
		dimension(nctrl,mvsiz)	n,
		dimension(3*nctrl,mvsiz)	b,
		dimension(mvsiz)	det )

Definition at line 28 of file ig3dderishap.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   G l o b a l   P a r a m e t e r s
C-----------------------------------------------
#include      "mvsiz_p.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER I, NCTRL
      my_real,
     .  DIMENSION(NCTRL,3) :: drdxi
      my_real,
     .  DIMENSION(NCTRL) :: r
      my_real
     .  detjac
      my_real,
     .  DIMENSION(NCTRL,MVSIZ) :: n 
      my_real,
     .  DIMENSION(3*NCTRL,MVSIZ) :: b
      my_real,
     .  DIMENSION(MVSIZ) :: det
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER J, INCTRL
C-----------------------------------------------
C     
C       i,j,k,...., NCTRL
C
C               NEL                               NEL
C          /-----------\                /--------------------\    ! 
C    [N] = [Ni  Nk  .. ] \  N     [B] = [dNi/dx  dNk/dx  ..  ] \  3
C          [Ni  Nl  .. ] |  C           [dNi/dy  dNk/dy  ..  ] |  *
C          [Nk  Nm  .. ] |  T           [dNi/dz  dNk/dz  ..  ] |  N
C          [..  ..  .. ] |  R           [dNj/dx  dNl/dx  ..  ] |  C
C          [..  ..  .. ] /  L           [dNj/dy  dNl/dy  ..  ] |  T
C                                       [dNj/dz  dNl/dz  ..  ] |  R
C                                       [  ..      ..    ..  ] /  L
 
      det(i)=detjac
 
      DO j=1,nctrl
       inctrl = (j-1)*3
       n(j,i)=r(j) 
       b(inctrl+1,i) = drdxi(j,1)
       b(inctrl+2,i) = drdxi(j,2)
       b(inctrl+3,i) = drdxi(j,3) 
      ENDDO
 
C-----------
      RETURN