nodm11.F File Reference

#include "implicit_f.inc"
#include "param_c.inc"
#include "com01_c.inc"
#include "com04_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	nodm11 (pm, ixs, ixq, icode)

Function/Subroutine Documentation

nodm11()

subroutine nodm11	(		pm,
		integer, dimension(nixs,*)	ixs,
		integer, dimension(nixq,*)	ixq,
		integer, dimension(*)	icode )

Definition at line 28 of file nodm11.F.

C     BLOCAGE DES NOEUDS MATERIAU 11 HORS DOMAINE CALCULE
C     FP 12/97
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   G l o b a l   P a r a m e t e r s
C-----------------------------------------------
#include      "param_c.inc"
#include      "com01_c.inc"
#include      "com04_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IXS(NIXS,*),IXQ(NIXQ,*),ICODE(*)
      my_real
     .   pm(npropm,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER, DIMENSION(NUMNOD) :: ITAG
      INTEGER I,J,IC,IC1,IC2,IC3,IC4
C=======================================================================
      DO i=1,numnod
       itag(i) =0
      ENDDO
C TAGAGE MAT 11
       IF(n2d==0)THEN
        DO i=1,numels
         IF(int(pm(19,ixs(1,i)))==11) THEN
          DO  j=2,9
           itag(ixs(j,i))=1
          ENDDO
         ENDIF
        ENDDO
       ELSE
        DO i=1,numelq
         IF(int(pm(19,ixq(1,i)))==11) THEN
          DO j=2,5
           itag(ixq(j,i))=1
          ENDDO
         ENDIF
        ENDDO
       ENDIF
       IF(n2d==0)THEN
        DO i=1,numels
         IF(int(pm(19,ixs(1,i)))/=11) THEN
          DO  j=2,9
           itag(ixs(j,i))=0
          ENDDO
         ENDIF
        ENDDO
       ELSE
        DO i=1,numelq
         IF(int(pm(19,ixq(1,i)))/=11) THEN
          DO j=2,5
           itag(ixq(j,i))=0
          ENDDO
         ENDIF
        ENDDO
       ENDIF
C
       DO i=1,numnod
        IF(itag(i)==1)THEN
         ic=icode(i)
         ic1=ic/512
         ic2=(ic-512*ic1)/64
         ic3=(ic-512*ic1-64*ic2)/8
         ic4=ic-512*ic1-64*ic2-8*ic3
         ic1=7
         icode(i)=ic1*512+ic2*64+ic3*8+ic4
        ENDIF
       ENDDO
      RETURN