C Rutina 013
C-------------------------------------------------------LEE_ELEMENTOS---
C
C SUBRUTINA LEE_ELEMENTOS
C
C LECTURA DE LOS ELEMENTOS DE UN CONTORNO
C
SUBROUTINE LEE_ELEMENTOS(CNT,ELTP,ELT)
C
IMPLICIT NONE
INCLUDE '_param.for'
INCLUDE '_struct.for'
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES EXTERNAS
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C CNT : Contorno
RECORD /S_contorno/ CNT
C ELTP : Tipo del elemento
RECORD /S_tipelt/ ELTP(N_MAX_TPEL)
C ELT : Elemento que se va a leer
RECORD /S_elemento/ ELT(N_MAX_ELT)
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES LOCALES
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C I,J,F,E,N1,N2 : Contadores
INTEGER i,j,e,f,sg,g1,g2,i1,i2,ng,ni
C TE : Tipo de elemento
INTEGER te
C PELT : Numero del primer elemento del contorno
INTEGER pelt
C AELT : Ultimo elmento generado
INTEGER aelt
C UELT : Ultimo elmento del contorno
INTEGER uelt
C CELT : Contador de los elementos
INTEGER celt
C NELT : Numero de elementos del contorno
INTEGER nelt
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VALORES INICIALES
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
nelt=CNT.NELT
pelt=0
uelt=1
aelt=0
celt=0
i=1
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C LECTURA VARIABLES ELEMENTO
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C BUCLE QUE RECORRE LOS ELEMENTOS DEL CONTORNO
C
DO WHILE (celt.LE.uelt)
C
READ(F_E,*) celt , ELT(celt).TIP , (ELT(celt).NG(j)
@ , j=1,ELTP(ELT(celt).TIP).NNG) , (ELT(celt).NI(j)
@ , j=1,ELTP(ELT(celt).TIP).NNI)
ELT(celt).NUM=celt
IF(pelt.EQ.0)THEN
pelt=celt
uelt=pelt-1+nelt
END IF
C
te=ELT(celt).TIP
IF(te.NE.ELTP(te).TIP)THEN
WRITE(F_W,01409) CNT.NUM , ELT(celt).NUM , ELT(celt).TIP
WRITE (*,01409) CNT.NUM , ELT(celt).NUM , ELT(celt).TIP
STOP
END IF
C
IF(celt.EQ.pelt) THEN
C
C PRIMER ELEMENTO DEL CONTORNO
C
CNT.ELT(i)=celt
aelt=celt
celt=celt+1
i=i+1
C
ELSE
C
C ELEMENTO CUALQUIERA
C
C Numero de nodos que hay que generar
e=ELT(celt).NUM-ELT(aelt).NUM
C
IF(e.eq.1) THEN
C UN SOLO ELEMENTO
CNT.ELT(i)=celt
aelt=celt
celt=celt+1
i=i+1
C
ELSE IF(e.GT.1.AND.(aelt+e).LE.N_MAX_NOD.
@ AND.(aelt+e).LE.uelt) THEN
C
C GENERACION AUTOMATICA
C
C Comprueba que la generacion automatica es posible
C
IF(ELT(aelt).TIP.NE.ELT(celt).TIP) THEN
WRITE(F_W,01404) ELT(celt-1).NUM , ELT(celt).NUM
WRITE (*,01404) ELT(celt-1).NUM , ELT(celt).NUM
STOP
END IF
sg=0
DO j = 2 , ELTP(te).NNG
g1=ELT(aelt).NG(j)-ELT(aelt).NG(j-1)
g2=ELT(celt).NG(j)-ELT(celt).NG(j-1)
IF(g1.NE.g2) THEN
WRITE(F_W,01404) ELT(aelt).NUM , ELT(celt).NUM
WRITE (*,01404) ELT(aelt).NUM , ELT(celt).NUM
STOP
END IF
sg=sg+g1
END DO
C
DO j = 2 , ELTP(te).NNI
i1=ELT(aelt).NI(j)-ELT(aelt).NI(j-1)
i2=ELT(celt).NI(j)-ELT(celt).NI(j-1)
IF(i1.NE.i2) THEN
WRITE(F_W,01404) ELT(aelt).NUM , ELT(celt).NUM
WRITE (*,01404) ELT(aelt).NUM , ELT(celt).NUM
STOP
END IF
END DO
C
ng=(ELT(celt).NG(1)-ELT(aelt).NG(1))/e
ni=(ELT(celt).NI(1)-ELT(aelt).NI(1))/e
IF(ng.NE.INT(ng).OR.ni.NE.INT(ni).OR.ng.NE.sg) THEN
WRITE(F_W,01404) ELT(aelt).NUM , ELT(celt).NUM
WRITE (*,01404) ELT(aelt).NUM , ELT(celt).NUM
STOP
ELSE
C
C Interpola los elementos
C
f=celt-1
DO WHILE (ELT(aelt).NUM.LE.f)
ELT(aelt+1).NUM=ELT(aelt).NUM+1
ELT(aelt+1).TIP=ELT(aelt).TIP
DO j = 1 , ELTP(te).NNG
ELT(aelt+1).NG(j)=ELT(aelt).NG(j)+ng
END DO
DO j = 1 , ELTP(te).NNI
ELT(aelt+1).NI(j)=ELT(aelt).NI(j)+ni
END DO
CNT.ELT(i)=aelt+1
aelt=aelt+1
i=i+1
END DO
celt=aelt+1
C
END IF
C
ELSE IF((aelt+e).GE.uelt) THEN
C DEMASIADOS ELEMENTOS
WRITE(F_W,01408) CNT.NUM , ELT(celt).NUM
WRITE (*,01408) CNT.NUM , ELT(celt).NUM
STOP
C
ELSE IF((aelt+e).GE.N_MAX_NOD) THEN
C DEMASIADOS ELEMENTOS
WRITE(F_W,01406) CNT.NUM , ELT(celt).NUM
WRITE (*,01406) CNT.NUM , ELT(celt).NUM
STOP
C
ELSE
C ELEMENTOS NO ORDENADOS DE FORMA CRECIENTE
WRITE(F_W,01405) CNT.NUM , ELT(celt).NUM
WRITE (*,01405) CNT.NUM , ELT(celt).NUM
STOP
END IF
C
END IF
C
C FIN BUCLE QUE RECORRE LOS ELEMENTOS DEL CONTORNO
C
END DO
C
C DEVUELVE EL ULTIMO ELEMENTO LEIDO
C
pelt=celt
C
C FIN DE LA RUTINA
C
RETURN
C
C FORMATOS DE ESCRITURA. SUBRUTINA ENTRADA
C
01400 FORMAT(/,2x,'ERROR 01400 - LEE_ELEMENTOS - Error Programa (#',
@ i1,'#)')
01401 FORMAT(/,2x,'ERROR 01401 - LEE_ELEMENTOS',/,' Nodo(G) ',i4,
@ ' del elemento ',i4,' fuera del rango permitido.')
01402 FORMAT(/,2x,'ERROR 01402 - LEE_ELEMENTOS',/,' Nodo(I) ',i4,
@ ' del elemento ',i4,' fuera del rango permitido.')
01403 FORMAT(/,2x,'ERROR 01403 - LEE_ELEMENTOS',/,' Contorno: ',i2,
@ '. Primer elemento no consecutivo.')
01404 FORMAT(/,2x,'ERROR 01404 - LEE_ELEMENTOS',/,' Elementos ',i4,
@ ' y',i4,' no compatibles para generacion automatica.')
01405 FORMAT(/,2x,'ERROR 01405 - LEE_ELEMENTOS - Contorno: ',i2,'.',/,
@ ' Elementos no ordenados de forma creciente (',i4,').')
01406 FORMAT(/,2x,'ERROR 01406 - LEE_ELEMENTOS - Contorno: ',i2,'.',/,
@ ' Elemento ',i4,' fuera del rango permitido.')
01407 FORMAT(/,2x,'ERROR 01407 - LEE_ELEMENTOS',/,
@ ' Los nodos del elemento ',i4,' deben ser consecutivos',
@ /,' para generacion automatica.')
01408 FORMAT(/,2x,'ERROR 01408 - LEE_ELEMENTOS - Contorno: ',i2,'.',/,
@ ' Elemento ',i4,' fuera del rango definido.')
01409 FORMAT(/,2x,'ERROR 01409 - LEE_ELEMENTOS - Contorno: ',i2,',',
@ ' elemento ',i4,'.',/,' Tipo de elmento (',i2,
@ ') no implementado.')
C
END
C
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
