anchor.diff:  
       IF (maxi) tranw = .TRUE.
   
  -    IF (step>1000) THEN
  +    IF (step>10000) THEN
           WRITE(ounit,*) 'TOO SMALL TIME STEPS. STOPPING'
           STOP
       END IF
  
positive.code: 
      IF (maxi) tranw = .TRUE.
  
  END DO
  
  ! Second Time Step
  
negative.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, &
                   aprad, apaxi, afrad
  USE InpOutp, ONLY: XS_updt, ounit, AsmPow, AxiPow, AsmFlux
  USE nodal, ONLY: nodal_coup4, outer4, outertf, outer4ad, PowTot, powdis
  USE th, ONLY: th_iter, th_trans, par_ave, par_max, par_ave_f
  
  