c_isph.F File Reference

#include "implicit_f.inc"
#include "com04_c.inc"
#include "sphcom.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	c_isph (kxsp, cepsp, proc, numsph_l, slonfsph_l, slprtsph, slprtsph_l, ipartsp, ssphveln_l, nsphsol_l, first_sphsol_l)

Function/Subroutine Documentation

c_isph()

subroutine c_isph	(	integer, dimension(nisp,*)	kxsp,
		integer, dimension(*)	cepsp,
		integer	proc,
		integer	numsph_l,
		integer	slonfsph_l,
		integer	slprtsph,
		integer	slprtsph_l,
		integer, dimension(*)	ipartsp,
		integer	ssphveln_l,
		integer	nsphsol_l,
		integer	first_sphsol_l )

Definition at line 29 of file c_isph.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      "sphcom.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER PROC, NUMSPH_L, NSPHSOL_L, FIRST_SPHSOL_L,
     .        KXSP(NISP,*), CEPSP(*)
 
c inlet outlet spmd
      INTEGER SLONFSPH_L, IPRT, N, IPARTSP(*),SLPRTSPH_L,SLPRTSPH,
     .SSPHVELN_L
 
C-----------------------------------------------
C   L o c a l  V a r i a b l e s
C-----------------------------------------------
      INTEGER I, NB_L, IUN
      my_real
     .        aux
      DATA iun/1/
C
      nb_l = 0
C
      DO i = 1, numsph
        IF(cepsp(i)==proc) THEN
          nb_l = nb_l + 1
        ENDIF
      ENDDO
C
      numsph_l = nb_l
C
c inlet spmd calcul SLONFSPH_L
      slonfsph_l = 0
 
      DO iprt=1,npart
       DO n=1,numsph
        IF(cepsp(n)==proc.and.ipartsp(n)==iprt.
     .     and.kxsp(2,n)/=0)THEN
         slonfsph_l=slonfsph_l+1
        ENDIF
       ENDDO
      ENDDO
 
      slprtsph_l=slprtsph
      ssphveln_l=nsphio*numsph_l*2
C
      nb_l = 0
      DO i = 1, first_sphsol-1
        IF(cepsp(i)==proc) THEN
          nb_l = nb_l + 1
        ENDIF
      ENDDO
      first_sphsol_l = nb_l+1
C
      nb_l = 0
      DO i = 1, nsphsol
        IF(cepsp(first_sphsol+i-1)==proc) THEN
          nb_l = nb_l + 1
        ENDIF
      ENDDO
      nsphsol_l = nb_l
 
      RETURN