extrapol.f

来自「二维热能方程的求解」· F 代码 · 共 33 行

!=======================================================================
      subroutine extrapol(nt,nx,ny,idsol,level,ht,tn,hx,hy,
     &               scale,x,y,z,ierr)
!=======================================================================
      implicit none
      integer nt,nx,ny,idsol,level,ierr
      real*8 ht,tn,hx,hy,scale
      real*8 x(nx,ny),y(nx,ny),z(nx,ny)

      integer ix,iy
      real*8 ftncu,t,u,fac

      if(level.ge.1) print'("EXTRAPLO: idsol=",i2)',idsol
      if(level.ge.2) then
         print*,"nt=",nt," nx=",nx," ny=",ny
         print*,"ht=",ht," tn=",tn," hx=",hx," hy=",hy
      endif


      t=tn+0.5d0*ht
      fac=scale*ht

      do iy=1,ny
      do ix=1,nx
         u=1.5d0*y(ix,iy)-0.5d0*z(ix,iy)
         x(ix,iy)=x(ix,iy)+fac*ftncu(t,u)
      enddo
      enddo


      return
      end