w_secbuf.F File Reference

#include "implicit_f.inc"
#include "com04_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	w_secbuf (secbuf, lsecbuf_l, nstrf, proc, len_am)

Function/Subroutine Documentation

w_secbuf()

subroutine w_secbuf	(		secbuf,
		integer	lsecbuf_l,
		integer, dimension(*)	nstrf,
		integer	proc,
		integer	len_am )

Definition at line 31 of file w_secbuf.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------  
 
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-----------------------------------------------
#include      "com04_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
       INTEGER  PROC, LSECBUF_L, LEN_AM,
     .          NSTRF(*)
      my_real
     .          secbuf(*)
C-----------------------------------------------
C   F u n c t i o n
C-----------------------------------------------
      INTEGER  NLOCAL
      EXTERNAL nlocal    
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER NNOD_S, NNOD_S_L, TYP_S, NSINT_S,
     .        N, IP, J, K0, K2, KR_L, KR, K
      my_real
     .        secbuf_l(lsecbuf_l)
C
      kr_l = 0
      kr = 0
C header commun
      DO j = 1, 10
        secbuf_l(kr_l+j) = secbuf(kr+j)
      END DO
      kr_l = kr_l + 10
      kr = kr + 10
C
      k0 = nstrf(25)
      DO n = 1, nsect
        typ_s  = nstrf(k0)
        nnod_s = nstrf(k0+6)
        nsint_s = nstrf(k0+14)
        k2=k0+30+nsint_s
C noeuds locaux
        nnod_s_l = 0
        DO j = 1, nnod_s
          k = nstrf(k2 + j -1)
          IF(nlocal(k,proc+1)==1)THEN
            nnod_s_l = nnod_s_l + 1
          END IF
        END DO
        DO j = 1, 10
          secbuf_l(kr_l+j) = secbuf(kr+j)
        END DO
        kr_l = kr_l + 10
        kr = kr + 10
        IF(typ_s>=100) THEN
          DO j = 1, 12*nnod_s_l
            secbuf_l(kr_l+j) = zero
          END DO
          kr_l = kr_l + 12*nnod_s_l
          kr = kr + 12*nnod_s
        END IF
        IF(typ_s>=101) THEN
          DO j = 1, 12*nnod_s_l
            secbuf_l(kr_l+j) = zero
          END DO
          kr_l = kr_l + 12*nnod_s_l
          kr = kr + 12*nnod_s
        END IF
        IF(typ_s>=102) THEN
          DO j = 1, 6*nnod_s_l
            secbuf_l(kr_l+j) = zero
          END DO
          kr_l = kr_l + 6*nnod_s_l
          kr = kr + 6*nnod_s
        END IF
        k0  = nstrf(k0+24)
      ENDDO
C
      CALL write_db(secbuf_l,lsecbuf_l)
      len_am = len_am + lsecbuf_l
      if(lsecbuf_l/=kr_l)print*,'error secbuf',lsecbuf_l,kr_l,nnod_s_l
C
      RETURN