Command for simple effects following 2 way anova with interaction

Dear All,

Suppose I do :-

head(warpbreaks)
model1<- aov(breaks ~ wool*tension,data = warpbreaks)
summary(model1)

There is significant interaction. So I need to test for simple effects of
wool at each level of  tension and vice versa. I can do a subset and then
do a one way anova for each level of tension and wool. My query is that is
there any command programmed which can do this automatically?

Best Regards,
Ashim