#include <cstdio>
#include <iostream>
#include <sstream>
#include <fstream>
#include <iomanip>
#include <algorithm>
#include <cmath>
#include <string>
#include <vector>
#include <list>
#include <queue>
#include <stack>
#include <set>
#include <map>
#include <bitset>
#include <numeric>
#include <limits>
#include <climits>
#include <cfloat>
#include <functional>
using namespace std;

const int MOD = 1000000007;

int main()
{
    int h, w;
    cin >> h >> w;
    vector<string> s(h);
    for(int i=0; i<h; ++i)
        cin >> s[i];

    int ans = 0;
    for(int i=0; i<2; ++i){
        vector<bool> a(w, false), b(w, false);
        for(int y=0; y<h; ++y){
            for(int x=0; x<w; ++x){
                if(s[y][x] != '?'){
                    if((s[y][x] == '0') ^ (y % 2 == 0))
                        a[x] = true;
                    else
                        b[x] = true;
                }
            }
        }

        int cnt = 1;
        for(int x=0; x<w; ++x){
            if(a[x] && b[x]){
                cnt *= 0;
            }
            else if(!a[x] && !b[x]){
                cnt *= 2;
                cnt %= MOD;
            }
        }
        ans += cnt;
        ans %= MOD;

        vector<string> t(w, string(h, ' '));
        for(int y=0; y<h; ++y){
            for(int x=0; x<w; ++x){
                t[x][y] = s[y][x];
            }
        }
        s.swap(t);
        swap(h, w);
    }

    bool a = true;
    bool b = true;
    for(int y=0; y<h; ++y){
        for(int x=0; x<w; ++x){
            if(s[y][x] != '?'){
                if((s[y][x] == '0') ^ ((y + x) % 2 == 0))
                    a = false;
                else
                    b = false;
            }
        }
    }
    if(a)
        ans = (ans - 1 + MOD) % MOD;
    if(b)
        ans = (ans - 1 + MOD) % MOD;

    cout << ans << endl;

    return 0;
}