// 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.49618984552235x):
#include<bits/stdc++.h>

#define int long long

#define ll long long

#define ld long double

#define mod ((ll)1e9+7)

#define maxx ((ll)1e6+5)

#define pb push_back

#define ins insert

#define usefile freopen("input.txt","r",stdin); freopen("output.txt","w",stdout);

#define fast ios::sync_with_stdio(0); cin.tie(0); cout.tie(0);

#define test ll t; cin>>t; while(t--)

ll exp(ll x,ll y,ll p){ ll res=1; while(y) {if(y%2) res=(res*x%p)%p; x=(x*x)%p; y/=2; } return res;}

ll expo(ll x,ll y){ ll res=1; while(y) {if(y%2) res=(res*x%mod)%mod; x=(x*x)%mod; y/=2; } return res;}

/*int sv[1000000]={0};

ll fact[1000007]={0};

void facto() {fact[0]=1;fact[1]=1;for(int i=2;i<100007;i++)fact[i]=(fact[i-1]*i)%mod;}

void sve() {sv[0]=1;sv[1]=1;for(int i=2;i*i<=maxx;i++){if(!sv[i]){for(int j=2*i;j<=maxx;j+=i)sv[j]=1;}}}//0 hai to prime

ll ncr(ll n,ll r) {ll res=1; res=fact[n]; res=(res*(expo(fact[r],mod-2)))%mod; res=(res*(expo(fact[n-r],mod-2)))%mod; return res;}

ll npr(ll n,ll r) {ll res=1; res=fact[n]; res=(res*(expo(fact[n-r],mod-2)))%mod; return res;}*/

int d[1003][1003];

using namespace std;

signed main()

{

    int r,c,i,j;

    cin >> r >> c;

    char a[r][c];

    for(i=0 ; i<r ; i++)

    	for(j=0 ; j<c ; j++)

    		cin >> a[i][j];

    d[0][0] = 1;

    for(i=0 ; i<r ; i++)

    	for(j=0 ; j<c ; j++)

    	{

    		if(a[i][j] == '#')

    			continue;

    		if(j>0) d[i][j] = (d[i][j] + d[i][j-1])%mod;

    		if(i>0) d[i][j] = (d[i][j] + d[i-1][j])%mod;

    	}

    cout << d[r-1][c-1] <<"\n";

    return 0;

}