preplyxfem.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	preplyxfem (ms_ply0, zi_ply0, iel, inod, ixc, ms_ply, zi_ply, addcne, msz20, msz2)
subroutine	fillcne_pxfem (iel, inod, ixc, cep, addcne, cne, cel)

Function/Subroutine Documentation

fillcne_pxfem()

subroutine fillcne_pxfem	(	integer, dimension(*)	iel,
		integer, dimension(*)	inod,
		integer, dimension(nixc,*)	ixc,
		integer, dimension(*)	cep,
		integer, dimension(0:*)	addcne,
		integer, dimension(*)	cne,
		integer, dimension(*)	cel )

Definition at line 102 of file preplyxfem.F.

      use element_mod , only : nixc
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      "com01_c.inc"
#include      "com04_c.inc"
#include      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IXC(NIXC,*),ADDCNE(0:*), CNE(*),IEL(*),INOD(*),
     .        CEP(*),CEL(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, J, K, N, ADSKY(0:NPLYXFE+1),N1,PROC,NG,
     .        INDX, INDEX(NUMELC),P,II,NIN,INDEXC(2*NUMELC),
     .        WORK(70000),ITRI(NUMELC)
C-----------------------------------------------
C   S o u r c e  L i n e s
C-----------------------------------------------
C   CALCULATION OF CNE ADDCNE
C-----------------------------------------------
      DO i = 0, nplyxfe + 1
        adsky(i) = addcne(i)
      ENDDO
C
C
C sorting the following local elements Num User
C
      DO i = 1, numelc
        itri(i) = ixc(7,i)
      ENDDO
      CALL my_orders(0,work,itri,indexc,numelc,1)
      
      indx = 0
      DO j=1,numelc
        i = indexc(j)
!!      DO   I = 1,NUMELC
       IF(iel(i) > 0 )THEN
         indx = indx + 1
         DO k=1,4
           n = ixc(k+1,i)
           n1 = inod(n)
           cne(adsky(n1)) = i
           adsky(n1) = adsky(n1) + 1
         ENDDO
         index(indx) = i
       ENDIF 
      ENDDO
C     
      DO proc = 1, nspmd
        nin = 0  
!!        DO I=1,INDX
          DO ii=1,numelc
!!           II = INDEX(I)
           ng = ii +  numels + numelq
           p = cep(ng) + 1 
           IF (p == proc) THEN
            j = iel(ii) 
            IF(j > 0) THEN 
              nin = nin + 1
              cel(j) = nin
            ENDIF  
           ENDIF
        ENDDO 
      ENDDO  
C      
      RETURN

preplyxfem()

subroutine preplyxfem	(		ms_ply0,
			zi_ply0,
		integer, dimension(*)	iel,
		integer, dimension(*)	inod,
		integer, dimension(nixc,*)	ixc,
			ms_ply,
			zi_ply,
		integer, dimension(0:*)	addcne,
			msz20,
			msz2 )

Definition at line 29 of file preplyxfem.F.

      use element_mod , only : nixc
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      "com04_c.inc"
#include      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER  IXC(NIXC,*), ADDCNE(0:*),IEL(*),INOD(*)
      my_real 
     .        ms_ply0(numnod,*),  zi_ply0(numnod,*),
     .        ms_ply(nplyxfe,*),  zi_ply(nplyxfe,*),
     .        msz20(*),msz2(*) 
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,II,K,J,N
      my_real
     .  msz
C-----------------------------------------------
C   S o u r c e  L i n e s
C-----------------------------------------------
         addcne(0)=0
         addcne(1)= 0
         DO i=1,numnod
            ii = inod(i)
            addcne(i + 1) = 0
            msz = zero
            IF( ii > 0 )THEN
              msz2(ii) = zero
              DO j=1,nplymax
                 ms_ply(ii,j) = ms_ply0(i,j)
                 zi_ply(ii,j) = zi_ply0(i,j)
C  ementary compute
cc                 MSZ2(II) = MSZ20(I)  
c nodale compute                 
                 msz2(ii) = msz2(ii) + 
     .                      ms_ply(ii,j)*zi_ply(ii,j)*zi_ply(ii,j)
              ENDDO
               IF(msz2(ii) == zero) msz2(ii) = ep30
            ENDIF
         ENDDO    
C
         DO  i=1,numelc
            IF(iel(i) > 0 ) THEN
              DO  k=2,5
                n = inod(ixc(k,i)) + 1
                addcne(n) = addcne(n) + 1
              ENDDO
            ENDIF
          ENDDO 
          addcne(1) = 1
          DO i=2,nplyxfe + 1
             addcne(i) = addcne(i) + addcne(i-1)
          END DO        
C          
          RETURN