anchor.diff:              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
  -            
  +            ! 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
  +               Names(indx_first) = trim(MeshName)//' '//Comp(j)//' moment, node '//trim(num2lstr(i))//', Nm'//UnitDesc
  +               indx_first = indx_first + 1
  +            end do      
  +         end do
  +      end if      
   
      END SUBROUTINE PackLoadMesh_Names
   !...............................................................................................................................
  
positive.code:          indx_first = indx_last + 1
        end do
  
        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 SUBROUTINE PackLoadMesh_dY
  !...............................................................................................................................
  !> This subroutine returns the names of rows/columns of motion meshes in the Jacobian matrices. It assumes all fields marked
  
negative.code:             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
  !...............................................................................................................................
  