added calculation of surface charge density (don't added Gorshkov fortran ported code)