#include<set>
#include<map>
#include<deque>
#include<queue>
#include<stack>
#include<cmath>
#include<ctime>
#include<bitset>
#include<string>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<climits>
#include<complex>
#include<iostream>
#include<algorithm>
#define ll long long
using namespace std;

const ll Mod = 1e9+7;
const int maxn = 2100;

int n,m;
ll c[maxn][maxn],pc[maxn][maxn];
ll f[maxn][maxn],g[maxn][maxn];

int main()
{
	scanf("%d%d",&n,&m);
	c[0][0]=1ll; pc[0][0]=1ll;
	for(int i=1;i<maxn;i++)
	{
		c[i][0]=1ll; pc[i][0]=1ll;
		for(int j=1;j<=i;j++) 
			c[i][j]=(c[i-1][j-1]+c[i-1][j])%Mod,
			pc[i][j]=(pc[i][j-1]+c[i][j])%Mod;
	}
	for(int i=0;i<=n;i++) f[i][0]=g[i][0]=1ll;
	for(int i=0;i<=n;i++)
	{
		for(int j=1;j<=m;j++)
			f[i][j]=pc[j-1][min(i,j-1)];
	}
	for(int j=1;j<=m;j++)
	{
		g[0][j]=f[0][j];
		for(int i=1;i<=n;i++) g[i][j]=(g[i-1][j]+f[i][j])%Mod;
	}
	ll ans=1;
	for(int i=n;i>=1;i--)
	{
		for(int y=0,u=min(m-1,i);y<=u;y++)
		{
			int x=i-y; //(i,m-1) -> (x,y)
			ans+=g[x][y]*c[i+m-1-x-y][i-x]%Mod;
		}
	}
	ans%=Mod;
	printf("%lld\n",ans);
	
	return 0;
}
