multi_ebcs.F File Reference

#include "implicit_f.inc"
#include "param_c.inc"
#include "tabsiz_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	multi_ebcs (itask, multi_fvm, ixs, ixq, ixtg, xgrid, wgrid, ipm, pm, func_value, ebcs_tab, npf, tf, fsavsurf, nsurf, timestep, numels, numelq, numeltg, numnod, ncycle, nummat, matparam, output, pred)

Function/Subroutine Documentation

multi_ebcs()

subroutine multi_ebcs	(	integer, intent(in)	itask,
		type(multi_fvm_struct), intent(inout)	multi_fvm,
		integer, dimension(nixs, *), intent(in)	ixs,
		integer, dimension(nixq, *), intent(in)	ixq,
		integer, dimension(nixtg, *), intent(in)	ixtg,
		dimension(3, *), intent(in)	xgrid,
		dimension(3, *), intent(in)	wgrid,
		integer, dimension(npropmi, *), intent(in)	ipm,
		dimension(npropm, *), intent(in)	pm,
		dimension(*), intent(in)	func_value,
		type(t_ebcs_tab), intent(in), target	ebcs_tab,
		integer, dimension(snpc), intent(in)	npf,
		dimension(stf), intent(in)	tf,
		dimension(th_surf_num_channel,nsurf), intent(inout)	fsavsurf,
		integer, intent(in)	nsurf,
		intent(in)	timestep,
		integer, intent(in)	numels,
		integer, intent(in)	numelq,
		integer, intent(in)	numeltg,
		integer, intent(in)	numnod,
		integer, intent(in)	ncycle,
		integer, intent(in)	nummat,
		type(matparam_struct_), dimension(nummat), intent(in)	matparam,
		type(output_), intent(inout)	output,
		logical, intent(in)	pred )

Parameters

[in]	nummat	array size
[in,out]	output	output structure
[in]	pred	predictor step if true, corrector step if false

Definition at line 40 of file multi_ebcs.F.

C-----------------------------------------------
C     M o d u l e s
C-----------------------------------------------
      USE multi_fvm_mod
      USE ebcs_mod
      USE th_surf_mod , only : th_surf_num_channel
      USE matparam_def_mod , ONLY : matparam_struct_
      USE output_mod , ONLY : output_
      USE multi_inlet_ebcs_mod , ONLY : multi_inlet_ebcs
      use element_mod , only : nixs,nixq,nixtg
      !USE multi_propellant_ebcs_mod , ONLY : MULTI_PROPELLANT_EBCS
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-----------------------------------------------
! NEBCS
! NPEBC, NIXS
#include      "param_c.inc"
! LISURF1
! SNPC,STF
#include      "tabsiz_c.inc"
C-----------------------------------------------
C     D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER,INTENT(IN) :: NUMMAT !< array size
      TYPE(MATPARAM_STRUCT_), DIMENSION(NUMMAT), INTENT(IN) :: MATPARAM !material data structure
      INTEGER,INTENT(IN) :: NCYCLE
      INTEGER,INTENT(IN) :: NSURF
      INTEGER,INTENT(IN) :: NUMELS, NUMELQ, NUMELTG, NUMNOD
      my_real,INTENT(INOUT) :: fsavsurf(th_surf_num_channel,nsurf)
      TYPE(MULTI_FVM_STRUCT), INTENT(INOUT) :: MULTI_FVM
      INTEGER, INTENT(IN) :: ITASK
      INTEGER, INTENT(IN) :: IXS(NIXS, *), IXQ(NIXQ, *), IXTG(NIXTG, *)
      my_real, INTENT(IN) :: xgrid(3, *), wgrid(3, *)
      INTEGER, INTENT(IN) :: IPM(NPROPMI, *)
      my_real, INTENT(IN) :: pm(npropm, *), func_value(*)
      TYPE(t_ebcs_tab), TARGET, INTENT(IN) :: EBCS_TAB
      INTEGER, INTENT(IN) :: NPF(SNPC)
      my_real, INTENT(IN) :: tf(stf), timestep
      TYPE(OUTPUT_), INTENT(INOUT) :: OUTPUT !< output structure
      LOGICAL, INTENT(IN) :: PRED !< predictor step if true, corrector step if false
C-----------------------------------------------
C     L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER :: I
      INTEGER :: TYPE, NELEM
      class(t_ebcs), POINTER :: ebcs
C-----------------------------------------------
C     B e g i n n i n g   o f   s u b r o u t i n e
C-----------------------------------------------
      DO i = 1, ebcs_tab%nebcs_fvm
         ebcs => ebcs_tab%tab(i)%poly
         TYPE = ebcs%type
         nelem = ebcs%nb_elem
         SELECT TYPE (twf => ebcs)
            TYPE IS (t_ebcs_inlet)
C     Inlet, outlet boundary condition
              CALL multi_inlet_ebcs(itask, i, multi_fvm, nelem,
     .           twf%ielem, twf%iface, twf%fvm_inlet_data,
     .           ixs, ixq, ixtg, xgrid, wgrid, ipm, pm, func_value,
     .           twf%surf_id,npf,tf,fsavsurf, timestep, matparam, output, pred)
            TYPE IS(t_ebcs_fluxout)
              CALL multi_inlet_ebcs(itask, i, multi_fvm, nelem,
     .           twf%ielem, twf%iface, twf%fvm_inlet_data,
     .           ixs, ixq, ixtg, xgrid, wgrid, ipm, pm, func_value,
     .           twf%surf_id,npf,tf,fsavsurf,timestep, matparam, output, pred)
            TYPE IS(t_ebcs_nrf)
              CALL multi_nrf_ebcs(itask, i, multi_fvm, nelem,
     .           twf%ielem, twf%iface, twf%fvm_inlet_data,
     .           ixs, ixq, ixtg, xgrid, wgrid, ipm, pm, func_value,
     .           twf,npf,tf,fsavsurf,timestep, matparam, output, pred)
            TYPE IS(t_ebcs_propellant)
!             work in progress
              write(6,*) 'MULTI_EBCS: Propellant EBCS not yet implemented'
              write(6,*) "NUMELS=",numels,"NUMELQ=",numelq,"NUMELTG=",numeltg
              write(6,*) "NUMNOD=",numnod,"NCYCLE=",ncycle,"NUMMAT=",nummat
              CALL arret(2)
!             CALL MULTI_PROPELLANT_EBCS(ITASK, I, MULTI_FVM, NELEM,
!    .           twf%ielem, twf%iface, twf%fvm_inlet_data,
!    .           IXS, IXQ, IXTG, XGRID, WGRID, IPM, PM, FUNC_VALUE,
!    .           twf,NPF,TF,FSAVSURF,TIMESTEP,NIXS, NIXQ, NIXTG, NPROPMI, NPROPM, NSURF, STF, SNPC,
!    .           NUMELS, NUMELQ, NUMELTG, NUMNOD, NCYCLE, NUMMAT ,MATPARAM)
            CLASS DEFAULT
         END SELECT
      ENDDO
C-----------------------------------------------
C     E n d   o f   s u b r o u t i n e 
C-----------------------------------------------