# Controlling False Discovery Rate in Stages

I have a three dimensional table of size $6\times6\times81$. Each cell of the table is a hypothesis test. Slicing the table on the third dimension produces $81$ sets of hypothesis tests which are independent between sets but dependent within sets. Originally I was thinking that I could just control the false discovery rate using the Benjamini-Hochberg procedure on all hypothesis tests simultaneously. Is that a reasonable way to attack this problem? My second thought is to control the false discovery rate within each slice along the third dimension of the table, then apply some other sort of correction after that. Does anyone have more information on this sort of procedure?