condrmat.F File Reference

#include "implicit_f.inc"
#include "com04_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	condrmat (icodt, icodr, irbym, lnrbym, icodrbym)

Function/Subroutine Documentation

condrmat()

subroutine condrmat	(	integer, dimension(*)	icodt,
		integer, dimension(*)	icodr,
		integer, dimension(nirbym,*)	irbym,
		integer, dimension(*)	lnrbym,
		integer, dimension(2,*)	icodrbym )

Definition at line 28 of file condrmat.F.

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  ICODT(*), ICODR(*), IRBYM(NIRBYM,*),LNRBYM(*),
     .         ICODRBYM(2,*)
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com04_c.inc"
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER N, M, NSLT, ITX,ITY,ITZ,IRX,IRY,IRZ,
     .        NSL,I,LCOD
C-----------------------------------------------
      nslt= 0
      DO  m=1,nrbym
C
          irx = 0
          iry = 0
          irz = 0 
          itx = 0
          ity = 0
          itz = 0        
          nsl = irbym(2,m)
          DO  i=1, nsl
            n = lnrbym(nslt + i)
C rotation bcs             
            lcod = icodr(n)
            IF(lcod > 0) THEN
               icodr(n)  = 0
               IF(lcod == 7) THEN
                icodrbym(1,m) = lcod
                irx = 4
                iry = 2
                irz = 1
               ELSE
                 IF(lcod == 4 .OR. lcod == 5 .OR. lcod == 6) irx = 4
                 IF(lcod == 2 .OR. lcod == 3 .OR. lcod == 6) iry = 2
                 IF(lcod == 1 .OR. lcod == 3 .OR. lcod == 5) irz = 1
               ENDIF  
            ENDIF          
C translation
            lcod = icodt(n)
            IF(lcod > 0) THEN
              icodt(n) = 0
               IF(lcod == 7) THEN
                icodrbym(2,m) = lcod
                itx = 4
                ity = 2
                itz = 1
               ELSE
                 IF(lcod == 4 .OR. lcod == 5 .OR. lcod == 6) itx = 4
                 IF(lcod == 2 .OR. lcod == 3 .OR. lcod == 6) ity = 2
                 IF(lcod == 1 .OR. lcod == 3 .OR. lcod == 5) itz = 1
               ENDIF 
            ENDIF 
          ENDDO 
C rotation       
           lcod = irx + iry + irz
           icodrbym(1,m) = lcod
C C translation          
           lcod = itx + ity + itz
           icodrbym(2,m) = lcod
          nslt =    nslt +  nsl   
      ENDDO
C
      RETURN