blero2.F File Reference

#include "implicit_f.inc"
#include "comlock.inc"
#include "mvsiz_p.inc"
#include "param_c.inc"
#include "com08_c.inc"
#include "vect01_c.inc"
#include "scr17_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	blero2 (alph, volo, rhon, flux, flu1, off, pm, ngl, voln, aire, d1, d2, d3, d4, d5, d6, vd2, dvol, vy1, vy2, vy3, vy4, vz1, vz2, vz3, vz4, mat)

Function/Subroutine Documentation

blero2()

subroutine blero2	(		alph,
			volo,
			rhon,
			flux,
			flu1,
			off,
			pm,
		integer, dimension(*)	ngl,
			voln,
			aire,
			d1,
			d2,
			d3,
			d4,
			d5,
			d6,
			vd2,
			dvol,
			vy1,
			vy2,
			vy3,
			vy4,
			vz1,
			vz2,
			vz3,
			vz4,
		integer, dimension(*)	mat )

Definition at line 33 of file blero2.F.

      USE message_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      "param_c.inc"
#include      "com08_c.inc"
#include      "vect01_c.inc"
#include      "scr17_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      my_real
     .   alph(*), volo(*), rhon(*), flux(4,*), flu1(*),
     .   off(*), pm(npropm,*),voln(mvsiz), 
     .   aire(*), d1(*), d2(*), d3(*), d4(*), d5(*), d6(*), vd2(*), dvol(*),
     .   vy1(*), vy2(*), vy3(*), vy4(*), 
     .   vz1(*), vz2(*), vz3(*), vz4(*)
      INTEGER MAT(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, II, NGL(*)
      my_real
     .   vy, vz, scal
C-----------------------------------------------
 
      DO i=lft,llt
        voln(i)=voln(i)*alph(i)
        aire(i)=aire(i)*alph(i)
        rhon(i) = rhon(i) / max(em30,voln(i))
        dvol(i) = voln(i)-volo(i)+dt1*half*(flu1(i)
     .          +flux(1,i)+flux(2,i)+flux(3,i)+flux(4,i))
        volo(i)=voln(i)
      ENDDO
 
      DO i=lft,llt
        IF(alph(i)==zero)THEN
          rhon(i) =pm( 1,mat(i))
          dvol(i) =zero
          off(i)  =zero
          idel7nok = 1
        ENDIF
      ENDDO
 
      DO i=lft,llt
        vy=vy1(i)+vy2(i)+vy3(i)+vy4(i)
        vz=vz1(i)+vz2(i)+vz3(i)+vz4(i)
        vd2(i)=(vy**2+vz**2)/sixteen
      ENDDO
 
      DO i=lft,llt
        IF((mtn==11).OR.(rhon(i)>zero))THEN
          IF(alph(i)<one)THEN
            scal=abs(dvol(i))/
     .           max(em15,voln(i)*dt1*abs(d1(i)+d2(i)+d3(i)))
            scal= min(one,scal)
            d1(i)=d1(i)*scal
            d2(i)=d2(i)*scal
            d3(i)=d3(i)*scal
            d4(i)=d4(i)*scal
            d5(i)=d5(i)*scal
            d6(i)=d6(i)*scal
          ENDIF
        ELSE
          ii = nft+i
          CALL ancmsg(msgid=12,anmode=aninfo,
     .                i1=ngl(i),r1=rhon(i))
          CALL arret(2)
        ENDIF
      ENDDO !next I
 
      RETURN