anchor.diff:                    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
  +                  EHYBRIDTH=21, EHYBRIDPL=22, EMAXJUTT=23
   
          type emis
            real(kind=8), dimension(:,:), allocatable :: j,K
  
positive.code: 
           integer :: i
           select  case(e%type)
             case(ehybridthpl)
                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, &
  
negative.code:              e%ncgs=ncgs; e%bcgs=bcgs; e%tcgs=tcgs; e%ncgsnth=ncgsnth 
             CASE (EHYBRIDPL)
               e%ncgs=ncgs; e%bcgs=bcgs; e%tcgs=tcgs; e%ncgsnth=ncgsnth 
             CASE (EMAXJUTT) !alwinnote 2015/03/05
               e%ncgs=ncgs; e%bcgs=bcgs; e%tcgs=tcgs
             CASE (EPOLSYNCHTH)
               e%ncgs=ncgs; e%bcgs=bcgs; e%tcgs=tcgs
             CASE (EPOLSYNCHPL)
  