m115_perturb.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	m115_perturb (pm, mat, rho, perturb, rnoise)

Function/Subroutine Documentation

m115_perturb()

subroutine m115_perturb	(		pm,
		integer, dimension(*)	mat,
			rho,
		integer, dimension(nperturb)	perturb,
			rnoise )

Definition at line 29 of file m115_perturb.F.

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"
#include      "vect01_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER PERTURB(NPERTURB),MAT(*)
      my_real  
     .   rho(2*llt),rnoise(nperturb,*),pm(npropm,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,J
      my_real
     .        rho0
C-----------------------------------------------
      ! Loop over elements
      DO i = lft,llt
        ! Initial constant density
        rho0   = pm(1,mat(i))
        rho(i) = rho0
        ! If /PERTURB/SOLID/PART card is set
        IF (nperturb /= 0) THEN
          DO j=1,nperturb
            IF (perturb(j) == 3 .AND. rnoise(j,numelc+numeltg+i+nft) /= zero) THEN
              rho(i)  = rho(i) * rnoise(j,numelc+numeltg+i+nft)
            ENDIF
          ENDDO
        ENDIF
        ! Copy of initial density on the second length of RHO
        rho(i+llt) = rho(i)
      ENDDO
C-----------------------------------------------