! pec.f90! ! Perfectly electric conducting boundaries! Perfekt elektrisch leitende Raender!!    Copyright (C) 2007  Paul Panserrieu, < peutetre@cs.tu-berlin.de >!!    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 CESTMODULE pecUSE fdtd_gitterIMPLICIT NONECONTAINSSUBROUTINE 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  ENDDOEND SUBROUTINE update_pec END MODULE pec