anchor.diff:            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
  
positive.code: c
  c.... end calculation of growth curve normal
  c
  c.... ---------------->   Re-position layered mesh   <-----------------
  c
  c.... loop over growth curves
  c
          flag = zero
          listcounter = 0
          ioffset = 1 ! the ID starts from 1 in phasta
          do ngc = 1, numgc
  
negative.code:         real*8, allocatable :: Estiff(:,:,:)
  c
  c.... ------------------->   layer base elements   <-------------------
  c
        if (numgc .gt. 0) then
  c
  c.... calculate the normal of each growth curve based on new boundary positions
  c
  