Dear leszec_malec,

You can bypass the error by using a functional parameterized for DFT-D3 (BJ), for example PBE, and then manually defining the dispersion parameters.

Please see the DFTD3 block below:

DFTD3
VERSION
4
FUNC
PBE
S6
1.0000
S8
1.6112
A1
0.0000
A2
7.3539
END

Best regards,

Lorenzo