link_with_unique_master_cell.F File Reference

#include "implicit_f.inc"
#include "comlock.inc"

Go to the source code of this file.

Functions/Subroutines
integer function	link_with_unique_main_cell (sharedface, centroid)

Function/Subroutine Documentation

link_with_unique_main_cell()

integer function link_with_unique_main_cell	(		sharedface,
			centroid )

Definition at line 33 of file link_with_unique_master_cell.F.

C-----------------------------------------------
C   D e s c r i p t i o n
C-----------------------------------------------
C Interface Type22 (/INTER/TYPE22) is an FSI coupling method based on cut cell method. 
C   This experimental cut cell method is not completed, abandoned, and is not an official option.
C
! With interface 22 Hexahedral cell is cut and
! main one is the biggest one which is greater
! then 50%. In case of multiple choice (2 are
! 50%-50% or several if criteria is reduced) we
! must ensure a unique choice which is repeatable
! with plateform and also different numbering.
!
! Here is chosen a sort by coordinates weighting
! X,Y,Z components. a score is computed
C-----------------------------------------------
C   M o d u l e s
C----------------------------------------------- 
      USE initbuf_mod
      USE i22bufbric_mod 
      USE i22tri_mod
      USE elbufdef_mod
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      "comlock.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER   :: LINK_WITH_UNIQUE_MAIN_CELL
      my_real   :: sharedface(6), centroid(3,6)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER   ::  IsMain, N, IPOS, LEVEL, IIP
      my_real   ::  score(6), max, nextmax, mysum
      my_real   ::  xyzmax(3)
C-----------------------------------------------
C   P r e - C o n d i t i o n
C----------------------------------------------- 
      !IF(INT22==0)RETURN   !already checked
      !IF(K<2)  RETURN   !already checked
C-----------------------------------------------
C   S o u r c e    L i n e s
C----------------------------------------------- 
      ipos              = 1
      level             = 3   !3:X, 2:Y, 1:Z
      score(1:6)        = zero
                     
      ipos      = maxloc(sharedface(1:6),1)
      
      IF(sharedface(ipos) == zero) THEN
        link_with_unique_main_cell = 0
        RETURN
      ENDIF
      
      mysum = sum(sharedface(1:ipos-1))
      IF (ipos/=6)THEN
        mysum = mysum + sum(sharedface(ipos+1:6))
      ENDIF
      IF (mysum == zero)THEN
        link_with_unique_main_cell = ipos
        RETURN
      ENDIF
 
      xyzmax(1) = maxval(centroid(1,1:6))
      xyzmax(2) = maxval(centroid(2,1:6))
      xyzmax(3) = maxval(centroid(3,1:6))            
           
      !find max
      max                          = sharedface(ipos)       !init with maximum value
      link_with_unique_main_cell = ipos
      !score it
      score(ipos) = 1000*(xyzmax(1)-centroid(1,ipos))+ 100*(xyzmax(2)-centroid(2,ipos))+ 10*(xyzmax(3)-centroid(3,ipos))      
      
      !find other ones equal to MAX
      IF (ipos == 6)RETURN
      iip = ipos
      DO WHILE(iip < 6)
        iip = maxloc(sharedface(iip+1:6),1)
        nextmax = sharedface(iip)
        IF (nextmax == max) THEN
          score(ipos) = 1000*(xyzmax(1)-centroid(1,iip))+100 * (xyzmax(2)-centroid(2,iip)) + 10*(xyzmax(3)-centroid(3,iip))
        ENDIF
      ENDDO
 
      ipos                         = maxloc(score(1:6),1)
      link_with_unique_main_cell = ipos
      
      RETURN