anchor.diff:            deallocate (if_kappa)
             nullify(if_kappa)
           endif
  +c... hacking, deallocation of the arrays used
  +        deallocate(w_normal_global)
  +        deallocate(length_temp)
  +        deallocate(w_normal_l0, w_normal_l1)
  +c... end of hacking        
   c
   c before the commu we need to rotate the residual vector for axisymmetric
   c boundary conditions (so that off processor periodicity is a dof add instead
  
positive.code: c        call calc_kappa_error(x,lcblkif(1,:),nelblif,nsd,nshg)
  c
          sum_vi_area = zero
          ifbc = zero
  c
          if_blocks: do iblk = 1, nelblif
  c
  
negative.code: c
          if_normal = zero
  c
  c... hacking, allocation and initialization of the weighted normal
             allocate(w_normal_global(nshg,nsd))
             allocate(length_temp(nshg))
  c
             w_normal_global = zero
             length_temp = zero
  c... end of hacking           
          if_blocks1: do iblk = 1, nelblif
  c
            iel     = lcblkif(1, iblk)
  