anchor.diff:     ! HD inputs:
      !...............
            
  -   !(Morison%LumpedMesh)
  -   do i=1,HD_M_LumpedMesh%NNodes
  +   !(Morison%Mesh)
  +   do i=1,HD_M_Mesh%NNodes
         do j=1,3
  -         MeshMapData%Jac_u_indx(index,1) =  9 !Module/Mesh/Field: u_HD%Morison%LumpedMesh%TranslationAcc = 9
  +         MeshMapData%Jac_u_indx(index,1) =  9 !Module/Mesh/Field: u_HD%Morison%Mesh%TranslationAcc = 9
            MeshMapData%Jac_u_indx(index,2) =  j !index:  j
            MeshMapData%Jac_u_indx(index,3) =  i !Node:   i
            index = index + 1
         end do !j      
      end do !i
      
  -   do i=1,HD_M_LumpedMesh%NNodes
  +   do i=1,HD_M_Mesh%NNodes
         do j=1,3
  -         MeshMapData%Jac_u_indx(index,1) = 10 !Module/Mesh/Field:  u_HD%Morison%LumpedMesh%RotationAcc = 10
  +         MeshMapData%Jac_u_indx(index,1) = 10 !Module/Mesh/Field:  u_HD%Morison%Mesh%RotationAcc = 10
            MeshMapData%Jac_u_indx(index,2) =  j !index:  j
            MeshMapData%Jac_u_indx(index,3) =  i !Node:   i
            index = index + 1
         end do !j      
      end do !i     
      
  -   
  -   !(Morison%DistribMesh)
  -   do i=1,HD_M_DistribMesh%NNodes
  -      do j=1,3
  -         MeshMapData%Jac_u_indx(index,1) = 11 !Module/Mesh/Field: u_HD%Morison%DistribMesh%TranslationAcc = 11
  -         MeshMapData%Jac_u_indx(index,2) =  j !index:  j
  -         MeshMapData%Jac_u_indx(index,3) =  i !Node:   i
  -         index = index + 1
  -      end do !j      
  -   end do !i
  -   
  -   do i=1,HD_M_DistribMesh%NNodes
  -      do j=1,3
  -         MeshMapData%Jac_u_indx(index,1) = 12 !Module/Mesh/Field:  u_HD%Morison%DistribMesh%RotationAcc = 12
  -         MeshMapData%Jac_u_indx(index,2) =  j !index:  j
  -         MeshMapData%Jac_u_indx(index,3) =  i !Node:   i
  -         index = index + 1
  -      end do !j      
  -   end do !i     
      
      !(Mesh)
      do i=1,HD_WAMIT_Mesh%NNodes
  
positive.code: !----------------------------------------------------------------------------------------------------------------------------------
  !> This routine transfers the HD outputs into inputs required for ED. Note that this *adds* to the values already in 
  !! u_SD_LMesh (so initialize it before calling this routine).
  SUBROUTINE Transfer_HD_to_SD( u_mapped, u_SD_LMesh, u_mapped_positions, y_HD, u_HD_M_LumpedMesh, u_HD_M_DistribMesh, MeshMapData, ErrStat, ErrMsg )
  !..................................................................................................................................
     TYPE(MeshType),                 INTENT(INOUT)  :: u_mapped                 !< temporary copy of SD mesh (an argument to avoid another temporary mesh copy)
     TYPE(MeshType),                 INTENT(INOUT)  :: u_SD_LMesh               !< SD Inputs on LMesh at t (separate so we can call from FullOpt1_InputOutputSolve with temp meshes)
     TYPE(MeshType),                 INTENT(IN   )  :: u_mapped_positions       !< Mesh sibling of u_mapped, with displaced positions
     TYPE(HydroDyn_OutputType),      INTENT(IN   )  :: y_HD                     !< HydroDyn outputs
     TYPE(MeshType),                 INTENT(IN   )  :: u_HD_M_LumpedMesh        !< HydroDyn input mesh (separate so we can call from FullOpt1_InputOutputSolve with temp meshes)
     TYPE(MeshType),                 INTENT(IN   )  :: u_HD_M_DistribMesh       !< HydroDyn input mesh (separate so we can call from FullOpt1_InputOutputSolve with temp meshes)
     
     TYPE(FAST_ModuleMapType),       INTENT(INOUT)  :: MeshMapData              !< Data for mapping between modules
     INTEGER(IntKi),                 INTENT(  OUT)  :: ErrStat                  !< Error status
  
negative.code: !----------------------------------------------------------------------------------------------------------------------------------
  !> This routine initializes the array that maps rows/columns of the Jacobian to specific mesh fields.
  !! Do not change the order of this packing without changing subroutine Create_FullOpt1_UVector()!
  SUBROUTINE Init_FullOpt1_Jacobian( p_FAST, MeshMapData, ED_PlatformPtMesh, SD_TPMesh, SD_LMesh, HD_M_Mesh,  &
                                     HD_WAMIT_Mesh, ED_HubPtLoad, u_BD, Orca_PtfmMesh, ExtPtfm_PtfmMesh, ErrStat, ErrMsg)
  
     TYPE(FAST_ParameterType)          , INTENT(INOUT) :: p_FAST                !< FAST parameters               
  