With this option the code computes the Berry phase polarization in a grid of geometries, the geometries specified as in what='mur_lc'. The results are written in different files in the directory elastic_constants. With a subsequent calculation that use what='mur_lc' or what='mur_lc_t' the code can read the file produced with this option and interpolate the polarization at the pressure (and temperature if what='mur_lc_t') dependent geometries using polynomials. The calculation of polarization is controled by the same variables used for the option what='scf_polarization'. The variables that control this option are:
fl_polar : the name of the file that contains the polarization.
Default : 'output_polar.dat'