ig3dmass3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	ig3dmass3 (rho, ms, partsav, x, y, z, vx, vy, vz, ipart, msig3d, volu, msnf, mssf, in, vr, ins, wma, rhocp, mcp, mcps, mssa, rhof, frac, nctrl, kxig3d, ixig3d, r, detjac, pgauss, i)

Function/Subroutine Documentation

ig3dmass3()

subroutine ig3dmass3	(		rho,
			ms,
			partsav,
			x,
			y,
			z,
			vx,
			vy,
			vz,
		integer, dimension(*)	ipart,
			msig3d,
			volu,
			msnf,
			mssf,
			in,
			vr,
			ins,
			wma,
			rhocp,
			mcp,
			mcps,
			mssa,
			rhof,
			frac,
		integer	nctrl,
		integer, dimension(nixig3d,*)	kxig3d,
		integer, dimension(*)	ixig3d,
			r,
			detjac,
			pgauss,
		integer	i )

Definition at line 29 of file ig3dmass3.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE ale_mod
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 IPART(*),NDDIM, NCTRL, I, KXIG3D(NIXIG3D,*),IXIG3D(*)
      my_real
     .   rho(*),ms(*),x(nctrl,*),
     .   y(nctrl,*),z(nctrl,*),
     .   vx(nctrl,mvsiz),vy(nctrl,mvsiz),vz(nctrl,mvsiz),
     .   partsav(20,*),volu(*),
     .   msig3d(numelig3d,*), mssf(8,*), msnf(*),
     .   in(*),vr(3,*),ins(8,*),wma(*),
     .   rhocp(*),mcp(*),mcps(8,*), mssa(*),rhof(*),
     .   frac(*),pgauss,detjac,r(*)
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"
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER J,IP,I1,I2,I3,I4,I5,I6,I7,I8, ITNCTRL, INCTRL
      my_real xx,yy,zz,xy,yz,zx,iner,rcp
      my_real mass(mvsiz),massf(mvsiz)
 
C=======================================================================
C   S o u r c e  L i n e s
C=======================================================================
C
      xx=zero
      xy=zero
      yy=zero
      yz=zero
      zz=zero
      zx=zero
C
      DO itnctrl=1,nctrl
        mass(i)= rho(i)*r(itnctrl)*detjac*pgauss 
        massf(i)= frac(i)*rhof(i)*r(itnctrl)*detjac*pgauss
        msig3d(i+nft,itnctrl)=msig3d(i+nft,itnctrl) + mass(i)
C
        ip=ipart(i)
        partsav(1,ip)=partsav(1,ip) + mass(i)
c
        partsav(2,ip)=partsav(2,ip) + mass(i)*x(itnctrl,i) 
        partsav(3,ip)=partsav(3,ip) + mass(i)*y(itnctrl,i)
        partsav(4,ip)=partsav(4,ip) + mass(i)*z(itnctrl,i)
c
        xx = xx + mass(i) * x(itnctrl,i)*x(itnctrl,i)  
        xy = xy + mass(i) * x(itnctrl,i)*y(itnctrl,i)
        yy = yy + mass(i) * y(itnctrl,i)*y(itnctrl,i)
        yz = yz + mass(i) * y(itnctrl,i)*z(itnctrl,i)
        zz = zz + mass(i) * z(itnctrl,i)*z(itnctrl,i)
        zx = zx + mass(i) * z(itnctrl,i)*x(itnctrl,i)
      ENDDO
C
      partsav(5,ip) =partsav(5,ip)  + (yy+zz)
      partsav(6,ip) =partsav(6,ip)  + (zz+xx)
      partsav(7,ip) =partsav(7,ip)  + (xx+yy)
      partsav(8,ip) =partsav(8,ip)  - xy
      partsav(9,ip) =partsav(9,ip)  - yz
      partsav(10,ip)=partsav(10,ip) - zx
C
      DO itnctrl=1,nctrl
      partsav(11,ip)=partsav(11,ip) + mass(i)*
     .     vx(itnctrl,i)
      partsav(12,ip)=partsav(12,ip) + mass(i)*
     .     vy(itnctrl,i)
      partsav(13,ip)=partsav(13,ip) + mass(i)*
     .     vz(itnctrl,i)
      partsav(14,ip)=partsav(14,ip) + half * mass(i) *
     .   (vx(itnctrl,i)*vx(itnctrl,i) +
     .    vy(itnctrl,i)*vy(itnctrl,i) +
     .    vz(itnctrl,i)*vz(itnctrl,i))
      ENDDO
C
      IF(irest_mselt /= 0)THEN
        mssa(nft+i)=mass(i)
      ENDIF
C       
      IF(jale == 3 .AND. jlag == 1)THEN
        DO itnctrl=1,nctrl
          mssf(itnctrl,i)=massf(i)
        ENDDO
      ELSEIF(jale+jeul>0)THEN
        DO itnctrl=1,nctrl
          mssf(itnctrl,i)=mass(i)
        ENDDO
      ENDIF
C
c       IF(JTHE < 0 ) THEN      
c           RCP=RHOCP(I)*VOLU(I)/NCTRL
c           DO ITNCTRL=1,NCTRL
c             MCPS(ITNCTRL,I) =RCP
c           ENDDO 
c       ENDIF
c       IF(ISROT==1)THEN
c        IF(IRODDL==0)THEN
C prov gw
c          STOP 1119
c        ENDIF
c         INER=(MASS(I)*VOLU(I)**TWO_THIRD)/SIX
c         DO ITNCTRL=1,NCTRL
c           INS(ITNCTRL,I)=INER
c         ENDDO
c       ENDIF
C
      IF(jale > 0 .AND. ale%GRID%NWALE == 4)THEN
         DO itnctrl=1,nctrl
          inctrl = ixig3d(kxig3d(4,i+nft)+itnctrl-1)         
          wma(inctrl)=wma(inctrl)+three_half
         ENDDO
      ENDIF
C-----------
      RETURN