anchor.diff:       $ddj(nlma,nnp1),
        $flmks(nlma,nnp1),dflmk(nlma,nnp1),
        $ds(nlma,nnp1),dds(nlma,nnp1)
  +      complex carg2	
   c
         dimension vr(nrp,ni),vt(nrp,ni),vp(nrp,ni),
        $dvrdr(nrp,ni),dvtdr(nrp,ni),dvpdr(nrp,ni),
  
positive.code:       call mei(enbp,enbt,apome,atome)
        enb=enbp+enbt
        ent=env+enb
  c urc: log-file with energy components written
  c
        topnuss=toplum/alum0
        botnuss=botlum/alum0
        if(mod(kstep,nlogstep).eq.0) write(15,'(f9.6,8f9.2,2f9.5)')
       $time/tscale,env/escale,envp/escale,enb/escale,enbp/escale,
       $adrke/escale,amcke/escale,apome/escale,atome/escale,
       $topnuss,botnuss
        if(nplog.gt.0) then
         if(mod(kstep,nplog).eq.0) write(17,'(f9.6,6(1x,f9.3))')
       $ time/tscale,
  
negative.code: c    -diagnostics
  c *********************************************************
  c
  c urc & plo modified
        if(mod(kstep,nlogstep).eq.0.or.kel.eq.2) then
           botlum=-4.*pi*y00*r(nn)**2*opr*
       $      (2.*real(ds(1,nn)))
  