id_list = [
      'ivcurve000',
      'ivcurve004',
      'ivcurve006',
      'ivcurve010',
      'ivcurve012',
      ]
for id_item in id_list:
   # Read the DeviceConfiguration
   configuration = nlread('ivcurve_selfconsistent_configurations_ndoped_e19.nc', object_id="%s" %id_item)[0]

   # -------------------------------------------------------------
   # Electrostatic Difference Potential
   # -------------------------------------------------------------
   electrostatic_difference_potential = ElectrostaticDifferencePotential(configuration)
   nlsave('l_ndoped_e19_EDP.nc', electrostatic_difference_potential)