// Slower C++ solution:
#include<bits/stdc++.h>

#define ll long long

using namespace std;

#define mod 1000000007

ll h,w;

char ar[1001][1001];

ll dp[1001][1001];



ll path(ll i,ll j){

    if(i==h-1&&j==w-1) return 1;

    if(dp[i][j]!=-1) return dp[i][j];

    ll x=0,y=0;

    if(i+1<h&&ar[i+1][j]!='#') x=path(i+1,j);

    if(j+1<w&&ar[i][j+1]!='#') y=path(i,j+1);

    return dp[i][j]=(x+y)%mod;

}



int main()

{

    ios_base::sync_with_stdio(false);

    cin.tie(NULL);

    //freopen("input.txt", "r", stdin);

    //freopen("output.txt", "w", stdout);



    cin>>h>>w;

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

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

            cin>>ar[i][j];

        }

    }

    memset(dp,-1,sizeof(dp));

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

    return 0;

}


// Faster C++ solution for the same problem (measured speedup: 9.366425119951955x):
#include <iostream>

#include <string>

#include<iomanip>

#include <algorithm>

#include <map>

#include <vector>

#include<queue>

#include<math.h>

#include<cmath>

#include<sstream>

#include<set>

#include <map>

#include<string.h>

#include<sstream>

#include<stack>

#include<unordered_map>

#include <iostream>

#include <stdio.h>

#include <stdlib.h>

#include <ctype.h>

#include <cstdio>

#include <cstdlib>

#include <complex>

#define ll long long

#define  point complex<double>

#define  pi acos(-1)

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

#define arr2dof vector  vector<vector<int>>dp(N,vector<int>(2, inf));

using namespace std;

const ll N = 2e5 + 5, mod = 1e9 + 7, inff = -1e12,inf=1e9+7;

const double ep = 0.0000000001;

//const double pi = 3.14159265359;



double dot(point a, point b) {

	return  (conj(a)*b).real();

}

double cross(point a, point b) {

	return  (conj(a)*b).imag();

}

double angle(point a, point b)

{

	double dt = dot(a, b);

	dt /= abs(a);

	dt /= abs(b);

	return acos(dt) * 180.0 / pi;

}



bool isrightangel(int ind) {//90 degree

	vector<point>v;

	return(angle(v[ind] - v[(ind + 1) % 3], v[ind] - v[(ind + 2) % 3]) == 90);//bshof hal eldot producty3ny fe

	//dl3en 3moden 3ashn yb2a right trinagle

}



ll fact[1000006], inv[1000006];



ll fp(ll base, ll exp)

{

	if (exp == 0)

		return 1;

	ll ans = fp(base, exp / 2);

	ans = (ans * ans) % mod;

	if (exp % 2 != 0)

		ans = (ans * (base % mod)) % mod;

	return ans;

}



void calcFacAndInv(ll n)

{

	fact[0] = inv[0] = 1;

	for (ll i = 1; i <= n; i++)

	{

		fact[i] = (i * fact[i - 1]) % mod;

		inv[i] = fp(fact[i], mod - 2);

	}

}



ll ncr(ll n, ll r)

{

	return ((fact[n] * inv[r]) % mod * inv[n - r]) % mod;

}



ll npr(ll n, ll r)

{

	return (fact[n] * inv[n - r]) % mod;

}

	







char arr[1005][1005];

int dp[1005][1005];

int h, w;





int main()

{

	cin >> h >> w;

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

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

			cin >> arr[i][l];

		}

	}

	dp[0][0] = 1;

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

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

			if (i == 0 && l == 0)

				continue;

			else if (arr[i][l] == '#')

				dp[i][l] = 0;

			else {

				if (i - 1 >= 0)

					dp[i][l] = dp[i - 1][l] ;

				else

					dp[i][l] = 0;

				if (l- 1 >= 0)

					dp[i][l] += dp[i][l-1] ;

				dp[i][l] %= mod;

					

			}

		}

	}

	cout << dp[h - 1][w - 1];

	return 0;

}