fdtd for inner grid pointd and pml calculation subroutines