#include <iostream>
#include <queue>
#include <set>
#include <stack>
#include <utility>
#include <tuple>
#include <iomanip>
#include <map>
#include <algorithm>
#include <math.h>
#include <string>
#include <vector>
#include <unordered_map>
#define ll long long
#define inf 10000000000007
#define mod 1000000007

using namespace std;

const int BLOCK = 450;


ll fastpow(ll n, ll x) {

	if (x == 0) {
		return 1;
	}
	else {
		ll ret = fastpow(n, x / 2);
		ret = ((ret % mod) * (ret % mod)) % mod;
		if (x % 2 == 0) {
			return ret;
		}
		else {
			return ((ret) * (n)) % mod;
		}
	}
}


ll gcd(ll a, ll b) {
	if (a == 0) {
		return b;
	}
	if (b == 0) {
		return a;
	}
	else {
		return gcd(b, a % b);
	}
}
ll lcm(ll a, ll b) {
	ll val = (a % mod * b % mod) % mod;
	val = (val * fastpow(gcd(a, b), mod - 2)) % mod;
	return val;
}
int Logk(ll n, ll k) {
	if (k == 1) {
		return 1;
	}
	if (n == 0) {
		return 0;
	}
	int count = -1;
	while (n > 0) {
		count++;
		n /= k;
	}
	return count;
}
struct Dsu {
	vector<int> par;

	void init(int n) {
		par.resize(n + 5, 0);
		for (int i = 1; i <= n; i++) par[i] = i;
	}

	int find(int u) {
		if (par[u] == u) return u;
		return par[u] = find(par[u]);
	}

	bool join(int u, int v) {
		u = find(u); v = find(v);
		if (u == v) return false;
		par[v] = u;
		return true;
	}
} dsu;

bool cmp(pair<vector<ll>, ll> a, pair<vector<ll>, ll> b) {
	if (a.first[0] == b.first[0]) {
		return a.second < b.second;
	}
	else {
		return a.first[0] < b.first[0];
	}
}
void solve() {
	int n, m, k;
	cin >> n >> m >> k;
	vector<pair<vector<ll>, ll>> a(n + 1);
	vector<pair<vector<ll>, ll>> b(n + 1);
	map<pair<ll,ll>, pair<ll, ll>> sum;
	for (int i = 1; i <= n; i++) {
		
		for (int j = 0; j <= m; j++) {
			ll x;
			cin >> x;
			a[i].second = i;
			a[i].first.push_back(x);
			if (j > 0) {
				if (x == -1) {
					sum[{a[i].first[0], i}].second += 1;
				}
				else {
					sum[{a[i].first[0], i}].first += x;
				}
				
			}
			
		}
	}
	b = a;
	sort(b.begin() + 1, b.end());
	ll prev = -10;
	ll prev_pos = 0;
	for (int i = 1; i <= n; i++) {
		ll cur_sum = sum[{b[i].first[0], b[i].second}].first;
		ll rate = b[i].first[0];
		ll id = b[i].second;
		ll slot = sum[{b[i].first[0], b[i].second}].second;
		
		if (rate == prev_pos) {
			if (cur_sum >= prev) {
				for (int j = 1; j <= m; j++) {
					if (a[id].first[j] == -1) {
						a[id].first[j] = 0;
					}
					
				}
				prev = cur_sum;
			}
			else {
				int need = prev - cur_sum;
				for (int j = 1; j <= m; j++) {
					if (a[id].first[j] == -1) {
						a[id].first[j] = min(k, need);
						need -= a[id].first[j];
						cur_sum += a[id].first[j];
					}

				}
				if (cur_sum < prev) {
					cout << "No\n";
					return;
				}
				else {
					prev = cur_sum;
				}
			

			}
		}
		else {
			if (cur_sum > prev) {
				for (int j = 1; j <= m; j++) {
					if (a[id].first[j] == -1) {
						a[id].first[j] = 0;
					}

				}
				prev = cur_sum;
			}
			else {
				int need = prev - cur_sum + 1;
				for (int j = 1; j <= m; j++) {
					if (a[id].first[j] == -1) {
						a[id].first[j] = min(k, need);
						need -= a[id].first[j];
						cur_sum += a[id].first[j];
					}

				}
				if (cur_sum <= prev) {
					cout << "No\n";
					return;
				}
				else {
					prev = cur_sum;
				}

			}
		}
		prev_pos = rate;
	}
	cout << "Yes\n";
	for (int i = 1; i <= n; i++) {
		for (int j = 1; j <= m; j++) {
			cout << a[i].first[j] << " ";
		}
		cout << endl;
	}
}



int main() {




	ios_base::sync_with_stdio(false);
	cin.tie(NULL);
	cout.tie(NULL);
	
	int t;
	cin >> t;
	while (t--) {
		solve();
	}






























}