i2curvv.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	i2curvv (nsn, nmn, ms, v, a, ar, vr, x, irect, nsv, msr, irtl, crst, weight)

Function/Subroutine Documentation

i2curvv()

subroutine i2curvv	(	integer	nsn,
		integer	nmn,
			ms,
			v,
			a,
			ar,
			vr,
			x,
		integer, dimension(4,*)	irect,
		integer, dimension(*)	nsv,
		integer, dimension(*)	msr,
		integer, dimension(*)	irtl,
			crst,
		integer, dimension(*)	weight )

Definition at line 32 of file i2curvv.F.

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 NSN,NMN,
     .   IRECT(4,*), MSR(*), NSV(*), IRTL(*),WEIGHT(*)
C     REAL
      my_real
     .   a(3,*),v(3,*),ar(3,*),vr(3,*),x(3,*),ms(*),crst(2,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER NIR,I,J,L,II,JJ,IS
      INTEGER INOD(4)
C     REAL
      my_real
     .   ls1,ls2,lt1,lt2,ls,lt,s,t
      my_real
     .   vs(6),as(6),vm(6,4),am(6,4),rot(9),roti(9),
     .   h(4),hh(4),hrs(4),hrt(4),hps(4),hpt(4),hprs(4),hprt(4),
     .   hxs(4),hxt(4)
C=======================================================================
      DO ii=1,nsn
        is=nsv(ii)
        IF (is > 0) THEN
          l=irtl(ii)
          nir = 4
          DO jj=1,nir
            inod(jj) = irect(jj,l)
          ENDDO
C---
          CALL i2curv_rep(inod    , x      ,v       ,ls1     ,ls2     ,
     .                    lt1     ,lt2     ,rot(1)  ,rot(4)  ,rot(7)  ,
     .                    rot(2)  ,rot(5)  ,rot(8)  ,rot(3)  ,rot(6)  ,
     .                    rot(9)  )                                    
          CALL inv3(rot,roti)
C---
          s=crst(1,ii)
          t=crst(2,ii)
          CALL i2_fform(
     .         nir,s,t,h,hh,hrs,hrt,hps,hpt,hprs,hprt,
     .         hxs,hxt,ls1,ls2,lt1,lt2,ls,lt)
C---
          DO jj=1,nir
            j=inod(jj)
            vm(1,jj) = rot(1)*v(1,j) + rot(4)*v(2,j) + rot(7)*v(3,j)
            vm(2,jj) = rot(2)*v(1,j) + rot(5)*v(2,j) + rot(8)*v(3,j)
            vm(3,jj) = rot(3)*v(1,j) + rot(6)*v(2,j) + rot(9)*v(3,j)
            am(1,jj) = rot(1)*a(1,j) + rot(4)*a(2,j) + rot(7)*a(3,j)
            am(2,jj) = rot(2)*a(1,j) + rot(5)*a(2,j) + rot(8)*a(3,j)
            am(3,jj) = rot(3)*a(1,j) + rot(6)*a(2,j) + rot(9)*a(3,j)
            vm(4,jj) = rot(1)*vr(1,j)+rot(4)*vr(2,j)+rot(7)*vr(3,j)  
            vm(5,jj) = rot(2)*vr(1,j)+rot(5)*vr(2,j)+rot(8)*vr(3,j)  
            vm(6,jj) = rot(3)*vr(1,j)+rot(6)*vr(2,j)+rot(9)*vr(3,j)  
            am(4,jj) = rot(1)*ar(1,j)+rot(4)*ar(2,j)+rot(7)*ar(3,j)  
            am(5,jj) = rot(2)*ar(1,j)+rot(5)*ar(2,j)+rot(8)*ar(3,j)  
            am(6,jj) = rot(3)*ar(1,j)+rot(6)*ar(2,j)+rot(9)*ar(3,j)  
          ENDDO
C---
          vs = zero                  
          as = zero                  
          DO jj=1,nir                                    
            vs(1) = vs(1) + h(jj) * vm(1,jj)                 
            vs(2) = vs(2) + h(jj) * vm(2,jj)                 
            vs(3) = vs(3) + hh(jj) *vm(3,jj)
     .                    + hrs(jj)*vm(4,jj) + hrt(jj) *vm(5,jj)
            vs(4) = vs(4) + hps(jj)*vm(3,jj) + hprs(jj)*vm(4,jj)
     .                                       + hxs(jj) *vm(5,jj)     
            vs(5) = vs(5) + hpt(jj)*vm(3,jj) + hprt(jj)*vm(5,jj)     
     .                                       + hxt(jj) *vm(4,jj)     
            vs(6) = vs(6) + h(jj) * vm(6,jj)                      
C
            as(1) = as(1) + h(jj) * am(1,jj)                 
            as(2) = as(2) + h(jj) * am(2,jj)                 
            as(3) = as(3) + hh(jj) *am(3,jj)
     .                    + hrs(jj)*am(4,jj)               
     .                    + hrt(jj)*am(5,jj)
            as(4) = as(4) + hps(jj)*am(3,jj) + hprs(jj)*am(4,jj)     
     .                                       + hxs(jj) *am(5,jj)     
            as(5) = as(5) + hpt(jj)*am(3,jj) + hprt(jj)*am(5,jj)     
     .                                       + hxt(jj) *am(4,jj)     
            as(6) = as(6) + h(jj) * am(6,jj)                      
          ENDDO    
C---
          a(1,is) = roti(1)*as(1)+roti(4)*as(2)+roti(7)*as(3) 
          a(2,is) = roti(2)*as(1)+roti(5)*as(2)+roti(8)*as(3)  
          a(3,is) = roti(3)*as(1)+roti(6)*as(2)+roti(9)*as(3) 
          v(1,is) = roti(1)*vs(1)+roti(4)*vs(2)+roti(7)*vs(3)
          v(2,is) = roti(2)*vs(1)+roti(5)*vs(2)+roti(8)*vs(3)
          v(3,is) = roti(3)*vs(1)+roti(6)*vs(2)+roti(9)*vs(3)
          ar(1,is)=roti(1)*as(4)+roti(4)*as(5)+roti(7)*as(6)   
          ar(2,is)=roti(2)*as(4)+roti(5)*as(5)+roti(8)*as(6)   
          ar(3,is)=roti(3)*as(4)+roti(6)*as(5)+roti(9)*as(6)   
          vr(1,is)=roti(1)*vs(4)+roti(4)*vs(5)+roti(7)*vs(6)   
          vr(2,is)=roti(2)*vs(4)+roti(5)*vs(5)+roti(8)*vs(6)   
          vr(3,is)=roti(3)*vs(4)+roti(6)*vs(5)+roti(9)*vs(6)   
C
        ENDIF
      ENDDO
C-----------
      RETURN