c_icnds10.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	c_icnds10 (icnds10, itagnd, proc, ns10e_l, nbddcndm, ms, ms_nd)

Function/Subroutine Documentation

c_icnds10()

subroutine c_icnds10	(	integer, dimension(3,*)	icnds10,
		integer, dimension(*)	itagnd,
		integer	proc,
		integer	ns10e_l,
		integer	nbddcndm,
			ms,
			ms_nd )

Definition at line 30 of file c_icnds10.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
 
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      "com01_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
       INTEGER  ICNDS10(3,*),PROC,NS10E_L,ITAGND(*),NBDDCNDM
      my_real
     .   ms(*),ms_nd
C-----------------------------------------------
C   F u n c t i o n
C-----------------------------------------------
      INTEGER  NLOCAL
      EXTERNAL nlocal 
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER N, NN,N1,N2,P,NF
      INTEGER, DIMENSION(:), ALLOCATABLE :: TAGI
C-----------------------------------------------
C------ COMPUTE MS_ND per PROC
        ms_nd = zero     
      DO n = 1, ns10e
       nn = icnds10(1,n)
        DO p = 1, proc
          IF(nlocal(nn,p)==1)GOTO 10
        ENDDO
       IF(nlocal(nn,proc+1)==1.AND.itagnd(nn)<=ns10e)THEN
        ms_nd = ms_nd + ms(nn)
       ENDIF
 10     CONTINUE
      ENDDO
C
      DO n = 1, ns10e
       nn = icnds10(1,n)
       IF(nlocal(nn,proc+1)==1.AND.itagnd(nn)<=ns10e)THEN
c        DO P = 1, PROC
c         IF(NLOCAL(NN,P)==1) GO TO 100           
c        ENDDO
        n1 = icnds10(2,n)
        n2 = icnds10(3,n)
        ns10e_l = ns10e_l + 1
 100    CONTINUE         
       ENDIF
      ENDDO
      IF (ns10e_l==0) RETURN
!     allocate 1d array
      ALLOCATE( tagi(numnod) )
      tagi(1:numnod) = 0 
C-------------pour main nodes---------
      DO n = 1, ns10e
        nn = icnds10(1,n)
        IF(nlocal(nn,proc+1)/=1.OR.itagnd(nn)>ns10e) cycle
        n1 = icnds10(2,n)
        n2 = icnds10(3,n)
C----- normally N1,N2 are local as  NN-----     
         nf = 0
         DO p = 1, nspmd
          nf = nf +nlocal(nn,p)           
         ENDDO
            IF(nlocal(n1,proc+1)==1) THEN
              IF(tagi(n1)==0.AND.nf > 1) THEN
C decompte des frontieres de main nodes
                DO p = 1, nspmd
                  IF(nlocal(n1,p)==1)THEN
                    nbddcndm = nbddcndm + 1
                  END IF
                END DO
C on ne se compte pas soi-meme
                nbddcndm = nbddcndm - 1
C pour ne pas prendre en compte 2 fois des noeuds main ds les frontieres
c                IF (NF > 1) NBDDCNDM = NBDDCNDM + 1
                tagi(n1) = 1
              END IF
            ENDIF
C----- N2 -------       
            IF(nlocal(n2,proc+1)==1) THEN
              IF(tagi(n2)==0.AND.nf > 1) THEN
C decompte des frontieres de main nodes
                DO p = 1, nspmd
                  IF(nlocal(n2,p)==1)THEN
                    nbddcndm = nbddcndm + 1
                  END IF
                END DO
C on ne se compte pas soi-meme
                nbddcndm = nbddcndm - 1
C pour ne pas prendre en compte 2 fois des noeuds main ds les frontieres
c                IF (NF > 1) NBDDCNDM = NBDDCNDM + 1
                tagi(n2) = 1
              END IF
            ENDIF
      END DO
C ----------------------------
      DEALLOCATE( tagi )
C ----------------------------
C
      RETURN