C Rutina 032
C--------------------------------------------------------------COND_CONT---
C
C SUBRUTINA COND_CONT
C
C CALCULO DE LAS GIJ (Condicion de contorno) EN EL NODO k
C (NODO 2D EXTREMO DE UN ELEMENTO LINEAL CONTINUO)
C
SUBROUTINE COND_CONT(NOD,Na,Nd,COD,CINC,MAT,TP)
C
IMPLICIT NONE
INCLUDE '_param.for'
INCLUDE '_struct.for'
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES EXTERNAS
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C NOD : Nodo actual (k)
RECORD /S_nodo/ NOD
C Na(),Nd() : Normal antes y despues del nodo
REAL*8 Na(MAX_DIM),Nd(MAX_DIM)
C COD : Cod=1 - Nodo normal
C Cod=2 - Nodo inicial (No existe k-1 => k nodo(1) del elto k)
C Cod=3 - Nodo final (No existe k+1 => k nodo(2) del elto k)
INTEGER COD
C CINC : Contador del numero de incognitas
INTEGER CINC
C MAT : Material
RECORD /S_material/ MAT
C TP : Tipo de problema
INTEGER TP
C
CCCCCCCCCCCCCCC
C RUTINA
CCCCCCCCCCCCCCC
C
SELECT CASE(TP)
CASE(1)
C PROBLEMA AXISIMETRICO
CALL COND_CONT1(NOD,Na,Nd,COD,CINC,MAT)
CASE(21,22,31,32,51,52)
C PROBLEMA PLANO
IF(COD.EQ.1)THEN
CALL COND_CONT2(NOD,Na,Nd,CINC,MAT)
ELSE
WRITE(F_W,04100) 2
WRITE( *,04100) 2
STOP
END IF
CASE DEFAULT
WRITE(F_W,04100) 1
WRITE( *,04100) 1
STOP
END SELECT
C
C FORMATOS DE ESCRITURA
C
04100 FORMAT(/,2x,'ERROR 04100 - COND_CONT - Error Programa (#',
@ i1,'#)')
C
END
C
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
