voln22.F File Reference

#include "implicit_f.inc"
#include "inter22.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	voln22 (voln, iad22, ixs, lft, llt, jale, jeul)

Function/Subroutine Documentation

voln22()

subroutine voln22	(	dimension(*), intent(inout)	voln,
		dimension(*), intent(in)	iad22,
		integer, dimension(nixs,*), intent(in)	ixs,
		integer, intent(inout)	lft,
		integer, intent(inout)	llt,
		integer, intent(inout)	jale,
		integer, intent(inout)	jeul )

Definition at line 32 of file voln22.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
C this subroutines update VOLN for cut cells.
C Main cells have extended values of volume due
C to secnd cell merging. This extension data is 
C stored in cut cell buffer and is now used here
C to update VOLN(*)
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      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      "inter22.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(INOUT) :: LFT
      INTEGER, INTENT(INOUT) :: LLT
      INTEGER, INTENT(INOUT) :: JALE
      INTEGER, INTENT(INOUT) :: JEUL
      my_real,INTENT(INOUT) :: voln(*)
      my_real,INTENT(IN)    :: iad22(*)
      INTEGER,INTENT(IN)    :: IXS(NIXS,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER :: I, II, IB, NCELL, IPOS, MCELL, NIN , NSecnds, IBV, ICELLv, NumSecnd,ISECND
      my_real :: var, aream, volm, volsecnds
C-----------------------------------------------
C   P r e - C o n d i t i o n s
C-----------------------------------------------
      IF(int22==0)RETURN
C-----------------------------------------------
C   S o u r c e   L i n e s
C-----------------------------------------------
      nin = 1
      
      IF(jeul+jale/=0)THEN    
        DO i=lft,llt
        
           ib        =  nint(iad22(i))
           IF(ib==0)cycle
           ncell     =  brick_list(nin,ib)%NBCUT      
           mcell     =  brick_list(nin,ib)%MainID
           voln(i)   =  em20*brick_list(nin,ib)%UncutVOL  !!!ZERO !in case of no main cell                               
           IF(mcell == 0)cycle          
          
           IF(brick_list(nin,ib)%POLY(mcell)%Vnew>zero)THEN
             voln(i)    =  brick_list(nin,ib)%Vnew_SCell
           ENDIF
                     
        ENDDO
      ENDIF
 
C-----------------------------------------------
      RETURN