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:              allocate(e%bcgs(npts))
               allocate(e%ncgsnth(npts)); allocate(e%p(npts))
               call initialize_polsynchpl(e%neq) !I guess this is right?
             CASE (EMAXJUTT) !alwinnote 2015/03/05
               allocate(e%tcgs(npts)); allocate(e%ncgs(npts))
               allocate(e%bcgs(npts))
             CASE (EPOLSYNCHTH)
               allocate(e%tcgs(npts)); allocate(e%ncgs(npts))
               allocate(e%bcgs(npts))
  