pnoise.F File Reference

#include "implicit_f.inc"
#include "scrnoi_c.inc"
#include "com01_c.inc"
#include "task_c.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	pnoise (elnoi, elg, noiadd, elbuf_tab, wa, iparg)

Function/Subroutine Documentation

pnoise()

subroutine pnoise	(	integer, dimension(*)	elnoi,
		integer, dimension(*)	elg,
		integer, dimension(*)	noiadd,
		type (elbuf_struct_), dimension(ngroup), target	elbuf_tab,
			wa,
		integer, dimension(nparg,*)	iparg )

Definition at line 33 of file pnoise.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE elbufdef_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      "scrnoi_c.inc"
#include      "com01_c.inc"
#include      "task_c.inc"
#include      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER ELNOI(*),ELG(*),NOIADD(*),IPARG(NPARG,*)
      my_real
     .   wa(*)
      TYPE (ELBUF_STRUCT_), DIMENSION(NGROUP), TARGET :: ELBUF_TAB
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,K,NG,IWA(NNOISE),NEL,IEL
      TYPE(G_BUFEL_)  ,POINTER :: GBUF     
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7--
C
      DO i=1,nnoise
       wa(i)=zero
      ENDDO
      DO i=1,nnoise
       DO k=noiadd(i),noiadd(i+1)-1
        ng = elg(k)
        gbuf => elbuf_tab(ng)%GBUF
        nel = iparg(2,ng)
        iel = elnoi(k)
        wa(i)=wa(i)+ gbuf%SIG(iel) + gbuf%SIG(iel+nel) + gbuf%SIG(iel+2*nel)
       ENDDO       
      ENDDO
      IF(nspmd==1) THEN
        DO i=1,nnoise
          k=noiadd(i+1)-noiadd(i)
          IF(k>0)wa(i)=wa(i)/three/float(k)
        ENDDO
      ELSE
        CALL spmd_glob_dsum9(wa,nnoise)
        DO i=1,nnoise
          iwa(i) = noiadd(i+1)-noiadd(i)
        ENDDO
        CALL spmd_glob_isum9(iwa,nnoise)
        IF(ispmd==0) THEN
          DO i=1,nnoise
            IF(iwa(i)>0)wa(i)=wa(i)/three/float(iwa(i))
          ENDDO
        ENDIF
      ENDIF
C-----------
      RETURN