anchor.diff:       
   
     ! The tower now
  -  nNodesTowerProps = SIZE(InitOut_AD%TwrElev)
     do k = p_OpFM%NumBl+1,p_OpFM%NMappings
  +     nNodesTowerProps = SIZE(InitOut_AD%TwrElev)
        ! Calculate the chord at the force nodes based on interpolation
        DO I=1,p_OpFM%NnodesForceTower
           Node = Node + 1
  
positive.code:          
        END DO !J = 1,p%BldNodes ! Loop through the blade nodes / elements
     END DO !K = 1,p%NumBl
     
     ! tower nodes
     DO J=1,u_AD%TowerMotion%nnodes
        Node = Node + 1
        OpFM%u%pxVel(Node) = u_AD%TowerMotion%TranslationDisp(1,J) + u_AD%TowerMotion%Position(1,J)
        OpFM%u%pyVel(Node) = u_AD%TowerMotion%TranslationDisp(2,J) + u_AD%TowerMotion%Position(2,J)
        OpFM%u%pzVel(Node) = u_AD%TowerMotion%TranslationDisp(3,J) + u_AD%TowerMotion%Position(3,J)
     END DO
     
     
     ! Do the Actuator Force nodes now
     Node = 1   ! displaced hub position 
     OpFM%u%pxForce(Node) = OpFM%u%pxVel(Node)
  
negative.code:       ! number of force nodes in the interface
     Opfm%p%NnodesForceBlade =  InitInp%NumActForcePtsBlade 
     OpFM%p%NnodesForceTower = InitInp%NumActForcePtsTower
     OpFM%p%NnodesForce = 1 + OpFM%p%NumBl * InitInp%NumActForcePtsBlade
     OpFM%p%BladeLength = InitInp%BladeLength
  
     if ( y_AD%TowerLoad%NNodes > 0 ) then
        OpFM%p%NMappings = OpFM%p%NumBl + 1
        OpFM%p%TowerHeight = InitInp%TowerHeight
        OpFM%p%NnodesForce = OpFM%p%NnodesForce + InitInp%NumActForcePtsTower
     else
        OpFM%p%NMappings = OpFM%p%NumBl
     end if
  
        ! air density, required for normalizing values sent to OpenFOAM:
     OpFM%p%AirDens = AirDens
  