liao.f90

Sfdtd Simple finite-difference time-domain

  TYPE(rand), INTENT(IN), DIMENSION(1:2)                   :: un
  INTEGER                                                  :: ix, iy, iz

  ! x Ebene
  DO iy = g%nyl, g%nyyh, 1
    DO iz = g%nzl+1, g%nzyh, 1
      ! E_y
      g%E(g%nxl, iy, iz, 2) = 2.0d0 * zero%faceXl(iy, iz, 1) - un(2)%faceXl(iy, iz, 1)
      g%E(g%nxgh, iy, iz, 2) = 2.0d0 * zero%faceXh(iy, iz, 1) - un(2)%faceXh(iy, iz, 1)
    ENDDO
  ENDDO
  DO iy = g%nyl+1, g%nyyh, 1
    DO iz = g%nzl, g%nzyh, 1
      ! E_z
      g%E(g%nxl, iy, iz, 3) = 2.0d0 * zero%faceXl(iy, iz, 2) - un(2)%faceXl(iy, iz, 2)
      g%E(g%nxgh, iy, iz, 3) = 2.0d0 * zero%faceXh(iy, iz, 2) - un(2)%faceXh(iy, iz, 2)
    ENDDO
  ENDDO
  ! y 
  DO ix = g%nxl, g%nxyh, 1
    DO iz = g%nzl+1, g%nzyh, 1
      ! E_x
      g%E(ix, g%nyl, iz, 1) = 2.0d0 * zero%faceYl(ix, iz, 1) - un(2)%faceYl(ix, iz, 1)
      g%E(ix, g%nygh, iz, 1) = 2.0d0 * zero%faceYh(ix, iz, 1) - un(2)%faceYh(ix, iz, 1)
    ENDDO
  ENDDO
  DO ix = g%nxl, g%nxgh, 1
    DO iz = g%nzl, g%nzyh, 1
      ! E_z
      g%E(ix, g%nyl, iz, 3) = 2.0d0 * zero%faceYl(ix, iz, 2) - un(2)%faceYl(ix, iz, 2)
      g%E(ix, g%nygh, iz, 3) = 2.0d0 * zero%faceYh(ix, iz, 2) - un(2)%faceYh(ix, iz, 2)
    ENDDO
  ENDDO
  ! und z 
  DO ix = g%nxl, g%nxyh, 1
    DO iy = g%nyl, g%nygh, 1
      ! E_x
      g%E(ix, iy, g%nzl, 1) = 2.0d0 * zero%faceZl(ix, iy, 1) - un(2)%faceZl(ix, iy, 1)
      g%E(ix, iy, g%nzgh, 1) = 2.0d0 * zero%faceZh(ix, iy, 1) - un(2)%faceZh(ix, iy, 1)
    ENDDO
  ENDDO
  DO ix = g%nxl, g%nxgh, 1
    DO iy = g%nyl, g%nyyh, 1
      ! E_y 
      g%E(ix, iy, g%nzl, 2) = 2.0d0 * zero%faceZl(ix, iy, 2) - un(2)%faceZl(ix, iy, 2)
      g%E(ix, iy, g%nzgh, 2) = 2.0d0 * zero%faceZh(ix, iy, 2) - un(2)%faceZh(ix, iy, 2)
    ENDDO
  ENDDO

END SUBROUTINE add_liao_second_order

