anchor.diff:    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
  +
  +!$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) &
  +!$    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
         RR(1) = +(-1.00d0/HY_m**2) * Ac_m(1, ix_m, iy_m-1) &
  
positive.code:   integer :: ix_m,iy_m
    ! calculate the Ohmic losses in the media
    ! (written by M.Uemoto on 2016-11-22)
    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
    return
  end subroutine calc_energy_joule
  !===========================================================
  
negative.code:       Bmag(:,ix_m,iy_m) = Rc(:) * c_light
      end do
    end do
  !$omp end parallel do
    return
  end subroutine calc_bmag_field_2d
  !===========================================================
  