anchor.diff:    integer :: ix_m,iy_m
     ! calculate the Ohmic losses in the media
     ! (written by M.Uemoto on 2016-11-22)
  +
  +!$omp parallel do collapse(2) default(none) &
  +!$    private(iy_m,ix_m) &
  +!$    shared(NYvacB_m,NYvacT_m,NXvacL_m,NXvacR_m,energy_joule,j_m,Elec,aLxyz) &
  +!$    firstprivate(dt)
     do iy_m=NYvacB_m,NYvacT_m
       do ix_m=NXvacL_m,NXvacR_m
         energy_joule(ix_m, iy_m) = energy_joule(ix_m, iy_m) &
             & + sum(-j_m(:,ix_m,iy_m) * Elec(:,ix_m,iy_m)) * dt * aLxyz
       end do
     end do
  +!$omp end parallel do
  +
     return
   end subroutine calc_energy_joule
   !===========================================================
  
positive.code:   integer :: ix_m,iy_m
    real(8) :: RR(3) ! rot rot Ac
    ! (written by M.Uemoto on 2016-11-22)
    Ac_old_m=Ac_m
    Ac_m=Ac_new_m
    do iy_m=NYvacB_m,NYvacT_m
      do ix_m=NXvacL_m,NXvacR_m
        RR(1) = +(-1.00d0/HY_m**2) * Ac_m(1, ix_m, iy_m-1) &
  
negative.code:   integer :: ix_m, iy_m
    real(8) :: Y, RR(3) ! rot rot Ac
    ! (written by M.Uemoto on 2016-11-22)
  
  !$omp parallel do collapse(2) default(none) &
  !$    private(iy_m,ix_m) &
  !$    shared(NYvacB_m,NYvacT_m,NXvacL_m,NXvacR_m,Ac_m,Ac_old_m,Ac_new_m)
    do iy_m=NYvacB_m-1,NYvacT_m+1
      do ix_m=NXvacL_m-1,NXvacR_m+1
        Ac_old_m(:,ix_m,iy_m) = Ac_m    (:,ix_m,iy_m)
        Ac_m    (:,ix_m,iy_m) = Ac_new_m(:,ix_m,iy_m)
      end do
    end do
  !$omp end parallel do
  
  !$omp parallel do collapse(2) default(none) &
  !$    private(iy_m,ix_m,RR,Y) &
  !$    shared(NYvacB_m,NYvacT_m,NXvacL_m,NXvacR_m,HX_m,HY_m,Ac_m,j_m,Ac_old_m,Ac_new_m) &
  !$    firstprivate(dt)
    do iy_m=NYvacB_m,NYvacT_m
      do ix_m=NXvacL_m,NXvacR_m
        Y = (iy_m - 0.50d0) * HY_m
        RR(1) = +(+0.50d0*(1.00d0/HY_m)*(1.00d0/Y)-(1.00d0/HY_m**2))*Ac_m(1,ix_m+0,iy_m-1) &
              & +2.00d0*(1.00d0/HY_m**2)*Ac_m(1,ix_m+0,iy_m+0) &
              & +(-0.50d0*(1.00d0/HY_m)*(1.00d0/Y)-(1.00d0/HY_m**2))*Ac_m(1,ix_m+0,iy_m+1) &
  