// Slower C++ solution:


/**

    And in the end, the love you take is equal to the love you make.



        I LOVE YOU LAURA.

**/



#include <bits/stdc++.h>

#include <ext/pb_ds/assoc_container.hpp>

#include <ext/pb_ds/tree_policy.hpp>

#include <ext/rope>



#define _CRT_SECURE_NO_DEPRECATE

#define _CRT_DISABLE_PERFCRIT_LOCKS



#define findOrder find_by_order

#define orderKey order_of_key



#define isOn(mask, i) (mask & (1 << i))

#define turn(mask, i) (mask |= (1 << i))

#define switch(mask, i) (mask &= ~(1 << i))

#define change(mask, i) (mask ^= (1 << i))

#define lowbit(mask) (mask & (-mask))

#define turnAll(mask, n) (mask = (1 << n) - 1)

#define switchLastBit(mask) (mask & (mask - 1))

#define turnLastZero(mask) (mask | (mask + 1))

#define modulo(x, mod) ((x) & (mod - 1))

#define log2(x) (31 - __builtin_clz(x))



#define mkp make_pair

#define pb push_back

#define pf push_front

#define eb emplace_back

#define ef emplace_front

#define ins insert

#define rsz resize

#define f first

#define s second

#define lb lower_bound

#define ub upper_bound



#define sz(x) ((int)x.size())

#define all(x) begin(x), end(x)

#define rall(x) rbegin(x), rend(x)



#define debug(x) (cerr << x << " ")



using namespace std;

using namespace __gnu_pbds;

using namespace __gnu_cxx;



typedef double db;

typedef long long ll;

typedef long double ld;

typedef unsigned long long ull;



typedef pair<int, int> ii;

typedef pair<ll,ll> pll;

typedef pair<ld,ld> pdb;

typedef pair<int, ii> tri;

typedef complex<ld> cd;



typedef vector<int> vi;

typedef vector<ii> vii;

typedef vector<ld> vdb;

typedef vector<cd> vcd;

typedef vector<ll> vll;

typedef vector<string> vs;

typedef vector<vll> matrix;



typedef tree<int, null_type, less<int>, rb_tree_tag, tree_order_statistics_node_update> indexed_tree;



const ll mod = 1e9+7;

const int maxn = 1005;

const int pi = 2 * acos(0.0);

const db eps = 1e-9;



int n, m;

char arr[maxn][maxn];



bool ok (int x, int y){

    return (x >= 1 && x <= n && y >= 1 && y <= m);

}



int dx[] = {0, 1};

int dy[] = {1, 0};



ll dp[maxn][maxn];



ll dfs (int x, int y){

    if (!ok(x, y) || arr[x][y] == '#') return 0LL;

    if (x == n && y == m){

        return 1LL;

    }



    ll &ans = dp[x][y];



    if (ans != -1) return ans % mod;



    ans = 0;

    for (int i = 0; i < 2; i++){

        int nx = x + dx[i];

        int ny = y + dy[i];

        ans = (ans + dfs(nx, ny)) % mod;

    }



    return ans % mod;

}



int main()

{

    ios_base::sync_with_stdio(0);

    cin.tie(0); cout.tie(0);



    memset(dp, -1, sizeof dp);



    cin >> n >> m;



    for (int i = 1; i <= n; i++){

        for (int j = 1; j <= m; j++){

            cin >> arr[i][j];

        }

    }



    cout << dfs(1, 1) % mod << "\n";



    return 0;

}


// Faster C++ solution for the same problem (measured speedup: 9.494822032690509x):
#include <bits/stdc++.h>

using namespace std;

int const mod = 1e9+7;



int h,w;

bool a[1001][1001];

int dp[1001][1001];



int solve(int x, int y){

    if (y == w-1 && x == h-1) return 1;

    else if (!a[x][y]) return dp[x][y] = 0;

    else if (dp[x][y] != -1) return dp[x][y];

    int aux = 0,aux2 = 0;

    if (x < h-1) aux = solve(x+1,y);

    if (y < w-1) aux2 = solve(x,y+1);

    return dp[x][y] = (aux+aux2)%mod;

}



int main(int argc, char const *argv[]) {

    cin >> h >> w;

    for (int i = 0; i < h; i++){

        for (int j = 0; j < w; j++){

            char c; cin >> c;

            if (c == '.') a[i][j] = 1;

            else a[i][j] = 0;

            dp[i][j] = -1;

        }

    }

    cout << solve(0,0) << endl;

    return 0;

}