anchor.diff:          & vclir,input%properties%exccplot%plot3d)
            call destroy_plotlablels(labels)
   
  -           labels=>create_plotlablels("Potential","VXC3d",3)
  +           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)
  -         Close (50)
  +        call destroy_plotlablels(labels)
            Write (*,*)
            Write (*, '("Info(potplot):")')
            Write (*, '(" 3D Coulomb potential plot written to VCL3d.xml")&
  
positive.code: ! write the potential plots to file
       If (associated(input%properties%exccplot%plot1d)) then
  
          labels=>create_plotlablels("Potential","VCL1D",1)
  		 call set_plotlabel_axis(labels,1,"Distance","a_0")
  		 call set_plotlabel_axis(labels,2,"Potential","E_h/(ea_0)")
           Call plot1d (labels, 1, input%groundstate%lmaxvr, lmmaxvr, &
          & vclmt, vclir,input%properties%exccplot%plot1d)
           call destroy_plotlablels(labels)
  
    		  labels=>create_plotlablels("Potential","VXC1D",1)
  		 call set_plotlabel_axis(labels,1,"Distance","a_0")
  		 call set_plotlabel_axis(labels,2,"Exchange Correlation Potential","E_h/(ea_0)")
           Call plot1d ( labels, 1, input%groundstate%lmaxvr, lmmaxvr, &
  
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")
  