ig3daverage.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	ig3daverage (sig, sigg, vol0, vol0g, rho, eint, eintg, rhog, vol, volg, eplas, eplasg, g_pla, epsd, epsdg, nel, israt)

Function/Subroutine Documentation

ig3daverage()

subroutine ig3daverage	(		sig,
			sigg,
			vol0,
			vol0g,
			rho,
			eint,
			eintg,
			rhog,
			vol,
			volg,
			eplas,
			eplasg,
		integer	g_pla,
			epsd,
			epsdg,
		integer	nel,
		integer, intent(in)	israt )

Definition at line 29 of file ig3daverage.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   C o m m o n   B l o c k s
C-----------------------------------------------
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: ISRAT
      INTEGER G_PLA,NEL
      my_real
     .   sig(nel,6), sigg(nel,6), vol(*), volg(*) , vol0(*),
     .   vol0g(*), eint(*)  , rho(*), eintg(*), rhog(*),
     .   eplas(*), eplasg(*), epsd(*),epsdg(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      my_real
     .   wip(mvsiz)
C-----------------------------------------------
c
      DO i=1,nel  
       wip(i) = vol(i)/volg(i)                                       
       sigg(i,1) = sigg(i,1) + wip(i) * sig(i,1)         
       sigg(i,2) = sigg(i,2) + wip(i) * sig(i,2)         
       sigg(i,3) = sigg(i,3) + wip(i) * sig(i,3)         
       sigg(i,4) = sigg(i,4) + wip(i) * sig(i,4)         
       sigg(i,5) = sigg(i,5) + wip(i) * sig(i,5)         
       sigg(i,6) = sigg(i,6) + wip(i) * sig(i,6)         
       rhog(i)   = rhog(i)   + wip(i) * rho(i)                       
       eintg(i)  = eintg(i)  + eint(i) * vol0(i)/vol0g(i)
      ENDDO
      IF (israt > 0) THEN                                                    
        DO i=1,nel                                                         
          epsdg(i) = epsdg(i) + wip(i) * epsd(i)                                
        ENDDO                                                                
      ENDIF 
      IF (g_pla > 0) THEN                                                    
        DO i=1,nel                                                          
          eplasg(i)  = eplasg(i)  + wip(i) * eplas(i)                              
        ENDDO                                                                 
      ENDIF
c
C-----------
      RETURN