function ps_HasLocalPotential(ps) result(p) type(ps_t), intent(in) :: ps logical :: p ! p = (initialized(ps%local%Vlocal%grid)) end function ps_HasLocalPotential