i24setnodes.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	i24setnodes (ipari, intbuf_tab, intercep, itab, i24maxnsne)
integer function	secnd_surface_on_domain (intercep, se, proc)

Function/Subroutine Documentation

i24setnodes()

subroutine i24setnodes	(	integer, dimension(npari,*)	ipari,
		type(intbuf_struct_), dimension(*)	intbuf_tab,
		type(intersurfp), dimension(3,ninter)	intercep,
		integer, dimension(*)	itab,
		integer	i24maxnsne )

Definition at line 37 of file i24setnodes.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------    
      USE message_mod
      USE i7i10splitmod 
      USE front_mod   
      USE intbufdef_mod    
C-----------------------------------------------
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      "param_c.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com04_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IPARI(NPARI,*),ITAB(*),I24MAXNSNE
      TYPE(INTERSURFP) :: INTERCEP(3,NINTER)
      TYPE(INTBUF_STRUCT_) INTBUF_TAB(*)
C-----------------------------------------------
C   F u n c t i o n
C-----------------------------------------------
      INTEGER  NLOCAL
      EXTERNAL nlocal      
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER ITY,NRTM,N,IED,N1,N2,N3,N4,NRTSE,NSNE
      INTEGER S1,S2,P1,P2,PP,NEXT,POS
      INTEGER NI,K,I,PROC,KM1(4),KM2(4),PROC0,IEDGE4
      DATA km1/1,2,3,4/
      DATA km2/2,3,4,1/
      INTEGER  INTERSURFL
      EXTERNAL intersurfl             
C--------------------------------------------------------------
      i24maxnsne = 0
      DO ni=1,ninter
 
        ity    = ipari(7,ni) 
        iedge4 = ipari(59,ni)
        IF (ity==24)THEN
 
C First set SECONDARY surface to avoid issues : Same SECONDARY surface not on same
C domain than MAIN surface
          nrtse = ipari(52,ni) 
          nsne  = ipari(55,ni) 
          i24maxnsne = max(i24maxnsne,nsne)
          IF (iedge4 > 0 ) THEN
            IF (.NOT.(ASSOCIATED(intercep(2,ni)%P)))THEN
C             In Type24 INTERCEP(2 gives Edge MPI Domain
              ALLOCATE(intercep(2,ni)%P(nrtse))
            ENDIF    
C 
            DO i=1,nrtse
              n1= intbuf_tab(ni)%IRTSE((i-1)*5+1)
              n2= intbuf_tab(ni)%IRTSE((i-1)*5+2)
              n3= intbuf_tab(ni)%IRTSE((i-1)*5+3)
              n4= intbuf_tab(ni)%IRTSE((i-1)*5+4)
 
                proc = intersurfl(n1,n2,n3,n4)
    
              intercep(2,ni)%P(i)=proc
           ENDDO
          ENDIF
 
 
          nrtm   = ipari(4,ni)
          DO k=1,nrtm    
             proc = intercep(1,ni)%P(k)
             IF(proc==0)THEN
                 CALL ancmsg(msgid=978,
     .                       msgtype=msgerror,
     .                       anmode=anstop,
     .                       i1=ipari(15,ni))
             ELSE
               DO i=1,8
                 n = iabs(intbuf_tab(ni)%NVOISIN((k-1)*8+i))
                 IF(n /=0) THEN              
                  CALL ifrontplus(n,proc)
                 ENDIF
               ENDDO
             END IF
          ENDDO
 
 
        ENDIF
      ENDDO 

secnd_surface_on_domain()

integer function secnd_surface_on_domain	(	type(intersurfp)	intercep,
		integer	se,
		integer	proc )

Definition at line 142 of file i24setnodes.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------    
      USE front_mod   
      USE intbufdef_mod    
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      TYPE(INTERSURFP)     :: INTERCEP
       INTEGER SE,PROC
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER NEXT,P
C-----------------------------------------------
      secnd_surface_on_domain = 0
 
      p= intercep%P(2+2*(se-1)+1)
      next=intercep%P(2+2*(se-1)+2) 
 
      IF (p==proc)THEN
         secnd_surface_on_domain = 1
         RETURN
      ENDIF
 
      DO WHILE(p/=proc.AND.next/=-1)
         p=intercep%P(2+2*(next-1)+1)
         next = intercep%P(2+2*(next-1)+2)
      ENDDO
 
      IF(p==proc)THEN
         secnd_surface_on_domain = 1
      ENDIF
      
      RETURN