alero2.F File Reference

#include "implicit_f.inc"
#include "comlock.inc"
#include "com08_c.inc"
#include "vect01_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	alero2 (off, volo, rhon, flux, flu1, w, vy1, vy2, vy3, vy4, vz1, vz2, vz3, vz4, voln, dvol, vd2, nc1, nc2, nc3, nc4, ngl)

Function/Subroutine Documentation

alero2()

subroutine alero2	(		off,
			volo,
			rhon,
			flux,
			flu1,
			w,
			vy1,
			vy2,
			vy3,
			vy4,
			vz1,
			vz2,
			vz3,
			vz4,
			voln,
			dvol,
			vd2,
		integer, dimension(*)	nc1,
		integer, dimension(*)	nc2,
		integer, dimension(*)	nc3,
		integer, dimension(*)	nc4,
		integer, dimension(*)	ngl )

Definition at line 33 of file alero2.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      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   D u m m y   A r g u m e n t s
C-----------------------------------------------
      my_real
     .   off(*), volo(*), rhon(*), flux(4,*), flu1(*), w(3,*),
     .   vy1(*),vy2(*),vy3(*),vy4(*),vz1(*),vz2(*),vz3(*),vz4(*),voln(*), dvol(*), vd2(*)
      INTEGER NC1(*), NC2(*), NC3(*), NC4(*), NGL(*)
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com08_c.inc"
#include      "vect01_c.inc"
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER  I
      my_real vy, vz
C-----------------------------------------------
C   S o u r c e   L i n e s
C-----------------------------------------------
 
      DO i=lft,llt
        rhon(i)=rhon(i)/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
C
      DO i=lft,llt
        vy=vy1(i)+vy2(i)+vy3(i)+vy4(i)-(w(2,nc1(i))+w(2,nc2(i))+w(2,nc3(i))+w(2,nc4(i)))
        vz=vz1(i)+vz2(i)+vz3(i)+vz4(i)-(w(3,nc1(i))+w(3,nc2(i))+w(3,nc3(i))+w(3,nc4(i)))
        vd2(i)=(vy**2+vz**2)/sixteen
      ENDDO
C
      DO i=lft,llt
        IF(off(i) /= 0.0 .AND. mtn /= 11)THEN
          IF(rhon(i) <= 0.0)THEN
            CALL ancmsg(msgid=12,anmode=aninfo,i1=ngl(i),r1=rhon(i))
            CALL arret(2)
          ENDIF
        ENDIF
      ENDDO
C-----------------------------------------------
      RETURN