dmasani6.F File Reference

#include "implicit_f.inc"
#include "vect01_c.inc"
#include "mvsiz_p.inc"
#include "com01_c.inc"
#include "com04_c.inc"
#include "sphcom.inc"
#include "scr17_c.inc"
#include "param_c.inc"
#include "task_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	dmasanis (elbuf_tab, x, d, geo, iparg, ixs, mas, pm, el2fa, nbf, ipart, ipartsp, isph3d)

Function/Subroutine Documentation

dmasanis()

subroutine dmasanis	(	type (elbuf_struct_), dimension(ngroup), target	elbuf_tab,
			x,
			d,
			geo,
		integer, dimension(nparg,*)	iparg,
		integer, dimension(nixs,*)	ixs,
			mas,
			pm,
		integer, dimension(*)	el2fa,
		integer	nbf,
		integer, dimension(lipart1,*)	ipart,
		integer, dimension(*)	ipartsp,
		integer	isph3d )

Definition at line 34 of file dmasani6.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE initbuf_mod
      USE elbufdef_mod
      use element_mod , only : nixs
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      "vect01_c.inc"
#include      "mvsiz_p.inc"
#include      "com01_c.inc"
#include      "com04_c.inc"
#include      "sphcom.inc"
#include      "scr17_c.inc"
#include      "param_c.inc"
#include      "task_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
C     REAL
      my_real
     .   mas(*) ,pm(npropm,*),geo(npropg,*),x(3,*),d(3,*)
      INTEGER IPARG(NPARG,*),IXS(NIXS,*),EL2FA(*),NBF,IPART(LIPART1,*),
     .        IPARTSP(*),ISPH3D
      TYPE (ELBUF_STRUCT_), DIMENSION(NGROUP), TARGET :: ELBUF_TAB
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
C     REAL
      my_real
     .   evar(mvsiz),
     .   off, p, vonm2, vonm, s1, s2, s12, s3, VALUE,thk0,a0,al0,
     .   rho0,xx1,xx2,xx3,yy1,yy2,yy3,zz1,zz2,zz3
      INTEGER I, NG, NEL, 
     .        IADD, N, J, MLW, 
     .        ISTRAIN,NN, K1, K2,JTURB,MT,IMID, IALEL,IPID,
     .        N1,N2,N3,N4,
     .        NN1,NN2,NN3,NN4,NN5,NN6,NN7,NN8,NN9,NN10,
     .        OFFSET,NEL_OLD,ITY_OLD,NFT_FA,N_FA,
     .        INOD, ISOLNOD, IPRT
      TYPE(G_BUFEL_)  ,POINTER :: GBUF     
C=======================================================================
      nn1 = 1
      nn2 = 1
      nn3 = nn2 + numels 
      nn4 = nn3 + isph3d*(numsph+maxpjet)
C-----------------------------------------------
C       face ext
C-----------------------------------------------
      nel_old = 0
      ity_old = 0
      DO 490 ng=1,ngroup
        CALL initbuf(iparg    ,ng      ,                    
     2          mlw     ,nel     ,nft     ,iad     ,ity     ,  
     3          npt     ,jale    ,ismstr  ,jeul    ,jtur    ,  
     4          jthe    ,jlag    ,jmult   ,jhbe    ,jivf    ,  
     5          nvaux   ,jpor    ,jcvt    ,jclose  ,jplasol ,  
     6          irep    ,iint    ,igtyp   ,israt   ,isrot   ,  
     7          icsen   ,isorth  ,isorthg ,ifailure,jsms    )
        isolnod = iparg(28,ng)
        IF (ispmd == 0) THEN
          IF (ity/=ity_old) THEN
            nel_old = 0
            ity_old= ity
          ENDIF
          nft_fa = nel_old
          nel_old = nel_old + nel
        ENDIF
        DO offset = 0,nel-1,nvsiz
          nft   =iparg(3,ng) + offset
          iad   =iparg(4,ng)
          lft=1
          llt=min(nvsiz,nel-offset)
          IF (ispmd == 0) THEN
            nft_fa = nel_old - nel + offset
          ELSE
            nft_fa = nft
          ENDIF
C-----------------------------------------------
C       SOLIDE
C-----------------------------------------------
          IF(ity == 1)THEN
           gbuf => elbuf_tab(ng)%GBUF
           ialel=iparg(7,ng)+iparg(11,ng)
           DO 130 i=lft,llt
             n = i + nft
             n_fa = i + nft_fa
             IF(el2fa(nn2+n_fa)/=0)THEN
               IF (mlw == 0 .or. mlw == 13 .or. igtyp == 0) THEN
                 VALUE = zero
               ELSEIF(ialel == 0)THEN
                 mt=ixs(1,n)
                 VALUE=pm(89,mt)*gbuf%VOL(i)
               ELSE
                 off = min(gbuf%OFF(i),one)
                 VALUE=gbuf%RHO(i)*gbuf%VOL(i)*off
               ENDIF
               IF( nfilsol /= 0 .AND. gbuf%G_FILL /= 0 ) 
     .           VALUE = VALUE * gbuf%FILL(i)
               IF (isolnod == 16) THEN
                 VALUE = fourth*VALUE
                 mas(el2fa(nn2+n_fa))   = VALUE
                 mas(el2fa(nn2+n_fa)+1) = VALUE
                 mas(el2fa(nn2+n_fa)+2) = VALUE
                 mas(el2fa(nn2+n_fa)+3) = VALUE
               ELSE
                 mas(el2fa(nn2+n_fa)) = VALUE
               ENDIF
             ENDIF
 130       CONTINUE
C
        ELSEIF(isph3d == 1 .AND. ity == 51)THEN
C-----------------------------------------------
C          TETRAS SPH.
C-----------------------------------------------
           gbuf => elbuf_tab(ng)%GBUF
           ialel=iparg(7,ng)+iparg(11,ng)
           DO 140 i=lft,llt
             n = i + nft
             n_fa = i + nft_fa
             IF(el2fa(nn3+n_fa)/=0)THEN
               IF(ialel == 0)THEN
                 iprt=ipartsp(n)
                 mt  =ipart(1,iprt)
                 VALUE=pm(89,mt)*gbuf%VOL(i)
               ELSE
                 off = min(gbuf%OFF(i),one)
                 VALUE=gbuf%RHO(i)*gbuf%VOL(i)*off
               ENDIF
               mas(el2fa(nn3+n_fa)) = VALUE
             ENDIF
 140       CONTINUE
C
        ELSE
        ENDIF
C-----------------------------------------------
C       end of loop over offsets
C-----------------------------------------------
       END DO
 490  CONTINUE
C-----------------------------------------------
C
      RETURN