anchor.diff:  
       endif
        If (associated(input%properties%exccplot%plot2d)) then
  -          labels=>create_plotlablels("Potential","VCL2d",2)
  +          labels=>create_plotlablels("Potential","VCL",2)
   		 call set_plotlabel_axis(labels,1,"a","lattice coordinate")
   		 call set_plotlabel_axis(labels,2,"b","lattice coordinate")
   		 call set_plotlabel_axis(labels,3,"Potential","E_h/(ea_0)")
  
positive.code:      endif
        If (associated(input%properties%exccplot%plot3d)) then
  
           labels=>create_plotlablels("Potential","VCL3d",3)
  		 call set_plotlabel_axis(labels,1,"a","lattice coordinate")
  		 call set_plotlabel_axis(labels,2,"b","lattice coordinate")
  		  call set_plotlabel_axis(labels,3,"c","lattice coordinate")
  
negative.code:         & vclir,input%properties%exccplot%plot3d)
           call destroy_plotlablels(labels)
  
             labels=>create_plotlablels("Potential","VXC",3)
  		 call set_plotlabel_axis(labels,1,"a","lattice coordinate")
  		 call set_plotlabel_axis(labels,2,"b","lattice coordinate")
  		  call set_plotlabel_axis(labels,3,"c","lattice coordinate")
  		 call set_plotlabel_axis(labels,4,"Potential","E_h/(ea_0)")
           Call plot3d (labels, 1, input%groundstate%lmaxvr, lmmaxvr, vxcmt, &
          & vxcir,input%properties%exccplot%plot3d)
          call destroy_plotlablels(labels)
           Write (*,*)
           Write (*, '("Info(potplot):")')
           Write (*, '(" 3D Coulomb potential plot written to VCL3d.xml")&
  