#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const ll N = 105;
ll s[N][N];
char t[N][N];
ll dx[] = {1, 1, 1, 0, 0, -1, -1, -1}, dy[] = {1, 0, -1, 1, -1, 1, 0, -1};
vector<ll>g[N * N];
bool st[N * N];
ll match[N * N];
ll find(ll x) {
	for (auto v : g[x]) {
		if (!st[v]) {
			st[v] = 1;
			if (match[v] == 0 || find(match[v])) {
				match[v] = x;
				return 1;
			}
		}
	}
	return 0;
}
ll calc(ll x, ll y) {
	ll res = 0;
	for (ll i = 0; i < 8; i++) {
		res += (t[x + dx[i]][y + dy[i]] == '#');
	}
	return res;
}
ll n, m;
ll get(ll x, ll y) {
	return (x - 1) * m + y;
}
void add(ll x, ll y, ll u, ll v) {
	// cout << "add " << u << " " << v << "\n";
	if ((x + y) & 1) {
		g[v].push_back(u);
	} else {
		g[u].push_back(v);

	}
}
int main() {
	cin >> n >> m;
	for (ll i = 0; i <= n * 2; i++) {
		for (ll j = 0; j <= m * 2; j++) {
			t[i][j] = '#';
		}
	}
	for (ll i = 1; i <= 2 * n - 1; i++) {
		for (ll j = 1; j <= 2 * m - 1; j++) {
			cin >> t[i][j];
		}
	}
	ll cur = 0, cnt = 0;
	for (ll i = 1; i <= n; i++) {
		for (ll j = 1; j <= m; j++) {
			s[i][j] = t[i * 2 - 1][j * 2 - 1] - '0' + calc(i * 2 - 1, j * 2 - 1);
		}
	}
	for (ll i = 1; i <= 2 * n - 1; i++) {
		for (ll j = 1; j <= 2 * m - 1; j++) {
			if (t[i][j] == '.') {
				t[i][j] = '#';
			}
		}
	}
	for (ll i = 1; i <= n; i++) {
		for (ll j = 1; j <= m; j++) {
			if (s[i][j] == 7) {
				cnt++;
			} else {
				continue;
			}
			if (j != 1 && j != m - 1) {
				if (s[i][j + 1] == 7) {
					add(i, j, get(i, j), get(i, j + 1));
				}
			}
			if (i != 1 && i != n - 1) {
				if (s[i + 1][j] == 7) {
					add(i, j, get(i, j), get(i + 1, j));
				}
			}
			if (i == 1 && j == 1) {
				if (s[i + 1][j] == 7) {
					add(i, j, get(i, j), get(i + 1, j));
				}
				if (s[i][j + 1] == 7) {
					add(i, j, get(i, j), get(i, j + 1));
				}
			}
			if (i == 1 && j == m) {
				if (s[i + 1][j] == 7) {
					add(i, j, get(i, j), get(i + 1, j ));
				}
				if (s[i][j - 1] == 7) {
					add(i, j, get(i, j), get(i, j - 1));
				}
			}
			if (i == n && j == 1) {
				if (s[i - 1][j] == 7) {
					add(i, j, get(i, j), get(i - 1, j));
				}
				if (s[i][j + 1] == 7) {
					add(i, j, get(i, j), get(i, j + 1));
				}
			}
			if (i == n && j == m) {
				if (s[i - 1][j] == 7) {
					add(i, j, get(i, j), get(i - 1, j ));
				}
				if (s[i][j - 1] == 7) {
					add(i, j, get(i, j), get(i, j - 1));
				}
			}
		}
	}
	cur = 0;
	ll res = 0;
	for (ll i = 1; i <= n; i++) {
		for (ll j = 1; j <= m; j++) {
			if (((i + j) % 2 == 0) && s[i][j] == 7) {
				memset(st, 0, sizeof(st));
				if (find(get(i, j))) {
					res++;
				}
				// cout << i << " " << j << " " << match[get(i, j)] << " " << get(i, j) << "\n";
			}
		}
	}
	if (res * 2 != cnt) {
		puts("NO");
		return 0;
	}
	puts("YES");
	cur = 0;
	for (ll i = 1; i <= n; i++) {
		for (ll j = 1; j <= m; j++) {
			cur++;
			if (((i + j) & 1) && s[i][j] == 7 && match[get(i, j)]) {
				if (match[get(i, j)] - get(i, j) == 1) {
					t[2 * i - 1][2 * j] = '.';
				} else if (match[get(i, j)] - get(i, j) == -1) {
					t[2 * i - 1][2 * j - 2] = '.';
				} else if (match[get(i, j)] - get(i, j) == m) {
					t[2 * i][2 * j - 1] = '.';
				} else if (match[get(i, j)] - get(i, j) == -m) {
					t[2 * i - 2][2 * j - 1] = '.';
				}
			}
		}
	}
	for (ll i = 1; i <= 2 * n - 1; i++) {
		for (ll j = 1; j <= 2 * m - 1; j++) {
			cout << t[i][j];
		}
		cout << '\n';
	}
}
/*
3 3
2.4.3
.....
5.8.5
.....
3.5.3
*/