cdkfint3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	cdkfint3 (jft, jlt, vol, thk0, for, mom, px2, py2, px3, py3, bz1, bz2, bz3, brx1, brx2, brx3, bry1, bry2, bry3, f11, f12, f13, f21, f22, f23, f32, f33, m11, m12, m13, m21, m22, m23, nel)

Function/Subroutine Documentation

cdkfint3()

subroutine cdkfint3	(	integer	jft,
		integer	jlt,
			vol,
			thk0,
			for,
			mom,
			px2,
			py2,
			px3,
			py3,
			bz1,
			bz2,
			bz3,
			brx1,
			brx2,
			brx3,
			bry1,
			bry2,
			bry3,
			f11,
			f12,
			f13,
			f21,
			f22,
			f23,
			f32,
			f33,
			m11,
			m12,
			m13,
			m21,
			m22,
			m23,
		integer	nel )

Definition at line 28 of file cdkfint3.F.

C--------------------------------------------------------------------------------------------------
C-----------------------------------------------
C   I M P L I C I T   T Y P E S
C-----------------------------------------------
#include      "implicit_f.inc"
#include      "mvsiz_p.inc"
C-----------------------------------------------
C   D U M M Y   A R G U M E N T S
C-----------------------------------------------
      INTEGER JFT,JLT,NG,NEL
      my_real 
     .   for(nel,5),mom(nel,3), 
     .   px2(*), py2(*), px3(*), py3(*), 
     .   brx1(mvsiz,3),brx2(mvsiz,3),brx3(mvsiz,3),bry1(mvsiz,3),bry2(mvsiz,3),bry3(mvsiz,3),
     .   bz1(mvsiz,2),bz2(mvsiz,2),bz3(mvsiz,2),
     .   f11(*), f12(*), f13(*),
     .   f21(*), f22(*), f23(*), f32(*), f33(*),
     .   m11(*), m12(*), m13(*),
     .   m21(*), m22(*), m23(*),
     .   thk0(*),vol(*)
C-----------------------------------------------
C   L O C A L   V A R I A B L E S
C-----------------------------------------------
      INTEGER I
      my_real 
     .    c2,fx2,fx3,fy2,fy3
C--------------------------------------------------------------------------------------------------
      DO i=jft,jlt
          c2=vol(i)
        fx2 =c2*(px2(i)*for(i,1)+py2(i)*for(i,3))
        fy2 =c2*(py2(i)*for(i,2)+px2(i)*for(i,3))
        fx3 =c2*(px3(i)*for(i,1)+py3(i)*for(i,3))
        fy3 =c2*(py3(i)*for(i,2)+px3(i)*for(i,3))
        f12(i)=f12(i)+fx2
        f22(i)=f22(i)+fy2
        f13(i)=f13(i)+fx3
        f23(i)=f23(i)+fy3
C        F11(I)=F11(I)-FX2-FX3
C        F21(I)=F21(I)-FY2-FY3
      ENDDO     
      DO i=jft,jlt
          c2=vol(i)*thk0(i)
        f32(i)=f32(i)+c2*(bz1(i,1)*mom(i,1)+bz2(i,1)*mom(i,2)+
     .                    bz3(i,1)*mom(i,3))
        f33(i)=f33(i)+c2*(bz1(i,2)*mom(i,1)+bz2(i,2)*mom(i,2)+
     .                    bz3(i,2)*mom(i,3))
        m11(i)=m11(i)+c2*(brx1(i,1)*mom(i,1)+brx2(i,1)*mom(i,2)+
     .                    brx3(i,1)*mom(i,3))
        m21(i)=m21(i)+c2*(bry1(i,1)*mom(i,1)+bry2(i,1)*mom(i,2)+
     .                    bry3(i,1)*mom(i,3))
        m12(i)=m12(i)+c2*(brx1(i,2)*mom(i,1)+brx2(i,2)*mom(i,2)+
     .                    brx3(i,2)*mom(i,3))
        m22(i)=m22(i)+c2*(bry1(i,2)*mom(i,1)+bry2(i,2)*mom(i,2)+
     .                    bry3(i,2)*mom(i,3))
        m13(i)=m13(i)+c2*(brx1(i,3)*mom(i,1)+brx2(i,3)*mom(i,2)+
     .                    brx3(i,3)*mom(i,3))
        m23(i)=m23(i)+c2*(bry1(i,3)*mom(i,1)+bry2(i,3)*mom(i,2)+
     .                    bry3(i,3)*mom(i,3))
      ENDDO     
C    
      RETURN