# include <bits/stdc++.h>
using namespace std;
using ll=long long;
constexpr int N=144,mod=1e9+7;
int pw2[150],fac[150],finv[150],f[2][80][5010];
char s[150][150];
inline int add(int x,int y){return (x+=y)<mod?x:x-mod;}
inline void upd(int &x,const int &y){if((x+=y)>=mod) x-=mod;}
inline void upd(int &x,const ll &y){x=(x+y)%mod;}
inline int P(int n,int m){return (ll)fac[n]*finv[n-m]%mod;}
int power(int a,int b)
{
	int ans=1;
	for(;b;b>>=1,a=(ll)a*a%mod)
		if(b&1) ans=(ll)ans*a%mod;
	return ans;
}
int main()
{
	int T,n,m,K;
	pw2[0]=fac[0]=fac[1]=1;
	for(int i=1;i<=N;i++) pw2[i]=add(pw2[i-1],pw2[i-1]);
	for(int i=2;i<=N;i++) fac[i]=(ll)fac[i-1]*i%mod;
	finv[N]=power(fac[N],mod-2);
	for(int i=N-1;i>=0;i--) finv[i]=(ll)finv[i+1]*(i+1)%mod;
	cin>>T;
	while(T--)
	{
		scanf("%d%d%d",&n,&m,&K);
		for(int i=1;i<=n;i++) scanf("%s",s[i]+1);
		if(n<m)
		{
			for(int i=1;i<=n;i++)
				for(int j=i+1;j<=m;j++)
					swap(s[i][j],s[j][i]);
			swap(n,m);
		}
		int cnt=0,cur=0,ans=0;
		for(int i=1;i<=n;i++)
			for(int j=1;j<=m;j++)
				cnt+=(s[i][j]=='.');
		f[cur][0][0]=1;
		for(int i=1;i<=n;i++)
			for(int j=1;j<=m;j++)
			{
				cur^=1;
				for(int k=0;k<=K;k++)
					for(int l=0;l<(1<<m);l++)
					{
						int &x=f[cur^1][k][l];
						if(!x) continue;
						if(s[i][j]=='X') upd(f[cur][k][l&~(1<<j-1)],x);
						else
						{
							upd(f[cur][k][l|(1<<j-1)],x);
							if(k<K && j>1 && (l>>j-2)&1) upd(f[cur][k+1][l&~(3<<j-2)],x);
							if(k<K && (l>>j-1)&1) upd(f[cur][k+1][l&~(1<<j-1)],x);
						}
						x=0;
					}
			}
		for(int i=0;i<=K;i++)
			for(int j=0;j<(1<<m);j++)
			{
				int &x=f[cur][i][j];
				upd(ans,(ll)x*(i&1?mod-pw2[i]:pw2[i])%mod*P(K,i)%mod*P(cnt-2*i,2*(K-i)));
				x=0;
			}
		printf("%d\n",ans);
	}
	return 0;
}