anchor.diff:  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)
  
positive.code:         use e3_solid_m
          use probe_m
          use ifbc_def_m
  c
          include "common.h"
          include "mpif.h"
  
negative.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
  