anchor.diff:                 torque = CROSS_PRODUCT( DisplacedPosition, Src%Force(:,i) / LoadsScaleFactor )
                  Dest%Moment(:,MeshMap%MapLoads(i)%OtherMesh_Element) = Dest%Moment(:,MeshMap%MapLoads(i)%OtherMesh_Element) + torque                                             
            enddo
  -      endif
  -      
  +      !endif      
  +      
  +   end if
  +   
  +      
  +   if (Src%FieldMask(MASKID_MOMENT) ) then
  +
  +      ! M_d += M_s      
  +      do i = 1, Src%NNodes
  +         !if ( MeshMap%MapLoads(i)%OtherMesh_Element < 1 )  CYCLE ! would only happen if we had non-point elements (or nodes not contained in an element)
  +
  +         Dest%Moment(:,MeshMap%MapLoads(i)%OtherMesh_Element) = Dest%Moment(:,MeshMap%MapLoads(i)%OtherMesh_Element) + (Src%Moment(:,i) / LoadsScaleFactor)
  +      end do
  +            
         Dest%Moment =  Dest%Moment * LoadsScaleFactor
      endif
         
  
positive.code:          END DO !loop through source elements
  
  
        end if ! Src has force, too      
        
                    
     end if !src has moment field
        
  #ifdef MESH_DEBUG     
  
negative.code: 
  
  !bjj note that we already checked that the following two conditions apply in this case:
  !   if Src%FieldMask(MASKID_FORCE),  Dest%FieldMask(MASKID_FORCE) and Dest%FieldMask(MASKID_MOMENT)
  !   if Src%FieldMask(MASKID_MOMENT), Dest%FieldMask(MASKID_MOMENT)
  
     if (Dest%FieldMask(MASKID_MOMENT) ) Dest%Moment = 0. ! whole array initialization; required to handle superposition of moments
     
     if (Src%FieldMask(MASKID_FORCE) ) THEN
        Dest%Force  = 0.     ! whole array initialization; required to handle superposition of forces
  