  We propose a new finite volume renormalization scheme. Our scheme is based on
the Gradient Flow applied to both fermion and gauge fields and, much like the
Schr\"odinger functional method, allows for a nonperturbative determination of
the scale dependence of operators using a step-scaling approach. We give some
preliminary results for the pseudo-scalar density in the quenched
approximation.