! Liao ABC 3. Ordnung ohne Interpolation (Speicherung 0)
SUBROUTINE store_liao_zero(g, zero)

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


  ! x Ebene
  DO iy = g%nyl, g%nyyh, 1
    DO iz = g%nzl+1, g%nzyh, 1
    ! E_y
      zero%faceXl(iy, iz, 1) = g%E(g%nxl+1, iy, iz, 2)
      zero%faceXh(iy, iz, 1) = g%E(g%nxyh, iy, iz, 2)
    ENDDO
  ENDDO
  DO iy = g%nyl+1, g%nyyh, 1
    DO iz = g%nzl, g%nzyh, 1
    ! E_z
      zero%faceXl(iy, iz, 2) = g%E(g%nxl+1, iy, iz, 3)
      zero%faceXh(iy, iz, 2) = g%E(g%nxyh, iy, iz, 3)
    ENDDO
  ENDDO
  ! y 
  DO ix = g%nxl, g%nxyh, 1
    DO iz = g%nzl+1, g%nzyh, 1
      ! E_x
      zero%faceYl(ix, iz, 1) = g%E(ix, g%nyl + 1, iz, 1)
      zero%faceYh(ix, iz, 1) = g%E(ix, g%nyyh, iz, 1)
    ENDDO
  ENDDO
  DO ix = g%nxl, g%nxgh, 1
    DO iz = g%nzl, g%nzyh, 1
      ! E_z
      zero%faceYl(ix, iz, 2) = g%E(ix, g%nyl + 1, iz, 3)
      zero%faceYh(ix, iz, 2) = g%E(ix, g%nyyh, iz, 3)
    ENDDO
  ENDDO
  ! und z 
  DO ix = g%nxl, g%nxyh, 1
    DO iy = g%nyl, g%nygh, 1
      ! E_x
      zero%faceZl(ix, iy, 1) = g%E(ix, iy, g%nzl + 1, 1)
      zero%faceZh(ix, iy, 1) = g%E(ix, iy, g%nzyh, 1)
    ENDDO
  ENDDO
  DO ix = g%nxl, g%nxgh, 1
    DO iy = g%nyl, g%nyyh, 1
      ! E_y
      zero%faceZl(ix, iy, 2) = g%E(ix, iy, g%nzl + 1, 2)
      zero%faceZh(ix, iy, 2) = g%E(ix, iy, g%nzyh, 2)
    ENDDO
  ENDDO

END SUBROUTINE store_liao_zero


! Liao ABC 3. Ordnung ohne Interpolation (Speicherung 1)
SUBROUTINE store_liao_one(g, un)

  TYPE(gitter), INTENT(IN)                                 :: g
  TYPE(rand), INTENT(INOUT), DIMENSION(1:2)                :: un
  INTEGER                                                  :: ix, iy, iz

  un(2)%faceXl(:,:,:) = un(1)%faceXl(:,:,:)
  un(2)%faceXh(:,:,:) = un(1)%faceXh(:,:,:)
  un(2)%faceYl(:,:,:) = un(1)%faceYl(:,:,:)
  un(2)%faceYh(:,:,:) = un(1)%faceYh(:,:,:)
  un(2)%faceZl(:,:,:) = un(1)%faceZl(:,:,:)
  un(2)%faceZh(:,:,:) = un(1)%faceZh(:,:,:)

  ! x Ebene
  DO iy = g%nyl, g%nyyh, 1
    DO iz = g%nzl+1, g%nzyh, 1
   ! E_y
      un(1)%faceXl(iy, iz, 1) = g%E(g%nxl+2, iy, iz, 2)
      un(1)%faceXh(iy, iz, 1) = g%E(g%nxyh-1, iy, iz, 2)
    ENDDO
  ENDDO
  DO iy = g%nyl+1, g%nyyh, 1
    DO iz = g%nzl, g%nzyh, 1
      ! E_z
      un(1)%faceXl(iy, iz, 2) = g%E(g%nxl+2, iy, iz, 3)
      un(1)%faceXh(iy, iz, 2) = g%E(g%nxyh-1, iy, iz, 3)
    ENDDO
  ENDDO
  ! y 
  DO ix = g%nxl, g%nxyh, 1
    DO iz = g%nzl+1, g%nzyh, 1
      ! E_x 
      un(1)%faceYl(ix, iz, 1) = g%E(ix, g%nyl + 2, iz, 1)
      un(1)%faceYh(ix, iz, 1) = g%E(ix, g%nyyh - 1, iz, 1)
    ENDDO
  ENDDO
  DO ix = g%nxl, g%nxgh, 1
    DO iz = g%nzl, g%nzyh, 1
      ! E_z
      un(1)%faceYl(ix, iz, 2) = g%E(ix, g%nyl + 2, iz, 3)
      un(1)%faceYh(ix, iz, 2) = g%E(ix, g%nyyh - 1, iz, 3)
    ENDDO
  ENDDO
  ! und z 
  DO ix = g%nxl, g%nxyh, 1
    DO iy = g%nyl, g%nygh, 1
      ! E_x 
      un(1)%faceZl(ix, iy, 1) = g%E(ix, iy, g%nzl + 2, 1)
      un(1)%faceZh(ix, iy, 1) = g%E(ix, iy, g%nzyh - 1, 1)
    ENDDO
  ENDDO
  DO ix = g%nxl, g%nxgh, 1
    DO iy = g%nyl, g%nyyh, 1
      ! E_y
      un(1)%faceZl(ix, iy, 2) = g%E(ix, iy, g%nzl + 2, 2)
      un(1)%faceZh(ix, iy, 2) = g%E(ix, iy, g%nzyh - 1, 2)
    ENDDO
  ENDDO

