#include <bits/stdc++.h>
using namespace std;
const long long mod = 1e9 + 7;
struct node {
  int x, y;
  node() {}
  node(int _x, int _y) {
    x = _x;
    y = _y;
  }
};
int n, m;
char mp[510][510];
int dp[2][510][510];
vector<node> vec[1010];
int num[510][510];
bool in(int x, int y) { return 1 <= x && x <= n && 1 <= y && y <= m; }
int main() {
  scanf("%d%d", &n, &m);
  for (int i = 1; i <= n; i++) {
    scanf("%s", mp[i] + 1);
  }
  for (int i = 1; i <= n; i++) {
    for (int j = 1; j <= m; j++) {
      vec[i + j - 1].push_back(node(i, j));
      num[i][j] = vec[i + j - 1].size() - 1;
    }
  }
  int now = 1;
  if (mp[1][1] == mp[n][m])
    dp[1][0][0] = 1;
  else
    dp[1][0][0] = 0;
  for (int i = 2; i <= (n + m - 1 + 1) / 2; i++) {
    now ^= 1;
    for (int j = 0; j < vec[i].size(); j++) {
      for (int k = 0; k < vec[n + m - i].size(); k++) {
        int x1 = vec[i][j].x;
        int y1 = vec[i][j].y;
        int x2 = vec[n + m - i][k].x;
        int y2 = vec[n + m - i][k].y;
        dp[now][j][k] = 0;
        if (mp[x1][y1] == mp[x2][y2]) {
          if (in(x1 - 1, y1) && in(x2 + 1, y2)) {
            dp[now][j][k] += dp[now ^ 1][num[x1 - 1][y1]][num[x2 + 1][y2]];
          }
          dp[now][j][k] %= mod;
          if (in(x1, y1 - 1) && in(x2 + 1, y2)) {
            dp[now][j][k] += dp[now ^ 1][num[x1][y1 - 1]][num[x2 + 1][y2]];
          }
          dp[now][j][k] %= mod;
          if (in(x1 - 1, y1) && in(x2, y2 + 1)) {
            dp[now][j][k] += dp[now ^ 1][num[x1 - 1][y1]][num[x2][y2 + 1]];
          }
          dp[now][j][k] %= mod;
          if (in(x1, y1 - 1) && in(x2, y2 + 1)) {
            dp[now][j][k] += dp[now ^ 1][num[x1][y1 - 1]][num[x2][y2 + 1]];
          }
          dp[now][j][k] %= mod;
        } else {
          dp[now][j][k] = 0;
        }
      }
    }
  }
  int ans = 0;
  if ((n + m - 1) % 2 == 0) {
    int kkk = ((n + m - 1) / 2) & 1;
    for (int i = 0; i < vec[(n + m - 1) / 2].size(); i++) {
      int x = vec[(n + m - 1) / 2][i].x;
      int y = vec[(n + m - 1) / 2][i].y;
      if (in(x + 1, y)) {
        ans += dp[kkk][i][num[x + 1][y]];
        ans %= mod;
      }
      if (in(x, y + 1)) {
        ans += dp[kkk][i][num[x][y + 1]];
        ans %= mod;
      }
    }
  } else {
    int kkk = ((n + m) / 2) & 1;
    for (int i = 0; i < vec[(n + m) / 2].size(); i++) {
      int x = vec[(n + m) / 2][i].x;
      int y = vec[(n + m) / 2][i].y;
      ans += dp[kkk][i][i];
      ans %= mod;
    }
  }
  printf("%d\n", ans);
  return 0;
}
