anchor.diff: +      module simmetrix_snap
  +      use iso_c_binding
  +
  +      interface
  +        subroutine sim_closest_pos ( dx, dy, dz, id )
  +     &    bind(C, NAME='sim_closest_pos')
  +        use iso_c_binding
  +          real(c_double), value :: dx, dy, dz
  +          integer(c_int), value :: id
  +        end subroutine
  +
  +c        subroutine get_model_velocity (v)
  +c     &    bind(C, NAME='get_model_velocity')
  +c        use iso_c_binding
  +c          real(c_double) :: v
  +c        end subroutine
  +      end interface
  +
  +      end module
  +c
  +c----------------------------------------------------------------------
  +c
  +c----------------------------------------------------------------------
  +c
           subroutine updateSnapSurfBC ( x,       disp_snap,
        &                                iBC,     BC)
   c
  
positive.code:         subroutine resetSnapBC ( x,       disp_snap,
       &                           iBC,     BC,     surfID_snap)
  c
          include "common.h"
  c
          real*8    x(numnp,nsd),  disp_snap(numnp,nsd)
  
negative.code:         real*8    normal(nshg, nsd)
  c
          if ( timeDepComp1Flag .eq. 1) then ! double check
  c
          write(*,*) "Modify comp1_elas BC on surfID",timeDepComp1ID
  c
          surfID_TDComp1 = zero
          normal = zero
  