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:          Call plot2d (labels, 1, input%groundstate%lmaxvr, lmmaxvr, vclmt, &
          & vclir,input%properties%exccplot%plot2d)
          call destroy_plotlablels(labels)
           labels=>create_plotlablels("Potential","VXC2d",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)")
  
negative.code:      endif
        If (associated(input%properties%exccplot%plot3d)) then
  
           labels=>create_plotlablels("Potential","VCL",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")
  