set_gcp_realspace_cutoff Subroutine

subroutine set_gcp_realspace_cutoff(verror, vgcp, bas, srb) &
      & bind(C, name=namespace//"set_gcp_realspace_cutoff")
   type(c_ptr), value :: verror
   type(vp_error), pointer :: error
   type(c_ptr), value :: vgcp
   type(vp_gcp), pointer :: gcp
   real(c_double), value, intent(in) :: bas
   real(c_double), value, intent(in) :: srb

   if (.not.c_associated(verror)) return
   call c_f_pointer(verror, error)

   if (.not.c_associated(vgcp)) then
      call fatal_error(error%ptr, "Counter-poise parameters are missing")
      return
   end if
   call c_f_pointer(vgcp, gcp)

   gcp%cutoff = realspace_cutoff(gcp=bas, srb=srb)
end subroutine set_gcp_realspace_cutoff