def main():
    from math import sqrt,sin,cos,tan,ceil,radians,floor,gcd,exp,log,log10,log2,factorial,fsum
    from heapq import heapify, heappop, heappush
    from bisect import bisect_left, bisect_right
    from copy import deepcopy
    import copy
    import random
    from collections import deque,Counter,defaultdict
    from itertools import permutations,combinations
    from decimal import Decimal,ROUND_HALF_UP
    #tmp = Decimal(mid).quantize(Decimal('0'), rounding=ROUND_HALF_UP)
    from functools import lru_cache, reduce
    #@lru_cache(maxsize=None)
    from operator import add,sub,mul,xor,and_,or_,itemgetter
    INF = 10**18
    mod1 = 10**9+7
    mod2 = 998244353
    
    #DecimalならPython
    #再帰ならPython!!!!!!!!!!!!!!!!!!!!!!!!!!
    
    
    '''
    または、か
    
    0が含まれてるパターンがえぐい
    
    二次元を1次元に変えてしまえばセグメントツリーが使えるのかな
    
    だから田なのか
    '''
    
    H,W = map(int, input().split())
    A = [list(map(int, input().split())) for _ in range(H)]
    
    acm1 = [[1]*(W+1) for _ in range(H+1)]
    acm2 = [[1]*(W+1) for _ in range(H+1)]
    acm3 = [[1]*(W+1) for _ in range(H+1)]
    acm4 = [[1]*(W+1) for _ in range(H+1)]
    
    for i in range(H):
        for j in range(W):
            acm1[i][j] = A[i][j]
            acm2[i][j] = A[i][j]
            acm3[i][j] = A[i][j]
            acm4[i][j] = A[i][j]
    
    for i in range(H):
        for j in range(W-1):
            acm1[i][j+1] *= acm1[i][j]
            acm1[i][j+1] %= mod1
    for j in range(W):
        for i in range(H-1):
            acm1[i+1][j] *= acm1[i][j]
            acm1[i+1][j] %= mod1
    
    for i in reversed(range(H)):
        for j in range(W-1):
            acm3[i][j+1] *= acm3[i][j]
            acm3[i][j+1] %= mod1
    for j in range(W):
        for i in reversed(range(1,H)):
            acm3[i-1][j] *= acm3[i][j]
            acm3[i-1][j] %= mod1
    
    for i in range(H):
        for j in reversed(range(1,W)):
            acm2[i][j-1] *= acm2[i][j]
            acm2[i][j-1] %= mod1
    for j in reversed(range(W)):
        for i in range(H-1):
            acm2[i+1][j] *= acm2[i][j]
            acm2[i+1][j] %= mod1
    
    for i in reversed(range(H)):
        for j in reversed(range(1,W)):
            acm4[i][j-1] *= acm4[i][j]
            acm4[i][j-1] %= mod1
    for j in reversed(range(W)):
        for i in reversed(range(1,H)):
            acm4[i-1][j] *= acm4[i][j]
            acm4[i-1][j] %= mod1
    
    Q = int(input())
    for _ in range(Q):
        r,c = map(int, input().split())
        
        print(acm1[r-2][c-2]%mod1*acm2[r-2][c]%mod1*acm3[r][c-2]%mod1*acm4[r][c]%mod1)
    
    
    
    
    
    
    
    
    
    
    
    
    
    
    
    
    
    
    
    
    
    
    
    
    
    
    
    
if __name__ == '__main__':
    main()