END SUBROUTINE store_liao_one

! Liao ABC 3. Ordnung ohne Interpolation (Speicherung 2)
SUBROUTINE store_liao_second(g, deux)

  TYPE(gitter), INTENT(IN)                                 :: g
  TYPE(rand), INTENT(INOUT), DIMENSION(1:3)                :: deux
  INTEGER                                                  :: ix, iy, iz, i

  DO i= 2, 1, -1
    deux(i+1)%faceXl(:,:,:) = deux(i)%faceXl(:,:,:)
    deux(i+1)%faceXh(:,:,:) = deux(i)%faceXh(:,:,:)
    deux(i+1)%faceYl(:,:,:) = deux(i)%faceYl(:,:,:)
    deux(i+1)%faceYh(:,:,:) = deux(i)%faceYh(:,:,:)
    deux(i+1)%faceZl(:,:,:) = deux(i)%faceZl(:,:,:)
    deux(i+1)%faceZh(:,:,:) = deux(i)%faceZh(:,:,:)
  ENDDO

  ! x Ebene
  DO iy = g%nyl, g%nyyh, 1
    DO iz = g%nzl+1, g%nzyh, 1
    ! E_y 
      deux(1)%faceXl(iy, iz, 1) = g%E(g%nxl+3, iy, iz, 2)
      deux(1)%faceXh(iy, iz, 1) = g%E(g%nxyh-2, iy, iz, 2)
    ENDDO
  ENDDO
  DO iy = g%nyl+1, g%nyyh, 1
    DO iz = g%nzl, g%nzyh, 1
      ! E_z
      deux(1)%faceXl(iy, iz, 2) = g%E(g%nxl+3, iy, iz, 3)
      deux(1)%faceXh(iy, iz, 2) = g%E(g%nxyh-2, iy, iz, 3)
    ENDDO
  ENDDO
  ! y 
  DO ix = g%nxl, g%nxyh, 1
    DO iz = g%nzl+1, g%nzyh, 1
     ! E_x 
      deux(1)%faceYl(ix, iz, 1) = g%E(ix, g%nyl + 3, iz, 1)
      deux(1)%faceYh(ix, iz, 1) = g%E(ix, g%nyyh - 2, iz, 1)
    ENDDO
  ENDDO
  DO ix = g%nxl, g%nxgh, 1
    DO iz = g%nzl, g%nzyh, 1
   ! E_z
      deux(1)%faceYl(ix, iz, 2) = g%E(ix, g%nyl + 3, iz, 3)
      deux(1)%faceYh(ix, iz, 2) = g%E(ix, g%nyyh - 2, iz, 3)
     ENDDO
  ENDDO
  ! und z 
  DO ix = g%nxl, g%nxyh, 1
    DO iy = g%nyl, g%nygh, 1
      ! E_x 
      deux(1)%faceZl(ix, iy, 1) = g%E(ix, iy, g%nzl + 3, 1)
      deux(1)%faceZh(ix, iy, 1) = g%E(ix, iy, g%nzyh - 2, 1)
    ENDDO
  ENDDO
  DO ix = g%nxl, g%nxgh, 1
    DO iy = g%nyl, g%nyyh, 1
     ! E_y
      deux(1)%faceZl(ix, iy, 2) = g%E(ix, iy, g%nzl + 3, 2)
      deux(1)%faceZh(ix, iy, 2) = g%E(ix, iy, g%nzyh - 2, 2)
    ENDDO
  ENDDO

END SUBROUTINE store_liao_second

