#include <bits/stdc++.h>
using namespace std;
#define rep(i, n) for (int i = 0; i < (int)(n); i++)
#define repr(i, n) for (int i = (int)(n); i >= 0; i--)
#define all(v) v.begin(), v.end()
#define mod1 1000000007
#define mod2 998244353
typedef long long ll;

int main(){
    ll H,W;
    cin >> H >> W;
    
    vector<string> S(H);
    rep(i,H){
        cin >> S[i];
    }

    ll len = H + W - 2;

    vector<vector<vector<ll>>> dp(len / 2 + 1, vector<vector<ll>>(H, vector<ll>(H, 0)));
    if (S[0][0] == S[H - 1][W - 1]){
        dp[0][0][0] = 1;
    }
    else{
        dp[0][0][0] = 0;
    }

    for (ll i = 0; i < len / 2; i++){
        //i + 1: 進んだ文字の数．j: 縦．(i + 1 - j)：横．H - 1 - k: 縦．
        for (ll j = 0; j <= i; j++){
            for (ll k = 0; k <= i; k++){
                if (j < H && i + 1 - j < W && H - 1 - k >= 0 && W - 2 - i + k >= 0){//横，横
                    if (S[j][i + 1 - j] == S[H - 1 - k][W - 1 - (i + 1 - k)]){
                        dp[i + 1][j][k] = (dp[i + 1][j][k] + dp[i][j][k]) % mod1;
                    }
                }
                if (j < H && i + 1 - j < W && H - 1 - k - 1 >= 0 && W - 1 - (i + 1 - k - 1) >= 0){//横，縦
                    if (S[j][i + 1 - j] == S[H - 1 - k - 1][W - 1 - i + k]){
                        dp[i + 1][j][k + 1] = (dp[i + 1][j][k + 1] + dp[i][j][k]) % mod1;
                    }
                }
                if (j + 1 < H && i - j < W && H - 1 - k >= 0 && W - 1 - (i + 1 - k) >= 0){//縦，横
                    if (S[j + 1][i - j] == S[H - 1 - k][W - 1 - (i + 1 - k)]){
                        dp[i + 1][j + 1][k] = (dp[i + 1][j + 1][k] + dp[i][j][k]) % mod1;
                    }
                }
                if (j + 1 < H && i - j < W && H - 1 - k - 1 >= 0 && W - 1 - (i + 1 - k - 1) >= 0){//縦，縦
                    if (S[j + 1][i - j] == S[H - 1 - k - 1][W - 1 - (i + 1 - k - 1)]){
                        dp[i + 1][j + 1][k + 1] = (dp[i + 1][j + 1][k + 1] + dp[i][j][k]) % mod1;
                    }
                }
            }
        }
    }
    /*for (ll i = 0; i <= len / 2; i++){
        for (ll j = 0; j < H; j++){
            for (ll k = 0; k < H; k++){
                cout << dp[i][j][k] << " ";
            }
            cout << endl;
        }
    }*/
    
    ll ans = 0;
    if (len % 2 == 0){
        for (ll j = 0; j < H; j++){
            for (ll k = 0; k < H; k++){
                if (j + k == H - 1) ans = (ans + dp[len / 2][j][k]) % mod1;
            }
        }
    }
    else{
        for (ll j = 0; j < H; j++){
            for (ll k = 0; k < H; k++){
                if (j + k == H - 1 || j + k == H - 2) ans = (ans + dp[len / 2][j][k]) % mod1;
            }
        }
    }
    cout << ans << endl;
}