anchor.diff:              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
  +         
  +!debug_print_unit = 77
  +!call WrNumAryFileNR(debug_print_unit,(/real(k,DbKi), this_time/), "1x,ES15.5E3", ErrStat, ErrMsg  ) 
  +!call WrNumAryFileNR(debug_print_unit,u_Orca%PtfmMesh%TranslationDisp(:,1), "1x,ES15.5E3", ErrStat, ErrMsg  ) 
  +!call WrNumAryFileNR(debug_print_unit,u_Orca%PtfmMesh%TranslationVel(:,1), "1x,ES15.5E3", ErrStat, ErrMsg  ) 
  +!call WrNumAryFileNR(debug_print_unit,u_Orca%PtfmMesh%TranslationAcc(:,1), "1x,ES15.5E3", ErrStat, ErrMsg  ) 
  +!call WrNumAryFileNR(debug_print_unit,u_Orca%PtfmMesh%Orientation(1,:,1), "1x,ES15.5E3", ErrStat, ErrMsg  ) 
  +!call WrNumAryFileNR(debug_print_unit,u_Orca%PtfmMesh%Orientation(2,:,1), "1x,ES15.5E3", ErrStat, ErrMsg  ) 
  +!call WrNumAryFileNR(debug_print_unit,u_Orca%PtfmMesh%Orientation(3,:,1), "1x,ES15.5E3", ErrStat, ErrMsg  ) 
  +!call WrNumAryFileNR(debug_print_unit,u_Orca%PtfmMesh%RotationVel(:,1), "1x,ES15.5E3", ErrStat, ErrMsg  ) 
  +!call WrNumAryFileNR(debug_print_unit,u_Orca%PtfmMesh%RotationAcc(:,1), "1x,ES15.5E3", ErrStat, ErrMsg  ) 
  +!write(debug_print_unit,'()')
  +!         
  +!debug_print_unit = 78
  +!call WrNumAryFileNR(debug_print_unit,(/real(k,DbKi), this_time/), "1x,ES15.5E3", ErrStat, ErrMsg  ) 
  +!call WrNumAryFileNR(debug_print_unit,y_Orca%PtfmMesh%Force(:,1), "1x,ES15.5E3", ErrStat, ErrMsg  ) 
  +!call WrNumAryFileNR(debug_print_unit,y_Orca%PtfmMesh%Moment(:,1), "1x,ES15.5E3", ErrStat, ErrMsg  ) 
  +!write(debug_print_unit,'()')
            
            
            IF ( ErrStat >= AbortErrLev ) THEN
  
positive.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)
                                                                             
        END IF
        
     !..................
  
negative.code:          END IF
           
           IF ( p_FAST%CompElast == Module_BD ) THEN 
              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  )
              end do            
           END IF
           
  