global_error.f90

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

! global_error.f90
! 
! Speicherung der E_z Feldkomponente fuer die Berechenung des globalen Fehlers
!
!    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: 26-07-2007 08:40:36 PM CEST

MODULE global_error

USE fdtd_gitter, ONLY: gitter
USE plot, ONLY: toch

CONTAINS

SUBROUTINE save_ez(do_global_error, g, time_step, limit)

  TYPE(gitter), INTENT(IN)                  :: g
  INTEGER, INTENT(IN)                       :: time_step
  INTEGER, INTENT(IN)                       :: limit
  INTEGER, INTENT(IN)                       :: do_global_error
  CHARACTER(5)                              :: zahl
  CHARACTER, DIMENSION(1:5)                 :: x
  INTEGER                                   :: i, j, k

  IF (do_global_error .EQ. 1) THEN
    CALL toch(x, time_step)
    zahl = x(1)//x(2)//x(3)//x(4)//x(5)

    OPEN(10,FILE='cube_'//g%id//'_'//zahl//'.data', ACTION='WRITE')
    WRITE(10,*) limit
    WRITE(10,*)
    DO k = -limit, limit-1, 1
      DO j = -limit, limit, 1
        DO i = -limit, limit, 1
          WRITE(10,*) i, j, k, g%E(i, j, k, 3)
        ENDDO
      ENDDO
    ENDDO
    CLOSE(10)
  ENDIF

END SUBROUTINE save_ez

END MODULE global_error