conv_binary.f

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

*--------------------------------------------------------------------*
      subroutine conv_binary(s)
      include 'xfdtmcom.f'
      character*1 s(nx,ny)
      drange = emax - emin
      if (idata.eq.1) then
         do i=1,nx
            do j=1,ny
                s(i,j)=char(int(255*(ezs(i,j) - emin)/drange))
            end do
         end do
      elseif (idata.eq.2) then
         do i=1,nx
            do j=1,ny
               s(i,j)=char(int(255*(hxs(i,j) - emin)/drange))
            end do
         end do
      else
         do i=1,nx
            do j=1,ny
               s(i,j)=char(int(255*(hys(i,j) - emin)/drange))
            end do
         end do
      endif
      return
      end