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:          SELECT CASE (e%type)
             CASE (EHYBRIDTHPL)
               e%ncgs=ncgs; e%bcgs=bcgs; e%tcgs=tcgs; e%ncgsnth=ncgsnth 
             CASE (EPOLSYNCHTH)
               e%ncgs=ncgs; e%bcgs=bcgs; e%tcgs=tcgs; e%ncgsnth=ncgsnth
             CASE (EPOLSYNCHPL)
  
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)
  