anchor.diff:  endif
   
   nullify(np)  
  -np=>getAttributeNode(thisnode,"broad")
  -getstructtddft%broad=0.01d0
  -if(associated(np)) then
  -       call extractDataAttribute(thisnode,"broad",getstructtddft%broad)
  -       call removeAttribute(thisnode,"broad")      
  -endif
  -
  -nullify(np)  
   np=>getAttributeNode(thisnode,"lindhard")
   getstructtddft%lindhard= .false.
   if(associated(np)) then
  
positive.code: #endif
        
  nullify(np)  
  np=>getAttributeNode(thisnode,"dfoffdiag")
  getstructtddft%dfoffdiag= .false.
  if(associated(np)) then
         call extractDataAttribute(thisnode,"dfoffdiag",getstructtddft%dfoffdiag)
         call removeAttribute(thisnode,"dfoffdiag")      
  endif
  
  nullify(np)  
  np=>getAttributeNode(thisnode,"emattype")
  getstructtddft%emattype=0
  if(associated(np)) then
         call extractDataAttribute(thisnode,"emattype",getstructtddft%emattype)
         call removeAttribute(thisnode,"emattype")      
  endif
  
  nullify(np)  
  np=>getAttributeNode(thisnode,"lmaxapwwf")
  getstructtddft%lmaxapwwf=-1
  if(associated(np)) then
         call extractDataAttribute(thisnode,"lmaxapwwf",getstructtddft%lmaxapwwf)
         call removeAttribute(thisnode,"lmaxapwwf")      
  endif
  
  nullify(np)  
  np=>getAttributeNode(thisnode,"lmaxemat")
  getstructtddft%lmaxemat=3
  if(associated(np)) then
         call extractDataAttribute(thisnode,"lmaxemat",getstructtddft%lmaxemat)
         call removeAttribute(thisnode,"lmaxemat")      
  endif
  
  nullify(np)  
  np=>getAttributeNode(thisnode,"scissor")
  getstructtddft%scissor=0
  if(associated(np)) then
         call extractDataAttribute(thisnode,"scissor",getstructtddft%scissor)
         call removeAttribute(thisnode,"scissor")      
  endif
  
  nullify(np)  
  np=>getAttributeNode(thisnode,"optswidth")
  getstructtddft%optswidth=0
  if(associated(np)) then
         call extractDataAttribute(thisnode,"optswidth",getstructtddft%optswidth)
         call removeAttribute(thisnode,"optswidth")      
  endif
  
  nullify(np)  
  np=>getAttributeNode(thisnode,"intraband")
  getstructtddft%intraband= .false.
  if(associated(np)) then
  
negative.code:   type(plot3d_type),pointer::plot3d
  end type
  type dos_type
   real(8)::sqados(3)
   logical::lmirep
   integer::nwdos
   integer::ngrdos
  