int18_alloc.F File Reference

#include "implicit_f.inc"
#include "comlock.inc"
#include "com04_c.inc"
#include "param_c.inc"
#include "task_c.inc"
#include "parit_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	int18_alloc (number_inter18, inter18_list, multi_fvm, ipari, xcell_remote, nspmd)

Function/Subroutine Documentation

int18_alloc()

subroutine int18_alloc	(	integer, intent(inout)	number_inter18,
		integer, dimension(ninter), intent(inout)	inter18_list,
		type(multi_fvm_struct)	multi_fvm,
		integer, dimension(npari,ninter), intent(in)	ipari,
		type(array_type), dimension(ninter), intent(inout)	xcell_remote,
		integer, intent(in)	nspmd )

Parameters

[in]	nspmd	number of spmd
[in,out]	number_inter18	number of interface 18
[in,out]	inter18_list	list of interface 18
[in,out]	xcell_remote	remote data structure for interface 18

Definition at line 33 of file int18_alloc.F.

!$COMMENT
!       INT18_ALLOC description
!       allocation of array for interface 18 & interface 18 combined 
!       with law151
!       
!       int18_alloc organization :
!       - check if /INT18 is used
!       - check if /INT18 + /LAW151 is used
!       - allocate the arrays
!$ENDCOMMENT
C-----------------------------------------------
C     M o d u l e s
C-----------------------------------------------           
        USE multi_fvm_mod
        USE groupdef_mod
        USE tri7box
        USE array_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
#include      "comlock.inc"
#include      "com04_c.inc"
#include      "param_c.inc"
#include      "task_c.inc" 
#include      "parit_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
        INTEGER, INTENT(in) :: NSPMD !< number of spmd
        INTEGER, INTENT(inout) :: NUMBER_INTER18 !< number of interface 18
        INTEGER, DIMENSION(NINTER), INTENT(inout) :: INTER18_LIST !< list of interface 18
        INTEGER, DIMENSION(NPARI,NINTER), INTENT(in) :: IPARI
        TYPE(MULTI_FVM_STRUCT) :: MULTI_FVM
        TYPE(array_type), DIMENSION(NINTER), INTENT(inout) :: XCELL_REMOTE !< remote data structure for interface 18 
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
        INTEGER :: N,NN,II,JJ,MY_SIZE,MY_SIZE_2,MY_SIZE_3
        INTEGER :: ISU1,NBRIC,NSN,NTY,INACTI,NODE_ID,IBRIC,NODFI
        INTEGER :: P
C-----------------------------------------------
        !   check if int18 + law151 is used in the model
        !   and create a list of int18
        multi_fvm%IS_INT18_LAW151 = .false.
        my_size = 0
        my_size_2 = 0
        number_inter18 = 0
        
        DO n=1,ninter
            nty   =ipari(7,n)
            inacti = ipari(22,n)           
            xcell_remote(n)%SIZE_MY_REAL_ARRAY_1D = 0
            !   int18 = int7 + inacti=7 (7+7=18) 
            IF( (nty==7).AND.(inacti ==7)) THEN
              IF(multi_fvm%IS_USED) THEN
                multi_fvm%IS_INT18_LAW151 = .true.
                my_size = numnod + numels
                my_size_2 = numels
            ENDIF      
            number_inter18 = number_inter18 + 1
            inter18_list(number_inter18) = n  !   list of interface int18
          ENDIF
        ENDDO
        !   number & list of interface 18 for MULTI_FVM solve
        multi_fvm%NUMBER_INT18 = 0
        IF(multi_fvm%IS_INT18_LAW151) multi_fvm%NUMBER_INT18 = number_inter18
        ALLOCATE( multi_fvm%INT18_LIST(multi_fvm%NUMBER_INT18) )
        !   allocation of X/V/MASS extended to NUMNOD+NUMELS
        !   1:NUMNOD --> classical x/v/mass
        !   NUMNOD+1:NUMNOD+NUMELS --> x/v/mass of phantom nodes (located to the center of 
        !                              the ALE elements)
        ALLOCATE( multi_fvm%X_APPEND(3*my_size) )
        ALLOCATE( multi_fvm%V_APPEND(3*my_size) )
        ALLOCATE( multi_fvm%MASS_APPEND(my_size) )
        ALLOCATE( multi_fvm%KINET_APPEND(my_size) )
        !   allocation of force array : size = NUMELS
        my_size_2 = my_size_2 * nthread
        
        IF(iparit/=0) THEN
            my_size_3 = 0 
        ELSE
            my_size_3 = my_size_2
        ENDIF
 
        multi_fvm%SIZE_FORCE_INT_1 = 3
        multi_fvm%SIZE_FORCE_INT_2 = my_size_3
        ALLOCATE( multi_fvm%FORCE_INT(3,my_size_3) )
        !   allocation of INT18_GLOBAL_LIST : marker for the interface /INT18+LAW151
        ALLOCATE( multi_fvm%INT18_GLOBAL_LIST(ninter) )
 
        ! --------------------------
        ! allocation of remote array for parith/on
        IF(multi_fvm%IS_USED) THEN
            ALLOCATE( multi_fvm%R_AFI(ninter) )  
            DO ii=1,multi_fvm%NUMBER_INT18
                n = inter18_list(ii) ! list of interface 18 + law151
                IF( ALLOCATED( multi_fvm%R_AFI(n)%R_FORCE_INT ) ) DEALLOCATE( multi_fvm%R_AFI(n)%R_FORCE_INT )
                nodfi = 0
                DO p = 1,nspmd
                  nodfi = nodfi + nsnfi(n)%P(p)
                ENDDO
                multi_fvm%R_AFI(n)%NODFI = nodfi
                ALLOCATE( multi_fvm%R_AFI(n)%R_FORCE_INT(3,6,nodfi*nthread) )
                multi_fvm%R_AFI(n)%R_FORCE_INT(1:3,1:6,1:nodfi*nthread) = 0d+00 
            ENDDO
        ELSE
            ALLOCATE( multi_fvm%R_AFI(0) )
        ENDIF 
        ! allocation of local array for parith/on
        IF(iparit/=0) THEN
            ALLOCATE( multi_fvm%FORCE_INT_PON(3,6,my_size_2) )
            multi_fvm%SIZE_FORCE_INT_PON = my_size_2
        ELSE
            ALLOCATE( multi_fvm%FORCE_INT_PON(0,0,0) )
            multi_fvm%SIZE_FORCE_INT_PON = 0
        ENDIF
        ! --------------------------
 
        IF( multi_fvm%IS_INT18_LAW151 ) THEN
            multi_fvm%INT18_LIST(1:multi_fvm%NUMBER_INT18) = inter18_list(1:multi_fvm%NUMBER_INT18)
        ENDIF
 
        RETURN