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:   !Do the blade first
    allocate(p_OpFM%forceBldRnodes(p_OpFM%NnodesForceBlade), stat=errStat2)
    dRforceNodes = p_OpFM%BladeLength/(p_OpFM%NnodesForceBlade-1)
    do i=1,p_OpFM%NnodesForceBlade
       p_OpFM%forceBldRnodes(i) = (i-1)*dRforceNodes 
    end do
  
  
    !Do the tower now
    allocate(p_OpFM%forceTwrHnodes(p_OpFM%NnodesForceTower), stat=errStat2)
    dRforceNodes = p_OpFM%TowerHeight/(p_OpFM%NnodesForceTower-1)
    do i=1,p_OpFM%NnodesForceTower
       p_OpFM%forceTwrHnodes(i) = (i-1)*dRforceNodes
    end do
    
    return
  
  END SUBROUTINE OpFM_CreateActForceBladeTowerNodes
  
negative.code:   ! Now calculate the positions of the force nodes based on interpolation
    DO I=1,p_OpFM%NnodesForceTower ! Calculate the position of the force nodes
       jLower=1
       do while ( ((hStructNodes(jLower) - p_OpFM%forceTwrHnodes(I))*(hStructNodes(jLower+1) - p_OpFM%forceTwrHnodes(I)) .gt. 0) .and. (jLower .lt. nStructNodes))
          jLower = jLower + 1
       end do
       hInterp =  (p_OpFM%forceTwrHnodes(I) - hStructNodes(jLower))/(hStructNodes(jLower+1)-hStructNodes(jLower)) ! The location of this force node in (0,1) co-ordinates between the jLower and jLower+1 nodes
  