anchor.diff:       
   
     ! The tower now
  +  nNodesTowerProps = SIZE(InitOut_AD%TwrElev)
     do k = p_OpFM%NumBl+1,p_OpFM%NMappings
        ! Calculate the chord at the force nodes based on interpolation
        DO I=1,p_OpFM%NnodesForceTower
           Node = Node + 1
           jLower=1
  -        do while ( (InitOut_AD%TwrElev(jLower) - p_OpFM%forceTwrHnodes(I))*(InitOut_AD%TwrElev(jLower+1) - p_OpFM%forceTwrHnodes(I)) .gt. 0 ) !Determine the closest two nodes at which the blade properties are specified
  +        do while ( ( (InitOut_AD%TwrElev(jLower) - p_OpFM%forceTwrHnodes(I))*(InitOut_AD%TwrElev(jLower+1) - p_OpFM%forceTwrHnodes(I)) .gt. 0) .and. (jLower .le. nNodesTowerProps) ) !Determine the closest two nodes at which the blade properties are specified
              jLower = jLower + 1
           end do
           rInterp =  (p_OpFM%forceTwrHnodes(I) - InitOut_AD%TwrElev(jLower))/(InitOut_AD%TwrElev(jLower+1)-InitOut_AD%TwrElev(jLower)) ! The location of this force node in (0,1) co-ordinates between the jLower and jLower+1 nodes
  
positive.code:   !Local variables
    INTEGER(IntKI)                         :: i,j,k,node  ! Loop variables
    INTEGER(IntKI)                         :: nNodesBladeProps ! Number of nodes in the blade properties for a given blade
    INTEGER(IntKI)                         :: jLower      ! Index of the blade properties node just smaller than the force node
    INTEGER(IntKi)                         :: ErrStat2    ! temporary Error status of the operation
    CHARACTER(ErrMsgLen)                   :: ErrMsg2     ! temporary Error message if ErrStat /= ErrID_None
  
negative.code:   ! Now calculate the positions of the force nodes based on interpolation
    DO I=1,p_OpFM%NnodesForceBlade ! Calculate the position of the force nodes
       jLower=1
       do while ( ( (rStructNodes(jLower) - p_OpFM%forceBldRnodes(I))*(rStructNodes(jLower+1) - p_OpFM%forceBldRnodes(I)) .gt. 0) .and. (jLower .lt. nStructNodes) )
          jLower = jLower + 1
       end do
       rInterp =  (p_OpFM%forceBldRnodes(I) - rStructNodes(jLower))/(rStructNodes(jLower+1)-rStructNodes(jLower)) ! The location of this force node in (0,1) co-ordinates between the jLower and jLower+1 nodes
  