anchor.diff:    !C*** hecmw_mpc_mat_init_explicit
     !C***
     !C
  -  subroutine hecmw_mpc_mat_init_explicit(hecMESH, hecMAT, hecMATmpc)
  -    implicit none
  -    type (hecmwST_local_mesh), intent(inout) :: hecMESH
  +  subroutine hecmw_mpc_mat_init_explicit(hecMESH, hecMAT, hecMESHmpc, hecMATmpc)
  +    implicit none
  +    type (hecmwST_local_mesh), intent(inout), target :: hecMESH
       type (hecmwST_matrix), intent(in), target :: hecMAT
  +    type (hecmwST_local_mesh), pointer :: hecMESHmpc
       type (hecmwST_matrix), pointer :: hecMATmpc
       integer(kind=kint) :: totalmpc, MPC_METHOD
   
  
positive.code:     endif
    end subroutine hecmw_tback_x
  
  end module hecmw_mpc_prepost
  
negative.code: 
      select case (MPC_METHOD)
      case (1)  ! penalty
        nullify(hecMESHmpc)
        nullify(hecMATmpc)
      case (2) ! MPCCG
        nullify(hecMESHmpc)
        nullify(hecMATmpc)
      case (3) ! elimination
        call hecmw_mpc_mesh_free(hecMESHmpc)
        deallocate(hecMESHmpc)
        nullify(hecMESHmpc)
        call hecmw_mat_finalize(hecMATmpc)
        deallocate(hecMATmpc)
        nullify(hecMATmpc)
  