Hi Evgeniy,
technically, this is not a problem.
Just set up your HF calculation and switch in dispersion in define (or add '$disp3 bj' to control).
In the output you can see that the HF parameter are in such lines:
"
DFT-D V3(BJ)
DF hf
parameters
s6 : 1.0000
s8 : 0.9171
a1 : 0.3385
a2 : 2.8830 "
Since TM 7.0, you can also use the HF-3c approach (parameterized in connection with the minix basis set) by Sure and Grimme, if you add this control:
$disp3 bj func hf-3c
Cheers,
Arnim