anchor.diff:  
       IF (maxi) tranw = .TRUE.
   
  -    IF (step>1000) THEN
  +    IF (step>10000) THEN
           WRITE(ounit,*) 'TOO SMALL TIME STEPS. STOPPING'
           STOP
       END IF
   
   END DO
  +
  +IF (aprad == 1 .OR. apaxi == 1) THEN
  +    CALL PowDis(npow)
  +END IF
  +
  +IF (aprad == 1) CALL AsmPow(npow)
  +
  +IF (apaxi == 1) CALL AxiPow(npow)
  +
  +IF (afrad == 1) CALL AsmFlux(f0, 1.e0_DP)
   
   
   END SUBROUTINE trod_eject
  
positive.code:                  f0, fx1, fy1, fz1, fx2, fy2, fz2, &
                   c0, cx1, cy1, cz1, cx2, cy2, cz2, tbeta, omeg, &
                   npow, pow, ppow, node_nf, ix, iy, iz, zdel, tranw, &
                   ft, ftx1, fty1, ftz1, ftx2, fty2, ftz2
  USE InpOutp, ONLY: XS_updt, ounit
  USE nodal, ONLY: nodal_coup4, outer4, outertf, outer4ad, PowTot, powdis
  USE th, ONLY: th_iter, th_trans, par_ave, par_max, par_ave_f
  
  
negative.code:     CALL PowTot(f0, tpow2)
  
      ! Calculate node power distribution
      CALL PowDis(npow)
  
      ! Power change
      xppow = ppow * tpow2/tpow1 * 0.01
  