C C-------------------------------------------------------------ESC_SAij--- C C SUBRUTINA ESC_SAij C C ESCRITURA DE LA SOLUCION EN EL FICHERO DE SALIDA(POTENCIAL) C SUBROUTINE ESC_SAij(PRB,SOL,NOD,SAij,SFi) 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 NOD : Nodos del problema RECORD /S_nodo/ NOD(N_MAX_NOD) C SFi() : Vector de incógnitas. REAL*8 SFi(N_MAX_INC) C SAIJ(,) : Matriz del sistema. REAL*8 SAij(N_MAX_INC,N_MAX_INC) C C CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C VARIABLES INTERNAS CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC C C M,N : Contadores INTEGER M,N,H,K C C K=0 DO m=1,PRB.NSOL DO N=SOL(M).NOD(1),SOL(M).NOD(SOL(M).NNOD) DO H=1,NOD(N).NINC K=K+1 WRITE(F_L,*) K,' NODO ',N,' - VARIABLE ' @ , PRB.NDTP.VTXT(NOD(N).VAR(H)) END DO END DO END DO WRITE (F_L,*) 'TOTAL ',K,' VARIABLES' C Bucle que recorre las FILAS DE SAij C WRITE(F_L,*) WRITE(F_L,*) ' ++++++++' WRITE(F_L,*) ' SAij ' WRITE(F_L,*) ' ++++++++' WRITE(F_L,*) C DO m = 1 , PRB.NAP+PRB.NBIE+PRB.NCE WRITE(F_L,15) m,(SAij(m,n),n=1 , PRB.NINC) C C Fin del bucle que recorre las FILAS DE SAij C END DO C C Bucle que recorre las FILAS DE SXj C WRITE(F_L,*) WRITE(F_L,*) ' ++++++++' WRITE(F_L,*) ' SFi ' WRITE(F_L,*) ' ++++++++' WRITE(F_L,*) C DO m = 1 , PRB.NAP+PRB.NBIE+PRB.NCE WRITE(F_L,10) m,SFi(m) C C Fin del bucle que recorre las FILAS DE SXj C END DO C C Fin del bucle que recorre los solidos C 10 FORMAT(4x,i4,4x,e22.16) 15 FORMAT(4x,i4,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16 @ ,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16,4x,e22.16) C C FIN DE LA RUTINA C RETURN C C FORMATOS DE ESCRITURA C END C