sphres44b.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "vect01_c.inc"
#include "com01_c.inc"
#include "param_c.inc"
#include "sphcom.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	sphres44b (kxsp, ixsp, nod2sp, iparg, spbuf)

Function/Subroutine Documentation

sphres44b()

subroutine sphres44b	(	integer, dimension(nisp,*)	kxsp,
		integer, dimension(kvoisph,*)	ixsp,
		integer, dimension(*)	nod2sp,
		integer, dimension(nparg,*)	iparg,
			spbuf )

Definition at line 28 of file sphres44b.F.

C----6---------------------------------------------------------------7---------8
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   G l o b a l   P a r a m e t e r s
C-----------------------------------------------
#include      "mvsiz_p.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "vect01_c.inc"
#include      "com01_c.inc"
#include      "param_c.inc"
#include      "sphcom.inc"
C-----------------------------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER 
     .    KXSP(NISP,*),IXSP(KVOISPH,*),NOD2SP(*),
     .    IPARG(NPARG,*)
C     REAL
      my_real
     .   spbuf(nspbuf,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER N,IG,NG,I,OFFSET,NEL
C-----------------------------------------------
C     Construit les Numeros d'appartenance aux groupes. 
C-----------------------------------------------
      DO ng=1,ngroup
      nel=iparg(2,ng)
      nft=iparg(3,ng)
      ity=iparg(5,ng)
      DO ig= 1,nel,mvsiz
        offset = ig - 1
        lft=1
        llt=min(mvsiz,nel-offset)  
        IF(ity==51)THEN
         DO i=lft,llt
          n =i+nft
          IF(spbuf(9,n)>zero)THEN
           kxsp(1,n) =0
           kxsp(2,n) =ng
           spbuf(9,n)=zero     
          ENDIF
         ENDDO
        ENDIF
      ENDDO
      ENDDO
C-----------------------------------------------
      RETURN