I have incorporated the dispersion effect for the OPBE functional by specifying the S6, SR6 and S8 parameters for GD3 in Gaussian G09 package. For this purpose, I have defined the environment variables in G09 as mentioned below:
export GAUSS_DFTD3_S6=1000000
export GAUSS_DFTD3_SR6=1217000
export GAUSS_DFTD3_S8=722000
This works for the Gaussian09 package. However, the same does not work for Gaussian16 even though the environment variables were set similar to that for Gaussian09 (as shown above). It is unable to read the S8 parameter in G16 package. Kindly help me with the use of environment variable for incorporating the GD3 dispersion in G16 for performing calculations using the OPBE functional.