anchor.diff:  c
           enddo
   c
  -c.... end re-position layered mesh
  -c
   c.... -------------------->     communication     <--------------------
   c
   c.... combine disp + flag into dispCommu
  
positive.code:         real*8, allocatable :: Estiff(:,:,:)
  c
  c.... ------------------->   layer base elements   <-------------------
  c
  c.... calculate the normal of each growth curve based on new boundary positions
  c
  
negative.code:           call dbl_to_iBCelas(flag(i), iBC(i))
            call setBLbc(disp(i,:), iBC(i), BC(i, ndof+2:ndof+5))
          enddo
  c
  c.... end if growth curve
  c
        endif ! end if numgc greater than 0
  c
  c.... end re-position layered mesh
  c
  c.... -------------------->   interior elements   <--------------------
  c
  