w_th_surf_loadp.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	w_th_surf_loadp (iloadp, nodlocal, proc, len_ia, cep, th_surf, numnod, sizloadp, nloadp, len_cep, off)
	Writing TH/SURF tabs in restart file.

Function/Subroutine Documentation

w_th_surf_loadp()

subroutine w_th_surf_loadp	(	integer, dimension(sizloadp,nloadp), intent(in)	iloadp,
		integer, dimension(numnod), intent(in)	nodlocal,
		integer, intent(in)	proc,
		integer, intent(inout)	len_ia,
		integer, dimension(len_cep), intent(in)	cep,
		type (th_surf_), intent(inout)	th_surf,
		integer, intent(in)	numnod,
		integer, intent(in)	sizloadp,
		integer, intent(in)	nloadp,
		integer, intent(in)	len_cep,
		integer, intent(in)	off )

Writing TH/SURF tabs in restart file.

Load pressures option

Parameters

[in]	len_cep	Parameters for size tables ( number of nodes,number of parameters in ILOADP tab, number of load pressure, size of CEP)
[in]	off	Index in table CEP
[in]	proc	Processor
[in,out]	len_ia	Length of integer tables written in restart
[in]	iloadp	Integer tabs for load pressures (/PFLUID, /PBLAST, /LOAD/PRESSURE )
[in]	nodlocal	table of local nodes
[in]	cep	Proc to which belong the element
[in,out]	th_surf	Type for /TH/SURF and load pressures output tabs

Definition at line 35 of file w_th_surf_loadp.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE my_alloc_mod
      USE th_surf_mod , ONLY : th_surf_
      USE message_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-----------------------------------------------
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER , INTENT(IN) :: NUMNOD   ,SIZLOADP ,NLOADP   ,LEN_CEP  !< Parameters for size tables ( number of nodes,number of parameters in ILOADP tab, number of load pressure, size of CEP)
      INTEGER , INTENT(IN) :: OFF               !< Index in table CEP
      INTEGER , INTENT(IN) :: PROC              !< Processor
      INTEGER , INTENT(INOUT) :: LEN_IA            !< Length of integer tables written in restart
      INTEGER , INTENT(IN) :: ILOADP(SIZLOADP,NLOADP)     !< Integer tabs for load pressures (/PFLUID, /PBLAST, /LOAD/PRESSURE )
      INTEGER , INTENT(IN) :: NODLOCAL(NUMNOD)   !< table of local nodes
      INTEGER , INTENT(IN) :: CEP(LEN_CEP)          !< Proc to which belong the element
      TYPE (TH_SURF_) , INTENT(INOUT) :: TH_SURF !< Type for /TH/SURF and load pressures output tabs
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, J, NSEGP, NSEGPL, NL, NS, K, ADCEP
      INTEGER, DIMENSION(:),ALLOCATABLE :: KSEGS_L, SEGS_L
C-----------------------------------------------
C   S o u r c e   L i n e s
C-----------------------------------------------
 
C-----------------------------------------------
C         Writing TH/SURF tabs in restart file
C          /PFLUID , /PBLAST, /LOADP_HYD option
C-----------------------------------------------
 
      CALL my_alloc(ksegs_l,th_surf%S_LOADP_KSEGS)
      CALL my_alloc(segs_l,th_surf%S_LOADP_SEGS)
      ksegs_l = 0
      segs_l = 0
      nsegp  = 0      
      nsegpl = 0
      adcep  = 0
      DO i = 1, nloadp
        nl = iloadp(1,i)
        DO j = 1, nl/4
          nsegp = nsegp + 1
          IF(cep(j+off+adcep)==proc) THEN
            nsegpl = nsegpl + 1
            ns = th_surf%LOADP_KSEGS (nsegp+1) - th_surf%LOADP_KSEGS (nsegp)
            ksegs_l(nsegpl+1) =ksegs_l(nsegpl)+ ns
            DO k=1,ns
               segs_l(ksegs_l(nsegpl)+k)= th_surf%LOADP_SEGS(th_surf%LOADP_KSEGS (nsegp)+k)
            ENDDO
          ENDIF
        ENDDO
        adcep = adcep + nl/4
      ENDDO
 
      CALL write_i_c(nsegpl+1,1)
      CALL write_i_c(ksegs_l,nsegpl+1)
      CALL write_i_c(ksegs_l(nsegpl+1),1)
      CALL write_i_c(segs_l,ksegs_l(nsegpl+1))
      len_ia = len_ia + nsegpl+3+ ksegs_l(nsegpl+1)
 
      DEALLOCATE(ksegs_l,segs_l)
 
      RETURN