Add option to scale AH via a biharmonic Reynolds #

This is done via parameter RE_AH: if nonzero, the biharmonic
coefficient is scaled so that the biharmonic Reynolds number
is equal to this.