#include <bits/stdc++.h>
#define IOS ios::sync_with_stdio(false);cin.tie(0);cout.tie(0);
#define i64 long long
#define endl '\n'
#define lb(x) ((x) & (-x))
using namespace std;
const i64 mod = 1e9+7;
const int maxn = 1e4+10;
int dir[4][2] = {-1,0,1,0,0,-1,0,1};
signed main()
{
    IOS;
    int n,m; cin >> n >> m;
    vector<vector<int> > v(n+1,vector<int> (m+1)),deg;
    deg = v;
    vector<vector<vector<i64>>> cnt(n+1,vector<vector<i64>> (m+1,vector<i64> (4)));
    for(int i = 1; i <= n; ++i)
    {
        for(int j = 1; j <= m; ++j) cin >> v[i][j];
    }
    for(int i = 1; i <= n; ++i)
    {
        for(int j = 1; j <= m; ++j)
        {
            bool ok = 1;
            for(int k = 0; k < 4; ++k)
            {
                int x = i + dir[k][0],y = j + dir[k][1];
                if(x < 1 || x > n || y < 1 || y > m) continue;
                if(v[x][y] - v[i][j] == 1) ok = 0;
            }
            if(ok) cnt[i][j][0] = 1;
        }
    }
    auto dfs = [&](auto self,int x,int y) -> vector<i64>
    {
        if(cnt[x][y][0] || cnt[x][y][1] || cnt[x][y][2] || cnt[x][y][3]) return cnt[x][y];
        for(int i = 0; i < 4; ++i)
        {
            int xx = x + dir[i][0],yy = y + dir[i][1];
            if(xx < 1 || xx > n || yy < 1 || yy > m || v[xx][yy] - v[x][y] != 1) continue;
            deg[xx][yy] = 1;
            vector<i64> s = self(self,xx,yy);
            cnt[x][y][1] = (cnt[x][y][1] + s[0]) % mod;
            cnt[x][y][2] = (cnt[x][y][2] + s[1]) % mod;
            cnt[x][y][3] = (cnt[x][y][3] + s[2] + s[3]) % mod;
        }
        return cnt[x][y];
    };
    for(int i = 1; i <= n; ++i)
    {
        for(int j = 1; j <= m; ++j)
        {
            dfs(dfs,i,j);
        }
    }
    i64 ans = 0;
    for(int i = 1; i <= n; ++i)
    {
        for(int j = 1; j <= m; ++j)
        {
            if(deg[i][j] == 0) ans = (ans + cnt[i][j][3]) % mod;
        }
    }
    cout << ans << endl;
    return 0;
}