! Trab1_TDMA2D_Transiente.f90 ! ! FUNCTIONS: ! Trab1_TDMA2D_Transiente - Entry point of console application. ! !**************************************************************************** ! ! PROGRAM: Trab1_TDMA2D_Transiente ! ! PURPOSE: Entry point for the console application. ! !**************************************************************************** program Trab1_TDMA2D_Transiente implicit none integer, parameter:: nvi=10 integer, parameter:: nvj=10 integer:: i, j, kk real:: L, W, qx, qy, T_ini, T_inf, cp, k, rho, h, dx, dy, tempo, tempof, dt real:: residuo, RMS, tol real, dimension(nvi, nvj):: Ap, Ap_0, Aw, Ae, An, As, Su, Sp real, dimension(nvi, nvj):: T, Tp_0 real, dimension(nvi):: x real, dimension(nvj):: y dt=1 L=1D0 W=1D0 T_ini=200.0D0 qx=0.0D0 qy=0.0D0 cp=4183D0 k=0.586D0 h=100D0 rho=998D0 T_inf=0.0D0 dx=L/nvi dy=W/nvj kk=0.0D0 residuo=0.0D0 RMS=0.0D0 tol=9.0D-03 write(*,*) 'Tempo final' read(*,*) tempof !Malha x(1)=0.5D0*dx do i=2,nvi x(i)=x(i-1)+dx end do y(1)=0.5D0*dy do j=2, nvj y(j)=y(j-1)+dy end do !Campo inicial de T do j=1, nvj do i=1, nvi T(i,j)=T_ini Tp_0(i,j)=T_ini end do end do tempo=0 do while (tempo1).and.(j1).and.(i1).and.(j1).and.(j1).and.(j1).and.(j1).and.(i1).and.(i1).and.(i1).and.(j1).and.(i1).and.(i