scal_vid.f

来自「a 2-d FDTD program using PML」· F 代码 · 共 20 行

*---------------------------------------------------------------------*
      subroutine scal_vid
      include 'fdtmcom.f'
      emax=ezs(1,1)
      emin=ezs(1,1)
      do i=1,nx
         do j=1,ny
            if(ezs(i,j).gt.emax)emax=ezs(i,j)
            if(ezs(i,j).lt.emin)emin=ezs(i,j)
         end do
      end do
      drange=emax-emin
      do i=1,nx
         do j=1,ny
            ezs(i,j)=255*(ezs(i,j) - emin)/drange
         end do
      end do
      return
      end