frho3.F File Reference

#include "implicit_f.inc"
#include "comlock.inc"
#include "mvsiz_p.inc"
#include "units_c.inc"
#include "scr06_c.inc"
#include "com08_c.inc"
#include "lockon.inc"
#include "lockoff.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	frho3 (uvar, rho0, dxx, dyy, dzz, flux, flu1, voln, ngl, mat, off, nel)

Function/Subroutine Documentation

frho3()

subroutine frho3	(	target	uvar,
			rho0,
			dxx,
			dyy,
			dzz,
			flux,
			flu1,
			voln,
		integer, dimension(*)	ngl,
		integer, dimension(*)	mat,
			off,
		integer	nel )

Definition at line 30 of file frho3.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"
#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"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "units_c.inc"
#include      "scr06_c.inc"
#include      "com08_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NGL(*),MAT(*),NEL
      my_real uvar(nel,*),flux(mvsiz,6), flu1(*), voln(*), dxx(*),dyy(*),dzz(*),off(*),rho0(*)
      TARGET :: uvar
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      my_real dvv, e0(nel),vavg,rv,rvp
      INTEGER I, J,COUNT,LIST(NEL),II
      my_real, DIMENSION(:), POINTER :: volo,rhon,eint,dvol
C-----------------------------------------------
      rhon => uvar(1:nel,1)   
      eint => uvar(1:nel,2)   
      volo => uvar(1:nel,3)   
      dvol => uvar(1:nel,4)   
C
C  incompressible flow avoir plus en cas de besoin
C      INCOMP = 1 not tested
      IF(ale%GLOBAL%INCOMP == 1)THEN
         e0 = zero ! 
         DO i=1,nel
          dvv=dt1*(dxx(i)+dyy(i)+dzz(i))
          rhon(i)=rhon(i)-dvv*rho0(i)
          vavg=half*(voln(i)+volo(i))
          dvol(i)=vavg*dvv
          eint(i)=eint(i)*volo(i)-e0(i)*dvv*vavg
         ENDDO
C compresssible flow
      ELSE
         DO i=1,nel
           rhon(i) = rhon(i)/voln(i)
          
           dvol(i) = voln(i) 
     .             - volo(i)+0.5*dt1*(flu1(i)+flux(i,1)+flux(i,2)
     .             +  flux(i,3)+flux(i,4)+flux(i,5)+flux(i,6) )
            volo(i)=voln(i)
           rhon(i) = max(rhon(i) , em30)
       
        ENDDO 
      ENDIF
C
        count=0
        DO  i=1,nel
          IF(rhon(i) <= zero .AND. off(i)/= zero)THEN
             count = count + 1
              list(count)=i
           ENDIF  
        ENDDO
        DO ii=1,count
         i=list(ii)
#include "lockon.inc"
        WRITE(iout,1000) ngl(i),rhon(i)
        WRITE(istdo,1000) ngl(i),rhon(i)
#include "lockoff.inc"
 1000   FORMAT('-STOP GAZ ',i8,' HAS NEGATIVE RHO',g12.5)
       ENDDO
C
      RETURN