In this line https://github.com/usnistgov/jarvis/blob/master/jarvis/io/vasp/outputs.py#L208 tmp_set should be reshaped as ng[2],ng[0],ng[1] to follow Fortran writing format and get correct average electrostatic potential etc.