get_size_BLACS Function

  function get_size_BLACS() result(numproc)
    integer :: numproc
    if(.not.blacs_status)stop "get_master_BLACS error: blacs_status=F. call init_BLACS"
    numproc = p_size
  end function get_size_BLACS