anchor.diff:         integer,parameter :: ELAMBDA=0,ESYNCHTH=1,EPOLSYNCHTH=2, &
                     ESYNCHPL=3,EPOLSYNCHPL=4,EBREMS=5,ELINE=6, &
                     EIRON=7,EBB=8,EBBPOL=9,EINTERP=10,EFBB=11,ERHO=12, &
  -                  ESYNCHTHAV=13, ESYNCHTHAVNOABS=14, EHYBRIDTHPL = 20
  +                  ESYNCHTHAV=13, ESYNCHTHAVNOABS=14, EHYBRIDTHPL = 20, &
  +                  EHYBRIDTH=21, EHYBRIDPL=22
   
          type emis
            double precision, dimension(:,:), allocatable :: j,K
  
positive.code:              allocate(e%bcgs(npts))
               allocate(e%ncgsnth(npts)); allocate(e%p(npts))
               call initialize_polsynchpl(e%neq) !I guess this is right?
             CASE (EPOLSYNCHTH)
               allocate(e%tcgs(npts)); allocate(e%ncgs(npts))
               allocate(e%bcgs(npts))
  
negative.code:              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
               K(:,5:11) = Kth(:,5:11) + Kpl(:,5:11)
  
             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
               K(:,5:11) = Kth(:,5:11) + Kpl(:,5:11)
  
             case(epolsynchth)
                write(6,*) 'polsynchth nu: ',nu
  