write_pcyl.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	write_pcyl (load_cyl, nload_cyl, nodlocal)

Function/Subroutine Documentation

write_pcyl()

subroutine write_pcyl	(	type (press_cyl_), dimension(nload_cyl), intent(in)	load_cyl,
		integer	nload_cyl,
		integer, dimension(*), intent(in)	nodlocal )

Definition at line 31 of file write_pcyl.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE message_mod
      USE pload_cyl_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 :: NLOAD_CYL
      TYPE (PRESS_CYL_) ,DIMENSION(NLOAD_CYL) ,INTENT(IN) :: LOAD_CYL
      INTEGER , DIMENSION(*) ,INTENT(IN) :: NODLOCAL
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER :: I,J,N1,N2,N3,N4,NSEG,LINT,IAD,OFF
      INTEGER, DIMENSION(:),ALLOCATABLE :: ITMP
      my_real, DIMENSION(3) :: rtmp
C-----------------------------------------------
c         -> PCYL_ID
c         -> NSEG
c         -> SEGNOD(NSEG,4)  (N1,N2,NB3,N4) by segment
c         -> FRAME_ID
c         -> SENS_ID
c         -> TABLE_ID
c         -> XSCALE_R
c         -> XSCALE_T
c         -> YSCALE_P
c=======================================================================
 
      ! --------------
      ! loop over the /load/pcyl
      DO i=1,nload_cyl        
        nseg = load_cyl(i)%NSEG ! get the number of segment 
        lint = 5 + nseg*4
        ALLOCATE(itmp(lint))
        rtmp(1) = load_cyl(i)%XSCALE_R
        rtmp(2) = load_cyl(i)%XSCALE_T
        rtmp(3) = load_cyl(i)%YSCALE
c
        itmp(1) = nseg
        itmp(2) = load_cyl(i)%ID
        itmp(3) = load_cyl(i)%ISENS
        itmp(4) = load_cyl(i)%ITABLE
        itmp(5) = load_cyl(i)%IFRAME
        iad     = 5
        ! --------------
        ! loop over the segment to convert global node id into local node id
        DO j=1,load_cyl(i)%NSEG
            n1 = load_cyl(i)%SEGNOD(j,1) ! get the global node id N1
            n2 = load_cyl(i)%SEGNOD(j,2) ! get the global node id N2
            n3 = load_cyl(i)%SEGNOD(j,3) ! get the global node id N3
            n4 = load_cyl(i)%SEGNOD(j,4) ! get the global node id N4
            itmp(iad+1) = nodlocal(n1) ! global --> local for N1
            itmp(iad+2) = nodlocal(n2) ! global --> local for N2
            itmp(iad+3) = nodlocal(n3) ! global --> local for N3
            ! global --> local for N4 : for triangle, N4 = 0, for shell N4/=0
            IF (n4 == 0) THEN
                itmp(iad+4) = 0
            ELSE
                itmp(iad+4) = nodlocal(n4)
            ENDIF
            iad = iad + 4
        ENDDO
        ! --------------
 
        ! --------------        
        ! write the data related to the load/pcyl
        CALL write_i_c(itmp,iad)
        CALL write_db(rtmp,3)
        ! /LOAD/PCYL option : adress for parith/on 
        CALL write_i_c(load_cyl(i)%SEGMENT_ADRESS,4*load_cyl(i)%S_SEGMENT_ADRESS)
        ! --------------
                
        DEALLOCATE(itmp)
      ENDDO
      ! --------------
 
 
      RETURN