reordr.F File Reference

#include "implicit_f.inc"
#include "param_c.inc"
#include "tablen_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	reordr (ix, nx, nel, pm, ipoint, iparts, ngrele, igrelem, nummat)

Function/Subroutine Documentation

reordr()

subroutine reordr	(	integer, dimension(nx,nel)	ix,
		integer	nx,
		integer	nel,
			pm,
		integer, dimension(2,*)	ipoint,
		integer, dimension(*)	iparts,
		integer	ngrele,
		type (group_), dimension(ngrele)	igrelem,
		integer, intent(in)	nummat )

Definition at line 30 of file reordr.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE groupdef_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "param_c.inc"
#include      "tablen_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER :: NX, NEL,NGRELE
      INTEGER :: IX(NX,NEL), IPOINT(2,*),IPARTS(*)
      INTEGER, INTENT(IN) :: NUMMAT
      my_real :: pm(npropm,nummat)
C-----------------------------------------------
      TYPE (GROUP_)  , DIMENSION(NGRELE) :: IGRELEM
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER IADR(0:MAXLAW+1), NL, M, I, MLAWP1, MT, MLAW, J,K,N,IAD
C-----------------------------------------------
C   S o u r c e   L i n e s
C-----------------------------------------------
      nl=maxlaw
      DO m=0,nl+1
        iadr(m)=0
      ENDDO
C---------------------------------------
C     NUMBER OF ELEMS FOR EACH MATERIAL LAW
C---------------------------------------
      DO i=1,nel
       IF(ix(1,i)<0) THEN
         mlawp1=17+1
       ELSE
         mt=ix(1,i)
         mlawp1=nint(pm(19,mt))+1
       ENDIF
       iadr(mlawp1)=iadr(mlawp1)+1
      ENDDO
C------------------------------------------
C     ELEM ADDRESSE FOR EACH MATERIAL LAW
C     SKIPPING LAW11
C------------------------------------------
      iadr(11+1)=0
      DO m=1,nl
        iadr(m)=iadr(m)+iadr(m-1)
      ENDDO
      iadr(11)=iadr(nl+1)+iadr(nl)
C--------------------------------------------------
C     CALL POINTER FOR EACH MLW
C--------------------------------------------------
      DO i=1,nel
        IF(ix(1,i)<0) THEN
          mlaw=17
        ELSE
          mt=ix(1,i)
          mlaw=nint(pm(19,mt))
        ENDIF
        iadr(mlaw)=iadr(mlaw)+1
        ipoint(1,i)=iadr(mlaw)
      END DO
C------------------------------------------
C     PERUMTATING ELEMS USING LAWS
C------------------------------------------
      DO j=1,nx
        DO i=1,nel
          ipoint(2,i)=ix(j,i)
        END DO
        DO i=1,nel
          ix(j,ipoint(1,i))=ipoint(2,i)
        END DO
      ENDDO
C------------------------------------------
C     PERMUTATING PARTS
C------------------------------------------
      DO i=1,nel
        ipoint(2,i)=iparts(i)
      ENDDO
      DO i=1,nel
        iparts(ipoint(1,i))=ipoint(2,i)
      ENDDO
C------------------------------------------
C     RRENUMBERING ELEMS INSIDE GROUPS
C------------------------------------------
      DO k=1,ngrele
        n = igrelem(k)%NENTITY
        DO i=1,n
          IF (igrelem(k)%ENTITY(i)<=nel)igrelem(k)%ENTITY(i) = ipoint(1,igrelem(k)%ENTITY(i))
        ENDDO
      ENDDO
C------------------------------------------
      RETURN