anchor.diff:                call polsynchpl(nu,e%ncgsnth,e%bcgs,e%incang,e%p,e%gmin, &
                  e%gmax,Kpl)
                 K = Kth + Kpl
  +           case(ehybridth)
  +              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 !get thermal emission
  +              K(:,5:11) = Kth(:,5:11) + Kpl(:,5:11) !and total absorption
  +           case(ehybridpl)
  +              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 = Kpl !get power law emission
  +              K(:,5:11) = Kth(:,5:11) + Kpl(:,5:11) !and total absorption
              case(epolsynchth)
   !              write(6,*) 'polsynchth ncgs: ',allocated(e%ncgs),size(e%ncgs)
   !              write(6,*) 'polsynchth bcgs: ',e%bcgs
  
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:            CASE (EHYBRIDTHPL)
               call emis_model_synchth(e,ep%mu)
               call emis_model_synchpl(e,ep%gmin,ep%gmax,ep%p1)
             CASE (EHYBRIDTH) !alwinnote
               call emis_model_synchth(e,ep%mu)
               call emis_model_synchpl(e,ep%gmin,ep%gmax,ep%p1)
             CASE (EHYBRIDPL) !alwinnote
               call emis_model_synchth(e,ep%mu)
               call emis_model_synchpl(e,ep%gmin,ep%gmax,ep%p1)
             CASE (EPOLSYNCHTH)
               call emis_model_synchth(e,ep%mu)
             CASE (EPOLSYNCHPL)
  