anchor.diff:     FieldMask(MASKID_TRANSLATIONVel) = .false.
      FieldMask(MASKID_RotationVel) = .true.
      call PackMotionMesh_Names(u%HubMotion, 'Hub', InitOut%LinNames_u, index, FieldMask=FieldMask)
  -   
  +
  +   index_last = index   
      FieldMask = .false.
      FieldMask(MASKID_Orientation) = .true.
      do k = 1,p%NumBlades
  
positive.code:    
     
        ! get the names of the linearized outputs:
     call AllocAry(InitOut%LinNames_y, p%Jac_ny,'',ErrStat2,ErrMsg2)
        call SetErrStat(ErrStat2,ErrMsg2,ErrStat,ErrMsg,RoutineName)
        if (ErrStat >= AbortErrLev) return
     
           
     indx_first = 1  
     call PackLoadMesh_Names(y%TowerLoad, 'Tower', InitOut%LinNames_y, indx_first)
     do k=1,p%NumBlades
        call PackLoadMesh_Names(y%BladeLoad(k), 'Blade '//trim(num2lstr(k)), InitOut%LinNames_y, indx_first)
     end do
  
     do i=1,p%NumOuts
        InitOut%LinNames_y(i+indx_first) = trim(p%OutParam(i)%Name)//', '//p%OutParam(i)%Units
     end do    
  
       
  END SUBROUTINE Init_Jacobian_y
  !----------------------------------------------------------------------------------------------------------------------------------
  !> This routine initializes the array that maps rows/columns of the Jacobian to specific mesh fields.
  
negative.code:    CHARACTER(*), PARAMETER                           :: RoutineName = 'Init_Jacobian'
     
        ! local variables:
     INTEGER(IntKi)                :: i, j, k, index, index_last, nu, i_meshField
     REAL(ReKi)                    :: perturb, perturb_t, perturb_b(MaxBl)
     LOGICAL                       :: FieldMask(FIELDMASK_SIZE)
     CHARACTER(1), PARAMETER       :: UVW(3) = (/'U','V','W'/)
  