i12nod3.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	i12nod3 (nrt, iseg, irect, segvar, jcodv, nodvar, fcount, nsn)

Function/Subroutine Documentation

i12nod3()

subroutine i12nod3	(	integer	nrt,
		integer, dimension(*)	iseg,
		integer, dimension(4,*)	irect,
		type(t_segvar), target	segvar,
		integer, dimension(*)	jcodv,
			nodvar,
		integer, dimension(*)	fcount,
		integer	nsn )

Definition at line 31 of file i12nod3.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE segvar_mod
      USE ale_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NRT,IRECT(4,*),JCODV(*),ISEG(*),FCOUNT(*),NSN
      my_real nodvar(*)
      TYPE(t_segvar),TARGET :: SEGVAR
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,J,KK,KVAR,NODAD,SEGAD
      my_real,DIMENSION(:),POINTER :: ptr
C-----------------------------------------------
      DO i=1,nsn*ale%GLOBAL%NVCONV
        nodvar(i)=zero
      ENDDO
 
 
      DO kvar=1,ale%GLOBAL%NVCONV
        DO i=1,nrt
          kk=iseg(i)                     
          SELECT CASE(kvar)              
            CASE(1)                      
              ptr(1:) => segvar%RHO(1:)  
            CASE(2)                      
              ptr(1:) => segvar%EINT(1:) 
            CASE(3)                      
              ptr(1:) => segvar%RK(1:)   
            CASE(4)                      
              ptr(1:) => segvar%RE(1:)   
            CASE(5)                      
              ptr(1:) => segvar%UVAR(1:) 
          END SELECT                     
               
          DO j=1,4
            nodad=ale%GLOBAL%NVCONV*(irect(j,i)-1)+kvar
            nodvar(nodad)=nodvar(nodad)+ptr(kk)
          ENDDO
        ENDDO
      ENDDO
 
      DO i=1,nsn
        DO kvar=1,ale%GLOBAL%NVCONV
          nodad=ale%GLOBAL%NVCONV*(i-1)+kvar
          nodvar(nodad)=nodvar(nodad)/float(abs(fcount(i)))
        ENDDO
      ENDDO
 
      RETURN