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:              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:            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)
  