anchor.diff:           type (emis), intent(out) :: e
            character(len=20), intent(in) :: ename
            if(ename=='POLSYNCHTH') then
  -           e%type=EPOLSYNCHTH
  -           e%neq=4
  +            e%type=EPOLSYNCHTH
  +            e%neq=4
  +         elseif(ename.eq.'HYBRIDTH') then
  +            e%type=EHYBRIDTH 
  +            e%neq=4
  +         elseif(ename.eq.'HYBRIDPL') then
  +            e%type=EHYBRIDPL 
  +            e%neq=4
            elseif(ename=='HYBRIDTHPL') then
  -            e%type=EHYBRIDTHPL
  +            e%type=EHYBRIDTHPL 
               e%neq=4
            elseif(ename=='SYNCHTHAV') then
               e%type=ESYNCHTHAV
  
positive.code:            CASE (EHYBRIDTHPL)
               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)
  
negative.code:        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, &
                    EHYBRIDTH=21, EHYBRIDPL=22
  
         type emis
           double precision, dimension(:,:), allocatable :: j,K
  