added [...,drop=F] for calculating Pressures with single region