initnoise.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	initnoise (in, noiadd, ixs, iwa, length, ixq)
subroutine	initnoise2 (in, elnoi, elg, noiadd, ixs, iwa, iparg, ixq)

Function/Subroutine Documentation

initnoise()

subroutine initnoise	(	integer, dimension(*)	in,
		integer, dimension(*)	noiadd,
		integer, dimension(nixs,*)	ixs,
		integer, dimension(0:*)	iwa,
		integer	length,
		integer, dimension(nixq,*)	ixq )

Definition at line 30 of file initnoise.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------        
      use element_mod , only : nixs,nixq
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      "scrnoi_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 IN(*),IXS(NIXS,*),IWA(0:*),IXQ(NIXQ,*)
      INTEGER NOIADD(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,J,LENGTH,K
C-----------------------------------------------
C
C ADDRESS DES NODES SAUVES DS NOISE
C
      DO i=0,numnod
       iwa(i)=0
      ENDDO
      DO i=1,nnoise
       iwa(in(i))=i                 ! In (i) = 0 if not present in SPMD
       noiadd(i)=0
      ENDDO
      noiadd(nnoise+1)=0
C
C STOCKE IN NOIADD(K) LE NUMBER d'ELEMENTS CONNECTES AU NODE K-1
C
      IF(n2d==0)THEN
        DO i=1,numels
         DO j=2,9
          k=iwa(ixs(j,i))
          IF(k/=0)
     .       noiadd(k+1)=noiadd(k+1)+1
         ENDDO
        ENDDO
      ELSE
        DO i=1,numelq
         DO j=2,5
          k=iwa(ixq(j,i))
          IF(k/=0)
     .       noiadd(k+1)=noiadd(k+1)+1
         ENDDO
        ENDDO
      ENDIF      
C
C address per node in the array of elements to be saved in nodes
C
      noiadd(1)=1
      DO i=2,nnoise+1
       noiadd(i)=noiadd(i)+noiadd(i-1)
      ENDDO
      length=noiadd(nnoise+1)-1
c-----------
      RETURN

initnoise2()

subroutine initnoise2	(	integer, dimension(*)	in,
		integer, dimension(*)	elnoi,
		integer, dimension(*)	elg,
		integer, dimension(*)	noiadd,
		integer, dimension(nixs,*)	ixs,
		integer, dimension(*)	iwa,
		integer, dimension(nparg,*)	iparg,
		integer, dimension(nixq,*)	ixq )

Definition at line 105 of file initnoise.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------        
      use element_mod , only : nixs,nixq
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      "scrnoi_c.inc"
#include      "com01_c.inc"
#include      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IN(*),IXS(NIXS,*),IXQ(NIXQ,*),IWA(*),IPARG(NPARG,*),
     .   ELNOI(*),ELG(*),NOIADD(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER , DIMENSION (:), ALLOCATABLE :: KN
      INTEGER I,J,K,II
      INTEGER NG,NEL,ITY,NB2,NFT
C-----------------------------------------------
      IF(nnoise>0)THEN
        ALLOCATE(kn(nnoise))
      ENDIF
      DO i=1,nnoise
        kn(i)=0
      ENDDO
C
C filling elnoi with the addresses from elbuf
C
      DO ng=1,ngroup
        ity   =iparg(5,ng)
        IF(ity==1)THEN
         nel   =iparg(2,ng)
         nb2   =iparg(4,ng) + nel
         nft   =iparg(3,ng)
         DO i=1,nel
          ii=i+nft
          DO j=2,9
            k=iwa(ixs(j,ii))
            IF(k /=0)THEN
              elnoi(noiadd(k)+kn(k))=i
              elg(noiadd(k)+kn(k))=ng
              kn(k)=kn(k)+1
            ENDIF
          ENDDO
         ENDDO
        ELSEIF(ity==2)THEN
         nel   =iparg(2,ng)
         nb2   =iparg(4,ng) + nel
         nft   =iparg(3,ng)
         DO i=1,nel
          ii=i+nft
          DO j=2,5
            k=iwa(ixq(j,ii))
            IF(k /=0)THEN
              elnoi(noiadd(k)+kn(k))=(i-1)*6
              elg(noiadd(k)+kn(k))=ng
              kn(k)=kn(k)+1
            ENDIF
          ENDDO
         ENDDO
        ENDIF
        
      ENDDO
C
      IF(nnoise>0)THEN
        DEALLOCATE(kn)
      ENDIF
C-----------
      RETURN