cevec3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	cevec3 (elbuf_str, dir_a, dir_b, jft, jlt, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, e1x, e1y, e1z, e2x, e2y, e2z, e3x, e3y, e3z, irep, nlay, nel)

Function/Subroutine Documentation

cevec3()

subroutine cevec3	(	type (elbuf_struct_)	elbuf_str,
			dir_a,
			dir_b,
		integer	jft,
		integer	jlt,
			x1,
			x2,
			x3,
			x4,
			y1,
			y2,
			y3,
			y4,
			z1,
			z2,
			z3,
			z4,
			e1x,
			e1y,
			e1z,
			e2x,
			e2y,
			e2z,
			e3x,
			e3y,
			e3z,
		integer	irep,
		integer	nlay,
		integer	nel )

Definition at line 32 of file cevec3.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE elbufdef_mod
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   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER JFT, JLT,J,IREP,NLAY,NEL
      my_real
     .   x1(*), x2(*), x3(*), x4(*), y1(*), y2(*), y3(*), y4(*),
     .   z1(*), z2(*), z3(*), z4(*), e1x(*), e1y(*), e1z(*), e2x(*),
     .   e2y(*), e2z(*), e3x(*), e3y(*), e3z(*), dir_a(*),dir_b(*)
      TYPE (ELBUF_STRUCT_) :: ELBUF_STR
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
C     REAL
      my_real
     .   x31(mvsiz), y31(mvsiz), z31(mvsiz), x42(mvsiz), y42(mvsiz),
     .   z42(mvsiz), x21(mvsiz), y21(mvsiz), z21(mvsiz),
     .   sum(mvsiz), suma
C=======================================================================
      DO i=jft,jlt                                      
        x21(i)=x2(i)-x1(i)                              
        y21(i)=y2(i)-y1(i)                              
        z21(i)=z2(i)-z1(i)                              
        x31(i)=x3(i)-x1(i)                              
        y31(i)=y3(i)-y1(i)                              
        z31(i)=z3(i)-z1(i)                              
        x42(i)=x4(i)-x2(i)                              
        y42(i)=y4(i)-y2(i)                              
        z42(i)=z4(i)-z2(i)                              
C
        e3x(i)=y31(i)*z42(i)-z31(i)*y42(i)              
        e3y(i)=z31(i)*x42(i)-x31(i)*z42(i)              
        e3z(i)=x31(i)*y42(i)-y31(i)*x42(i)              
        suma=e3x(i)*e3x(i)+e3y(i)*e3y(i)+e3z(i)*e3z(i)  
        suma=max(sqrt(suma),em20)                       
        e3x(i)=e3x(i)/suma                              
        e3y(i)=e3y(i)/suma                              
        e3z(i)=e3z(i)/suma                              
      ENDDO                                             
C
      DO i=jft,jlt                                       
        suma= x21(i)*e3x(i)+y21(i)*e3y(i)+z21(i)*e3z(i)  
        e1x(i)= x21(i)-e3x(i)*suma                       
        e1y(i)= y21(i)-e3y(i)*suma                       
        e1z(i)= z21(i)-e3z(i)*suma                       
      ENDDO                                              
C
      DO i=jft,jlt                                      
        suma=e1x(i)*e1x(i)+e1y(i)*e1y(i)+e1z(i)*e1z(i)  
        suma=max(sqrt(suma),em20)                       
        e1x(i)=e1x(i)/suma                              
        e1y(i)=e1y(i)/suma                              
        e1z(i)=e1z(i)/suma                              
      ENDDO                                             
C
      DO i=jft,jlt                                        
        e2x(i)=e3y(i)*e1z(i)-e3z(i)*e1y(i)                
        e2y(i)=e3z(i)*e1x(i)-e3x(i)*e1z(i)                
        e2z(i)=e3x(i)*e1y(i)-e3y(i)*e1x(i)                
        suma  =e2x(i)*e2x(i)+e2y(i)*e2y(i)+e2z(i)*e2z(i)  
        suma  =max(sqrt(suma),em20)                       
        e2x(i)=e2x(i)/suma                                
        e2y(i)=e2y(i)/suma                                
        e2z(i)=e2z(i)/suma                                
      ENDDO                                               
C
C---  directions orthotropie / anisotropie convectes
C
      CALL cortdir3(elbuf_str,dir_a,dir_b  ,jft    ,jlt    ,
     .              nlay   ,irep   ,x21    ,y21    ,z21    ,
     .              x31    ,y31    ,z31    ,e1x    ,e1y    ,   
     .              e1z    ,e2x    ,e2y    ,e2z    ,nel    )
C-----------
      RETURN