anchor.diff:  getstructgroundstate%mixernumber=stringtonumbermixer(getstructgroundstate%mixer)
   
   nullify(np)  
  +np=>getAttributeNode(thisnode,"beta0")
  +getstructgroundstate%beta0=0.4
  +if(associated(np)) then
  +       call extractDataAttribute(thisnode,"beta0",getstructgroundstate%beta0)
  +       call removeAttribute(thisnode,"beta0")      
  +endif
  +
  +nullify(np)  
  +np=>getAttributeNode(thisnode,"betainc")
  +getstructgroundstate%betainc=1.1
  +if(associated(np)) then
  +       call extractDataAttribute(thisnode,"betainc",getstructgroundstate%betainc)
  +       call removeAttribute(thisnode,"betainc")      
  +endif
  +
  +nullify(np)  
  +np=>getAttributeNode(thisnode,"betadec")
  +getstructgroundstate%betadec=0.6
  +if(associated(np)) then
  +       call extractDataAttribute(thisnode,"betadec",getstructgroundstate%betadec)
  +       call removeAttribute(thisnode,"betadec")      
  +endif
  +
  +nullify(np)  
   np=>getAttributeNode(thisnode,"fromscratch")
   getstructgroundstate%fromscratch= .true.
   if(associated(np)) then
  
positive.code: implicit none
  type(Node),pointer::thisnode
  type(eliashberg_type),pointer::getstructeliashberg
  		
  integer::len=1,i=0
  allocate(getstructeliashberg)  
  #ifdef INPUTDEBUG      
        write(*,*)"we are at eliashberg"
  #endif
        
        i=0
        len=0
        call  handleunknownnodes(thisnode)
  
negative.code:        call removeAttribute(thisnode,"epsarpack")      
  endif
  
  nullify(np)  
  np=>getAttributeNode(thisnode,"evaltol")
  getstructsolver%evaltol=1e-8
  if(associated(np)) then
         call extractDataAttribute(thisnode,"evaltol",getstructsolver%evaltol)
         call removeAttribute(thisnode,"evaltol")      
  endif
  
        i=0
        len=0
        call  handleunknownnodes(thisnode)
  end function
  
  function getstructOEP(thisnode)
  
  implicit none
  type(Node),pointer::thisnode
  type(OEP_type),pointer::getstructOEP
  		type(Node),pointer::np
  
  
  integer::len=1,i=0
  allocate(getstructOEP)  
  #ifdef INPUTDEBUG      
        write(*,*)"we are at OEP"
  #endif
        
  nullify(np)  
  np=>getAttributeNode(thisnode,"maxitoep")
  getstructOEP%maxitoep=120
  if(associated(np)) then
         call extractDataAttribute(thisnode,"maxitoep",getstructOEP%maxitoep)
         call removeAttribute(thisnode,"maxitoep")      
  endif
  
  nullify(np)  
  np=>getAttributeNode(thisnode,"tauoep")
  getstructOEP%tauoep=(/1.0,0.2,1.5/)
  if(associated(np)) then
         call extractDataAttribute(thisnode,"tauoep",getstructOEP%tauoep)
         call removeAttribute(thisnode,"tauoep")      
  endif
  
        i=0
        len=0
        call  handleunknownnodes(thisnode)
  end function
  
  function getstructRDMFT(thisnode)
  
  implicit none
  type(Node),pointer::thisnode
  type(RDMFT_type),pointer::getstructRDMFT
  		type(Node),pointer::np
  
  
  integer::len=1,i=0
  allocate(getstructRDMFT)  
  #ifdef INPUTDEBUG      
        write(*,*)"we are at RDMFT"
  #endif
        
  nullify(np)  
  np=>getAttributeNode(thisnode,"rdmxctype")
  getstructRDMFT%rdmxctype=2
  if(associated(np)) then
         call extractDataAttribute(thisnode,"rdmxctype",getstructRDMFT%rdmxctype)
         call removeAttribute(thisnode,"rdmxctype")      
  endif
  
  nullify(np)  
  np=>getAttributeNode(thisnode,"rdmmaxscl")
  getstructRDMFT%rdmmaxscl=1
  if(associated(np)) then
         call extractDataAttribute(thisnode,"rdmmaxscl",getstructRDMFT%rdmmaxscl)
         call removeAttribute(thisnode,"rdmmaxscl")      
  endif
  
  nullify(np)  
  np=>getAttributeNode(thisnode,"maxitn")
  getstructRDMFT%maxitn=250
  if(associated(np)) then
         call extractDataAttribute(thisnode,"maxitn",getstructRDMFT%maxitn)
         call removeAttribute(thisnode,"maxitn")      
  endif
  
  nullify(np)  
  np=>getAttributeNode(thisnode,"maxitc")
  getstructRDMFT%maxitc=10
  if(associated(np)) then
         call extractDataAttribute(thisnode,"maxitc",getstructRDMFT%maxitc)
         call removeAttribute(thisnode,"maxitc")      
  endif
  
  nullify(np)  
  np=>getAttributeNode(thisnode,"taurdmn")
  getstructRDMFT%taurdmn=1.0
  if(associated(np)) then
         call extractDataAttribute(thisnode,"taurdmn",getstructRDMFT%taurdmn)
         call removeAttribute(thisnode,"taurdmn")      
  endif
  
  nullify(np)  
  np=>getAttributeNode(thisnode,"taurdmc")
  getstructRDMFT%taurdmc=0.5
  if(associated(np)) then
         call extractDataAttribute(thisnode,"taurdmc",getstructRDMFT%taurdmc)
         call removeAttribute(thisnode,"taurdmc")      
  endif
  
  nullify(np)  
  np=>getAttributeNode(thisnode,"rdmalpha")
  getstructRDMFT%rdmalpha=0.7
  if(associated(np)) then
         call extractDataAttribute(thisnode,"rdmalpha",getstructRDMFT%rdmalpha)
         call removeAttribute(thisnode,"rdmalpha")      
  endif
  
  nullify(np)  
  np=>getAttributeNode(thisnode,"rdmtemp")
  getstructRDMFT%rdmtemp=0.0
  if(associated(np)) then
         call extractDataAttribute(thisnode,"rdmtemp",getstructRDMFT%rdmtemp)
         call removeAttribute(thisnode,"rdmtemp")      
  endif
  
        i=0
        len=0
        call  handleunknownnodes(thisnode)
  