! Third order Liao ABC without interpolation (Zufuegung)
SUBROUTINE add_liao_no_interpolation(g, zero, un, deux)

  TYPE(gitter), INTENT(INOUT)                              :: g
  TYPE(rand), INTENT(IN)                                   :: zero
  TYPE(rand), INTENT(IN), DIMENSION(1:2)                   :: un
  TYPE(rand), INTENT(IN), DIMENSION(1:3)                   :: deux
  INTEGER                                                  :: ix, iy, iz

  ! x Ebene
  DO iy = g%nyl, g%nyyh, 1
    DO iz = g%nzl+1, g%nzyh, 1
      ! E_y 
      g%E(g%nxl, iy, iz, 2) = 3.0d0 * (zero%faceXl(iy, iz, 1) - un(2)%faceXl(iy, iz, 1)) &
                                + deux(3)%faceXl(iy, iz, 1)
      g%E(g%nxgh, iy, iz, 2) = 3.0d0 * (zero%faceXh(iy, iz, 1) - un(2)%faceXh(iy, iz, 1)) &
                                + deux(3)%faceXh(iy, iz, 1)
    ENDDO
  ENDDO
  DO iy = g%nyl+1, g%nyyh, 1
    DO iz = g%nzl, g%nzyh, 1
      ! E_z
      g%E(g%nxl, iy, iz, 3) = 3.0d0 * (zero%faceXl(iy, iz, 2) - un(2)%faceXl(iy, iz, 2)) &
                                + deux(3)%faceXl(iy, iz, 2)
      g%E(g%nxgh, iy, iz, 3) = 3.0d0 * (zero%faceXh(iy, iz, 2) - un(2)%faceXh(iy, iz, 2)) &
                                + deux(3)%faceXh(iy, iz, 2)
    ENDDO
  ENDDO
  ! y 
  DO ix = g%nxl, g%nxyh, 1
    DO iz = g%nzl+1, g%nzyh, 1
      ! E_x 
      g%E(ix, g%nyl, iz, 1) = 3.0d0 * (zero%faceYl(ix, iz, 1) - un(2)%faceYl(ix, iz, 1)) &
                                + deux(3)%faceYl(ix, iz, 1)
      g%E(ix, g%nygh, iz, 1) = 3.0d0 * (zero%faceYh(ix, iz, 1) - un(2)%faceYh(ix, iz, 1)) &
                                + deux(3)%faceYh(ix, iz, 1)
    ENDDO
  ENDDO
  DO ix = g%nxl, g%nxgh, 1
    DO iz = g%nzl, g%nzyh, 1
      ! E_z
      g%E(ix, g%nyl, iz, 3) = 3.0d0 * (zero%faceYl(ix, iz, 2) - un(2)%faceYl(ix, iz, 2)) &
                                + deux(3)%faceYl(ix, iz, 2)
      g%E(ix, g%nygh, iz, 3) = 3.0d0 * (zero%faceYh(ix, iz, 2) - un(2)%faceYh(ix, iz, 2)) &
                                + deux(3)%faceYh(ix, iz, 2)

    ENDDO
  ENDDO
  ! und z 
  DO ix = g%nxl, g%nxyh, 1
    DO iy = g%nyl, g%nygh, 1
      ! E_x 
      g%E(ix, iy, g%nzl, 1) = 3.0d0 * (zero%faceZl(ix, iy, 1) - un(2)%faceZl(ix, iy, 1)) &
                                + deux(3)%faceZl(ix, iy, 1)
      g%E(ix, iy, g%nzgh, 1) = 3.0d0 * (zero%faceZh(ix, iy, 1) - un(2)%faceZh(ix, iy, 1)) &
                                + deux(3)%faceZh(ix, iy, 1)

    ENDDO
  ENDDO
  DO ix = g%nxl, g%nxgh, 1
    DO iy = g%nyl, g%nyyh, 1
     !  E_y 
      g%E(ix, iy, g%nzl, 2) = 3.0d0 * (zero%faceZl(ix, iy, 2) - un(2)%faceZl(ix, iy, 2)) &
                                + deux(3)%faceZl(ix, iy, 2)
      g%E(ix, iy, g%nzgh, 2) = 3.0d0 * (zero%faceZh(ix, iy, 2) - un(2)%faceZh(ix, iy, 2)) &
                                + deux(3)%faceZh(ix, iy, 2)

    ENDDO
  ENDDO
  
END SUBROUTINE add_liao_no_interpolation

END MODULE liao