C Rutina 031
C-------------------------------------------------------------NOR_ELT---
C
C SUBRUTINA NOR_ELT
C
C CALCULO DE LA NORMAL A UN ELEMENTO EN UNA POSICION DADA
C
SUBROUTINE NOR_ELT(N,ELT,NOD)
C
IMPLICIT NONE
INCLUDE '_param.for'
INCLUDE '_struct.for'
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C NOD.VARIABLES EXTERNAS
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C N : Normal al elemento
REAL*8 N(MAX_DIM)
C ELT : Elemento
RECORD /S_elemento/ ELT
C NOD : Nodos del problema
RECORD /S_nodo/ NOD(N_MAX_NOD)
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C NOD.VARIABLES INTERNAS
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C MOD : Modulo del vector normal (antes de normalizar)
REAL*8 MOD
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VALORES INICIALES
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C CALCULO DEL VECTOR NORMAL
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
SELECT CASE (ELT.TIP)
C
CASE(1)
C
C ELEMENTO LINEAL CONTINUO
C
N(1)= NOD(ELT.NG(2)).XYZ(2)-NOD(ELT.NG(1)).XYZ(2)
N(2)=-NOD(ELT.NG(2)).XYZ(1)+NOD(ELT.NG(1)).XYZ(1)
MOD=DSQRT(N(1)**2.d0+N(2)**2.d0)
N(1)=N(1)/MOD
N(2)=N(2)/MOD
C
CASE DEFAULT
WRITE(F_W,03201) ELT.TIP
WRITE (*,03201) ELT.TIP
STOP
END SELECT
C
C FIN DE LA RUTINA
C
RETURN
C
C FORMATOS DE ESCRITURA
C
03200 FORMAT(/,2x,'ERROR 03200 - NOR_ELT - Error Programa (#',
@ i1,'#)')
03201 FORMAT(/,2x,'ERROR 03201 - NOR_ELT - Tipo de elemento ',i2
@ ,' no implementado.')
C
END
C
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
