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:               call polsynchpl(nu,e%ncgsnth,e%bcgs,e%incang,e%p,e%gmin, &
                 e%gmax,Kpl)
                K = Kth + Kpl
             case(epolsynchth)
  !              write(6,*) 'polsynchth ncgs: ',allocated(e%ncgs),size(e%ncgs)
  !              write(6,*) 'polsynchth bcgs: ',e%bcgs
  
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)
  