anchor.diff:           END IF
            
            IF ( p_FAST%CompElast == Module_BD ) THEN 
  -!if (k==p_FAST%KMax) debug_print = .true.
  -!debug_print_unit = 70
               do nb=1,p_FAST%nBeams
                  CALL BD_CalcOutput( this_time, u_BD(nb), p_BD(nb), x_BD(nb), xd_BD(nb), z_BD(nb), OtherSt_BD(nb), y_BD(nb), ErrStat2, ErrMsg2 )
                     CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName  )
  -!debug_print = .false.
               end do            
            END IF
            
  
positive.code:             CALL Orca_CalcOutput( this_time, u_Orca, p_Orca, x_Orca, xd_Orca, z_Orca, OtherSt_Orca, y_Orca, ErrStat2, ErrMsg2 )
                 CALL SetErrStat( ErrStat2, ErrMsg2, ErrStat, ErrMsg, RoutineName  )
           END IF
           
           
           IF ( ErrStat >= AbortErrLev ) THEN
  
negative.code:             ! we're mapping loads, so we also need the sibling meshes' displacements:
           CALL Transfer_Point_to_Point( y_HD2%AllHdroOrigin, MeshMapData%u_ED_PlatformPtMesh, MeshMapData%HD_W_P_2_ED_P, ErrStat2, ErrMsg2, MeshMapData%u_HD_Mesh, y_ED2%PlatformPtMesh) !u_HD and u_mapped_positions contain the displaced positions for load calculations
              CALL SetErrStat(ErrStat2,ErrMsg2, ErrStat, ErrMsg, RoutineName)
                                        
        ELSE
           
              ! When using OrcaFlex, we need to zero this out
           MeshMapData%u_ED_PlatformPtMesh%Force  = 0.0_ReKi         
           MeshMapData%u_ED_PlatformPtMesh%Moment = 0.0_ReKi
           
        END IF
        
     !..................
  