Fix dangerous code
Description
The use of an explicit unit is dangerous. This is fixed using Fortran 2008 newunit feature.
This might explain some strange error like a "End of file" occurring while reading the D3 coefficients.
Checklist
-
I have checked that my code follows the Octopus coding standards -
I have added tests for all the new features added in this request.