spmd_cell_size_exchange_init.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	spmd_cell_size_exchange_init (ircvfrom, isendto, ipari, nin, inter_struct, sort_comm)

Function/Subroutine Documentation

spmd_cell_size_exchange_init()

subroutine spmd_cell_size_exchange_init	(	integer, dimension(ninter+1,nspmd+1), intent(in)	ircvfrom,
		integer, dimension(ninter+1,nspmd+1), intent(in)	isendto,
		integer, dimension(npari,ninter), intent(in)	ipari,
		integer, intent(in)	nin,
		type(inter_struct_type), dimension(ninter), intent(inout)	inter_struct,
		type(sorting_comm_type), dimension(ninter), intent(inout)	sort_comm )

Definition at line 32 of file spmd_cell_size_exchange_init.F.

!$COMMENT
!       SPMD_CELL_SIZE_EXCHANGE_INIT description :
!       initialization of buffer size + allocation
!       SPMD_CELL_SIZE_EXCHANGE_INIT organization :
!$ENDCOMMENT
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
        USE inter_sorting_mod
        USE inter_struct_mod
        USE tri7box
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      "task_c.inc"
#include      "sms_c.inc"
#include      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
        INTEGER, INTENT(in) :: NIN
        INTEGER, DIMENSION(NPARI,NINTER), INTENT(in) ::  IPARI
        INTEGER, DIMENSION(NINTER+1,NSPMD+1), INTENT(in) :: ISENDTO,IRCVFROM
        TYPE(inter_struct_type), DIMENSION(NINTER), INTENT(inout) :: INTER_STRUCT   !   structure for interface
        TYPE(sorting_comm_type), DIMENSION(NINTER), INTENT(inout) :: SORT_COMM   ! structure for interface sorting comm
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C----------------------------------------------- 
#ifdef MPI        
        INTEGER :: I,J,NOD,L,L2,KK,IJK,KJI
        INTEGER :: P,P_LOC
        INTEGER :: ADRESS,SHIFT_
        INTEGER :: ISIZ,RSIZ,IDEB,JDEB
        INTEGER :: NSN,NMN,IGAP,INTTH,INTFRIC,ITYP,ITIED
        INTEGER :: IFQ,INACTI
        INTEGER :: LOC_PROC
!   --------------------------------------------------------------------     
 
        loc_proc = ispmd + 1
        sort_comm(nin)%RSIZ = 0
        sort_comm(nin)%ISIZ = 0
 
        igap = ipari(21,nin)
        intth = ipari(47,nin)
        intfric = ipari(72,nin)
        ityp = ipari(7,nin)
        itied = ipari(85,nin)
        nmn = ipari(6,nin)
        nsn = ipari(5,nin)
        inacti = ipari(22,nin)
        ifq =ipari(31,nin)
 
        IF(inacti==5.OR.inacti==6.OR.inacti==7.OR.ifq>0.OR.itied/=0) THEN
            IF(.NOT.ALLOCATED(inter_struct(nin)%NSNFIOLD) ) THEN
                ALLOCATE(inter_struct(nin)%NSNFIOLD(nspmd))
                inter_struct(nin)%NSNFIOLD(1:nspmd) = 0
            ENDIF
            DO p = 1, nspmd
                inter_struct(nin)%NSNFIOLD(p) = nsnfi(nin)%P(p)
            END DO
        END IF
        sort_comm(nin)%NBSEND_NB = 0
        sort_comm(nin)%NBRECV_NB = 0
 
        IF(ircvfrom(nin,loc_proc)/=0.OR.isendto(nin,loc_proc)/=0) THEN
 
            IF(.NOT.ALLOCATED(sort_comm(nin)%REQUEST_NB_R)) THEN
                ALLOCATE(sort_comm(nin)%REQUEST_NB_R(nspmd) )
            ENDIF
            IF(.NOT.ALLOCATED(sort_comm(nin)%REQUEST_NB_S)) THEN
                ALLOCATE(sort_comm(nin)%REQUEST_NB_S(nspmd) )
            ENDIF
            IF(.NOT.ALLOCATED(sort_comm(nin)%NB)) THEN
                ALLOCATE(sort_comm(nin)%NB(nspmd) )
                sort_comm(nin)%NB(1:nspmd) = 0
            ENDIF
            IF(.NOT.ALLOCATED(sort_comm(nin)%INDEX_RCV)) THEN
                ALLOCATE( sort_comm(nin)%INDEX_RCV(sort_comm(nin)%PROC_NUMBER) )
            ENDIF
        
            rsiz = 8    
            isiz = 6
 
c specific cases 
c IGAP=1 or IGAP=2
            IF(igap==1 .OR. igap==2)THEN
                rsiz = rsiz + 1
c IGAP=3
            ELSEIF(igap==3)THEN
                rsiz = rsiz + 2
            ENDIF
 
C thermic      
            IF(intth > 0 ) THEN    
                rsiz = rsiz + 2
                      isiz = isiz + 1
            ENDIF
C Friction      
            IF(intfric > 0 ) THEN    
                      isiz = isiz + 1
            ENDIF
 
C -- IDTMINS==2      
            IF(idtmins == 2)THEN     
                isiz = isiz + 2
C -- IDTMINS_INT /= 0
            ELSEIF(idtmins_int/=0)THEN    
                isiz = isiz + 1
            END IF
 
            sort_comm(nin)%RSIZ = rsiz
            sort_comm(nin)%ISIZ = isiz
 
            IF(isendto(nin,loc_proc)>0) THEN
                IF(.NOT.ALLOCATED(sort_comm(nin)%SEND_NB)) ALLOCATE(sort_comm(nin)%SEND_NB(nspmd))
            ENDIF
            IF(ircvfrom(nin,loc_proc)>0) THEN
                IF(.NOT.ALLOCATED(sort_comm(nin)%RECV_NB)) ALLOCATE(sort_comm(nin)%RECV_NB(nspmd))
            ENDIF
        ENDIF
#endif
        RETURN