#include <bits/stdc++.h>
using namespace std;
char a[502][502];
int n, m;
int ds[502][502][251];
int dfs(int i, int j, int x, int y, int st) {
  if (i > x || j > y) {
    return 0;
  }
  if (a[i][j] == '*' || a[x][y] == '*' || a[i][j] != a[x][y]) {
    return 0;
  }
  if (ds[i][x][st] != -1) {
    return ds[i][x][st];
  }
  if ((i == x && j == y - 1) || (i == x - 1 && j == y)) {
    return 1;
  }
  if (i == x && j == y) {
    return 1;
  }
  long long int ans = 0;
  ans += dfs(i + 1, j, x - 1, y, st + 1);
  ans %= 1000000007;
  ans += dfs(i, j + 1, x - 1, y, st + 1);
  ans %= 1000000007;
  ans += dfs(i + 1, j, x, y - 1, st + 1);
  ans %= 1000000007;
  ans += dfs(i, j + 1, x, y - 1, st + 1);
  ans %= 1000000007;
  ds[i][x][st] = ans;
  return ans;
}
int main() {
  cin >> n >> m;
  for (int i = 0; i < 502; i++) {
    for (int j = 0; j < 502; j++) {
      a[i][j] = '*';
      for (int k = 0; k < 251; k++) {
        ds[i][j][k] = -1;
      }
    }
  }
  for (int i = 1; i <= n; i++) {
    for (int j = 1; j <= m; j++) {
      cin >> a[i][j];
    }
  }
  int ans = dfs(1, 1, n, m, 0);
  ans %= 1000000007;
  cout << ans << endl;
  return 0;
}
