rbilan33.F File Reference

#include "implicit_f.inc"
#include "comlock.inc"
#include "mvsiz_p.inc"
#include "param_c.inc"
#include "com06_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	rbilan33 (jft, jlt, eint, partsav, umas, v, ipartr, gresav, grth, igrth, nc1, nc2, igre)

Function/Subroutine Documentation

rbilan33()

subroutine rbilan33	(	integer	jft,
		integer	jlt,
			eint,
			partsav,
			umas,
			v,
		integer, dimension(*)	ipartr,
			gresav,
		integer, dimension(*)	grth,
		integer, dimension(*)	igrth,
		integer, dimension(*)	nc1,
		integer, dimension(*)	nc2,
		integer, intent(in)	igre )

Definition at line 30 of file rbilan33.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
#include "comlock.inc"
C-----------------------------------------------
C   G l o b a l   P a r a m e t e r s
C-----------------------------------------------
#include      "mvsiz_p.inc"
#include      "param_c.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com06_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: IGRE
      INTEGER JFT,JLT, IPARTR(*),GRTH(*),IGRTH(*),NC1(*),NC2(*)
      my_real umas(*),eint(*),partsav(npsav,*),v(3,*),gresav(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,MX, FLAG
      my_real
     .   vxa, vya, vza, xmas2, va2, reintt, 
     .   ek(mvsiz), xm(mvsiz), ym(mvsiz), zm(mvsiz),
     .   rbidon(1),off(mvsiz)
C-----------------------------------------------
      flag = 0
      rbidon = zero
      off    = zero
      reintt = zero
      DO i=jft,jlt
        reintt=reintt + eint(i)
      ENDDO
C
!$OMP ATOMIC
        reint = reint + reintt
C
      DO i=jft,jlt
        xmas2=umas(i)*half
        vxa=v(1,nc1(i))+v(1,nc2(i))
        vya=v(2,nc1(i))+v(2,nc2(i))
        vza=v(3,nc1(i))+v(3,nc2(i))
        xm(i)= xmas2*vxa
        ym(i)= xmas2*vya
        zm(i)= xmas2*vza
        va2   =v(1,nc1(i))*v(1,nc1(i))+v(1,nc2(i))*v(1,nc2(i))
     .        +v(2,nc1(i))*v(2,nc1(i))+v(2,nc2(i))*v(2,nc2(i))
     .        +v(3,nc1(i))*v(3,nc1(i))+v(3,nc2(i))*v(3,nc2(i))
        ek(i)= xmas2*va2*half
      ENDDO
C
      IF (igre /= 0) THEN
        flag = 0
        CALL grelem_sav(jft   ,jlt   ,gresav,igrth ,grth  ,
     2                  off   ,eint  ,ek    ,xm    ,ym    ,
     3                  zm    ,umas  ,rbidon,rbidon,rbidon,
     4                  rbidon,rbidon,rbidon,rbidon,rbidon,
     5                  rbidon,rbidon,rbidon,rbidon,rbidon,
     6                  rbidon,flag)
      ENDIF
C
      DO i=jft,jlt
         mx = ipartr(i)
         partsav(1,mx)= partsav(1,mx) + eint(i)
         partsav(2,mx)= partsav(2,mx) + ek(i)
         partsav(3,mx)= partsav(3,mx) + xm(i)
         partsav(4,mx)= partsav(4,mx) + ym(i)
         partsav(5,mx)= partsav(5,mx) + zm(i)
         partsav(6,mx)= partsav(6,mx) + umas(i)
      ENDDO
C
      RETURN