anchor.diff:           character(len=20), intent(in) :: ename
            if(ename=='POLSYNCHTH') then
               e%type=EPOLSYNCHTH
  +            e%neq=4
  +         if(ename=='MAXJUTT') then
  +            e%type=EMAXJUTT
               e%neq=4
            elseif(ename.eq.'HYBRIDTH') then
               e%type=EHYBRIDTH 
  
positive.code:              deallocate(e%bcgs); deallocate(e%incang)
               deallocate(e%ncgsnth); deallocate(e%p)
               call del_polsynchpl(e%neq)              
             CASE (EPOLSYNCHTH)
               deallocate(e%tcgs); deallocate(e%ncgs)
               deallocate(e%bcgs); deallocate(e%incang)
  
negative.code: 
           integer :: i
           select  case(e%type)
             case(emaxjutt)
  !alwinnote 2015/03/05 Below is just copied from ehybridthpl not currently functional
  !need to ensure unit conversions are safe before continuing
  !convert theta_min in units of theta = kT/m_e c^2 to tcgs units
  !nu is just a frequency
  !ncgs just gets multiplied by dimensionless weights
  !bcgs is just an input
                call polsynchth(nu,e%ncgs,e%bcgs,e%tcgs,e%incang,Kth)
                call polsynchpl(nu,e%ncgsnth,e%bcgs,e%incang,e%p,e%gmin, &
                 e%gmax,Kpl)
                K = Kth + Kpl
             case(ehybridthpl)
                call polsynchth(nu,e%ncgs,e%bcgs,e%tcgs,e%incang,Kth)
                call polsynchpl(nu,e%ncgsnth,e%bcgs,e%incang,e%p,e%gmin, &
  