input paramters for the nonlinear solvation read from the INCAR file: LNSOL, P0_K, EINFTY_K, NSOL_K, T_K
 LSOL and LNSOL must be set to TRUE for nonlinear solvation. Just LSOL=.TRUE. defaults to the linear solvation

 some other minor cleanup and verifications in pot_nln_k.F: use of RGRID type for work arrays, setunb and scaling by the number of plane waves after the transform to t

 successful compilation