!=====================================================================================
!                                 program   main
!=====================================================================================
program main

use fdtd
implicit none

integer i,j



call init 
call model 
call test_model

time=0.0d0

do i=1,time_start
  call renew
end do

do i=1,(time_end-time_start),time_step
      do j=1,time_step
         call renew                         
      end do

call output
write(*,*) 'TIME STEP=',time_start+i

end do
 
end program

!------------------------------------------------------------------------------------------------
subroutine renew

use fdtd

call electric_field
!call electric_boundary_condition
call plane_wave_source

time=time+dt/2.0d0

call magnetic_field
call absorbing_boundary_condition
time=time+dt/2.0d0

return
end subroutine