#include <bits/stdc++.h>
char s[1001][1001];
int n, m, cnt, tot, T, lx, ly, x[1001], y[1001];
void Check(int x1, int y1, int x2, int y2, int x3, int y3, int x4, int y4) {
  if (x4 - x3 < 2 || y4 - y3 < 2) return;
  T++;
  cnt = 0;
  for (int j = y1; j <= y2; j++) {
    if (s[x1][j] < T) s[x1][j] = T, cnt++;
    if (s[x2][j] < T) s[x2][j] = T, cnt++;
  }
  for (int i = x1; i <= x2; i++) {
    if (s[i][y1] < T) s[i][y1] = T, cnt++;
    if (s[i][y2] < T) s[i][y2] = T, cnt++;
  }
  for (int j = y3; j <= y4; j++) {
    if ((s[x3][j] * s[x4][j]) == 0) return;
    if (s[x3][j] < T) s[x3][j] = T, cnt++;
    if (s[x4][j] < T) s[x4][j] = T, cnt++;
  }
  for (int i = x3; i <= x4; i++) {
    if ((s[i][y3] * s[i][y4]) == 0) return;
    if (s[i][y3] < T) s[i][y3] = T, cnt++;
    if (s[i][y4] < T) s[i][y4] = T, cnt++;
  }
  if (cnt == tot) {
    puts("YES");
    printf("%d %d %d %d\n%d %d %d %d", x1 + 1, y1 + 1, x2 + 1, y2 + 1, x3 + 1,
           y3 + 1, x4 + 1, y4 + 1);
    exit(0);
  }
}
void Calc(int x1, int y1, int x2, int y2) {
  if (x2 - x1 < 2 || y2 - y1 < 2) return;
  for (int i = x1; i <= x2; i++)
    if ((s[i][y1] * s[i][y2]) == 0) return;
  for (int j = y1; j <= y2; j++)
    if ((s[x1][j] * s[x2][j]) == 0) return;
  for (int x3 = 0; x3 < 4; x3++)
    for (int x4 = 0; x4 < 4; x4++)
      for (int y3 = 0; y3 < 4; y3++)
        for (int y4 = 0; y4 < 4; y4++)
          Check(x1, y1, x2, y2, x[x3], y[y3], x[x4], y[y4]);
}
int main() {
  scanf("%d%d\n", &n, &m);
  for (int i = 0; i < n; i++) gets(s[i]);
  for (int i = 0; i < n; i++)
    for (int j = 0; j < m; j++) tot += s[i][j] = s[i][j] == '#';
  for (int i = 0; i < n; i++)
    for (int j = 0; j < m - 2; j++)
      if (s[i][j] & s[i][j + 1] & s[i][j + 2]) {
        x[lx++] = i;
        break;
      }
  for (int j = 0; j < m; j++)
    for (int i = 0; i < n - 2; i++)
      if (s[i][j] & s[i + 1][j] & s[i + 2][j]) {
        y[ly++] = j;
        break;
      }
  if (lx > 3) x[2] = x[lx - 2], x[3] = x[lx - 1];
  if (ly > 3) y[2] = y[ly - 2], y[3] = y[ly - 1];
  T = 1;
  for (int x1 = 0; x1 < 4; x1++)
    for (int x2 = 0; x2 < 4; x2++)
      for (int y1 = 0; y1 < 4; y1++)
        for (int y2 = 0; y2 < 4; y2++) Calc(x[x1], y[y1], x[x2], y[y2]);
  puts("NO");
  return 0;
}
