C Rutina 013
C-----------------------------------------------------------LEE_NODOS---
C
C SUBRUTINA LEE_NODOS
C
C LECTURA DE DATOS DE LOS NODOS
C
SUBROUTINE LEE_NODOS(ETAP,SOL,NOD,DIM)
C
IMPLICIT NONE
INCLUDE '_param.for'
INCLUDE '_struct.for'
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES EXTERNAS
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C ETAP : Etapa de carga
INTEGER ETAP
C SOL : Solido
RECORD /S_solido/ SOL
C NOD : Nodos
RECORD /S_nodo/ NOD(N_MAX_NOD)
C DIM : Numero de coordenadas de los nodos
INTEGER DIM
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES LOCALES
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C PNOD : Primer nodo del solido
INTEGER PNOD
C UNOD : Ultimo nodo del solido
INTEGER UNOD
C I : Contador
INTEGER I
C N,CNOD : Numero de nodo local/global
INTEGER N,CNOD
C XYZ() : Coordenadas
REAL*8 XYZ(MAX_DIM)
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VALORES INICIALES
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
i=0
pnod=0
cnod=0
unod=1
n=1
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C LECTURA VARIABLES NODALES
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C BUCLE QUE RECORRE LOS NODOS DEL SOLIDO
C
DO WHILE (cnod.LT.UNOD)
C
C LECTURA DE DATOS DE UN NODO
C
READ(F_E,*) cnod , (XYZ(i) , i=1,DIM)
C
IF(n.EQ.1) THEN
C Primer nodo del sólido
pnod=cnod
unod=SOL.NNOD+pnod-1
C
NOD(cnod).NUM=cnod
SOL.NOD(n)=cnod
DO i = 1 , DIM
NOD(cnod).XYZo(i)=XYZ(i)
IF(ETAP.EQ.1)THEN
NOD(cnod).XYZ(i)=NOD(cnod).XYZo(i)
END IF
END DO
n=n+1
pnod=pnod+1
C
ELSE IF(pnod.LE.cnod.AND.cnod.LE.(unod+1)) THEN
C Genera los nodos entre pnod y cnod (ambos inclusive)
NOD(cnod).NUM=cnod
DO i = 1 , DIM
NOD(cnod).XYZo(i)=XYZ(i)
IF(ETAP.EQ.1)THEN
NOD(cnod).XYZ(i)=NOD(cnod).XYZo(i)
END IF
XYZ(i)=(XYZ(i)-NOD(pnod-1).XYZo(i))/(cnod-pnod+1)
END DO
DO WHILE (pnod.LT.cnod)
NOD(pnod).NUM=pnod
DO i= 1,DIM
NOD(pnod).XYZo(i)=NOD(cnod).XYZo(i)-(cnod-pnod)*XYZ(i)
IF(ETAP.EQ.1)THEN
NOD(pnod).XYZ(i)=NOD(pnod).XYZo(i)
END IF
SOL.NOD(n)=pnod
END DO
pnod=pnod+1
n=n+1
END DO
SOL.NOD(n)=cnod
pnod=pnod+1
n=n+1
xC
xC Borra el ultimo nodo cuando es NNOD+1
xC
x IF(n.EQ.SOL.NNOD+2) THEN
x NOD(n-1).NUM=0
x END IF
ELSE
C
C El nodo leido no es valido
C
IF(N.LT.cnod) THEN
WRITE(F_W,01202) cnod,SOL.NUM
WRITE (*,01202) cnod,SOL.NUM
STOP
ELSE IF(cnod.GT.unod+1) THEN
WRITE(F_W,01203) SOL.NUM
WRITE (*,01203) SOL.NUM
STOP
ELSE
WRITE(F_W,01200) 1
WRITE(*, 01200) 1
STOP
END IF
END IF
C
C FIN DEL BUCLE QUE RECORRE LOS NODOS DEL SOLIDO
C
END DO
C
C FIN DE LA RUTINA
C
RETURN
C
C FORMATOS DE ESCRITURA
C
01200 FORMAT(/,2x,'ERROR 01200 - LEE_NODOS - Error Programa (#',i1,'#)')
01201 FORMAT(/,2x,'ERROR 01201 - LEE_NODOS',/,
* ' Primer nodo del solido ',i2,' erroneo.')
01202 FORMAT(/,2x,'ERROR 01202 - LEE_NODOS',/,
* ' Nodo',i4,' del solido ',i2,' repetido.')
01203 FORMAT(/,2x,'ERROR 01203 - LEE_NODOS',/,
* ' El solido ',i2,' tiene mas nodos de los esperados')
C
END
C
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
