write_th.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	write_th (n, nspmd, nn, nvar, ittyp, eltype_struct, wa_eltype_p0)

Function/Subroutine Documentation

write_th()

subroutine write_th	(	integer, intent(in)	n,
		integer, intent(in)	nspmd,
		integer, intent(in)	nn,
		integer, intent(in)	nvar,
		integer, intent(in)	ittyp,
		type(th_proc_type), intent(in)	eltype_struct,
		type(th_wa_real), intent(in)	wa_eltype_p0 )

Definition at line 32 of file write_th.F.

        USE th_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
        INTEGER, INTENT(in) :: N,NSPMD,NN,NVAR,ITTYP
        TYPE(TH_WA_REAL), INTENT(in) ::WA_ELTYPE_P0
        TYPE(TH_PROC_TYPE), INTENT(in) :: ELTYPE_STRUCT
!       -*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-**-*-*-*-*-*-*-*-*-*-*-*-*
!       N : integer ; current TH group 
!       NSPMD : integer ; number of MPI domains
!       NN : integer, number of element group
!       NVAR : integer ; number of value per element group
!       ITTYP : integer ; type of TH group
!       -*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-*-**-*-*-*-*-*-*-*-*-*-*-*-*
C-----------------------------------------------
C   L o c a l  V a r i a b l e s
C-----------------------------------------------
        LOGICAL :: BOOL
        INTEGER :: I,J,K,II,IJK,LOCAL_SIZE
        INTEGER :: NEXT,CURRENT,SIZE_N
        INTEGER :: CURRENT_J,REST
        my_real, DIMENSION(:), ALLOCATABLE :: wa_local
!       ----------------------------------------
!$COMMENT
!       WRITE_TH_COQ description
!       write all the values for shell element
!       and for a given group N by the PROC0
!       
!       WRITE_TH_COQ organization :
!        * loop over the NSPMD processor and:
!           - check if a processor must write its values (SIZE_N>0)
!           - compute the number of group for a given processor (REST)
!           - compute the position in the TH file (CURRENT_J)
!           - initialization of the written values WA_LOCAL with the position CURRENT_J + J
!        * write the WA_LOCAL values 
!$ENDCOMMENT
 
!       allocation + initialization of local array
        ALLOCATE( wa_local(nn*nvar) )
        wa_local(1:nn*nvar) = zero
 
 
 
        DO i=1,nspmd
            local_size = eltype_struct%TH_PROC(i)%TH_ELM_SIZE
            bool=.true.
            DO k=1,local_size
                IF(bool.EQV..true.) THEN
                    IF(eltype_struct%TH_PROC(i)%TH_ELM(k,2)==n) THEN
                        bool=.false.
                        ijk=k
                    ENDIF
                ENDIF
            ENDDO
            IF(bool.EQV..false.) THEN
                current = eltype_struct%TH_PROC(i)%TH_ELM(ijk,1) !   index of the current proc and N
                next = eltype_struct%TH_PROC(i)%TH_ELM(ijk+1,1)  !   index of the next proc and n
                size_n = next-current   !   nbr of WA element of the current proc and for the current N
!                IF( SIZE_N>0 ) THEN     !   if SIZE_N>0, must write some data
                rest = size_n / (nvar+1)    ! nbr of th group for the current proc
                ii = 0
                DO k=1,rest
                    current_j = nint(wa_eltype_p0%WA_REAL( current+k*(nvar+1) ) ) ! the position in WA_LOCAL for the current TH 
                    DO j=1,nvar   
                        ii = ii + 1   
                        wa_local(current_j+j) = wa_eltype_p0%WA_REAL(current+ii)                
                    ENDDO
                    ii = ii + 1  ! skip the NVAR+1 value = the position in WA_LOCAL
                ENDDO
!                ENDIF
            ENDIF
        ENDDO
 
!       write the data
        CALL wrtdes(wa_local,wa_local,nn*nvar,ittyp,1)
!       deallocation
        DEALLOCATE( wa_local )
 
!       ----------------------------------------
        RETURN