mg_get_face_coords Subroutine

  subroutine mg_get_face_coords(box, nb, nc, x)
    type(mg_box_t), intent(in) :: box
    integer, intent(in)        :: nb
    integer, intent(in)        :: nc
    real(dp), intent(out)      :: x(1)
    integer                    :: nb_dim
    real(dp)                   :: rmin(1)

    nb_dim = mg_neighb_dim(nb)


    rmin = box%r_min
    if (.not. mg_neighb_low(nb)) then
       rmin(nb_dim) = rmin(nb_dim) + box%dr(nb_dim) * nc
    end if

    x(1) = rmin(1)
  end subroutine mg_get_face_coords