anchor.diff:              indx_first = indx_first + 1
            end do      
         end do
  -                           
  -      do i=1,M%NNodes
  -         do j=1,3
  -            Ary(indx_first) = M%Moment(j,i)
  -            indx_first = indx_first + 1
  -         end do      
  -      end do
  -            
  +      
  +         ! This is needed for MAP meshes because it only contains the Force field not the Moment field
  +      if ( M%fieldmask(MASKID_Moment) .AND. ALLOCATED(M%Moment)) then     
  +         do i=1,M%NNodes
  +            do j=1,3
  +               Ary(indx_first) = M%Moment(j,i)
  +               indx_first = indx_first + 1
  +            end do      
  +         end do
  +      end if     
   
      END SUBROUTINE PackLoadMesh
   !...............................................................................................................................
  
positive.code:             indx_first = indx_first + 1
           end do      
        end do
                             
        do i=1,M%NNodes
           do j=1,3
              Names(indx_first) = trim(MeshName)//' '//Comp(j)//' moment, node '//trim(num2lstr(i))//', Nm'//UnitDesc
              indx_first = indx_first + 1
           end do      
        end do
              
  
     END SUBROUTINE PackLoadMesh_Names
  !...............................................................................................................................
  
negative.code:          indx_first = indx_last + 1
        end do
  
           ! This is needed for MAP meshes because it only contains the Force field not the Moment field
        if ( M_p%fieldmask(MASKID_Moment) .AND. ALLOCATED(M_p%Moment)) then     
           do i=1,M_p%NNodes
              indx_last  = indx_first + 2 
              dY(indx_first:indx_last) = M_p%Moment(:,i) - M_m%Moment(:,i)
              indx_first = indx_last + 1
           end do
        end if
        
     END SUBROUTINE PackLoadMesh_dY
  !...............................................................................................................................
  !> This subroutine returns the names of rows/columns of motion meshes in the Jacobian matrices. It assumes all fields marked
  