import sys

input=sys.stdin.readline

mod=10**9+7

H,W=map(int,input().split())
A=[list(map(int,input().split())) for i in range(H)]

UL=[[A[i][j] for j in range(W)] for i in range(H)]
for i in range(H):
    for j in range(1,W):
        UL[i][j]=(UL[i][j-1]*UL[i][j])%mod
for j in range(W):
    for i in range(1,H):
        UL[i][j]=(UL[i][j]*UL[i-1][j])%mod

UR=[[A[i][j] for j in range(W)] for i in range(H)]
for i in range(H):
    for j in range(W-2,-1,-1):
        UR[i][j]=(UR[i][j+1]*UR[i][j])%mod
for j in range(W):
    for i in range(1,H):
        UR[i][j]=(UR[i][j]*UR[i-1][j])%mod

DL=[[A[i][j] for j in range(W)] for i in range(H)]
for i in range(H):
    for j in range(1,W):
        DL[i][j]=(DL[i][j]*DL[i][j-1])%mod
for i in range(W):
    for j in range(H-2,-1,-1):
        DL[j][i]=(DL[j][i]*DL[j+1][i])%mod

DR=[[A[i][j] for j in range(W)] for i in range(H)]
for i in range(H):
    for j in range(W-2,-1,-1):
        DR[i][j]=(DR[i][j]*DR[i][j+1])%mod
for i in range(W):
    for j in range(H-2,-1,-1):
        DR[j][i]=(DR[j][i]*DR[j+1][i])%mod

for _ in range(int(input())):
    r,c=map(int,input().split())
    r-=1;c-=1
    res=1
    if r!=0 and c!=0:
        res=(res*UL[r-1][c-1])%mod
    if r!=0 and c!=W-1:
        res=(res*UR[r-1][c+1])%mod
    if r!=H-1 and c!=0:
        res=(res*DL[r+1][c-1])%mod
    if r!=H-1 and c!=W-1:
        res=(res*DR[r+1][c+1])%mod
    print(res)