C Rutina 101
C------------------------------------------------------NODO_LIBRE_DEL---
C
C SUBRUTINA NODO_LIBRE_DEL
C
C MONTAJE DE LAS CONDICIONES DE CONTORNO CORRESPONDIENTES A
C UN NODO LIBRE DE LA ZONA DE CONTACTO
C
C CONTACTO DEBIL
C
SUBROUTINE NODO_LIBRE_DEL(NOD,SAij,ceq,COD)
C
IMPLICIT NONE
C
INCLUDE '_param.for'
INCLUDE '_struct.for'
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES EXTERNAS
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C NOD : Nodo
RECORD /S_nodo/ NOD
C SAIJ(,) : Matriz del sistema.
REAL*8 SAij(N_MAX_INC,N_MAX_INC)
C CEQ : Contador de ecuacuones de contacto
INTEGER CEQ
C COD : Localizacion del nodo
C = 1 , nodo origen
C = 2 , nodo intermedio
C = 3 , nodo extremo
INTEGER COD
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES INTERNAS
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C I : Contador
INTEGER I
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VALORES INICIALES
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C ECUACIONES DE COMPATIBILIDAD
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
SELECT CASE(COD)
CASE(1)
DO i=1,NOD.NINC
IF(NOD.VAR(i).EQ.13)THEN
SAij(ceq,NOD.INC(i))=1.d0
END IF
IF(NOD.VAR(i).EQ.14)THEN
SAij(ceq+1,NOD.INC(i))=1.d0
END IF
END DO
CASE(2,3)
DO i=1,NOD.NINC
IF(NOD.VAR(i).EQ.11)THEN
SAij(ceq,NOD.INC(i))=1.d0
END IF
IF(NOD.VAR(i).EQ.12)THEN
SAij(ceq+1,NOD.INC(i))=1.d0
END IF
END DO
CASE DEFAULT
WRITE(F_W,10300) 1
WRITE( *,10300) 1
STOP
END SELECT
C
C FIN DE LA RUTINA
C
RETURN
C
C FORMATOS DE ESCRITURA.
C
10300 FORMAT(/,2x,'ERROR 10300 - NODO_LIBRE_DEL - Error Programa (#',
@ i1,'#)')
C
END
C
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
