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
          if_normal = zero
  c
          if_blocks1: do iblk = 1, nelblif
  c
            iel     = lcblkif(1, iblk)
  
negative.code:         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"
  