import numpy as np

MOD = 1000000007
mul = np.frompyfunc(lambda x, y: x * y % MOD, 2, 1)
int1 = lambda x: int(x) - 1

H, W = map(int, input().split())
A = np.array([tuple(map(int, input().split())) for _ in range(H)], dtype=np.object)

Au = mul.accumulate(A, axis=0)
Ad = np.flipud(mul.accumulate(np.flipud(A), axis=0))

Aul = mul.accumulate(Au, axis=1).tolist()
Aur = np.fliplr(mul.accumulate(np.fliplr(Au), axis=1)).tolist()
Adl = mul.accumulate(Ad, axis=1).tolist()
Adr = np.fliplr(mul.accumulate(np.fliplr(Ad), axis=1)).tolist()

Q = int(input())
for _ in range(Q):
    ans = 1
    r, c = map(int1, input().split())
    if 0 < r and 0 < c:
        ans = ans * Aul[r - 1][c - 1] % MOD
    if 0 < r and c < W - 1:
        ans = ans * Aur[r - 1][c + 1] % MOD
    if r < H - 1 and 0 < c:
        ans = ans * Adl[r + 1][c - 1] % MOD
    if r < H - 1 and c < W - 1:
        ans = ans * Adr[r + 1][c + 1] % MOD
    print(ans)