i2vit10.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

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

Function/Subroutine Documentation

i2rot10()

subroutine i2rot10	(	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(*)	irupt,
		integer, dimension(*)	weight )

Definition at line 102 of file i2vit10.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(*),IRUPT(*),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, II, JJ, L
C     REAL
      my_real
     . h(4), ss, tt, amx, amy, amz, vmx, vmy, vmz,sp,sm,tp,tm,
     . v1x2, v2x1, v2x3, v3x2, v3x1, v1x3,xc0,yc0,zc0
C=======================================================================
      nir=4
C
      DO ii=1,nsn
        i = nsv(ii)
        l = irtl(ii)
C
        IF(i > 0 .AND. irupt(ii) == 0) THEN
          ss=crst(1,ii)                             
          tt=crst(2,ii)                             
          sp=one + ss                               
          sm=one - ss                               
          tp=fourth*(one + tt)                        
          tm=fourth*(one - tt)                        
          h(1)=tm*sm                                
          h(2)=tm*sp                                
          h(3)=tp*sp                                
          h(4)=tp*sm                                
C
          xc0 = x(1,i)                              
          yc0 = x(2,i)                              
          zc0 = x(3,i)                              
C
          amx=zero                                    
          amy=zero                                    
          amz=zero                                    
          vmx=zero                                    
          vmy=zero                                    
          vmz=zero                                    
          DO jj=1,4                                 
            j=irect(jj,l)                           
            amx=amx+ar(1,j)*h(jj)                   
            amy=amy+ar(2,j)*h(jj)                   
            amz=amz+ar(3,j)*h(jj)                   
            vmx=vmx+vr(1,j)*h(jj)                   
            vmy=vmy+vr(2,j)*h(jj)                   
            vmz=vmz+vr(3,j)*h(jj)                   
            xc0=xc0 - x(1,j) * h(jj)                
            yc0=yc0 - x(2,j) * h(jj)                
            zc0=zc0 - x(3,j) * h(jj)                
          ENDDO                                     
          ar(1,i)=amx                               
          ar(2,i)=amy                               
          ar(3,i)=amz                               
          vr(1,i)=vmx                               
          vr(2,i)=vmy                               
          vr(3,i)=vmz                               
          v1x2=vmx * yc0                            
          v2x1=vmy * xc0                            
          v2x3=vmy * zc0                            
          v3x2=vmz * yc0                            
          v3x1=vmz * xc0                            
          v1x3=vmx * zc0                            
 
          a(1,i)= a(1,i) + amy * zc0 - amz * yc0    
     .          - half*(vmy*v2x1+vmz*v3x1)           
          a(2,i)= a(2,i) + amz * xc0 - amx * zc0    
     .          - half*(vmz*v3x2+vmx*v1x2)           
          a(3,i)= a(3,i) + amx * yc0 - amy * xc0    
     .          - half*(vmx*v1x3+vmy*v2x3)           
          v(1,i)= v(1,i) + v2x3 - v3x2              
          v(2,i)= v(2,i) + v3x1 - v1x3              
          v(3,i)= v(3,i) + v1x2 - v2x1              
        ENDIF                                       
      ENDDO
C----
      RETURN

i2vit10()

subroutine i2vit10	(	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(*)	irupt,
		integer, dimension(*)	weight )

Definition at line 28 of file i2vit10.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(*),IRUPT(*),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, II, JJ, L
C     REAL
      my_real
     . h(4), ss, tt, amx, amy, amz, vmx, vmy, vmz,sp,sm,tp,tm
C=======================================================================
      nir=4
C
      DO ii=1,nsn
        i = nsv(ii)
        l = irtl(ii)
C
        IF(i > 0 .AND. irupt(ii) == 0) THEN
           ss=crst(1,ii)
           tt=crst(2,ii)
           sp=one + ss
           sm=one - ss
           tp=fourth*(one + tt)
           tm=fourth*(one - tt)
           h(1)=tm*sm
           h(2)=tm*sp
           h(3)=tp*sp
           h(4)=tp*sm
           amx=zero
           amy=zero
           amz=zero
           vmx=zero
           vmy=zero
           vmz=zero    
C
           DO jj=1,nir
             j=irect(jj,l)
             amx=amx+a(1,j)*h(jj)
             amy=amy+a(2,j)*h(jj)
             amz=amz+a(3,j)*h(jj)
             vmx=vmx+v(1,j)*h(jj)
             vmy=vmy+v(2,j)*h(jj)
             vmz=vmz+v(3,j)*h(jj)
           ENDDO
           a(1,i)=amx
           a(2,i)=amy
           a(3,i)=amz
           v(1,i)=vmx
           v(2,i)=vmy
           v(3,i)=vmz
C
        ENDIF
      ENDDO
C----
      RETURN