C Rutina 080
C-------------------------------------------------------------ESC_ENT---
C
C SUBRUTINA ESC_ENT
C
C ESCRITURA DE LOS DATOS DE ENTRADA EN EL FICHERO DE SALIDA
C
SUBROUTINE ESC_ENT(PRB,SOL,MAT,CNT,ELT,NOD,PINT)
C
IMPLICIT NONE
INCLUDE '_param.for'
INCLUDE '_struct.for'
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES EXTERNAS
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C PRB : Variables globales del problema
RECORD /S_problema/ PRB
C SOL : Solidos del problema
RECORD /S_solido/ SOL(N_MAX_SOL)
C MAT : Materiales del problema
RECORD /S_material/ MAT(N_MAX_SOL)
C CNT : Contornos del problema
RECORD /S_contorno/ CNT(N_MAX_CNT)
C ELT : Elementos del problema
RECORD /S_elemento/ ELT(N_MAX_ELT)
C NOD : Nodos del problema
RECORD /S_nodo/ NOD(N_MAX_NOD)
C PINT : Puntos internos del problema
RECORD /S_pntint/ PINT(N_MAX_PINT)
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VARIABLES INTERNAS
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
C I,J,K,M,N : Contadores
INTEGER I,J,K,M,N
C CSOL,CCNT,CELT,PNOD,ANOD,CNOD,CPNT : Contadores
INTEGER CSOL,CCNT,CELT,PNOD,ANOD,CNOD,CPNT
C CMATP,CELTP,CNDTP : Contador
INTEGER CMATP,CELTP,CNDTP
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C VALORES INICIALES
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
ANOD=0
PNOD=0
C
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C ESCRITURA FICHERO DE SALIDA
CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC
C
WRITE(F_S,*)
WRITE(F_S,*) ' Titulo: ',PRB.TITULO
C
WRITE(F_S,*)
WRITE(F_S,10) 'Tipo de problema : ',PRB.TP.TIP,' - ',
@ PRB.TP.TXT
WRITE(F_S,10) 'Tipo de formulacion : ',PRB.TF.TIP,' - ',
@ PRB.TF.TXT
WRITE(F_S,10) 'Tipo de integracion : ',PRB.TI.TIP,' - ',
@ PRB.TI.TXT
WRITE(F_S,10) ' Puntos de Gauss : ',PRB.TI.NPG
C
WRITE(F_S,10) 'Numero de solidos....: ',PRB.NSOL
C
WRITE(F_S,*)
C
WRITE(F_S,*) ' *****************'
WRITE(F_S,10) ' TIPOS DE NODO '
WRITE(F_S,*) ' *****************'
DO i = 1 , PRB.NNDTP
WRITE(F_S,10) ' Tipo : ',PRB.NDTP(i).TIP,' - ',
@ PRB.NDTP(i).TXT
END DO
C
WRITE(F_S,*)
WRITE(F_S,*) ' *********************'
WRITE(F_S,10) ' TIPOS DE ELEMENTO '
WRITE(F_S,*) ' *********************'
DO i = 1 , PRB.NELTP
WRITE(F_S,10) ' Tipo : ',PRB.ELTP(i).TIP,' - ',
@ PRB.ELTP(i).TXT
END DO
C
WRITE(F_S,*)
C
WRITE(F_S,*) ' ****************************'
WRITE(F_S,10) ' TIPOS DE PUNTOS INTERNOS '
WRITE(F_S,*) ' ****************************'
DO i = 1 , PRB.NPINTP
WRITE(F_S,10) ' Tipo : ',PRB.PINTP(i).TIP,' - ',
@ PRB.PINTP(i).TXT
END DO
C
C Bucle que recorre los solidos
C
DO i = 1 , PRB.NSOL
csol=SOL(i).NUM
WRITE(F_S,*)
WRITE(F_S,*) ' ***********************'
WRITE(F_S,10) ' DATOS DEL SOLIDO ',csol
WRITE(F_S,*) ' ***********************'
WRITE(F_S,*)
cmatp=MAT(i).TIP
WRITE(F_S,10) ' Tipo de material :',
@ cmatp,' - ',PRB.MATP(cmatp).TXT
DO j = 1 , PRB.MATP(cmatp).NPRP
WRITE(F_S,12) PRB.MATP(cmatp).PTXT(j) , MAT(i).PRP(j)
END DO
WRITE(F_S,*)
WRITE(F_S,*) ' ++++++++++++++++++++'
WRITE(F_S,*) ' NODOS DEL SOLIDO '
WRITE(F_S,*) ' ++++++++++++++++++++'
WRITE(F_S,*)
WRITE(F_S,*) ' Nodo'
@ ,(' Coordenada ',PRB.TP.DTXt(k),' ',k=1,PRB.TP.DIM)
WRITE(F_S,*)
DO j = 1 , SOL(csol).NNOD
cnod=SOL(csol).NOD(j)
WRITE(F_S,14) NOD(cnod).NUM,(NOD(cnod).XYZo(k),k=1,
@ PRB.TP.DIM)
END DO
WRITE(F_S,*)
WRITE(F_S,*) ' ++++++++++++++++++++++++'
WRITE(F_S,*) ' CONTORNOS DEL SOLIDO '
WRITE(F_S,*) ' ++++++++++++++++++++++++'
WRITE(F_S,*)
DO j = 1 , SOL(csol).NCNT
ccnt=SOL(csol).CNT(j)
WRITE(F_S,*) ' ---------------'
WRITE(F_S,10) ' CONTORNO ',ccnt
WRITE(F_S,*) ' ---------------'
WRITE(F_S,*)
WRITE(F_S,10) ' Numero de elementos: ',CNT(ccnt).NELT
WRITE(F_S,*)
WRITE(F_S,*) ' Elto Tipo Nodos'
WRITE(F_S,*)
DO k = 1 , CNT(ccnt).NELT
celt=CNT(ccnt).ELT(k)
celtp=ELT(celt).TIP
WRITE(F_S,16) ELT(celt).NUM,celtp
@ ,(' - NG(',m,'): ',ELT(celt).NG(m),m=1,PRB.ELTP(celtp).NNG)
@ ,(' - NI(',m,'): ',ELT(celt).NI(m),m=1,PRB.ELTP(celtp).NNI)
END DO
END DO
WRITE(F_S,*)
WRITE(F_S,*) ' +++++++++++++++++++++++++++'
WRITE(F_S,*) ' CONDICIONES DE CONTORNO '
WRITE(F_S,*) ' +++++++++++++++++++++++++++'
WRITE(F_S,*)
WRITE(F_S,*) ' Nodo Cdc. Valores'
WRITE(F_S,*)
cndtp=PRB.ELTP(ELT(CNT(SOL(csol).CNT(1)).ELT(1)).TIP).TNI(1)
DO j = 1 , SOL(csol).NCNT
ccnt=SOL(csol).CNT(j)
DO k = 1 , CNT(ccnt).NELT
celt=CNT(ccnt).ELT(k)
celtp=ELT(celt).TIP
DO m = 1 , PRB.ELTP(celtp).NNI
cnod=ELT(celt).NI(m)
cndtp=PRB.ELTP(celtp).TNI(m)
IF(cnod.NE.anod) THEN
WRITE(F_S,18) cnod,(NOD(cnod).CBC(n),n=1,N_MAX_CODBC)
@ ,(PRB.NDTP(cndtp).VTXT(NOD(cnod).BC(n))
@ ,NOD(cnod).VBC(n),n=1,NOD(cnod).NBC)
anod=cnod
END IF
END DO
END DO
END DO
x pnod=ELT(CNT(SOL(csol).CNT(1)).ELT(1)).NI(1)
x cndtp=PRB.ELTP(ELT(CNT(SOL(csol).CNT(1)).ELT(1)).TIP).TNI(1)
x WRITE(F_S,18) pnod,NOD(pnod).CBC,
x @ (' - ',PRB.NDTP(cndtp).VTXT(NOD(pnod).BC(n)),
x @ NOD(pnod).VBC(n),n=1,NOD(pnod).NBC)
x DO j = 1 , SOL(csol).NCNT
x ccnt=SOL(csol).CNT(j)
x DO k = 1 , CNT(ccnt).NELT
x celt=CNT(ccnt).ELT(k)
x celtp=ELT(celt).TIP
x DO m = 1 , PRB.ELTP(celtp).NNI
x cnod=ELT(celt).NI(m)
x cndtp=PRB.ELTP(celtp).TNI(m)
x IF(cnod.NE.anod.AND.cnod.NE.pnod) THEN
x WRITE(F_S,18) cnod,NOD(cnod).CBC,
x @ (' - ',PRB.NDTP(cndtp).VTXT(NOD(cnod).BC(n)),
x @ NOD(cnod).VBC(n),n=1,NOD(cnod).NBC)
x anod=cnod
x END IF
x END DO
x END DO
x END DO
IF(SOL(csol).NPINT.GT.0)THEN
WRITE(F_S,*)
WRITE(F_S,*) ' ++++++++++++++++++++++++++++++'
WRITE(F_S,*) ' PUNTOS INTERNOS DEL SOLIDO '
WRITE(F_S,*) ' ++++++++++++++++++++++++++++++'
WRITE(F_S,*)
WRITE(F_S,*) ' Puntos internos:',
@ (' Coordenada ',PRB.TP.DTXT(k),' ',k=1,PRB.TP.DIM)
WRITE(F_S,*)
DO j = 1 , SOL(csol).NPINT
cpnt=SOL(csol).PINT(j)
WRITE(F_S,14) PINT(cpnt).NUM,
@ (PINT(cpnt).XYZ(k),k=1,PRB.TP.DIM)
END DO
END IF
C
C Fin bucle que recorre los solidos
C
END DO
C
C FIN DE LA RUTINA
C
RETURN
C
C FORMATOS DE ESCRITURA
C
10 FORMAT(2x,a,i2,a,a)
12 FORMAT(4x,a,': ',e22.16)
14 FORMAT(4x,i4,2x,e22.16,2x,e22.16,2x,e22.16,2x)
16 FORMAT(6x,i4,4x,i2,a,i1,a,i4,a,i1,a,i4,a,i1,a,i4,a,i1,a,i4)
18 FORMAT(6x,i4,4x,3(i2,1x),4('- ',a3,': ',e12.6,1x))
C
END
C
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
