C
C---------------------------------------------------------CONV_T_TERM---
C
C SUBRUTINA CONV_T_TERM
C
C CONVERGENCIA TERMICA(CALCULO DE EPSILON)
C
SUBROUTINE CONV_T_TERM(NSOL,PRBT,SOLT,MATT,CNTT,ELTT,NODT,
@ PINTT,ZCTOT,NPGT,CONV,NODTU,FIN)
C
IMPLICIT NONE
C
INCLUDE '_param.for'
INCLUDE '_struct.for'
C
CCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C VARIABLES EXTERNAS
C
CCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C PRBT : Variables globales del problema
RECORD /S_problema/ PRBT
C SOLT : Solidos del problema
RECORD /S_solido/ SOLT(N_MAX_SOL)
C MAT : Materiales
RECORD /S_material/ MATT(N_MAX_SOL)
C CNTT : Contornos del problema
RECORD /S_contorno/ CNTT(N_MAX_CNT)
C ELTT : Elementos del problema
RECORD /S_elemento/ ELTT(N_MAX_ELT)
C NODT,NODTU : Nodos del problema
RECORD /S_nodo/ NODT(N_MAX_NOD)
RECORD /S_nodo/ NODTU(N_MAX_NOD)
C PINTT : Puntos internos del problema
RECORD /S_pntint/ PINTT(N_MAX_PINT)
C ZCTOT : Zonas de contacto del problema
RECORD /S_zncto/ ZCTOT(N_MAX_ZCTO)
C NPGT : Numero de puntos de Gauss
INTEGER NPGT
C
CCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C VARIABLES EXTERNAS
C
CCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C CONV,FIN,NSOL : Indicadores
INTEGER CONV,FIN,NSOL
C
C
CCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C VARIABLES INTERNAS
C
CCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C CSOL : Indicadores
INTEGER CSOL
C EPSILON,NUM,SUMA1,SUMA2,SUM1,SUM2 : Acumuladores,sumadores
REAL*8 EPSILON,NUM,SUMA1,SUMA2,SUM1,SUM2
C
C
CCCCCCCCCCCCCCC
C
C RUTINA
C
CCCCCCCCCCCCCCC
C
C INICIALIZAMOS LAS SUMAS
C
SUMA1=0
SUMA2=0
C
C CALCULAMOS LA SUMA PARA TODOS LOS NODOS
C
C Bucle que recorra todos los solidos
c
DO csol = 1, nsol
C
CALL CONV1_T_TERM(SOLT(csol),NODT,NODTU,SUM1,SUM2)
C
SUMA1=SUMA1+SUM1
SUMA2=SUMA2+SUM2
C
END DO
C
C CALCULO DE NUM(comparar con epsilon)
C
C Modificado por: Diógenes Morales Sierra Fecha:15/07/05
C Objeto de la modificación:
C En el caso de temperaturas nulas, aparece un float divide by cero
C para corregirlo uso un bucle if y dentro un select con la convergencia
C
IF(SUMA2.EQ.0.d0)THEN
SELECT CASE(CONV)
CASE(0)
CONV=1
CASE(1)
FIN=1
END SELECT
ELSE
NUM=SUMA1/SUMA2
C
EPSILON=0.0002
C
C COMPARAMOS EPSILON Y NUM
C
IF(NUM.LE.EPSILON)THEN
C
SELECT CASE(CONV)
C
CASE(0)
C
CONV=1
C
CASE(1)
C
FIN=1
C
CASE DEFAULT
WRITE(F_W,01000)
WRITE(*,01000)
END SELECT
C Modificado por: Diógenes Morales Sierra Fecha:20/09/05
C Objeto de la modificación:
C QUE EN EL CASO DE HABERSE ACERCADO A LA CONVERGENCIA, SI EN LA SIGUIENTE
C ITERACION SE ALEJA, QUE DEJE DE INDICARSE LA CONVERGENCIA
C
ELSE
CONV=0
END IF
END IF
C
C FIN DE LA RUTINA LECTURA
C
RETURN
C
C FORMATOS DE ESCRITURA.
C
01000 FORMAT(/,2x,'ERROR 01000 - CONV_T_TERM - Error Programa')
C
END
C
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
