pec.f90

来自「Sfdtd Simple finite-difference time-doma」· F90 代码 · 共 64 行

! pec.f90
! 
! Perfectly electric conducting boundaries
! Perfekt elektrisch leitende Raender
!
!    Copyright (C) 2007  Paul Panserrieu, < [e-mail] >
!
!    This program is free software: you can redistribute it and/or modify
!    it under the terms of the GNU General Public License as published by
!    the Free Software Foundation, either version 3 of the License.
! 
! last modified: 14-09-2007 06:03:18 PM CEST

MODULE pec

USE fdtd_gitter

IMPLICIT NONE

CONTAINS

SUBROUTINE update_pec(g)

  TYPE(gitter), INTENT(INOUT)    :: g
  INTEGER                        :: ix, iy, iz

  ! x Ebene
  DO iy = g%nyl, g%nyyh, 1
    DO iz = g%nzl+1, g%nzyh, 1
      g%E(g%nxl, iy, iz, 2) = 0.0d0; g%E(g%nxgh, iy, iz, 2) = 0.0d0
    ENDDO
  ENDDO
  DO iy = g%nyl+1, g%nygh-1, 1
    DO iz = g%nzl, g%nzyh, 1
      g%E(g%nxl, iy, iz, 3) = 0.0d0; g%E(g%nxgh, iy, iz, 3) = 0.0d0
    ENDDO
  ENDDO
  ! y 
  DO ix = g%nxl, g%nxyh, 1
    DO iz = g%nzl+1, g%nzyh, 1
      g%E(ix, g%nyl, iz, 1) = 0.0d0; g%E(ix, g%nygh, iz, 1) = 0.0d0
    ENDDO
  ENDDO
  DO ix = g%nxl, g%nxgh, 1
    DO iz = g%nzl, g%nzyh, 1
      g%E(ix, g%nyl, iz, 3) = 0.0d0; g%E(ix, g%nygh, iz, 3) = 0.0d0
    ENDDO
  ENDDO
  ! und z 
  DO ix = g%nxl, g%nxyh, 1
    DO iy = g%nyl, g%nygh, 1
      g%E(ix, iy, g%nzl, 1) = 0.0d0; g%E(ix, iy, g%nzgh, 1) = 0.0d0
    ENDDO
  ENDDO
  DO ix = g%nxl, g%nxgh, 1
    DO iy = g%nyl, g%nyyh, 1
      g%E(ix, iy, g%nzl, 2) = 0.0d0; g%E(ix, iy, g%nzgh, 2) = 0.0d0
    ENDDO
  ENDDO

END SUBROUTINE update_pec 

END MODULE pec