bhol2.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	bhol2 (itask, iparg, pm, bhole, ixq)

Function/Subroutine Documentation

bhol2()

subroutine bhol2	(	integer	itask,
		integer, dimension(nparg,ngroup)	iparg,
			pm,
		integer, dimension(nmult,*)	bhole,
		integer, dimension(nixq,numelq)	ixq )

Definition at line 28 of file bhol2.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      "com01_c.inc"
#include      "com04_c.inc"
#include      "vect01_c.inc"
#include      "param_c.inc"
#include      "task_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER ITASK,IPARG(NPARG,NGROUP),BHOLE(NMULT,*),IXQ(NIXQ,NUMELQ)
      my_real pm(npropm,nummat)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER NM, NG, JMUL, I, J, MX
C-----------------------------------------------
C   S o u r c e   L i n e s
C------------------------------------------------
      DO nm=1,nmult
       DO ng=itask+1,ngroup,nthread
C       ALE ON / OFF
        IF (iparg(76, ng)  ==  1) cycle ! --> OFF
        jale=iparg(7,ng)
        jeul=iparg(11,ng)
        IF(jale+jeul == 0)   cycle
        IF(iparg(8,ng) == 1) cycle
        jmul=iparg(20,ng)
        llt=iparg(2,ng)
        nft=iparg(3,ng)
        lft=1
        IF(jmul == 0)THEN
          DO i=lft,llt
           j=i+nft
           mx=ixq(1,j)
           bhole(nm,j)=mx
          ENDDO
C       ELSEIF(JMUL<NM)THEN
C         DO I=LFT,LLT
C          J=I+NFT
C          MX=IXQ(1,J)
C          BHOLE(NM,J)=0
C         ENDDO
        ELSE
          DO i=lft,llt
           j=i+nft
           mx=ixq(1,j)
           bhole(nm,j)=nint(pm( 20+nm,mx))
           IF(nint(pm(185+nm,mx)) /= 0)bhole(nm,j)=-bhole(nm,j)
          ENDDO
        ENDIF
       enddo! next NG
      ENDDO !next NM=1,NMULT
C-----------------------------------------------      
      RETURN