s11fx3.F File Reference

#include "implicit_f.inc"
#include "comlock.inc"
#include "task_c.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	s11fx3 (pm, flux, ale_connect, ixs, ipm, bufmat, nel, nft, jale, jeul, nummat, s_ixs)

Function/Subroutine Documentation

s11fx3()

subroutine s11fx3	(		pm,
			flux,
		type(t_ale_connectivity), intent(in)	ale_connect,
		integer, dimension(nixs,s_ixs), intent(in)	ixs,
		integer, dimension(npropmi,nummat), intent(in)	ipm,
			bufmat,
		integer, intent(in)	nel,
		integer, intent(in)	nft,
		integer, intent(in)	jale,
		integer, intent(in)	jeul,
		integer, intent(in)	nummat,
		integer, intent(in)	s_ixs )

Definition at line 32 of file s11fx3.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE ale_connectivity_mod
      use element_mod , only : nixs
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
#include      "comlock.inc"
#include      "task_c.inc"
C-----------------------------------------------
C   G l o b a l   P a r a m e t e r s
C-----------------------------------------------
#include      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NEL
      INTEGER, INTENT(IN) :: NFT
      INTEGER, INTENT(IN) :: JALE
      INTEGER, INTENT(IN) :: JEUL
      INTEGER, INTENT(IN) :: NUMMAT
      INTEGER, INTENT(IN) :: S_IXS
      my_real pm(npropm,nummat),flux(6,*),bufmat(*)
      INTEGER, INTENT(IN) :: IXS(NIXS,S_IXS),IPM(NPROPMI,NUMMAT)
      TYPE(t_ale_connectivity), INTENT(IN) :: ALE_CONNECT
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, J,JJ,K,IFLG,IADBUF,IAD2
C-----------------------------------------------
C     Incoming Enthalpy material law 11
C-----------------------------------------------
      IF(jeul+jale==1)THEN
          DO i=1,nel
             iad2 = ale_connect%ee_connect%iad_connect(i + nft)
            DO k=1,6
            jj = ale_connect%ee_connect%connected(iad2 + k - 1)
            IF(jj>0)THEN         ! negatif si ebcs
              j = ixs(1,jj)
              IF(nint(pm(50,j))==8)THEN
                !TEST LAW11
                IF(nint(pm(19,j))==11)THEN
                    flux(k,i) = max(flux(k,i),zero)
                !TEST LAW51 IFORM=3
                ELSEIF(nint(pm(19,j))==51)THEN
                  iadbuf = ipm(27,j)
                  iflg   = nint(bufmat(31+iadbuf-1))
                  IF (iflg == 6)THEN
                    flux(k,i) = max(flux(k,i),zero)
                  ENDIF
                endif!LAW TEST
              endif!ISOLNOD==8             
            endif!JJ>0
          ENDDO
        ENDDO
      ENDIF
C-----------------------------------------------
      RETURN