#include<bits/stdc++.h>
using namespace std;
#define int long long
#define endl '\n'
int go[4][2]={{0,1},{1,0},{-1,0},{0,-1}};
const int MOD=1e9+7;
const int N=1005;
int dp[5][N][N];
int n,m;
struct node{
	int x,y;
};
bool in(int x,int y){
	if(x>=0&&x<n&&y>=0&&y<m) return true;
	return false;
}
void solve(){
	cin>>n>>m;
	int a[n][m],c[n][m];
	for(int i=0;i<n;i++){
		for(int j=0;j<m;j++){
			for(int k=0;k<5;k++) dp[k][i][j]=0;
			cin>>a[i][j];
			c[i][j]=0;
		}
	}
	queue<node>q;
	for(int i=0;i<n;i++){
		for(int j=0;j<m;j++){
			for(int k=0;k<4;k++){
				int x=i+go[k][0],y=j+go[k][1];
				if(!in(x,y)) continue;
				if(a[x][y]==a[i][j]-1) c[i][j]++;
			}
			if(c[i][j]==0){
				node tmp;
				tmp.x=i,tmp.y=j;
				dp[1][i][j]=1;
				q.push(tmp);
			} 
		}
	}
	while(!q.empty()){
		node tmp=q.front();
		q.pop();
		int x=tmp.x,y=tmp.y;
		for(int i=0;i<4;i++){
			int x1=x+go[i][0],y1=y+go[i][1];
			if(!in(x1,y1)||a[x1][y1]!=a[x][y]+1) continue;
			c[x1][y1]--;
			for(int j=1;j<=4;j++){
				dp[j][x1][y1]=(dp[j][x1][y1]+dp[j-1][x][y])%MOD;
			}
			dp[4][x1][y1]=(dp[4][x1][y1]+dp[4][x][y])%MOD;
			if(c[x1][y1]==0){
				node tmp;
				tmp.x=x1,tmp.y=y1;
				q.push(tmp);
			}
		}
	}
	int ans=0;
	for(int i=0;i<n;i++){
		for(int j=0;j<m;j++){
			bool ok=true;
			for(int k=0;k<4;k++){
				int x=i+go[k][0],y=j+go[k][1];
				if(!in(x,y)) continue;
				if(a[x][y]==a[i][j]+1) ok=false;
			}
			if(ok) ans=(ans+dp[4][i][j])%MOD;
		}
	}
	cout<<ans<<endl;
}
signed main(){
	ios::sync_with_stdio(false);
	cin.tie(0),cout.tie(0);
	int t=1;
//	cin>>t;
	while(t--){
		solve();
	}
	return 0;
}