i25cor3_e2s.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	i25cor3_e2s (jlt, ledge, irect, x, cand_s, cand_m, ex, ey, ez, xxs1, xxs2, xys1, xys2, xzs1, xzs2, xxm1, xxm2, xym1, xym2, xzm1, xzm2, n1, n2, m1, m2, nrts, gape, gapve, fx, fy, fz, iedge, admsr, lbound, edg_bisector, vtx_bisector, itab)

Function/Subroutine Documentation

i25cor3_e2s()

subroutine i25cor3_e2s	(	integer	jlt,
		integer, dimension(nledge,*)	ledge,
		integer, dimension(4,*)	irect,
			x,
		integer, dimension(*)	cand_s,
		integer, dimension(*)	cand_m,
			ex,
			ey,
			ez,
			xxs1,
			xxs2,
			xys1,
			xys2,
			xzs1,
			xzs2,
			xxm1,
			xxm2,
			xym1,
			xym2,
			xzm1,
			xzm2,
		integer, dimension(mvsiz)	n1,
		integer, dimension(mvsiz)	n2,
		integer, dimension(4,mvsiz)	m1,
		integer, dimension(4,mvsiz)	m2,
		integer	nrts,
			gape,
			gapve,
			fx,
			fy,
			fz,
		integer	iedge,
		integer, dimension(4,*)	admsr,
		integer, dimension(*)	lbound,
		real*4, dimension(3,4,*)	edg_bisector,
		real*4, dimension(3,2,*)	vtx_bisector,
		integer, dimension(*)	itab )

Definition at line 30 of file i25cor3_e2s.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE tri7box
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   G l o b a l   P a r a m e t e r s
C-----------------------------------------------
#include      "mvsiz_p.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER LEDGE(NLEDGE,*), IRECT(4,*), CAND_M(*), CAND_S(*), ADMSR(4,*),
     .        LBOUND(*), JLT, NRTS, IEDGE, ITAB(*),
     .        N1(MVSIZ), N2(MVSIZ), 
     .        M1(4,MVSIZ), M2(4,MVSIZ)
C     REAL
      my_real
     .        x(3,*), 
     .        xxs1(mvsiz), xxs2(mvsiz), xys1(mvsiz), xys2(mvsiz),
     .        xzs1(mvsiz), xzs2(mvsiz), xxm1(4,mvsiz), xxm2(4,mvsiz),
     .        xym1(4,mvsiz), xym2(4,mvsiz), xzm1(4,mvsiz), xzm2(4,mvsiz),
     .        gape(*),gapve(mvsiz),
     .        ex(4,mvsiz), ey(4,mvsiz), ez(4,mvsiz), fx(mvsiz), fy(mvsiz), fz(mvsiz)
      real*4 edg_bisector(3,4,*), vtx_bisector(3,2,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I ,NN, J, JRM, K, KRM, I1, J1, I2, J2, EJ, ES,
     .         IE, JE, SOL_EDGE, SH_EDGE
      INTEGER IAM(MVSIZ),JAM(4,MVSIZ),IAS(MVSIZ),JAS(MVSIZ)
      my_real
     .   aaa, dx, dy, dz, dd, nni, ni2, invcos, gape_m(mvsiz), gape_s(mvsiz)
C-----------------------------------------------
      DO i=1,jlt
        IF(cand_s(i)<=nrts) THEN
 
          es    =cand_s(i)
          ias(i)=ledge(1,es)
          jas(i)=ledge(2,es)
          n1(i)=ledge(5,es)
          n2(i)=ledge(6,es)
 
          xxs1(i) = x(1,n1(i))
          xys1(i) = x(2,n1(i))
          xzs1(i) = x(3,n1(i))
          xxs2(i) = x(1,n2(i))
          xys2(i) = x(2,n2(i))
          xzs2(i) = x(3,n2(i))
 
          iam(i)=cand_m(i)
          DO ej=1,4
            jam(ej,i)=ej
            m1(ej,i)=irect(ej,iam(i))
            m2(ej,i)=irect(mod(ej,4)+1,iam(i))
 
            xxm1(ej,i) = x(1,m1(ej,i))
            xym1(ej,i) = x(2,m1(ej,i))
            xzm1(ej,i) = x(3,m1(ej,i))
            xxm2(ej,i) = x(1,m2(ej,i))
            xym2(ej,i) = x(2,m2(ej,i))
            xzm2(ej,i) = x(3,m2(ej,i))
          END DO
        END IF
      END DO
 
      DO i=1,jlt
        gape_m(i)=zero ! Solids 
                       ! If edge is shared by solid and shell : edge is considered as a shell edge
          IF(cand_s(i)<=nrts) THEN
          gape_s(i)=gape(cand_s(i))
          END IF
        gapve(i)=zero
      END DO
 
C
      sol_edge=iedge/10 ! solids
      sh_edge =iedge-10*sol_edge ! shells
 
      DO i=1,jlt
        DO ej=1,4
          ex(ej,i)=edg_bisector(1,ej,iam(i))
          ey(ej,i)=edg_bisector(2,ej,iam(i))
          ez(ej,i)=edg_bisector(3,ej,iam(i))
        END DO
      END DO
C
      DO i=1,jlt
        fx(i) = edg_bisector(1,jas(i),ias(i))
        fy(i) = edg_bisector(2,jas(i),ias(i))
        fz(i) = edg_bisector(3,jas(i),ias(i))
      END DO
C
      RETURN