parsor0.F File Reference

#include "implicit_f.inc"
#include "com01_c.inc"
#include "com04_c.inc"
#include "sphcom.inc"
#include "param_c.inc"
#include "task_c.inc"
#include "spmd_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	parsor0 (iadd, iparg, mater, el2fa, dd_iad, iadg, kxsp, ipartsp, nodglob)

Function/Subroutine Documentation

parsor0()

subroutine parsor0	(	integer, dimension(*)	iadd,
		integer, dimension(nparg,*)	iparg,
		integer, dimension(*)	mater,
		integer, dimension(*)	el2fa,
		integer, dimension(nspmd+1,*)	dd_iad,
		integer, dimension(nspmd,*)	iadg,
		integer, dimension(nisp,*)	kxsp,
		integer, dimension(*)	ipartsp,
		integer, dimension(*)	nodglob )

Definition at line 32 of file parsor0.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      "sphcom.inc"
#include      "param_c.inc"
#include      "task_c.inc"
#include      "spmd_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
C     REAL
      INTEGER IADD(*),IPARG(NPARG,*),
     .        MATER(*),EL2FA(*),
     .        IADG(NSPMD,*),
     .        DD_IAD(NSPMD+1,*),
     .        KXSP(NISP,*),IPARTSP(*),
     .        NODGLOB(*)
C-----------------------------------------------
C     REAL
      INTEGER II, IE, NG, ITY, LFT, LLT, N, I, 
     .        NEL, IAD, NPAR, NFT, IPRT,
     . NN1,NN2,
     .        JJ, INOD,
     .        NP(NUMSPH),BUF,BUFP
C-----------------------------------------------
C
      nn1 = 1
      nn2 = nn1 + numsph+maxpjet
      ie = 0
C-----------------------------------------------
C     MID
C-----------------------------------------------
      npar = 0
      jj = 0
C
      DO 100 iprt = 1 , npart
C
       IF(mater(iprt)/=4) GOTO 100
       npar = npar + 1
       DO 90 ng=1,ngroup
        nel   =iparg(2,ng)
        nft   =iparg(3,ng)
        iad   =iparg(4,ng)
        ity   =iparg(5,ng)
        lft=1
        llt=nel
        IF(ity==51)THEN
C-----------------------------------------------
C         PARTICULES SPH.
C-----------------------------------------------
          DO 20 i=lft,llt
           n = i + nft
           IF(ipartsp(n)/=iprt) GOTO 20
             inod=kxsp(3,n)
             IF (nspmd == 1) THEN 
               ii = inod-1
               CALL write_i_c(ii,1)
             ELSE
               np(jj+1) = nodglob(inod)-1
             END IF
             jj = jj+1
             ie = ie + 1
             el2fa(nn1+n) = ie
 20      CONTINUE
        ELSE
        ENDIF
 90    CONTINUE
C-----------------------------------------------
C       PART ADRESS
C-----------------------------------------------
       iadd(npar) = ie
 100  CONTINUE
C
      IF (nspmd > 1) THEN
        IF (ispmd==0) THEN
          bufp = npart
          buf = numsphg
        ELSE
          bufp = 1
          buf=1
        END IF
 
          CALL spmd_iglob_partn(iadd,npar,iadg,bufp)
          CALL spmd_iget_partn(1,jj,np,npar,iadg,buf,1)
      ELSE ! IADG filling for mono/multi compatibility
        DO i = 1, npart
          iadg(1,i) = iadd(i)
        END DO
      ENDIF
C-----------------------------------------------
      RETURN