import math,string,itertools,fractions,heapq,collections,re,array,bisect,sys,random,time
sys.setrecursionlimit(10**7)
inf = 10**20
mod = 10**9 + 7
def LI(): return list(map(int, input().split()))
def LF(): return list(map(float, input().split()))
def LS(): return input().split()
def I(): return int(input())
def F(): return float(input())
def S(): return input()
def main():
h,w = LI()
a = [LI() for _ in range(h)]
r = 0
ff = [[False]*w for _ in range(h)]
dy = [0,0,1,-1]
dx = [1,-1,0,0]
m = {}
def f(y,x):
if (y,x) in m:
return m[(y,x)]
r = 1
k = a[y][x]
for i in range(4):
ty,tx = y+dy[i], x+dx[i]
if ty < 0 or ty >= h or tx < 0 or tx >= w:
continue
if a[ty][tx] >= k:
continue
r += f(ty,tx)
m[(y,x)] = r
return r
rr = 0
for i in range(h):
for j in range(w):
rr += f(i,j)
rr %= mod
return rr
print(main())