w_th_surf_pload.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	w_th_surf_pload (ibcl, nodlocal, proc, len_ia, cep, th_surf, numnod, nibcld, nconld, len_cep, off)
	Writing TH/SURF tabs in restart file.

Function/Subroutine Documentation

w_th_surf_pload()

subroutine w_th_surf_pload	(	integer, dimension(nibcld,nconld), intent(in)	ibcl,
		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)	nibcld,
		integer, intent(in)	nconld,
		integer, intent(in)	len_cep,
		integer, intent(in)	off )

Writing TH/SURF tabs in restart file.

/PLOAD option

Parameters

[in]	len_cep	Parameters for size tables ( number of nodes,number of parameters in IBCL tab, number of pload, size of CEP)
[in]	off	Index in table CEP
[in]	proc	Processor
[in,out]	len_ia	Length of integer tables written in restart
[in]	ibcl	list of segments where pressure /PLOAD is applied
[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_pload.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   ,NIBCLD ,NCONLD   ,LEN_CEP  !< Parameters for size tables ( number of nodes,number of parameters in IBCL tab, number of pload, 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) :: IBCL(NIBCLD,NCONLD)   !< list of segments where pressure /PLOAD is applied 
      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, N4
      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          /PLOAD option
C-----------------------------------------------
 
      CALL my_alloc(ksegs_l,th_surf%S_PLOAD_KSEGS)
      CALL my_alloc(segs_l,th_surf%S_PLOAD_SEGS)
      ksegs_l = 0
      segs_l = 0
      nsegp  = 0      
      nsegpl = 0
      DO i = 1, nconld
        n4 = ibcl(4,i)
        IF(n4/=-1) THEN
          nsegp = nsegp + 1
          IF(cep(i+off)==proc) THEN
            nsegpl = nsegpl + 1
            ns = th_surf%PLOAD_KSEGS (nsegp+1) - th_surf%PLOAD_KSEGS (nsegp)
            ksegs_l(nsegpl+1) =ksegs_l(nsegpl)+ ns
            DO j=1,ns
               segs_l(ksegs_l(nsegpl)+j)= th_surf%PLOAD_SEGS(th_surf%PLOAD_KSEGS (nsegp)+j)
            ENDDO
          ENDIF
        ENDIF
      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