save_binary.f

来自「xy时域有限差分的fortran数值模拟 xy时域有限差分的fortran数值」· F 代码 · 共 19 行

*--------------------------------------------------------------------*
      subroutine save_binary
*
*  Write E-Field to a binary file
*
      include 'fdtmcom.f'
      character*1 saveout(nx,ny)
      call scal_vid
      do i=1, nx
         do j=1,ny
            saveout(i,j)=char(int(ezs(i,j)))
         end do
      end do
      open(8,file=file_out,status='unknown',form='unformatted')
      write(8)saveout
      close(8)
      return
      end