anchor.diff:          use e3_solid_m
           use probe_m
           use ifbc_def_m
  +        use weighted_normal_m !hacking for weighted normal
  +        use interfaceflag ! hacking for weighted normal        
   c
           include "common.h"
           include "mpif.h"
  
positive.code: c
          if_normal = zero
  c
          if_blocks1: do iblk = 1, nelblif
  c
            iel     = lcblkif(1, iblk)
  
negative.code:           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
  