! global_error.f90! ! Speicherung der E_z Feldkomponente fuer die Berechenung des globalen Fehlers!!    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: 26-07-2007 08:40:36 PM CESTMODULE global_errorUSE fdtd_gitter, ONLY: gitterUSE plot, ONLY: tochCONTAINSSUBROUTINE 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)  ENDIFEND SUBROUTINE save_ezEND MODULE global_error